我正在尝试复制我在SAS中使用R的伽马回归。不幸的是,我得到了完全不同的参数估计,我不知道是什么原因导致的。我一直在寻找我的问题的答案,但我没有找到任何覆盖它的东西。
我的SAS代码是:
proc genmod data=data;
class var1 factor1 ;
Make 'Obstats' Out=Outdata_G;
model answer = var2 factor1 var3/
offset= offset dist=Gamma link=log obstats;
run;
我看到var1在类声明中使用,但未在模型语句中使用 - 它能改变结果吗? (我不能重新运行代码)。
我在R中所做的是:
glm(answer~var2+factor1+var3, data =(data), offset = offset, family = Gamma(link = "log"))
什么可能导致这两个代码的结果出现差异?
提前致谢!