求高手指点写 ell,optim code

HongYu_柳

mu<-0
0 W1 B, k0 a0 I0 ? sigma2<-1
+ ]9 y, F7 p$ V- ~" c; S5 y, ]2 Y  n<-20
  L<--0.25
% w3 @4 K7 _" x  y<-rnorm(n,mu,sqrt(sigma2))
  x<-y
2 X" d' p' y% X- z0 _5 X! \  c<-(y<=L)
x[c]<-L
8 l, q3 t0 N) H# S: c" Q 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)
if(theta[2]<=0)
  N5 c" W- p: w$ Q return(NA)
3 N' l  y" ?" V" [5 ~) P* |0 b else(c<-(x==L))
1 C1 i9 p" _( L/ I& Psum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}
: V- Z  I, D* N1 [% }5 {
我的求最大值试子是init<-c(mu=0,sigma2=1)
( s4 \) v8 x$ Y- \! X6 I( b1 C6 T
0 h# |$ o# l; `3 R$ ?. U; g) E+ |opt<-optim(par=init,fn=ell,method="L-BFGS-B",
6 X: `; q( Q$ k: `) N5 k  lower=c(-Inf,0.01),
  upper=c(Inf,Inf),hessian=TRUE,
$ h) c4 e0 ~# H/ L, t  control=list(fnscale=-1))
谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧