求高手指点写 ell,optim code

HongYu_柳

mu<-0
sigma2<-1
9 W; d# e5 X, x8 Y  n<-20
  L<--0.25
  y<-rnorm(n,mu,sqrt(sigma2))
# }! S. w& B4 E1 }  x<-y
  c<-(y<=L)
; P  I+ O4 O# j4 [- {* n# |# m7 F x[c]<-L
: J2 m1 B/ h3 n0 ?7 e" ] x
Q1:写一个ell function compute log-likelihood od x.
  v! |+ G2 w9 NQ2用 R optim来maximize这个likelihood.

我现在写了一个ell的试子了，但是怎么也求不出最大值。
theta<-c(mu,sigma2)
ell<-function(theta,x,L=-0.25){
7 _# B/ e" J+ }+ q! f7 f$ B, l stopifnot(is.vector(theta),length(theta)==2)
) o2 S6 e2 o( T. o7 R if(theta[2]<=0)
4 A5 u8 f2 M  \8 h return(NA)
else(c<-(x==L))
% S+ ^9 O( v5 V+ M0 \, T/ `sum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}

  r' `3 h8 _" l/ B& ~; L. N7 r我的求最大值试子是init<-c(mu=0,sigma2=1)
6 Y% n; h5 R- e
opt<-optim(par=init,fn=ell,method="L-BFGS-B",
  lower=c(-Inf,0.01),
6 P5 @/ G8 _3 k$ _) u# {$ [: N  upper=c(Inf,Inf),hessian=TRUE,
, K  ?; k, R: x& f  control=list(fnscale=-1))
谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧