Displaying 1 result from an estimated 1 matches for "zi0".
Did you mean:
i0
2008 Sep 07
1
an error to call 'gee' function in R
...s problem? I appreciate your help.
In this project I will use the 'gee' or 'geeglm' and 'glmer' to fit the simulated multivariate count responses. I generated the data like this:
Set ?0 = ?1 = 1, ?0 = 3, and n = 50.
For each 1 ? i ? n,
Simulate xi from N (1, 1).
Simulate zi0 and zit from
zi0 follows i.d. Poisson (?0) ,
zit | xi follows i.d. Poisson (?it) , 1 ? t ? 3,
log (?it) = log(E (zit | xi)) = ?0t + xi?1t = 1+xi.
Let yit = zi0 + zit, 1 ? t ? 3.
So my data frame, let me call it 'simdata', the first 10 rows look like this:
id y.1 y.2 y.3 x
1...