求高手指点写 ell,optim code

HongYu_柳

mu<-0
sigma2<-1
  n<-20
  L<--0.25
  y<-rnorm(n,mu,sqrt(sigma2))
  x<-y
  c<-(y<=L)
x[c]<-L
x
- }2 @, d; w( d8 w+ DQ1:写一个ell function compute log-likelihood od x.
Q2用 R optim来maximize这个likelihood.

) x4 j; j, P$ q我现在写了一个ell的试子了，但是怎么也求不出最大值。
) F, ~" q8 F+ s/ x/ x& J3 Ttheta<-c(mu,sigma2)
: t5 |& m0 E$ Y* W/ X. a9 U ell<-function(theta,x,L=-0.25){
stopifnot(is.vector(theta),length(theta)==2)
if(theta[2]<=0)
return(NA)
else(c<-(x==L))
sum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
+ F; a7 x. Q* d( H# V3 t1 i  v sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}
. _$ }/ a1 k; k: @9 g: R
我的求最大值试子是init<-c(mu=0,sigma2=1)
0 @1 c4 W) x8 N
9 i: {/ ~7 e9 Y* ~% Vopt<-optim(par=init,fn=ell,method="L-BFGS-B",
7 y2 \3 Z) F2 Z5 k3 {  lower=c(-Inf,0.01),
  upper=c(Inf,Inf),hessian=TRUE,
  control=list(fnscale=-1))
. Z9 {. e5 a: U: Z谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧