代码: [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不同，请问上述的代码有没有问题？谢谢