前两天用SAS/IML写了一个用EM算法估计generalized linear model中的参数的问题,
结果总是有问题,error在不停的变大,effect 的估计值也在不停的变大,后来调整了
初始值,有了一点点进展,但还是没完全解决问题。estimation 还是不理想。模型中
是有interaction的,但是调整初始值之后,只能估计出main effects,估计不出
interaction。后来把程序翻译成R语言,结果,这下无论怎么调初始值,程序运行到一
定地步后,总是有不能取倒数的值出现。而且前几轮算出的值也没有特别大的值出现。
彻底晕了,SAS, R which one should I trust?