求高手指点写 ell,optim code

HongYu_柳

mu<-0
sigma2<-1
9 \) a2 ~3 a" I6 x  ^  n<-20
( i$ e  q' u* ]  L<--0.25
& E& Z* W$ x# K* P  y<-rnorm(n,mu,sqrt(sigma2))
1 V3 p- k3 M- F+ h9 k' S' f4 D  x<-y
  c<-(y<=L)
x[c]<-L
- a8 D$ }5 q- I7 [2 G x
Q1:写一个ell function compute log-likelihood od x.
Q2用 R optim来maximize这个likelihood.
! s) D" m" d: g' u& N& L6 J
我现在写了一个ell的试子了，但是怎么也求不出最大值。
7 t# M5 J" B$ k7 F% Xtheta<-c(mu,sigma2)
ell<-function(theta,x,L=-0.25){
stopifnot(is.vector(theta),length(theta)==2)
# ^5 @, }: G; h! l9 O( n if(theta[2]<=0)
' \) L9 H9 t9 m return(NA)
! E7 V! I. c, u& ~0 K2 \" m; x else(c<-(x==L))
sum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
: v: j! j$ g# J4 T& X/ i3 s sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}
/ u# ~- ^+ |. l2 J
我的求最大值试子是init<-c(mu=0,sigma2=1)

opt<-optim(par=init,fn=ell,method="L-BFGS-B",
% R: R6 d; _1 A8 M  lower=c(-Inf,0.01),
  upper=c(Inf,Inf),hessian=TRUE,
2 D) j, P. t2 ~, ?+ ]$ P3 C  control=list(fnscale=-1))
谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧