求高手指点写 ell,optim code

HongYu_柳

mu<-0
: H, I% g1 H" M) H, C sigma2<-1
8 s& {  I6 y8 i  n<-20
  L<--0.25
$ ^3 L) \& ~( m! A- f# Y$ f! |  y<-rnorm(n,mu,sqrt(sigma2))
  x<-y
; d3 V( S7 E1 E/ _* h# h  c<-(y<=L)
x[c]<-L
x
4 I+ j/ \7 o4 X& S% d1 tQ1:写一个ell function compute log-likelihood od x.
Q2用 R optim来maximize这个likelihood.

0 V! I& Z0 l' v4 x* v  v我现在写了一个ell的试子了，但是怎么也求不出最大值。
/ B; S! r* u& [2 b3 P, I( k5 ?theta<-c(mu,sigma2)
ell<-function(theta,x,L=-0.25){
; q2 m; b! U- C7 p  C7 P. K stopifnot(is.vector(theta),length(theta)==2)
/ C% x5 \7 f* w# U if(theta[2]<=0)
; [& F1 {6 m3 r2 `: P! w return(NA)
1 U' `% g) ^4 w1 ] else(c<-(x==L))
sum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}
3 X2 l5 Y# q, W/ t5 n
我的求最大值试子是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,
4 x7 o4 s7 m5 B0 H& v% G) M  control=list(fnscale=-1))
谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧