求高手指点写 ell,optim code

HongYu_柳

mu<-0
sigma2<-1
  n<-20
  L<--0.25
! _4 M4 f3 {4 {) D  y<-rnorm(n,mu,sqrt(sigma2))
/ Y" k& F5 y# ?% m$ g6 I% E. }  x<-y
  c<-(y<=L)
x[c]<-L
8 Z5 p5 |$ v1 B  V0 y  z1 A x
" Y( Z7 T& p" a! jQ1:写一个ell function compute log-likelihood od x.
: f# `( g" z" J0 Y2 ]0 Y- YQ2用 R optim来maximize这个likelihood.
: g$ q$ }; p/ A, F$ J
3 a/ v; b2 {  M& ~; {我现在写了一个ell的试子了，但是怎么也求不出最大值。
7 t4 p1 E; e+ n  X' {theta<-c(mu,sigma2)
ell<-function(theta,x,L=-0.25){
0 F7 M' a! o% f& g. r  q* V+ e2 r stopifnot(is.vector(theta),length(theta)==2)
3 K' ]" e  h( Y( r0 M if(theta[2]<=0)
return(NA)
* f6 d; s% ~2 Y$ u  P( l3 c' j6 a8 G else(c<-(x==L))
7 F. H5 k5 V4 B* K* l7 k3 r( _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)
9 o) P1 f* H. d. E9 _6 {3 i
2 A: h0 m+ U, Y" k+ |& G; j8 e* f) qopt<-optim(par=init,fn=ell,method="L-BFGS-B",
% _% P" B9 M8 r6 D  lower=c(-Inf,0.01),
  upper=c(Inf,Inf),hessian=TRUE,
8 l+ a6 e. ^+ d( e" S6 d! x  control=list(fnscale=-1))
谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧