Long time listener, first time caller...

I am still not completely comfortable that I have this specified correctly:

mix<-lmer(ERPamp~TV_avg_week * Category * QuestionnaireScore + (1|subject)
+(0+TV_avg_week|subject),ERPdata,REML=TRUE)

ERPamp: continuous, Gaussian
TV_avg_week: continuous, positive skew
Category: 4-level categorical predictor, same number of levels for every
subject, repeated within subject
QuestionnaireScore: continuous, Gaussian
subject: random, uncorrelated with other variables, every subject has a
value for all other variables

Category is well-replicated to predict ERPamp, so we are h(x):
Does Category interact with QuestionnaireScore in predicting ERPamp even
when TV_avg_week is controlled?

Looking forward to excellent advice, but please feel free to ask questions
about my less-than-excellent specification.

	[[alternative HTML version deleted]]


