poisson回归拟合度问题 - 经管之家

› 论坛 › 数据科学与人工智能 › 数据分析与数据科学 › SAS专版

poisson回归拟合度问题

7646

26

收藏 2015-05-13

      我做出了poisson回归模型到拟合度的检验这里，SAS的计算过程不给出P值，需要自己查表，请问我查卡方表的时候要查VALUE,还是查value/df呢？

                                          Criteria For Assessing Goodness Of Fit

                              Criterion                      DF          Value          Value/DF
                              Deviance                      2       5.9023       2.9511
                              Scaled Deviance          2       5.9023       2.9511
                              Pearson Chi-Square       2       6.4180       3.2090
                              Scaled Pearson X2       2       6.4180       3.2090
                              Log Likelihood                            838.6918

二维码

扫码加我拉你入群

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

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

全部回复

2015-5-13 13:56:47

（1）用Value计算卡方值；
（2）Value/DF是用来直观判断模型的拟合度，它的值越接近1越好。Deviance可以认为是拟合模型与饱和模型相比，Pearson 可以看成是观测与预测相比，在两者相近的情况下Value/DF的值才约等于1。
用SAS的函数可以求出卡方值：
data m;
x=1-probchi(5.9023,2);
put x;
run;
日志里得出结果：0.052279549

二维码

扫码加我拉你入群

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

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

2015-5-13 16:05:05

就喜欢马发表于 2015-5-13 13:56
（1）用Value计算卡方值；
（2）Value/DF是用来直观判断模型的拟合度，它的值越接近1越好。Deviance可以认 ...

谢谢，明白了，那么请问poisson回归模型中是不是还可以计算RR值

二维码

扫码加我拉你入群

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

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

2015-5-14 13:11:22

就喜欢马发表于 2015-5-13 13:56
（1）用Value计算卡方值；
（2）Value/DF是用来直观判断模型的拟合度，它的值越接近1越好。Deviance可以认 ...

请问你当时用poisson回归画图了吗

二维码

扫码加我拉你入群

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

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

2015-5-14 14:46:56

llkknnllyytt 发表于 2015-5-13 16:05
谢谢，明白了，那么请问poisson回归模型中是不是还可以计算RR值

（1）可以。比如
proc genmod data=m;
class age ; /*假设年龄已分4组*/
model y=age / offset=log dist=poisson link=log;
estimate 'age' age -1 1 0 0/exp;/*以第1组为参照，计算第2组相对第一组的RR*/
estimate 'age' age -1 0 1 0/exp;/*以第1组为参照，计算第3组相对第一组的RR*/
estimate 'age' age -1 1 0 1/exp;/*以第1组为参照，计算第4组相对第一组的RR*/
run;
（2）没画过poisson图，帮不了你

二维码

扫码加我拉你入群

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

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

2015-5-14 15:25:17

就喜欢马发表于 2015-5-14 14:46
（1）可以。比如
proc genmod data=m;
class age ; /*假设年龄已分4组*/

我想请问下这个程序对不对
data example1;
input hp path c n @@;
ln=log(n);
cards;
0 1 1 981
0 2 5 230
0 3 6 97
1 1 6 1010
1 2 20 604
1 3 21 287
;
proc genmod;
class hp path;
model c=hp path/dist=poisson
link=log offset=ln type1 type3;
run;

class那里是用hp 和path这两个因素吗，还是只写一个hp因素？这里hp代表的是年龄，path代表的是时间因素的影响

二维码

扫码加我拉你入群

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

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

点击查看更多内容…

2015-5-14 23:25:53

llkknnllyytt 发表于 2015-5-14 15:25
我想请问下这个程序对不对
data example1;
input hp path c n @@;

两个都写

二维码

扫码加我拉你入群

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

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

2015-5-16 22:12:08

就喜欢马发表于 2015-5-14 23:25
两个都写

你好，我根据你的程序算不出RR值，能不能就我的程序给编一下。还有请问你得出的是预测率吗？预测的率是哪个变量，是estimate吗？

二维码

扫码加我拉你入群

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

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

2015-5-16 23:04:39

就喜欢马发表于 2015-5-14 23:25
两个都写

最大似然参数估计值的分析
参数             自由度估计值    标准误差    Wald 95% 置信限       Wald          卡方 Pr > 卡方
Intercept          1 -3.6453    0.5990    -4.8193 -2.4714          37.04                   <.0001
hp          1       1 -0.6028    0.3513    -1.2914 0.0858          2.94                         0.0862
hp          2       1 -1.0021    0.3537    -1.6954 -0.3088          8.02                      0.0046
hp          3    1 -1.4095    0.3833    -2.1608 -0.6582          13.52                      0.0002
hp          4    1 -0.5991    0.4025    -1.3881 0.1898             2.22                         0.1367
hp          5       0 0.0000       0.0000    0.0000 0.0000 . .
path       1       1 3.8995       0.5428    2.8356 4.9633          51.61                      <.0001
path          2       1 3.0110       0.5236    1.9848 4.0371          33.07                      <.0001
path          3       1 2.6790       0.5220    1.6558 3.7022          26.34                         <.0001
path          4       1 2.2589       0.5858    1.1106 3.4071          14.87                         0.0001
path          5       1 0.0502       0.6175    -1.1600 1.2604          0.01                         0.9352
帮忙解释下这几个参数的意思呗，比如估计值是什么的估计值？还有最终公式里面的参数是不是还得运行别的命令才能出来？比如β？

二维码

扫码加我拉你入群

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

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

2015-5-17 14:54:09

llkknnllyytt 发表于 2015-5-16 23:04
最大似然参数估计值的分析
参数自由度估计值标准误差 Wald 95% 置信限 ...

data example1;
input hp path c n @@;
ln=log(n);
cards;
0 1 1 981
0 2 5 230
0 3 6 97
1 1 6 1010
1 2 20 604
1 3 21 287
;
proc genmod;
class hp(param=ref ref=first) path(param=ref ref=first);
model c=hp path/dist=poisson link=log offset=ln type1 type3;
output out=fit pred=fitmean l=lowerci u=upperci resraw=r;/*数据集中变量fitmean为预测数，在通过对数转换后除以n就可以得到率*/
estimate 'hp' hp -1 1 /exp;/*得到hp=1比hp=0增加的风险倍数*/
estimate 'path' path -1 1 0 /exp;/*得到path=2比path=1增加的风险倍数*/
estimate 'path' path -1 0 1/exp;/*得到path=3比path=1增加的风险倍数*/
run;
/*估计值是最后拟合的回归方程log(c)=log(n)+(-3.6453+(-0.6028)hp1+(-1.0021))hp2+...的每个自变量的系数，系数再经过转换得到OR值*/

二维码

扫码加我拉你入群

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

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

2015-5-17 22:59:21

就喜欢马发表于 2015-5-17 14:54
data example1;
input hp path c n @@;
ln=log(n);

(param=ref ref=first)这个话的意思是以第一组为参照吗

二维码

扫码加我拉你入群

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

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

2015-5-17 23:05:16

就喜欢马发表于 2015-5-17 14:54
data example1;
input hp path c n @@;
ln=log(n);

请问是ln(fitmean)/n为最终的预测率吗

二维码

扫码加我拉你入群

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

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

2015-5-17 23:11:31

就喜欢马发表于 2015-5-17 14:54
data example1;
input hp path c n @@;
ln=log(n);

还有亲，estimate语句运行不出来，显示错误语句无效或未按正确顺序使用

二维码

扫码加我拉你入群

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

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

2015-5-17 23:14:32

就喜欢马发表于 2015-5-17 14:54
data example1;
input hp path c n @@;
ln=log(n);

poisson回归是不是算RR值就行，还要算OR值吗

二维码

扫码加我拉你入群

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

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

2015-5-18 13:42:37

llkknnllyytt 发表于 2015-5-17 23:14
poisson回归是不是算RR值就行，还要算OR值吗

param=ref ref=first是设定第一组为参照组。
poisson回归可以根据总数c预测发病数n，预测率你有没有相关的文献，传上来看看。
我把程序复制到我的电脑可以运行，你可以把日志和程序传上来看看。
RR=暴露组的事件发生率/非暴露组事件的发生率，poisson中从c/n来看是RR，当c很小时，c/n近似等于c/(n-c)，即RR约等于OR。

二维码

扫码加我拉你入群

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

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

2015-5-18 14:04:48

就喜欢马发表于 2015-5-18 13:42
param=ref ref=first是设定第一组为参照组。
poisson回归可以根据总数c预测发病数n，预测率你有没有相关 ...

不知道怎么上传，我另外发了一个帖子传的率。我看到有人是预测率的。你这个意思是单纯从发病的总人数N来预测未来发病人数吗？那C就没用了呗。

二维码

扫码加我拉你入群

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

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

2015-5-18 14:05:34

就喜欢马发表于 2015-5-18 13:42
param=ref ref=first是设定第一组为参照组。
poisson回归可以根据总数c预测发病数n，预测率你有没有相关 ...

C/N应该等于发病率呀，RR值等于发病率？

二维码

扫码加我拉你入群

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

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

2015-5-18 14:09:35

就喜欢马发表于 2015-5-17 14:54
data example1;
input hp path c n @@;
ln=log(n);

poisson模型能得到预测的人数，模型拟合度好不好。还可以得出别的吗

二维码

扫码加我拉你入群

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

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

2015-5-18 14:12:58

就喜欢马发表于 2015-5-17 14:54
data example1;
input hp path c n @@;
ln=log(n);

Contrast Estimate Results

                                       Standard                                     Chi-
               Label       Estimate    Error    Alpha    Confidence Limits    Square Pr > ChiSq

               hp          -0.4870    0.3275    0.05    -1.1288    0.1549    2.21       0.1370
               Exp(hp)       0.6145    0.2012    0.05    0.3234    1.1675
               path          0.8426    0.2776    0.05    0.2985    1.3867    9.21       0.0024
               Exp(path)    2.3224    0.6447    0.05    1.3478    4.0017
               path       -2.0449    0.4318    0.05    -2.8913    -1.1985    22.42       <.0001
               Exp(path)    0.1294    0.0559    0.05    0.0555    0.3017

结果出来了，该利用哪个变量来解释呢？还有问下estimat的语句怎么使用，因为我有很多组，改怎么编这个程序呢

二维码

扫码加我拉你入群

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

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

2015-5-18 14:29:34

llkknnllyytt 发表于 2015-5-18 14:04
不知道怎么上传，我另外发了一个帖子传的率。我看到有人是预测率的。你这个意思是单纯从发病的总人数N来预 ...

c应该是没用的。
比如你的模型是Ln(c)=log(n)+age;年龄分为两组（50以下=0和50以上=1）。c/n是发病率的话，RR代表50以上的年龄组比50以下增加的患病危险。
你把很多组的那个变量和组数说一下。

二维码

扫码加我拉你入群

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

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

2015-5-18 17:42:10

就喜欢马发表于 2015-5-18 14:29
c应该是没用的。
比如你的模型是Ln(c)=log(n)+age;年龄分为两组（50以下=0和50以上=1）。c/n是发病率的话 ...

data example1;
input hp path c n @@;
ln=log(n);
cards;
1 1 23 158 1 2 41 655 1 3 31 1048 1 4 8 308 1 5 8 483 1 6 7 1003 1 7 1 605 1 8 1 158
2 1 14 78 2 2 32 338 2 3 19 462 2 4 4 85 2 5 1 157 2 6 2 262 2 7 7 438 2 8 2 178
3 1 7 23 3 2 13 86 3 3 15 143 3 4 0 35 3 5 0 204 3 6 0 40 3 7 0 47 3 8 1 67
4 1 4 4 4 2 7 14 4 3 4 30 4 4 5 16 4 5 0 239 4 6 1 44 4 7 1 15 4 8 1 30
5 1 1 3 5 2 2 3 5 3 4 8 5 4 1 4 5 5 1 52 5 6 0 68 5 7 1 7 5 8 2 18
;

proc genmod data=example1 order=data;
class hp path;
model c=hp path/dist=poisson
link=log offset=ln type1 type3;
run;
proc genmod data=example1 order=data;
class hp path;
make 'obstats' out=d1;
model c=path/dist=poi link=log offset=ln type1 type3 obstats;
run;
set example1;

if hp=1;
rate=exp(xbeta);
run;
proc print;
run;
data g1;
set d1;

rr=pred/exp(ln);
proc print data=g1;
run;

这是老师给的程序，我加的自己的数，但是总觉得有点运行不出来。我QQ：469691061 方便的话QQ聊下，万分感谢

二维码

扫码加我拉你入群

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

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

2015-5-19 11:37:49

就喜欢马发表于 2015-5-18 14:29
c应该是没用的。
比如你的模型是Ln(c)=log(n)+age;年龄分为两组（50以下=0和50以上=1）。c/n是发病率的话 ...

还有问下RR值的可信区间能求吗

二维码

扫码加我拉你入群

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

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

2015-5-20 17:12:46

llkknnllyytt 发表于 2015-5-19 11:37
还有问下RR值的可信区间能求吗

proc genmod data=example1 order=data;
class hp(param=ref ref=first) path(param=ref ref=first);
model c=hp path/dist=poisson
link=log offset=ln type1 type3;
/*hp2、3、4、5与1的比较*/
estimate 'hp_rr21' hp -1 1 0 0 0/exp;
estimate 'hp_rr31' hp -1 0 1 0 0/exp;
estimate 'hp_rr41' hp -1 0 0 1 0/exp;
estimate 'hp_rr51' hp -1 0 0 0 1/exp;
/*path2、3、4、5、6、7、8与1的比较*/
estimate 'path_rr21' path -1 1 0 0 0 0 0 0/exp;/*path_rr21名字代表2与1对比*/
estimate 'path_rr31' path -1 0 1 0 0 0 0 0/exp;
estimate 'path_rr41' path -1 0 0 1 0 0 0 0/exp;
estimate 'path_rr51' path -1 0 0 0 1 0 0 0/exp;
estimate 'path_rr61' path -1 0 0 0 0 1 0 0/exp;
estimate 'path_rr71' path -1 0 0 0 0 0 1 0/exp;
estimate 'path_rr81' path -1 0 0 0 0 0 0 1/exp;
run;

以hp=1为对照，hp=2的RR=1.5491，置信区间为

1.0445

2.2975

二维码

扫码加我拉你入群

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

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

2015-5-20 18:00:20

就喜欢马发表于 2015-5-20 17:12
proc genmod data=example1 order=data;
class hp(param=ref ref=first) path(param=ref ref=first);
...

我想问下β值在哪里显示

二维码

扫码加我拉你入群

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

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

2015-5-20 18:05:45

就喜欢马发表于 2015-5-20 17:12
proc genmod data=example1 order=data;
class hp(param=ref ref=first) path(param=ref ref=first);
...

均值估计值就指的是RR呗

二维码

扫码加我拉你入群

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

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

2015-5-20 18:17:38

就喜欢马发表于 2015-5-20 17:12
proc genmod data=example1 order=data;
class hp(param=ref ref=first) path(param=ref ref=first);
...

我现在很迷茫，有本书上说poisson用estimate语句算OR值，和95%可信区间。迷茫了

附件列表

原图尺寸 292.13 KB

二维码

扫码加我拉你入群

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

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

2017-2-25 20:14:35

就喜欢马发表于 2015-5-13 13:56
（1）用Value计算卡方值；
（2）Value/DF是用来直观判断模型的拟合度，它的值越接近1越好。Deviance可以认 ...

得出的值0.052279要怎么看呢？

二维码

扫码加我拉你入群

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

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

相关推荐

栏目导航

热门文章

推荐文章

扫码加好友，拉您进群

各岗位、行业、专业交流群