similar to: A shorter version of ".Last.value"?

Hi, This is really strange. Can anyone help explain what's going on here (on 3 and 7)? > targets <- seq(from=.1, to=.9, by=.1) > targets[1]==.1 [1] TRUE > targets[2]==.2 [1] TRUE > targets[3]==.3 [1] FALSE > targets[4]==.4 [1] TRUE > targets[5]==.5 [1] TRUE > targets[6]==.6 [1] TRUE > targets[7]==.7 [1] FALSE > targets[8]==.8 [1] TRUE >

Hi, I'm trying to randomize a sequence of trials for an experimental design. The trials consist of values for each of two factors. As is there are 30 combinations of the two factors, and I want them to be ordered randomly but with the requirement that for one of the factors, the value can never be the same as the previous value. I'm currently randomizing my dataframe by using:

Hi, I want to add error bars to a plot generated with xyplot. I've tried both errbar() and plotCI(), but in both cases the points are not in the same place. It's as if the two functions are using a different frame of reference for the plotting area. for example: means <- c(92.5, 92.25, 90.9, 91.0, 94.15, 90.05) #means time <- c(1,1,2,2,3,3) #occasion variable group <-

Hi, I often use plotmeans() from the gplots package to quickly visualize a pattern of change. I would like to be able to plot separate lines for different groups, but the function gives an error when a grouping variable is included in the formula argument. For instance, > require(gplots) > x <- data.frame(Score=rnorm(100), Time=rep(1:10, 10),

Hi, I've looked all over for a solution to this, but haven't had much look in specifying what I want to do with appropriate search terms. Thus I'm turning to R-help. In the process of trying to write a simple function to rename individual column names in a data frame, I ran into the following problem: When I rename the columns within my function, I can't seem to get it to

Hello, I make frequent use of the *step()* and, for plotting, *all.effects() *functions for *lm()* objects. I am now doing more with *lme()* random effects models, and haven't been able to find functions similar to *step()* or *all.effects()* which will accept *lme()* objects. Do similar functions for either exist? If someone could point me in the right direction it would be much

Hi, I'm having some problems with the aggregate() function in the {stats} package, and the documentation doesn't address them. 1) Why would the first line work, but the second not? According to the help file, it accepts a "data=" argument. > with(tsrc, aggregate(x=DistRatio, by=list(Condition), FUN=mean)) Group.1 x 1 Congruent 1.741789 2 Mismatch 1.771425

Hello, I am confused about the order of nesting in mixed models using functions like aov(), lme(), lmer(). I have the following data: n subjects in either condition A or B each subject tested at each of 3 numerical values ("distance" = 40,50,60), repeated 4 times for each of the 3 numerical values ("trial" = 1,2,3,4) Variable summary: Condition: 2 level factor Distance:

Hi, I've done my best to search for a solution to this, but had no luck. How can I create a lattice plot (I'm using xyplot() ) that does not have a border on the top and right side, but keeps the bottom/left axes? So far all I've found is this, which inserted into the xyplot call removes all 4 borders: / par.settings = list(axis.line = list(col = 0))/ xyplot( Sepal.Length ~

Is there a way to estimate a nominal response model? To be more specific let's say I want to calibrate: \pi_{v}(\theta_j)=\frac{e^{\xi_{v}+\lambda_{v}\theta_j}}{\sum_{h=1}^m e^{\xi_{h}+\lambda_{h}\theta_j}} Where $\theta_j$ is a the dependent variable and I need to estimate $\xi_{h}$ and $\lambda_{h}$ for $h \in {1...,m}$. Thank you, Mauricio Romero Quantil S.A.S. Cel: 3112231150

Hello, FileListBox with VisualBasic 5.0 (32Bit) and NT4.0SP3 displays only Filenames shorter than 8 character. I am using Samba1.9.18.p3. The problem doesn't occur if the connected services is on an NT-Plattform. Does anybody know a solution to this problem. Or hasn?t anybody had this problem yet and am I just goofed ? regards Michael Emrich

This is a continuation from a previous posting of mine: The following algorithm below is what I want to accomplish: Z(xk) = Average(Yi, i belongs to Ik), where Ik contains all i such that for each j, |Xi,j - xkj?? 2. Here, j = 1, 2 and i corresponds to the elements in each X and/or xk >data x1 x2 y 1 1 2 2 2 6 3 3 12 Now, consider a second data frame or matrix (xk):

Hi, How to cbind or rbind different lengths vectors/arrays without repeating the elements of the shorter vectors/arrays ? > cbind(1:2, 1:10) [,1] [,2] [1,] 1 1 [2,] 2 2 [3,] 1 3 [4,] 2 4 [5,] 1 5 [6,] 2 6 [7,] 1 7 [8,] 2 8 [9,] 1 9 [10,] 2 10 [[alternative HTML version deleted]]

Hello, I am putting a legend with lines in a line plot and I would like to make the lines in the legend shorter. Does anybody knows how to do this? Thank you Felipe Parra [[alternative HTML version deleted]]

Hello, Sorry, there was an mistake in my previous mail. FileListBox with VisualBasic 5.0 (32Bit) and NT4.0SP3 displays only files if directoryname is shorter than 9 character. I am using Samba1.9.18.p3. The problem doesn't occur if the connected services is on an NT-Plattform. Does anybody know a solution to this problem. Or hasn?t anybody had this problem yet and am I just goofed ?

Hi folks, I'm following An Introduction to R http://cran.r-project.org/doc/manuals/R-intro.html#R-and-statistics to learn R. Coming to; 2.2 Vector arithmetic > v <- 2*x + y + 1 Warning message: In 2 * x + y : longer object length is not a multiple of shorter object length What does it mean? How to rectify it? Please help. TIA B.R. Stephen L

Hi! I suspect there must be an easy way to access components of a data frame by name, i.e. the input should look like "name1 name2 name3 ..." and the output be a data frame of those components with the corresponding names. I ´ve been trying for hours, but only found the long way to do it (which is not feasible, since I have lots of components to select):

Hi, I've done some searching where others have had trouble with this error (or "warning" actually), but I'm unable to solve my problem. I have a data sheet with 13 columns and 36 rows. Each column has exactly the same number of rows. I've created glms and now want to do cross-validation on 2 of them. Please be gentle-- I'm new to R (and statistics, too, for that

In datamatrix[, "y"] == datamatrix[, "y"][-1] : longer object length is not a multiple of shorter object length out = c(FALSE,datamatrix[,'y'] == datamatrix[,'y'][-1]) and I do not know why I get that error, the resulting out matrix is somehow one row larger than datamatrix... all I try to do is filter matrix by dropping rows where [,'y'][-1] ==

Dear R users, I am using the code below to generate a fitted value of b. I have about 300 different values for for y (y1, y2, ...y300) which means I will have to write the code below 300 times to generate the 300 different fitted values for y. Is there a short way of doing that ? Many thanks in advance Mintewab library(mgcv) dat <- read.table("e:/minti's