Гамма регрессия R GLM против SAS GENMOD
Я пытаюсь повторить гамма-регрессию, которую я имею в SAS с R. К сожалению, я получаю совершенно другую оценку параметров, и у меня нет никакой идеи, что может вызвать это. Я искал ответ на мою проблему на SO, но я не нашел ничего, что могло бы ее охватить.
Мой код 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"))
Что может вызвать различия в результатах этих двух кусков кода?
Заранее спасибо!