求高手指点写 ell,optim code

HongYu_柳

mu<-0
sigma2<-1
% v7 h1 p5 m1 K# O3 K  n<-20
+ x6 ^/ `; f- }4 I% d$ H) x1 F  L<--0.25
: M+ V- b4 e$ C& k& @  y<-rnorm(n,mu,sqrt(sigma2))
1 D, m/ L, }4 }9 _$ I( t  x<-y
% ~' k: C4 `- h, \/ u  c<-(y<=L)
; ~5 w2 q$ q& e) g5 ? x[c]<-L
x
: {9 B! c! b( eQ1:写一个ell function compute log-likelihood od x.
" }' a7 v2 H2 k7 L% [8 l2 D; ?Q2用 R optim来maximize这个likelihood.
9 i# W0 B, c, b5 o( s# p5 K
我现在写了一个ell的试子了，但是怎么也求不出最大值。
theta<-c(mu,sigma2)
ell<-function(theta,x,L=-0.25){
6 [* q+ B  e# ]4 R# Z stopifnot(is.vector(theta),length(theta)==2)
1 Q# R: ^5 }6 C. O8 O# T2 E  l+ x8 N if(theta[2]<=0)
return(NA)
* F+ J8 T/ k2 e  x% d0 k else(c<-(x==L))
sum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
$ l# x' T( Y* H  w2 g+ e sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}

我的求最大值试子是init<-c(mu=0,sigma2=1)

opt<-optim(par=init,fn=ell,method="L-BFGS-B",
  lower=c(-Inf,0.01),
  upper=c(Inf,Inf),hessian=TRUE,
  control=list(fnscale=-1))
# Y7 h. W( N4 ^( T5 R/ U0 h2 d谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧