求解 不明白

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

程序1：

data dd;

input f e n d@@;

lgn=log(n);

cards;

1 1 240 216  1  0  100  63  2  1  60  48  2  0  210  84

;

proc genmod data=dd  order=data;

class f;

model d=e f/dist=poisson

link=log   type1 type3;

run;

程序2：

data a;

input hp path c n@@;

ln=log(n);

cards;

1 1 216 240 0 1 36 100 1 2 48 60 0 2 84 210

;

proc genmod;

class path;

model c=hp path/dist=poisson link=log offset=ln type1 type3;

run;

这两个程序说的是不是一回事 怎么结果不同呢

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：不明白 poisson genmod Offset 是不是一回事 程序

看你这几天都在弄这个poisson分布 真是辛苦

sushe1527 发表于 2015-2-11 12:58
看你这几天都在弄这个poisson分布 真是辛苦

是呀，整的我都混乱了，不知道哪个是对的程序了，都能跑出来结果

1、你的dd数据中63 应为36，要与a数据集要相同；
2、
proc genmod data=dd  order=data;
class f;
model d=e f/dist=poisson
link=log   type1 type3;
run;
这一poisson回归没有设offset位移变量，就是你已经假设所有的抽样单元大小相同，这种情况下就不用在model语句中设offset关键字。反之，如果抽样单元大小不同，需要加入offset位移变量。这里的抽样单元可以理解为事件发生数的分母，比较值分母不同的情况下需要加入offset位移变量
3、你的第一个程序假设分母相同，而第二个假设分母不同并进行了校正，所以结果不同。该选择那一个，请问你的分母是什么？

天穹下 发表于 2015-2-12 20:31
1、你的dd数据中63 应为36，要与a数据集要相同；
2、
proc genmod data=dd  order=data;

data dd;
input f e n d@@;
lgn=log(n);
cards;
1 1 240 216  1  0  100  63  2  1  60  48  2  0  210  84

;
proc genmod data=dd  order=data;
class f;
model d=e f/dist=poisson
link=log   type1 type3;

run;

proc genmod data=dd;
class f;
make 'obstats' out=outdata;
model d=f/dist=poisson link=log offset=lgn
type1 type3  obstats;
run;

proc print data=outdata;
run;

data k; set outdata;
set dd;
rate=exp(xbeta);run;proc print;
run;
这是完整 的程序  这个程序有问题吗

天穹下 发表于 2015-2-12 20:31
1、你的dd数据中63 应为36，要与a数据集要相同；
2、
proc genmod data=dd  order=data;

data dd;
input f e n d@@;
lgn=log(n);
cards;
1 1 240 216  1  0  100  63  2  1  60  48  2  0  210  84

;
proc genmod data=dd  order=data;
class f;
model d=e f/dist=poisson
link=log   type1 type3;

run;

proc genmod data=dd;
class f;
make 'obstats' out=outdata;
model d=f/dist=poisson link=log offset=lgn
type1 type3  obstats;
run;

proc print data=outdata;
run;

data k; set outdata;
set dd;
rate=exp(xbeta);run;proc print;
run;
这是完整 的程序  这个程序有问题吗

llkknnllyytt 发表于 2015-2-14 16:34
data dd;
input f e n d@@;
lgn=log(n);

data k; set outdata;
set dd;
rate=exp(xbeta);run;proc print;
run;
其他的没错，最后这部分有问题，应改为
proc sql;
create table k as
select *, exp(xbeta) as rate
from dd left join outdata on  dd.d=outdata.d;