高手求解

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

data e;
input year month$ n rem@@;
ln=log(n);
cards;
1983  jan  0    0  1983  feb  10   0  1983  mar  26   0
1983  apr  192  0  1983  may  238  0  1983  jun  374  1
1983  jul  356  0  1983  aug  358  4  1983  sep  594  0
1983  oct  786  1  1983  nov  885  2  1983  dec  981  2
1984  jan  1044  0  1984  feb  1266  1  1984  mar  1533  1
1984  apr  1510  2  1984  may  1818  2  1984  jun  2210  7
1984  jul  2003  2  1984  aug  2489  4  1984  sep  2841  2
1984  oct  3236  3  1984  nov  3104  5  1984  dec  2919  3
1985  jan  2849  3  1985  feb  3119  2  1985  mar  3596  7
1985  apr  4003  5  1985  may  4067  5  1985  jun  3690  2
1985  jul  3509  2  1985  aug  3653  9  1985  sep  4186  3
1985  oct  4562  2  1985  nov  4277  2  1985  dec  3838  2

1986  jan  4253  3  1986  feb  4242  5  1986  mar  5119  5
1986  apr  5281  7  1986  may  5163  4  1986  jun  4977  2
1986  jul  4663  4  1986  aug  5465  4  1986  sep  5314  5
1986  oct  5485  4  1986  nov  5688  6  1986  dec  5403  0
1987  jan  5749  2  1987  feb  5682  4  1987  mar  6460  3
1987  apr  6681  3  1987  may  7215  3  1987  jun  6650  8
1987  jul  7094  2  1987  aug  6600  6  1987  sep  7649  3
1987  oct  7615  9  1987  nov  6733  4  1987  dec  6540  10
1988  jan  6680  4  1988  feb  7646  6  1988  mar  8139  4
1988  apr  7829  4  1988  may  8220  3  1988  jun  7671  5
1988  jul  7120  3  1988  aug  7293  4  1988  sep  8045  5
1988  oct  8567  3  1988  nov  7682  6  1988  dec  7048  3
1989  jan  7369  2  1989  feb  7270  6  1989  mar  8124  1
1989  apr  7636  5  1989  may  7512  5  1989  jun  7049  4
1989  jul  7286  2  1989  aug  7624  2  1989  sep  7623  2
1989  oct  7970  5  1989  nov  7569  1  1989  dec  7156  10

1990  jan  7404  3  1990  feb  7447  8  1990  mar  7951  12
1990  apr  8056  7  1990  may  7742  3  1990  jun  7109  2
1990  jul  7229  4  1990  aug  7249  3  1990  sep  7366  0
1990  oct  7955  6  1990  nov  7044  6  1990  dec  3929  3
;
proc genmod data=e order=data;
class month year;
model rem=month year/dist=poi link=log offset=ln type1 type3;
run;
proc genmod data=e order=data;
class month year;
make 'obstats' out=d;
model rem=year/dist=poi link=log offset=ln type1 type3 obstats;
run;

data rates;
keep year rate month rem ln;
set d;
set e;
if month='dec';rate=exp(xbeta);
run;
proc print;
run;
data g; set d;
rr=pred/exp(ln);
proc print data=g;
run;


proc genmod data=e order=data;
class month year;
model rem/n=month year/dist=poi link=log  type1 type3;
run;

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：genmod Offset OBSTAT month model

llkknnllyytt 发表于 2015-1-15 19:14
data e;
input year month$ n rem@@;
ln=log(n);

我读书少，真心看不懂

support .......

这里面的内容是进行Poisson过程和对数Possible过程的拟合，然后再看看其预测率和按月计算的相关指标。

ermutuxia 发表于 2015-1-15 20:27
这里面的内容是进行Poisson过程和对数Possible过程的拟合，然后再看看其预测率和按月计算的相关指标。

你好，SAS帮助在哪里

ermutuxia 发表于 2015-1-15 20:27
这里面的内容是进行Poisson过程和对数Possible过程的拟合，然后再看看其预测率和按月计算的相关指标。

我想计算好几年的，不想要月的，不知道哪个可以去掉

ermutuxia 发表于 2015-1-15 20:27
这里面的内容是进行Poisson过程和对数Possible过程的拟合，然后再看看其预测率和按月计算的相关指标。

请问那个type1和type3是什么

llkknnllyytt 发表于 2015-1-15 23:16
你好，SAS帮助在哪里

我现在有几个相关的程序，不过想换成我自己的数据，能否给指点下，我的Q469691061 要开题了，着急呀，万分感谢

llkknnllyytt 发表于 2015-1-15 23:22
请问那个type1和type3是什么

你只要知道type1和type3的区别就可以了，
type1 与自变量进入方程的顺序有关，例如y= a  b c和y=b a c这两个在type1 的结果中a,b的p值是不同的，y=a b c中a因素是在没有校正b,c 情况下计算的统计量，y=a b c 中的b 是在只调整之前的a的情况下（没有校正c）得出的统计量。
而
type3与自变量的顺序无关，在y=a b c中，a的统计量是在校正了b c 后得出的，一般情况下看这个结果就行。

单独按年算的话
proc sql;
create table m as
select year,sum(n) as sumn ,sum(rem) as sumrem
from e
group by year;
quit;
data m;
set m;
ln=log(sumn);
proc genmod data=m;
class year;
make 'obstats' out=m2;
model sumrem=year/dist=poisson link=log offset=ln type3 obstats;
run;
这个语句只分析1983-1989这7条观测，你会看到结果里预测的准确性无偏，因为观测太少了。你要算年的话把12个月的估计值加起来就行了，还是以月为单位，以年为单位数据太少，意义不大