ranic
我用Genmod过程来拟合广义线性模型,连接函数用log,分布用负二项分布,程序如下:
proc genmod data=sasuser.ship;
class style year circle;
a=_MEAN_;y=_RESP_;
r=0.42;
c=2*(y*log(a/y)+(y+r)*log((y+r)/(a+r)));
d=a*(1+a/r);
variance var=d;
deviance dev=c;
model times=style year circle/ link=log offset=ln;
run;
结果老是报错: ERROR: Error in computing deviance function.想知道偏差函数错在哪里啊?请教达人...谢谢了
