求高手指点写 ell,optim code

HongYu_柳

mu<-0
sigma2<-1
* Y, x; ?5 e4 i" m" n$ v  n<-20
  L<--0.25
# n  f% z7 r' {0 W9 W) S/ \5 a: o  y<-rnorm(n,mu,sqrt(sigma2))
  x<-y
  c<-(y<=L)
; a1 \# f) s: k8 y2 l x[c]<-L
x
. T4 R5 X1 X7 q" y% b1 GQ1:写一个ell function compute log-likelihood od x.
Q2用 R optim来maximize这个likelihood.

1 h1 U# O* W8 N1 y9 A& T2 j我现在写了一个ell的试子了，但是怎么也求不出最大值。
theta<-c(mu,sigma2)
" T2 F( i6 }' H! P( u( c ell<-function(theta,x,L=-0.25){
. F: ~* X; ]5 V/ v, X stopifnot(is.vector(theta),length(theta)==2)
if(theta[2]<=0)
return(NA)
. V9 }" @; Q  [+ j4 L/ ^9 w2 V6 R+ v9 D else(c<-(x==L))
4 Y# Q+ y* E+ a1 ^9 I2 {# Lsum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}

* X7 F$ l; a7 ^( R' L+ E我的求最大值试子是init<-c(mu=0,sigma2=1)
8 G  _5 Y; V) q: h/ t" v
8 p9 s* {, {9 T* A. J6 _* Bopt<-optim(par=init,fn=ell,method="L-BFGS-B",
  lower=c(-Inf,0.01),
  upper=c(Inf,Inf),hessian=TRUE,
  control=list(fnscale=-1))
谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧