各位高手,如何在R软件中求零膨胀poisson模型的参数呢?我R软件不懂,请帮帮我,谢谢!
如:在d:\data.xls有如下数据:
| x1 | x2 | x3 | x4 | x5 | x6 | y |
| 23.6 | 78 | 37 | 5 | 1.4 | 64.5 | 1 |
| 17.6 | 85 | 56 | 6 | 1 | 49.3 | 0 |
| 30.5 | 84 | 24 | 6 | 1 | 91.1 | 3 |
| 29.8 | 85 | 39 | 4 | 0.6 | 97.3 | 0 |
| 30.4 | 77 | 22 | 6 | 1.3 | 109.6 | 1 |
| 32.4 | 80 | 28 | 6 | 1.1 | 97.9 | 5 |
| 235 | 81 | 33 | 5 | 0.9 | 50.1 | 0 |
| 27.6 | 81 | 24 | 6 | 1.1 | 88.7 | 2 |
| 33 | 75 | 24 | 5 | 1.6 | 165.7 | 7 |
log(p/(1-p))=a1+a2*x1+a3*x2+a4*x3+a5*x4+z6*x5
λ=exp(b1+b2*x1+b3*x2+b4*x3+b5*x4)
请高手帮个忙,写个代码帮我求出模型的参数,并算出误差!谢谢了!