求高手指点写 ell,optim code

HongYu_柳

mu<-0
sigma2<-1
/ E! r; }) A" L: O4 c7 |  n<-20
  L<--0.25
  y<-rnorm(n,mu,sqrt(sigma2))
; G0 e) m( [6 E5 G0 x  x<-y
  c<-(y<=L)
x[c]<-L
$ ~" S5 c; [. u8 Q/ { x
0 g' a4 F1 c8 x5 \8 K+ s, IQ1:写一个ell function compute log-likelihood od x.
Q2用 R optim来maximize这个likelihood.

我现在写了一个ell的试子了，但是怎么也求不出最大值。
5 T$ Y0 [0 h  f, Mtheta<-c(mu,sigma2)
$ C& e+ T- S7 W8 t/ L* M ell<-function(theta,x,L=-0.25){
6 n5 M3 ?" h+ T* G7 E6 G. H/ N stopifnot(is.vector(theta),length(theta)==2)
if(theta[2]<=0)
return(NA)
. n$ b% I( a+ e! H/ U1 m else(c<-(x==L))
sum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
! G5 e. z. w2 |9 [" @ sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}
8 w" ?7 \* l7 Q1 I6 `& `1 L
+ y+ S* K% L8 ?7 x) @我的求最大值试子是init<-c(mu=0,sigma2=1)
; h2 V/ w/ p2 o
opt<-optim(par=init,fn=ell,method="L-BFGS-B",
2 P! R/ Y- ]7 g% b  lower=c(-Inf,0.01),
  upper=c(Inf,Inf),hessian=TRUE,
  control=list(fnscale=-1))
谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧