求高手指点写 ell,optim code

HongYu_柳

mu<-0
sigma2<-1
6 I5 w5 ]* e, Y9 ~  n<-20
  L<--0.25
4 o! C: ^2 C2 m5 D2 |# K3 }* i  y<-rnorm(n,mu,sqrt(sigma2))
# H( `% h: Y) q! e1 d) I$ C  x<-y
8 W" [5 f1 D; V9 A9 X" l/ ~4 {  c<-(y<=L)
3 p2 m# n  j3 Q4 ^, c x[c]<-L
x
9 V. \0 _+ x3 l% HQ1:写一个ell function compute log-likelihood od x.
0 |+ k& x3 ^3 q$ ]: e3 s/ @Q2用 R optim来maximize这个likelihood.

我现在写了一个ell的试子了，但是怎么也求不出最大值。
theta<-c(mu,sigma2)
4 z0 l7 u& e- W( I+ v( B- w% ? ell<-function(theta,x,L=-0.25){
/ t; _7 W+ k3 f% { stopifnot(is.vector(theta),length(theta)==2)
if(theta[2]<=0)
% T: |# k. E- z! K, Z* g; @/ D. G return(NA)
( N% v" C5 w: a  w0 a6 |" s else(c<-(x==L))
5 J0 k5 Y; `% C% X2 ]sum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
5 A3 f; d' a3 J$ J# O sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}
% |2 O" A/ X0 l  n; s( I
我的求最大值试子是init<-c(mu=0,sigma2=1)

' p2 N4 r+ r. b6 W. Vopt<-optim(par=init,fn=ell,method="L-BFGS-B",
  lower=c(-Inf,0.01),
  upper=c(Inf,Inf),hessian=TRUE,
! O7 j5 |# [( P6 g6 `2 r  control=list(fnscale=-1))
谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧