求高手指点写 ell,optim code

HongYu_柳

mu<-0
sigma2<-1
  n<-20
' A% p+ X- t+ \+ ]2 T  L<--0.25
  y<-rnorm(n,mu,sqrt(sigma2))
+ @: `  G$ K: ~$ O- ]  x<-y
' V" ?; a. O5 i% F) Z  c<-(y<=L)
' J% A$ x4 d8 j0 x' x2 |/ K7 q x[c]<-L
* C# ~0 d5 u% V7 O x
- i6 v) o( k8 D7 q" X+ J- QQ1:写一个ell function compute log-likelihood od x.
! e+ I* B, N6 f. \& `" o5 bQ2用 R optim来maximize这个likelihood.

  d/ W2 w/ D7 S$ ]6 Y我现在写了一个ell的试子了，但是怎么也求不出最大值。
; O! A1 v" B) e9 g# ntheta<-c(mu,sigma2)
ell<-function(theta,x,L=-0.25){
stopifnot(is.vector(theta),length(theta)==2)
if(theta[2]<=0)
% o* N& x  i* H4 R+ ~ return(NA)
+ g) E, J  n6 ]+ P0 _$ F6 a8 D0 w 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))}

$ S; Z" e! ]0 P我的求最大值试子是init<-c(mu=0,sigma2=1)
4 f) }% y" c5 J( {
opt<-optim(par=init,fn=ell,method="L-BFGS-B",
  lower=c(-Inf,0.01),
* B# j5 q' i7 G! L  upper=c(Inf,Inf),hessian=TRUE,
  control=list(fnscale=-1))
谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧