ziyenano 发表于 2013-1-17 20:03 
能不能把你写的程序贴上来~
下面是我本来的程序，但是只能给出表格，不能给出直观的图形。我想加点什么，做出像最上面提问中的那样直观的累计残差检验图，但是完全不会编。。还请您指点指点啊，我都没找到类似的程序小样儿。。
data accident;
input tar vot x2 x5 x6 x8;
ln=log(vot);
datalines;
0.18 4.11 0.21 -1.25 0.856 0.81
1.00 4.11 0.32 -2.70 1.325 0.81
1.18 4.11 0.37 -2.34 1.291 0.81
2.72 4.11 0.96 -3.34 1.355 0.81
5.28 4.11 0.53 -3.85 1.538 0.81
1.64 4.11 0.46 -2.48 2.277 0.81
0.54 4.11 0.66 -1.70 2.838 0.81
3.36 4.11 0.36 -4.02 4.026 0.81
5.54 4.11 0.54 -3.42 5.026 0.81
8.36 4.11 0.91 -3.16 5.054 0.81
12.34 4.11 0.06 -3.95 6.695 0.81
28.76 4.11 0.94 -3.31 5.614 0.81
45.28 4.11 0.54 -2.53 7.944 0.81
3.36 4.20 0.21 -1.25 0.856 0.74
4.46 4.20 0.32 -2.70 1.325 0.74
3.08 4.20 0.37 -2.34 1.291 0.74
8.00 4.20 0.96 -3.34 1.355 0.74
7.10 4.20 0.53 -3.85 1.538 0.74
2.82 4.20 0.46 -2.48 2.277 0.74
2.36 4.20 0.66 -1.70 2.838 0.74
5.08 4.20 0.36 -4.02 4.026 0.74
13.06 4.20 0.54 -3.42 5.026 0.74
28.08 4.20 0.91 -3.16 5.054 0.74
37.32 4.20 0.06 -3.95 6.695 0.74
55.20 4.20 0.94 -3.31 5.614 0.74
39.62 4.20 0.54 -2.53 7.944 0.74
3.82 4.24 0.21 -1.25 0.856 0.73
3.54 4.24 0.32 -2.70 1.325 0.73
7.08 4.24 0.37 -2.34 1.291 0.73
13.26 4.24 0.96 -3.34 1.355 0.73
16.48 4.24 0.53 -3.85 1.538 0.73
4.28 4.24 0.46 -2.48 2.277 0.73
1.90 4.24 0.66 -1.70 2.838 0.73
6.72 4.24 0.36 -4.02 4.026 0.73
10.36 4.24 0.54 -3.42 5.026 0.73
13.26 4.24 0.91 -3.16 5.054 0.73
19.18 4.24 0.06 -3.95 6.695 0.73
20.18 4.24 0.94 -3.31 5.614 0.73
17.76 4.24 0.54 -2.53 7.944 0.73
;
proc genmod data=accident;
model tar=x2 x5 x6 x8/dist=negbin
offset=ln;
proc genmod data=accident;
k = .14;
a = _MEAN_;
y = _RESP_;
variance var = a+k*a*a;
if (y>0) then d = 2 * (y*log(y/a)-(1+k*y)/k *
log((1+k*y)/(1+k*a)));
else if (y=0) then d = 2 * log(1+k*a)/k;
deviance dev = d;
model tar=x2 x5 x6 x8 / offset = ln noscale link = log
type1
type3;
run;
扫码加好友，拉您进群
全部版块
我的主页
收藏
