关于最大似然估计

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

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)
       }
      
       library(maxLik)
       res=maxLik(loglik,start=c(0.1,1,1))
       summary(res)

Error in qr.solve(H[!fixed, !fixed, drop = FALSE], G0[!fixed], tol = slot(control,  :
  singular matrix 'a' in solve

请问到底哪错了？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：最大似然估计 似然估计 最大似然 function singular function library control return matrix

没法迭代了，有奇异阵出现