急！SAS 的 IRR 怎么求？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

现在在做负二项回归，看别的大侠都用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)  
           
请教各位大侠
拜谢！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：IRR p-value genmod value 负二项回归 杂志

我想很多人没有研究负二项回归,还有一些人不知道IRR是什么意思(google除外,包括本人).所以回答你的人不会很多.
除非你能在具体下你的问题.

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衡量的是自变量每变动一个单位时，发生率的增长倍数！在考查其中一个自变量的效应时，必须固定其它的变量！

tj0412ymy 发表于 2011-8-10 11:04
proc genmod data = y;
  class x2 (param=ref ref=first);
  model y = x1 x2 / type3 dist=negbin;

非常非常非常非常感谢。 这个正是我需要的。
就是这个问题解决不了， 论文没法往下写。
我现在马上试试看

实在太感谢了

我的命令是这么写的：

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.
这个是怎么回事啊？

yugao1986 发表于 2011-8-10 11:02
我想很多人没有研究负二项回归,还有一些人不知道IRR是什么意思(google除外,包括本人).所以回答你的人不会很 ...

谢谢您的建议
我马上把我的问题具体化

yanmeilin 发表于 2011-8-10 14:07
我的命令是这么写的：

proc genmod data=a;

  estimate 'label1' intercept 1 exch 0 0 gdp 2 / exp;
  estimate 'label2' intercept 1 exch 0 0 gdp 3 / exp;

我觉得你对GDP值的设定有问题，2和3不一定在你所处理数据的合理范围内，故而造成hessian矩阵非正定！

tj0412ymy 发表于 2011-8-10 15:08
estimate 'label1' intercept 1 exch 0 0 gdp 2 / exp;
  estimate 'label2' intercept 1 exch 0 0 gd ...

我没有对gdp的值做过设定啊，直接从excel导入sas的
那么这个问题应该怎么解决呢？

不知道能否告知您的电话号码及合适的通话时间，我想电话请教您

yanmeilin 发表于 2011-8-10 15:33
我没有对gdp的值做过设定啊，直接从excel导入sas的
那么这个问题应该怎么解决呢？

请查看“消息”

IRR=incident rate ratio.
In negative binomial model, Y=log(P), so the estimate of IRR is the exponential of the estimated parameter.

BTW,  Sometimes you may run into convergence problems when you use NB model...