代码: [b,dev,stats] =glmfit(X, Y, 'binomial', 'link', 'logit') ; % X自变量,Y应变量(1和0)
r=stats.resid;
......
L = -(T/2)*(1+log(2*pi)+log((r'*r)/T)); %求对数似然值:T为样本容量,pi=3.1415...,r为残差向量,k为方程参数个数;
aic = -2*(L)+2*(k+1);
L的计算结果同SPSS和Stata不同,请问上述的代码有没有问题?谢谢


雷达卡


京公网安备 11010802022788号







