求高手指点写 ell,optim code

HongYu_柳

mu<-0
sigma2<-1
/ e, {6 l* r8 `  n<-20
5 S) o  P4 m; z7 |! o+ p2 T- _  L<--0.25
/ P. `+ M* z/ H3 B9 t( m  y<-rnorm(n,mu,sqrt(sigma2))
  x<-y
  c<-(y<=L)
x[c]<-L
! ?8 O. C: z7 a x
; U% T- O- o2 l  w9 xQ1:写一个ell function compute log-likelihood od x.
Q2用 R optim来maximize这个likelihood.

我现在写了一个ell的试子了，但是怎么也求不出最大值。
# \# A& r  g3 ntheta<-c(mu,sigma2)
ell<-function(theta,x,L=-0.25){
stopifnot(is.vector(theta),length(theta)==2)
if(theta[2]<=0)
( R$ ^/ S- g( ?3 m. S: P6 j return(NA)
else(c<-(x==L))
sum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}

我的求最大值试子是init<-c(mu=0,sigma2=1)
" A9 m' q$ C  X6 M
) q" w* e2 ]0 Q3 r& ]opt<-optim(par=init,fn=ell,method="L-BFGS-B",
- z1 m( r  f3 p5 s$ e- p3 q* y  lower=c(-Inf,0.01),
  upper=c(Inf,Inf),hessian=TRUE,
  control=list(fnscale=-1))
谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧