求高手指点写 ell,optim code

HongYu_柳

mu<-0
sigma2<-1
' e) @: j( x$ E7 t; i  n<-20
  L<--0.25
5 l; \) y1 @' c3 w  y<-rnorm(n,mu,sqrt(sigma2))
  x<-y
; ^& ~8 M% Z" S# C* D  c<-(y<=L)
! T( L/ ]$ k# j4 s3 W x[c]<-L
x
- ^  q! t1 ?# v0 t! {9 M7 eQ1:写一个ell function compute log-likelihood od x.
/ i2 A8 e% B! u* J9 PQ2用 R optim来maximize这个likelihood.

我现在写了一个ell的试子了，但是怎么也求不出最大值。
theta<-c(mu,sigma2)
, p1 w: `; o0 V7 F' f ell<-function(theta,x,L=-0.25){
stopifnot(is.vector(theta),length(theta)==2)
if(theta[2]<=0)
1 r6 w: j, T5 W3 { return(NA)
else(c<-(x==L))
9 A& Z2 u1 F; R! \, Y2 I# @sum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
  z7 J. k8 d( ? sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}

我的求最大值试子是init<-c(mu=0,sigma2=1)
, ]. ^- x; z9 [
opt<-optim(par=init,fn=ell,method="L-BFGS-B",
  lower=c(-Inf,0.01),
/ Z' O4 F  y" A# r7 F& {$ D8 q3 o  upper=c(Inf,Inf),hessian=TRUE,
# ~, i% i( Z% Q) E- t* b  control=list(fnscale=-1))
谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧