[请教达人]:关于偏差函数的问题

shiyiming

我用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.想知道偏差函数错在哪里啊?请教达人...谢谢了