求高手指点写 ell,optim code

HongYu_柳

mu<-0
1 [( m. @8 O1 E+ e sigma2<-1
  n<-20
  L<--0.25
  y<-rnorm(n,mu,sqrt(sigma2))
  x<-y
  c<-(y<=L)
x[c]<-L
* r% i# c" P$ V' e  u x
# [" V, ?1 `5 h1 {  D$ LQ1:写一个ell function compute log-likelihood od x.
  n( b7 _8 f; _, }' _Q2用 R optim来maximize这个likelihood.

我现在写了一个ell的试子了，但是怎么也求不出最大值。
( ?' }; v/ j5 k0 C$ R/ S+ |5 B. ]theta<-c(mu,sigma2)
  H$ G) ~' j/ T9 f6 a4 R ell<-function(theta,x,L=-0.25){
4 [% e# G; B5 C- o1 z. ~, p6 ^ stopifnot(is.vector(theta),length(theta)==2)
' p" l6 n* y1 I$ Q. N if(theta[2]<=0)
return(NA)
7 v' P1 f% F; F% u0 N( D else(c<-(x==L))
sum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
  M- J5 t( d  F4 K; [- v sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}
" f- F7 ?3 r. `9 m& H
我的求最大值试子是init<-c(mu=0,sigma2=1)
( W* X2 h# F9 k7 N$ a
opt<-optim(par=init,fn=ell,method="L-BFGS-B",
  lower=c(-Inf,0.01),
+ |. J) D* c# c  upper=c(Inf,Inf),hessian=TRUE,
4 W( |$ d' u4 r6 b  control=list(fnscale=-1))
谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧