使用 Poisson 回归计算死亡率

现需要使用 Poisson regression 计算某疾病不同暴露水平的死亡率，并且调整其它因素，如性别、年龄和年代。有朋友熟悉不？以下是我的程序：

1. * prepare data for Poisson regression;
   
2. data poisson;
   
3.         set poisson;
   
4.         logpy=log(py); * the offset;
   
5. run;
   
6. 
   
7. * Poisson regression;
   
8. proc genmod data=poisson;
   
9.         class period age;
   
10.         model ncase=gender period age/dist=poisson link=log offset=logpy covb type3 wald ;
    
11.         output out=out xbeta=xb stdxbeta=std;
    
12.         by exposure;
    
13. run;
    
14. 
    
15. * calculate predicted rates;
    
16. data predrates;
    
17.         set out;
    
18.         obsrate=ncase/py;
    
19.         lograte=xb-logpy;
    
20.         prate=exp(lograte); * predicted rate;
    
21.         lcl=exp(lograte-1.96*std);
    
22.         ucl=exp(lograte+1.96*std);
    
23. run;

复制代码

以上程序貌似只能计算某接触水平各个调整因素不同水平的死亡率，没办法计算某接触水平（exposure）的“总体”死亡率，谢谢。

我再次阅读一下你的整个帖子，大概如此：A问题的答案是什么？我的B答案没有错啊。原来C的答案是如此啊。只是跟你抬个杠而已。京剧

"总体“死亡率是什么呢？INTERCEPT-ONLY MODEL或许可以给出“总体”的概念。/京剧

那怎么调整其他因素呢？

比如说，调整age， 你假设不同的年龄有不同的死亡率。所以总体的概念是什么呢？如果你需要平均age的死亡估测，lsmeans age或许可以满足要求。其实有很多不同的‘raw’死亡率的表示法，一个是你的obsrate，还有比如#death in total/total follow-up, etc.\jingju

jingju11 发表于 2012-7-16 20:37
比如说，调整age， 你假设不同的年龄有不同的死亡率。所以总体的概念是什么呢？如果你需要平均age的死亡估测 ...

我这里要计算的其实就是人年死亡率：No deaths/person-years；obsrate 的算法，也是基于此。

不是很明白楼主的“总体”含义？

YueweiLiu 发表于 2012-7-16 20:41
我这里要计算的其实就是人年死亡率：No deaths/person-years；obsrate 的算法，也是基于此。

呵呵。说白了，你并不是很清楚自己想要什么。为什么不问问你的老板呢？京剧

好吧。 可能是我误导你们了。 问题已经解决，所谓调整的死亡率，其实是通过计算不同暴露组相对于对照组的RR，和对照组的粗死亡率来计算的。
我可能受到其它一些算法的干扰，在考虑这个问题的时候，有些不对路，实际上很简单。

jingju11 发表于 2012-7-17 01:17
呵呵。说白了，你并不是很清楚自己想要什么。为什么不问问你的老板呢？京剧

呵呵，我很清楚自己想要做什么。