[R] anova leads to an error

Dear R-list, the following code had been running well over the last months: exam <- matrix(rnorm(100,0,1), 10, 10) gg <- factor(c(rep("A", 5), rep("B", 5))) mlmfit <- lm(exam ~ 1); mlmfitG <- lm(exam ~ gg) result <- anova(mlmfitG, mlmfit, X=~0, M=~1) Until, all of a sudd

Re: [R] Anova and unbalanced designs

fore, this problem is so simple, that I find it hard to understand > where there's room for error, but I wanted to check against SAS to test my > sanity (a procedure that will likely get a rise out of some list members). > > Maybe you should send a message to the SPSS help list. >

Re: [R] Anova and unbalanced designs

- I hope this helps, John -- John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario, Canada web: socserv.mcmaster.ca/jfox -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]

[R] Anova and unbalanced designs

Dear R-list! My question is related to an Anova including within and between subject factors and unequal group sizes. Here is a minimal example of what I did: library(car) within1 <- c(1,2,3,4,5,6,4,5,3,2); within2 <- c(3,4,3,4,3,4,3,4,5,4) values <- data.frame(w1 = within1, w2 = within2) valu

Re: [R] assign a list using expression?

Thank you Patrick and Gabor! Sorry, I think I have not explainend it well. The purpose is as follows: names <- letters[1:3] values <- data.frame(a = 1:3, b = 4:6, c = 7:9) With more complicated objects similar to 'names' and 'values' I wrote the following line to assign the elements of the

[R] assign a list using expression?

Dear R-users, I would like to assign elements to a list in the following manner: mylist <- list(a = a, b = b, c = c) To do this I tried myexpr <- expression(a = a, b = b, c = c) mylist <- list( eval(myexpr) ) It ends up by overwriting a when b is assigned and b when c is assigned. Additionally

Re: [R] How to get Greenhouse-Geisser epsilons from anova?

Dear John, thank you for the kind offer! Sorry, I just made a mistake anywhere I can not trace back, now it works as you described it. Thank you again! Dear Peter, thank you for the information, I did not know about the quotation marks. It indeed works using "G-G Pr"! The SPSS and R output for

Re: [R] How to get Greenhouse-Geisser epsilons from anova?

Dear John and Peter, thank you both very much for your help! Everything works fine now! John, Anova also works very fine. Thank you very much! However, if I had more than 2 levels for the between factor the same thing as mentioned occured. The degrees of freedom showed that Anova calculated it

Re: [R] How to get Greenhouse-Geisser epsilons from anova?

group but for example 24 from 2 groups, the output treats it as if all subjects came from the same group, for example for main effect A the dfs are 1 and 35. SPSS puts out 1 and 33 which is what I would have expected.. .. Peter Dalgaard schrieb: Nils Skotara wrote: Thank you, this helped me

Re: [R] How to get Greenhouse-Geisser epsilons from anova?

but whatever I tried afterwards did not seem logical to me. I am afraid I do not understand how to include the between factor. I cannot include ~D into M or X because it has length 24 whereas the other factors have 28... Zitat von Peter Dalgaard <[EMAIL PROTECTED]>: > Skotara wrote:

Re: [R] How to get Greenhouse-Geisser epsilons from anova?

4 I also don't know how I can calculate the various interactions.. My read is I should change the second argument mlmfit0, too, but I can't figure out how... Do you know what to do? Thank you very much! Peter Dalgaard schrieb: Skotara wrote: Dear all, I apologize for my basic q

[R] How to get Greenhouse-Geisser epsilons from anova?

Dear all, I apologize for my basic question. I try to calculate an anova for repeated measurements with 3 factors (A,B,C) having 2, 2, and 7 levels. or with an additional fourth between subjects factor D. Everything works fine using aov(val ~ A*B*C + Error(subject/ (A*B*C) ) ) or aov(val ~ (D