求高手指点写 ell,optim code

HongYu_柳

mu<-0
  f5 b2 b1 B$ N" G sigma2<-1
3 h# @8 v& U5 H, B+ ?  N  n<-20
  L<--0.25
  y<-rnorm(n,mu,sqrt(sigma2))
  x<-y
  c<-(y<=L)
& I: k1 K6 p7 G x[c]<-L
3 a6 I% x, F# L7 H- U2 r" A x
Q1:写一个ell function compute log-likelihood od x.
Q2用 R optim来maximize这个likelihood.

& q5 W& j5 c4 R3 \2 v7 l4 W/ M0 |我现在写了一个ell的试子了，但是怎么也求不出最大值。
theta<-c(mu,sigma2)
; j4 Q, o) Z, X/ D. G ell<-function(theta,x,L=-0.25){
stopifnot(is.vector(theta),length(theta)==2)
4 U! T& o7 O4 @9 _8 G) k8 I if(theta[2]<=0)
* X/ M. R  f4 G3 l! m' C3 h# _ return(NA)
0 `: K0 Z4 r( p1 o' e: G/ [0 f% z else(c<-(x==L))
sum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
3 m$ E3 G, B  w* z( r; y: E sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}
& g0 u6 r. A  G" r
我的求最大值试子是init<-c(mu=0,sigma2=1)

7 K( B8 o+ m3 x2 ]0 Q. ~4 K  c. ~; }opt<-optim(par=init,fn=ell,method="L-BFGS-B",
  lower=c(-Inf,0.01),
  O- _/ ^% t- _& @; j. J  upper=c(Inf,Inf),hessian=TRUE,
  control=list(fnscale=-1))
谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧