求高手指点写 ell,optim code

HongYu_柳

mu<-0
sigma2<-1
4 e# b3 m! V: [, v/ ?# s  n<-20
- w" C6 P( `: e2 q  c  L<--0.25
  y<-rnorm(n,mu,sqrt(sigma2))
. Z* |% w5 j) F  x<-y
  c<-(y<=L)
6 v' K/ @4 `- W) a8 b6 s% T x[c]<-L
3 u! d, V) d+ z" g# U' t x
1 R9 ?+ ~% u" s8 G1 f8 bQ1:写一个ell function compute log-likelihood od x.
Q2用 R optim来maximize这个likelihood.

我现在写了一个ell的试子了，但是怎么也求不出最大值。
* K% L# y+ o. w  [theta<-c(mu,sigma2)
/ r' _# X% H3 Z ell<-function(theta,x,L=-0.25){
+ K9 }( E4 p2 @! Y, R9 x stopifnot(is.vector(theta),length(theta)==2)
7 d* o/ w1 f3 `* x4 l+ T if(theta[2]<=0)
return(NA)
else(c<-(x==L))
sum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
* x% f' _2 I! q/ e sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}

我的求最大值试子是init<-c(mu=0,sigma2=1)
# {: l! p/ @6 x5 q" b
opt<-optim(par=init,fn=ell,method="L-BFGS-B",
6 q/ h; M7 L% Z0 M8 ^- A  Y/ u2 ]  lower=c(-Inf,0.01),
1 T& B$ M0 m9 ~/ ^: D! ]' s* V' s3 x  upper=c(Inf,Inf),hessian=TRUE,
. f3 }5 [$ W' z9 |4 v5 m" c/ p$ a% ]  control=list(fnscale=-1))
谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧