R：计算最大似然估计答案

【问题标题】：R : Computing maximum likelihood estimatorR：计算最大似然估计
【发布时间】：2023-03-07 03:54:01
【问题描述】：

我想计算方差 = 36 的正态分布的均值 mu 的最大似然估计量。

但以下程序显示错误：

set.seed(36)

x <- rnorm(50,mean=20,sd=6)

loglikelikelihood<- function(mu,x){
 L = (-50)*log(sqrt(2*pi)*6)-((1/72)*(sum(x^2)-2*mu*sum(x)+50*mu^2))
 }

optimize(f=loglik,X=x,interval=c(0,100),maximum=T)

【问题讨论】：

R 区分大小写，并坚持完全匹配：试试optimize(f=loglikelihood,x=x,interval=c(0,100),maximum=T)

标签： r mathematical-optimization inference

【解决方案1】：

你拼错了两个参数。尝试写loglikelihood 代替loglik 和x 代替X，它会工作。

> optimize(f=loglikelihood,x=x,interval=c(0,100),maximum=T)
$maximum
[1] 20.86679

$objective
[1] -157.6814

【讨论】：