I am following the wiki page for LME analysis and I have a quick question. The Mass-univariate example near the bottom of the page proposes an initial model that contains intercept, linear and quadratic terms as random effects. However, the examples just below for lme_mass_fit_EM_init(), lme_mass_fit_EM_Rgw() only have [1 2] as selected random effects. Should the vector Zcols contain [1 2 3] as selected random effects in order to test the proposed model?
Thanks, Francois
François Lalonde, Ph.D. Child Psychiatry Branch NIMH / NIH 10 Center Drive, Room 3N202 Bethesda, MD 20892