全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 SAS专版
6942 10
2011-08-10
现在在做负二项回归,看别的大侠都用IRR来解释结果, 我的输出里怎么看不到有IRR啊?
需要什么特别的命令吗?
我在sas里是这么写我的命令的:

proc genmod data=a;
model ad=gdp exch impt uept rtr / dist=negbin;
run;

仅仅得到系数的测定和对应的p-value。 不知道学术杂志上的 IRR 是怎么求出如下这个表格来的,
              (1)          (2)          (3)
GDP     1.00
            (0.268)

EXCH    1.0303    1.0296   1.0336
           (2.749)    (2.486)    (2.897)
IMPT                     1.6015
                             (1.981)  
           
请教各位大侠
拜谢!
二维码

扫码加我 拉你入群

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

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

全部回复
2011-8-10 11:02:08
我想很多人没有研究负二项回归,还有一些人不知道IRR是什么意思(google除外,包括本人).所以回答你的人不会很多.
除非你能在具体下你的问题.
二维码

扫码加我 拉你入群

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

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

2011-8-10 11:04:41
proc genmod data = y;
  class x2 (param=ref ref=first);
  model y = x1 x2 / type3 dist=negbin;
  estimate 'label1' intercept 1 x2 0 0 x1 2 / exp;
  estimate 'label2' intercept 1 x2 0 0 x1 3 / exp;
run;
x2是一分类变量,x1是一连续变量,y是因变量;在固定x2的前提下,计算x1的一阶IRR。
对计算得到的Mean Estimate做对数运算即可得到IRR:
IRR=log(u2/u1)。若是计算n阶IRR,则为IRR=logn(u2/u1)。

注:IRR衡量的是自变量每变动一个单位时,发生率的增长倍数!在考查其中一个自变量的效应时,必须固定其它的变量!
二维码

扫码加我 拉你入群

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

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

2011-8-10 13:34:32
tj0412ymy 发表于 2011-8-10 11:04
proc genmod data = y;
  class x2 (param=ref ref=first);
  model y = x1 x2 / type3 dist=negbin;
非常非常非常非常感谢。 这个正是我需要的。
就是这个问题解决不了, 论文没法往下写。
我现在马上试试看

实在太感谢了
二维码

扫码加我 拉你入群

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

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

2011-8-10 14:07:17
我的命令是这么写的:

proc genmod data=a;
class exch;
model ad = gdp exch / type3 dist=negbin;
  estimate 'label1' intercept 1 exch 0 0 gdp 2 / exp;
  estimate 'label2' intercept 1 exch 0 0 gdp 3 / exp;
run;

我的观察值是336个, 输出显示 df 291
出来很多 estimate (每个exch 对应的值),不知道该怎么对您所说的 “Mean Estimate“做对数运算。

还有一个警告如下
Warning: negative of hessian not positive definite.
这个是怎么回事啊?
二维码

扫码加我 拉你入群

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

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

2011-8-10 14:08:22
yugao1986 发表于 2011-8-10 11:02
我想很多人没有研究负二项回归,还有一些人不知道IRR是什么意思(google除外,包括本人).所以回答你的人不会很 ...
谢谢您的建议
我马上把我的问题具体化
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群