求高手指点写 ell,optim code

HongYu_柳

mu<-0
  Y; [, T3 T8 v; x' ~  e sigma2<-1
; Q0 d9 [) d7 @7 V0 l3 X; a  n<-20
  L<--0.25
2 u4 G' x! ?( P; e  y<-rnorm(n,mu,sqrt(sigma2))
7 U3 Z+ l  O0 }: _  x<-y
  c<-(y<=L)
x[c]<-L
3 Z. E3 C$ i+ Q  ?+ w x
; I# [: a% ]0 \' TQ1:写一个ell function compute log-likelihood od x.
Q2用 R optim来maximize这个likelihood.

. g" X9 h% Y2 s9 ^) j我现在写了一个ell的试子了，但是怎么也求不出最大值。
theta<-c(mu,sigma2)
ell<-function(theta,x,L=-0.25){
stopifnot(is.vector(theta),length(theta)==2)
if(theta[2]<=0)
return(NA)
* u2 r3 v0 i7 }8 p3 Q' @' s5 O! a7 y% V9 U 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))}
" H3 o8 \0 }6 b# c0 `7 r
我的求最大值试子是init<-c(mu=0,sigma2=1)
1 `  q4 N( X" C& f% ?$ j
opt<-optim(par=init,fn=ell,method="L-BFGS-B",
4 {6 U5 Y: [- p6 R% M  lower=c(-Inf,0.01),
  upper=c(Inf,Inf),hessian=TRUE,
  control=list(fnscale=-1))
- t' ?% q) ]. Q. K谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧