求高手指点写 ell,optim code

HongYu_柳

mu<-0
sigma2<-1
  n<-20
  L<--0.25
6 w3 f8 [+ h% [  y<-rnorm(n,mu,sqrt(sigma2))
  x<-y
9 I! l& W0 W3 U3 ^  c<-(y<=L)
x[c]<-L
x
Q1:写一个ell function compute log-likelihood od x.
Q2用 R optim来maximize这个likelihood.
) w7 f% i  Q* j& e, r
7 W' Y8 L3 j7 F3 R; w0 W9 M& ?我现在写了一个ell的试子了，但是怎么也求不出最大值。
; S" J3 s* h( y, ^" F; [theta<-c(mu,sigma2)
9 _' o) E9 b# U  D- J ell<-function(theta,x,L=-0.25){
8 E" P' ^4 @2 [7 i stopifnot(is.vector(theta),length(theta)==2)
, \: j: l) o! Z3 e if(theta[2]<=0)
return(NA)
else(c<-(x==L))
* ~! Y1 i' F  }sum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
% o0 q# v4 }) h: U) o sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}

我的求最大值试子是init<-c(mu=0,sigma2=1)

opt<-optim(par=init,fn=ell,method="L-BFGS-B",
  lower=c(-Inf,0.01),
$ z3 A7 A1 d5 r% e0 X& n$ ], u1 G  upper=c(Inf,Inf),hessian=TRUE,
  control=list(fnscale=-1))
$ \0 Y2 \, L7 x1 p谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧