求高手指点写 ell,optim code

HongYu_柳

mu<-0
sigma2<-1
! G+ C) W1 L+ L  n<-20
; p: `8 F8 w- E/ n% O- z  L<--0.25
% V! t; f+ V/ o' e- j  y<-rnorm(n,mu,sqrt(sigma2))
  x<-y
  c<-(y<=L)
x[c]<-L
8 [9 j. x6 T9 x. ^3 j$ Y x
Q1:写一个ell function compute log-likelihood od x.
7 p- d' {# J' r6 u! C. [6 E1 P: XQ2用 R optim来maximize这个likelihood.
6 Q7 s# L+ m6 k
我现在写了一个ell的试子了，但是怎么也求不出最大值。
theta<-c(mu,sigma2)
ell<-function(theta,x,L=-0.25){
( r  I+ G) T/ a stopifnot(is.vector(theta),length(theta)==2)
+ a' z9 d. ]3 }, Q3 b if(theta[2]<=0)
return(NA)
else(c<-(x==L))
) i+ y$ }! m- x6 n8 Psum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}

8 R6 C. }/ M+ r' M7 r# s) r' L) n我的求最大值试子是init<-c(mu=0,sigma2=1)
% g3 D0 _& Z8 \8 B8 |! G
opt<-optim(par=init,fn=ell,method="L-BFGS-B",
0 h  H  \% I- b1 v4 _  lower=c(-Inf,0.01),
+ T' \+ v. L" ~! v  upper=c(Inf,Inf),hessian=TRUE,
/ k- n; t  h! B1 a/ i1 C  control=list(fnscale=-1))
2 |0 w0 N3 K7 \7 Y' A  c/ W谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧