similar to: anova

I would like to do a one way anova and then a tukeyHSD. I have three vectors A,B and C. In a previous help message, I was told to do the following for the anova: y = c(A,B,C) group = factor(rep(a:3,c(7,9,13))) #provided there a 7 elements in A,9 in B and 13 in C and then anova(lm(y~group)) Looking at the tukeyHSD method it looks like it wants the aov method which I don't understand.

Does anyone know of some methods already programmed up ( ie freeware : ) ) in R for the post hoc methods in anova. Particularly Scheffe's method or Tukey's? Anna

Please excuse me if this is a very basic question and I am just misunderstanding the basic structure of R. Right now, I am writing everything in a buffer (emacs) and then copying it into R when it is correct. Is there a way to save a file of commands and then just type the file name or something in R and have the commands executed within R? Anna

When I am trying to put a legend on a plot where I am using "lines", R just ignores it. I can do it with boxplot or plot, but just not with lines. Am I doing something wrong? Maybe I am just making a mistake? Anna

Is there a modulus operator in R? Anna

Is there a simpler way then the solution to the one that was posted here? I'm not very proficient with legend, and I don't understand this solution. All I have is two or more lines on one plot that I want to put a legend on and I can't figure out how to do it from the examples. Can you give a very simple example? It does not have to be fancy!! I have never worked with a

I didn't explain myself well. You are right about not needing a list on the right hand side of the equation below. What I end up with is a list of arrays. Now what i would like to do is to access the individual elements in the arrays in the list. For example, when I type, y[[1]], I just get the whole first array. How do I get the first element in the first array for instance? Is

I know how to print my own labels along the x axis: plot(1:10,xlab = "My label", axes = FALSE) axis(1,at=seq(1,10,by=2),labels=c("first","second","third","fourth","fifth")) Is there a way to make the labels print vertically along the axis rather than horizontally? (I'm not sure if this is even going to look all right, but I

I've tried looking at ?colors and ?palette and if I'm understanding it correctly, I'm supposed to type in (for example) palette(rainbow(13)) before I type in my plot (of 13 lines) if I want 13 different colors. But this does not work. Other things that i have tried besides rainbow give me errors. Am I just doing something completely wrong? Anna

I have a plot with 13 lines and R repeats colors. Is there any way to make R not repeat colors, or is this just a fundamental limitation on R's part to say have only 8 colors that it can use? Anna

I have two separate questions. The first is when I use xlab, the labels on my tic marks disappear. Is there a way to make them reappear? Also, is there a way to make my own labels for the tick marks? My second question is, after I have spent a lot of time getting my plot to look just right, I want to save it to a file, but the only way I seem to know how to do it is to first use

Is there a function that calculates the power of a repeated-measure ANOVA design, e.g., 2 groups, 4 within-subject factors, an average 0.40 correlation between the 4 within factors, etc. I don't think I can use power.anova.test() because it does not consider corr=0.40. I am hoping that someone has already implemented the method by Muller & Barton (1989, JASA 84: 549-555) on how to

I am trying to make an array of lists. I don't think that I am doing it right however because I cannot access the individual elements of the lists once I have created the array of lists. Can someone help? for(i in 1:3){ y[[i]] = list(name[((i-1)*index+1):(i*index)]) } Anna

I've been programming in one directory and recently switched to another directory. It appears that in doing so I've uncovered a problem. My environment was saving something so that my code would work and now I don't know how to fix it. I have the following bit of code: for(i in 1:index){ indexList = lst[i] for(j in 2:jobs-1){ indexList =

Dear R people: I wish to switch from SPSS to R, but there is one particular type of ANOVA design that cannot be done in R. Or more likely, it can be done, but it is nowhere documented. The problem is typical for psychologists: You have a repeated measures design with different groups of subjects. Now, this can be done with the aov command, but the number of subjects in both groups must be

I don't know how to retain my row labels. I read in a table with row labels, ie: Price Floor Area Rooms 52. 111. 830 5 54. 128. 710 5 47.4 101 1000 4 ... but then I need to transform the data raggedRidzeros <- function(x){ y = list(ridzeros(x[[1]])) for(i in 2:length(x)){ y = c(y,list(ridzeros(x[[i]]))) } y } (where ridzeros is a function

When I am trying to use the source function to read in some lines of code, I get an error. The code is simply a boxplot: boxplot(s$fitness2weightedSum,s$OffNader10xWeightweightedSum,s$OffNader5xWeightweightedSum,main="weightedSum",col=8) axis(1,at=seq(1,3,by=1),las=3,labels=c("fitness2","OffNader10xWeight","OffNader5xWeight)) The error I get is: Error in

Hi R users, I google the website and I found that there are three ways to perform repeated measure ANOVA: aov, lme and lmer. http://www.mail-archive.com/r-help at stat.math.ethz.ch/msg58502.html But the questions are which one is good to use and how to do post-hoc test? I use the example that is provided in the above link and I try > tt<-aov(p.pa~group*time+Error(subject/time),data=P.PA)

Hi All, I have a question which does not pertain directly to the use of R but comes from my use of R! I have data which can be described as 3-dimensional e.g. (x,y,z), with no negative component. The suggested way to analyze this data is via multivariate techniques or by calculating what amounts to a levene's test on the data and then an ANOVA on the three components if the first test is

Has anyone written a function that returns all possible subsets of r elements out of a vector that contains n elements? I have not been able to find an answer in the FAQ or the archives. Say I have a vector: letters[1:10]. I want to get a matrix of 120 rows that contains all possible subsets of 3 letters out of the 10. And I would like to be able to generalize to any r out of n, like