求高手指点写 ell,optim code

HongYu_柳

mu<-0
sigma2<-1
  n<-20
1 O3 Z% Z! t0 ]1 T2 ?. w  L<--0.25
) H3 b* Q+ K2 l. `& b+ |/ |$ N$ j* Q  y<-rnorm(n,mu,sqrt(sigma2))
  x<-y
  c<-(y<=L)
x[c]<-L
; {7 W% _+ Q" H! @ x
Q1:写一个ell function compute log-likelihood od x.
9 M8 ]6 x; _% O( ^$ \Q2用 R optim来maximize这个likelihood.
5 K( X+ p' C1 i$ w  `+ s
8 [0 f( K+ [! z( C我现在写了一个ell的试子了，但是怎么也求不出最大值。
! z9 T& P! J! h1 L+ xtheta<-c(mu,sigma2)
( r* B3 `( h* p, p8 m( y ell<-function(theta,x,L=-0.25){
stopifnot(is.vector(theta),length(theta)==2)
if(theta[2]<=0)
$ t& I: w/ F) [" Z return(NA)
* ~$ t, u" T: a, O" Q( v5 Q; r else(c<-(x==L))
! \6 M, u$ d% Wsum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}

3 D  B# T+ C/ N, t' q$ r我的求最大值试子是init<-c(mu=0,sigma2=1)
" |  F$ U7 M' |) q1 k. V2 V
opt<-optim(par=init,fn=ell,method="L-BFGS-B",
  lower=c(-Inf,0.01),
  [/ a* I1 \, p, ]( ^  upper=c(Inf,Inf),hessian=TRUE,
2 e9 t$ V2 ?" N' o; y4 I  control=list(fnscale=-1))
1 V9 x& V* D: ?! x: J谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧