Hi,
 
I specified the following model in zelig() :
 
z.out <- zelig(qasbat ~ group + t + t2 + t3 + tag(1 + t + t2 | ID), data=matched.1.mtch.long,   model="ls.mixed")
 
But when it comes to simulate the quantities of interest with sim(), I get the following message :
 
Error in .local(object, ...) : Code not written yet

Theres is no error message from sim() when I delete t and t2 from the random part in tag(). But I'd like to let the first two polynomial time (t and t2) to be random.
 
Am I suppose to know how to use lmer() from lme4 package ? Or is it a problem with the way I write the zelig() formula ?
 
To help analyse my problem I write down my model in two forms :
 
LEVEL FORM :
 

Level 1 :                           qasbatij = π0i + π1itimeij + π2itime2ij + π3itime3ij + εij

 

Level 2 :                           π0i = γ00 + γ01groupi + ζ0i

                                        π 1i = γ10 + ζ1i

                                        π 2i = γ20 + ζ2i

                                        π 3i = γ30

 
 
ONE FORMULA FORM : Y = [Fixe part] + [Random part]
 
qasbatij = [γ00 + γ01groupi + γ10timeij + γ20time2ij + γ30time3ij] + [ζ0i +  ζ1itimeij + ζ2itime2ij + εij]
 
 
Merci beaucoup,

François Maurice, B. Sc., A. Stat.
Candidat à la maîtrise
Département de sociologie
Université de Montréal