Thank you @stijn for your constructive comments. The command is working now. Let me show you some values of Y
and mu[1]
. Do you think that the estimated Poisson rate parameter is health (seems underestimated?)?
y mu[1]
1 104 9.17
2 83 10.77
3 78 9.07
4 60 7.39
5 45 4.82
6 45 8.69
7 41 4.15
8 2 0.70
9 2 0.73
10 2 0.41
11 2 0.84
12 2 0.62
13 2 0.96
14 2 0.73
15 2 0.74
16 2 1.22
Therefore, would you please give me some suggestions over it.