求助

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

请问哪位高手可以帮我解决一下：用R程序求最大似然估计的问题（正太分布的），多谢了！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：最大似然估计 似然估计 最大似然 R程序 求助

mlFun <- function(par) {
  xbeta <- X %*% par[1:K]
  Sig <- par[(K+1)]
  ll <- -(1/2) * log(2*pi) - (1/2) * log(Sig^2) - (1/(2*Sig^2))*(Y-xbeta)^2
  return(sum(ll))
}
n <- 100
K <- 1                                  #解释变量个数
paras <- c(1,1)                         #初值（斜率系数，扰动项方差)
X <- matrix(rnorm(n),ncol=K)
Y <- matrix(2*X+rnorm(n),ncol=1)
mlfit <- optim(par=paras,fn=mlFun)
mlfit$par
[1] 7.906692e-01 1.884238e-52