R:计算最大似然估计量

R:计算最大似然估计量,r,mathematical-optimization,inference,R,Mathematical Optimization,Inference,我想计算方差=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) se

我想计算方差=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)
set.seed(36)

你拼错了两个参数。试着写
loglike
而不是
loglik
x
而不是
x
,它会工作的

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

$objective
[1] -157.6814

你拼错了两个论点。试着写
loglike
而不是
loglik
x
而不是
x
,它会工作的

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

$objective
[1] -157.6814

R区分大小写,并坚持精确匹配:try
optimize(f=loglikelibility,x=x,interval=c(0100),max=T)
R区分大小写,并坚持精确匹配:try
optimize(f=loglikelibility,x=x,interval=c(0100),max=T)