求高手指点写 ell,optim code

HongYu_柳

mu<-0
sigma2<-1
  n<-20
' d5 O; @2 l0 F% g# N  L<--0.25
% h$ N/ E+ v- J* n  y<-rnorm(n,mu,sqrt(sigma2))
  x<-y
' f7 O# s/ A1 g0 N* Y+ N# Y  c<-(y<=L)
x[c]<-L
x
Q1:写一个ell function compute log-likelihood od x.
Q2用 R optim来maximize这个likelihood.

我现在写了一个ell的试子了，但是怎么也求不出最大值。
theta<-c(mu,sigma2)
ell<-function(theta,x,L=-0.25){
stopifnot(is.vector(theta),length(theta)==2)
' w1 A: Q& T+ E- v) r5 @- T$ m if(theta[2]<=0)
return(NA)
  k1 @! O  C) T! C# l* t0 S else(c<-(x==L))
! m! \1 I. H0 m; hsum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
5 T9 `- s$ b2 G- _) z, u sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}

  @- f4 v+ G  L& G+ y2 q( H7 G我的求最大值试子是init<-c(mu=0,sigma2=1)
, O9 |+ e+ ]2 Z1 E2 P
& K; D- X( P' ]opt<-optim(par=init,fn=ell,method="L-BFGS-B",
  lower=c(-Inf,0.01),
/ n7 Q5 g- n* P  upper=c(Inf,Inf),hessian=TRUE,
( ]2 `* [' L  a( e  ]" W  control=list(fnscale=-1))
谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧