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
Show replies by date