求高手指点写 ell,optim code

HongYu_柳

mu<-0
$ ~6 T2 F" _: u& I; f( G# O; {) v sigma2<-1
  n<-20
  L<--0.25
3 U9 Z5 z4 M: R3 ^4 F  y<-rnorm(n,mu,sqrt(sigma2))
  x<-y
  c<-(y<=L)
; E' g  a% X( b% G* L: b) V8 [ x[c]<-L
# H, P9 g' Z4 V0 j1 A x
Q1:写一个ell function compute log-likelihood od x.
' N: }' b! d5 S' RQ2用 R optim来maximize这个likelihood.

8 \( T8 l3 H1 C1 T% w我现在写了一个ell的试子了，但是怎么也求不出最大值。
theta<-c(mu,sigma2)
ell<-function(theta,x,L=-0.25){
) P, Q& f6 Q2 j stopifnot(is.vector(theta),length(theta)==2)
+ k% y6 r3 y/ x8 }1 y if(theta[2]<=0)
return(NA)
8 v( p" T' W  R* }6 T& C" U else(c<-(x==L))
sum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
9 ?2 y2 u$ k" J$ T sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}

8 y: P$ K7 q0 h1 @我的求最大值试子是init<-c(mu=0,sigma2=1)
* Y& `) c& P! h
opt<-optim(par=init,fn=ell,method="L-BFGS-B",
  lower=c(-Inf,0.01),
  upper=c(Inf,Inf),hessian=TRUE,
  control=list(fnscale=-1))
谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧