小编use*_*254的帖子

我正在尝试优化正态分布的对数似然.对数似然的函数有效,它识别传入的数据集,但优化无法识别数据集是否存在？如果我们设置data_x=rnorm(100,0,1)此代码返回正确答案,但我需要能够传入不同的数据集.

x <- rnorm ( 100, 0, 1)
loglike <- function( pars, data_x=x) {
  mu <- pars[1]
  sigma2 <- pars[2]^2
  numobs <- length( data_x )
  sumsq <- sum( ( data_x-mu )^2 )
  val.log.like <- -numobs / 2 * log( sigma2 ) - ( 1 / (2*sigma2) ) * sumsq
  return( val.log.like )
}

optimization <- optim( c( 0, 1), loglike)
answer <- matrix( optimization$par, 2, 1) 
answer