search for: mlmfit1

...come ~ group * time + baseline + Error(subj+subj:time)) (from: http://cran.r-project.org/doc/contrib/Baron-rpsych.pdf, p. 41) C. "Repeated measures" MANOVA as follows (using data in wide format): response <- cbind(time1,time2,time3,time4,time5,time6) mlmfit <- lm(response ~ group) mlmfit1 <- lm(response ~ 1) mlmfit0 <- lm(response ~ 0) # Test time*group effect anova.mlm(mlmfit, mlmfit1, X=~1, test="Spherical") # Test overall group effect anova.mlm(mlmfit, mlmfit1, M=~1) # Test overall time effect anova.mlm(mlmfit1, mlmfit0, X=~1, test="Spherical") (taken f...

Hi, yet another anova.mlm problem - it doesn't seem to end. This time, I have a setup with a few within-subject factors and a between-subject factor (SGROUP). Consider the most simple case with only one within-factor (apo): > mlmfit0 <- lm(data.n ~ 0 + SGROUP) > mlmfit1 <- lm(data.n ~ 1 + SGROUP) > anova(mlmfit1,mlmfit0,test="Spherical",M=~hemi,X=~1) Analysis of Variance Table Model 1: data.n ~ 1 + SGROUP Model 2: data.n ~ 0 + SGROUP Contrasts orthogonal to ~1 Contrasts spanned by ~hemi Greenhouse-Geisser epsilon: 0.6829 Huynh-Feldt epsilon:...

Dear R folk, First off I want to thank those of you who responded with comments for my R quick and dirty stats tutorial. They've been quite helpful, and I'm in the process of revising them. When it comes to repeated measures MANOVA, I'm in a bit of a bind, however. I'm beginning to see that all of the documentation is written for psychologists, who have a slightly