母鲎的例子，泊松对数回归中参数的标准差问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

各位高手，小弟现在在学习属性数据分析。书中有个关于母鲎追随者的例子，响应变量是追随者数量，解释变量是母鲎的宽度。对其建立泊松对数模型，请问结果中width前的回归系数的标准差0.0200是如何计算的？万分感谢！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：标准差 回归系数 解释变量 对数模型 数据分析 标准差 追随者 模型 如何

Maximum Likelihood --->Hessian matrix H(the final one)--->covariance matrix for parameters =- inverse of H--->std = square root of diagonal of COV matrix
Poisson has fixed scale =1, no std at all. Of note: sometimes forcing scale =1 is not realistic at all, due to overdispersion in the real data. seen here:
http://blog.sina.com.cn/s/blog_a3a92636010128im.html
you can look at the covariance matrix by:

1. proc genmod;
   
2. model y = x1 x2/covb;
   
3. run;

复制代码

jingju

jingju11 发表于 2012-4-15 23:56
Maximum Likelihood --->Hessian matrix H(the final one)--->covariance matrix for parameters =- invers ...

Thank you for your help!
I notice that the Hessian matrix is the expectation of the d2(L)/d(b1)d(b2), whether could I use the information matrix J , and covariance matrix cov(b)=J^-1 , where J * J^-1 = I ? Then I could just find the variance of the i th parameter bi from the main diagonal element.  

wh_wing 发表于 2012-4-16 00:17
Thank you for your help!
I notice that the Hessian matrix is the expectation of the d2(L)/d(b1)d ...

agree. same for using observed information matrix.

jingju11 发表于 2012-4-16 00:28
agree. same for using observed information matrix.

有一个关于芯片瑕疵数量的模型：
说A、B工艺分别生产出10个芯片，其瑕疵数量为：
A: 8, 7, 6, 6, 3, 4, 7, 2, 3, 4
B: 9, 9, 8, 14, 8, 13, 11, 5, 7, 6
假设瑕疵数服从独立的泊松分布muA和muB。
模型是log（mu）=a+bx，A对应x=0，B对应x=1。由题意，muA、muB的估计值分别为5和9。
通过log（5）=a ；log（9）=a+b 解得a=1.6094，b=0.5878
我用信息矩阵的方法：J=x'*w*x
x=[1 0; 1 1], w=[5 0; 0 9],
得到J=[14 9; 9 9]
J^-1=  [0.2000  -0.2000; -0.2000 0.3111]
那么var（a）=0.2 ， var（b）=0.311？ 可是与sas的数据不一样啊？

having trouble to understand what you have done. you might be correct anyway.

我查了一下公式，从asymptotic 理论来讲，是不是应该除以N啊？
se1 =sqrt(0.2/10), etc.
jingju

jingju11 发表于 2012-4-17 11:45
我查了一下公式，从asymptotic 理论来讲，是不是应该除以N啊？
se1 =sqrt(0.2/10), etc.
jingju

哎，好像是这样的！太感谢了！！！！！！！！！
请问你参考的是那本书呀？渐进性理论是不是主要用在MLE上面的？真心觉得统计没好好学。。。