求高手指点写 ell,optim code

HongYu_柳

mu<-0
2 k8 p* |4 h  R& [ sigma2<-1
  d2 m$ F6 J; O4 ~  n<-20
  L<--0.25
  y<-rnorm(n,mu,sqrt(sigma2))
  x<-y
  c<-(y<=L)
x[c]<-L
x
8 D/ R. K3 o; D% q" }Q1:写一个ell function compute log-likelihood od x.
) Z5 Y6 V, H' h; a: T$ AQ2用 R optim来maximize这个likelihood.
. L0 ?0 p, v8 u8 P" r
我现在写了一个ell的试子了，但是怎么也求不出最大值。
/ o6 U- K2 J% m) Jtheta<-c(mu,sigma2)
* B5 r0 v& L* I& V% g4 T ell<-function(theta,x,L=-0.25){
stopifnot(is.vector(theta),length(theta)==2)
if(theta[2]<=0)
return(NA)
& J" v, {" P+ q! w* b' L) M else(c<-(x==L))
sum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
- }- S. r% A/ t sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}

我的求最大值试子是init<-c(mu=0,sigma2=1)

opt<-optim(par=init,fn=ell,method="L-BFGS-B",
6 }- w0 ]3 }& z  h6 M9 C+ L8 f5 A  lower=c(-Inf,0.01),
  upper=c(Inf,Inf),hessian=TRUE,
  control=list(fnscale=-1))
谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧