求高手指点写 ell,optim code

HongYu_柳

mu<-0
) v  w% c/ J1 R' D+ T2 n sigma2<-1
9 N4 f& F9 V- |- v$ f  n<-20
8 A2 x  f; F8 G0 ]  L<--0.25
# z5 G+ R( [! Z7 \  y<-rnorm(n,mu,sqrt(sigma2))
  x<-y
8 V: W& q* t: s/ v7 v* O  c<-(y<=L)
x[c]<-L
. J/ P8 ^5 u; y: j x
7 R/ G  N9 k$ [Q1:写一个ell function compute log-likelihood od x.
Q2用 R optim来maximize这个likelihood.

  T5 c; T8 y$ v- E" m我现在写了一个ell的试子了，但是怎么也求不出最大值。
7 Z; F  a+ B) b+ T5 b0 N- ?theta<-c(mu,sigma2)
6 H6 B# p5 d* L0 u$ U ell<-function(theta,x,L=-0.25){
stopifnot(is.vector(theta),length(theta)==2)
if(theta[2]<=0)
return(NA)
' f, H5 ^4 G2 G6 V1 m$ g  B3 [- u else(c<-(x==L))
/ B+ K( h' B+ n+ Qsum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}
( G5 R7 Q$ I4 k: E+ v4 s0 ^
我的求最大值试子是init<-c(mu=0,sigma2=1)
: ]2 C: @: l1 o: g: d/ f$ q4 m2 t
opt<-optim(par=init,fn=ell,method="L-BFGS-B",
  lower=c(-Inf,0.01),
  upper=c(Inf,Inf),hessian=TRUE,
9 ~* f$ K3 Y3 }: y, f/ C8 ^: p+ i  control=list(fnscale=-1))
$ R! S# e# O1 \8 d  ]- `  X) k谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧