方老师,
在高级版视频里面第一讲第一部分,关于极大似然估计有如下函数
loglik=function (para){
N=length(totalexp)
e=foodexp-para[1]-para[2]*totalexp
ll=-0.5*N*log(2*pi)-0.5*N*log(para[3]^2)-0.5*sum(e^2/para[3]^2)
return(ll)
}
res=maxLik(loglik,start=c(0.1,1,1))
summary(res)
其中para[3]是标准差,para[1]是beta1对应截距,para[2]是beta2对应自变量的系数,那么为什么在res=maxLik(loglik,start=c(0.1,1,1))的调用中将初始值设为c(0.1,1,1),原因何在?
望回复