r - 伽玛回归 R GLM 与 SAS GENMOD

Question

我正在尝试用 R 复制我在 SAS 中的伽马回归。不幸的是，我得到了完全不同的参数估计，我不知道是什么原因造成的。我一直在寻找关于我的问题的答案，但我没有找到任何涵盖它的东西。

我的 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"))

什么可能导致这两个代码的结果不同？

提前致谢！