求高手指点写 ell,optim code

HongYu_柳

mu<-0
! c1 ]$ a5 `* Q1 u( s2 @: T/ ~ sigma2<-1
& k/ T% j$ C7 n9 U6 m/ N  n<-20
6 u* Z5 z& w4 Q- c+ e3 X$ l6 I  L<--0.25
7 }6 O/ x% {# j- Y) J6 q5 s  y<-rnorm(n,mu,sqrt(sigma2))
  x<-y
; d+ j% r$ q3 j6 \) I" F3 K7 q  c<-(y<=L)
x[c]<-L
x
  G: ]9 x! c+ S! g$ `* ZQ1:写一个ell function compute log-likelihood od x.
0 C) `6 A/ s1 g4 I6 z( e( vQ2用 R optim来maximize这个likelihood.

我现在写了一个ell的试子了，但是怎么也求不出最大值。
theta<-c(mu,sigma2)
ell<-function(theta,x,L=-0.25){
stopifnot(is.vector(theta),length(theta)==2)
! e; u& M2 V4 X, P& U2 s if(theta[2]<=0)
return(NA)
4 y" G7 q% l: K9 ~( y& y 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)
( ^6 n8 b, {! X, ~) ], d8 U
5 E2 s" P+ i5 i8 m* l3 [6 }opt<-optim(par=init,fn=ell,method="L-BFGS-B",
9 q$ \8 N4 B) B5 F6 h' a* [  lower=c(-Inf,0.01),
, \3 q: c" b8 x" ?  upper=c(Inf,Inf),hessian=TRUE,
  control=list(fnscale=-1))
谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧