求高手指点写 ell,optim code

HongYu_柳

mu<-0
: F" Q3 z  B. }! b sigma2<-1
# n8 P: w6 ]/ D: S" z8 t% R7 c  n<-20
  L<--0.25
. V2 d6 {8 g" ?5 d' Y  y<-rnorm(n,mu,sqrt(sigma2))
  x<-y
5 _; ~4 V, y2 I. U( z  c<-(y<=L)
; T% a( }3 s8 d8 m, B) y: e% e* \5 v0 ? x[c]<-L
x
Q1:写一个ell function compute log-likelihood od x.
9 Y$ K( p6 }! A/ D0 Z1 \1 d$ KQ2用 R optim来maximize这个likelihood.
5 k" \$ L' s# o; B- d
# h# t, z+ ]; |. D我现在写了一个ell的试子了，但是怎么也求不出最大值。
theta<-c(mu,sigma2)
# {/ d6 m! H$ o4 L% T# f ell<-function(theta,x,L=-0.25){
stopifnot(is.vector(theta),length(theta)==2)
if(theta[2]<=0)
return(NA)
, x5 z$ b$ @. c( z$ N else(c<-(x==L))
7 |5 K) z+ ]  H0 p, @sum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
" }7 y1 H* O4 o- O sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}

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

0 N0 z2 X7 M& E! n0 L" Jopt<-optim(par=init,fn=ell,method="L-BFGS-B",
) i$ Q  i6 h( I  lower=c(-Inf,0.01),
  upper=c(Inf,Inf),hessian=TRUE,
9 m6 d7 H5 C/ F6 Z* z- F5 M* ~) N  control=list(fnscale=-1))
: F2 F% ~$ ?3 ~8 D谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧