求高手指点写 ell,optim code

HongYu_柳

mu<-0
: R! g+ S) o% x( L sigma2<-1
; p- H; D- p6 |/ z- ^+ h  n<-20
2 e0 M8 g5 R( K' a+ o% `7 I! N4 O  L<--0.25
  y<-rnorm(n,mu,sqrt(sigma2))
, L- h( ^" ~, ~& Y5 H  x<-y
3 `( ]/ C" [( {2 a! r& f- S  c<-(y<=L)
x[c]<-L
x
2 C0 \9 f5 P7 g! R' F  zQ1:写一个ell function compute log-likelihood od x.
- x. t( e9 @' p& z8 [8 tQ2用 R optim来maximize这个likelihood.

我现在写了一个ell的试子了，但是怎么也求不出最大值。
theta<-c(mu,sigma2)
5 Q8 |+ [* ]& J/ v4 S ell<-function(theta,x,L=-0.25){
* E6 U; _) t  [ stopifnot(is.vector(theta),length(theta)==2)
if(theta[2]<=0)
4 E$ e% I# I1 y) N5 u2 U- F return(NA)
4 O4 T, _9 }. K, M else(c<-(x==L))
sum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}
7 `" o* n; R2 }1 j
我的求最大值试子是init<-c(mu=0,sigma2=1)

( Q6 J* Y4 Y1 c2 x) u% [  U2 \opt<-optim(par=init,fn=ell,method="L-BFGS-B",
- N; n* i1 z! w  lower=c(-Inf,0.01),
  upper=c(Inf,Inf),hessian=TRUE,
  control=list(fnscale=-1))
' I( h% K1 R. r1 V谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧