[R] lme() direction

[Dieter Menne]

Mike Lawrence <mike <at> thatmike.com> writes:

Thanks for the excellent reproducible sample set!

> I'm most interested in the interaction between color and type, but I
> know that there is likely an effect of word. Yet since word is not
> completely crossed with type, simply adding it to an aov() won't work.
> A colleague recommended I look into lme() but so far I can't figure
> out the proper call.

Without word, it would be

summary(lme(rt~type*color, data=a,random=~1|id))

With the interaction, the extreme would be 
summary(lme(rt~type*color*word, data=a,random=~1|id))

or, less extreme

summary(lme(rt~type*color+color:word, data=a,random=~1|id))

but all these fail because of the rather degenerate structure 
of you data set. While lmer in package lme4 allows for a wider
set of solutions, I currently do not see how it could help,
but I might be wrong with p=0.5.


               word happy joy sad grumpy
type     color                          
positive white         93  90   0      0
         red           90  88   0      0
         green         88  87   0      0
negative white          0   0  88     95
         red            0   0  91     85
         green          0   0  88     88
>
 
> Another issue is whether to collapse across repetition before running
> the stats, particularly since errors will leave unequal numbers of
> observations per cell if it's left in.

That's one of the points where you have little to bother with the lme
approach. Collapsing would give equal weights to unequal numbers of
repeat, and might of minor importance when not too extreme, though.

Dieter


set.seed(1)
a=rbind(
	cbind(
		type='positive'
		,expand.grid(
			id=1:10
			,color=c('white','red','green')
			,word=c('happy','joy')
			,repetition = 1:10
		)
	)
	,cbind(
		type='negative'
		,expand.grid(
			id=1:10
			,color=c('white','red','green')
			,word=c('sad','grumpy')
			,repetition = 1:10
		)
	)
)

#add some fake rt data
a$rt=rnorm(length(a[,1]))

#And because people make errors sometimes:
a$error = rbinom(length(a[,1]),1,.1)

#remove error trials because they're not psychologically interesting:
a=a[a$error==0,]

library(nlme)
ftable(a[,c(1,3,4)])
summary(lme(rt~type*color, data=a,random=~1|id))