求高手指点写 ell,optim code

HongYu_柳

mu<-0
sigma2<-1
3 @1 y  H/ r! }  M0 L. G  n<-20
4 d0 I  g, P' T9 S& L4 a# r: {4 X  L<--0.25
* H" n; e# c0 D, ]  y<-rnorm(n,mu,sqrt(sigma2))
9 X9 K: u; D1 a6 x$ s' h* v" P  x<-y
4 @4 x7 j8 ?2 `5 i, e  c<-(y<=L)
: p/ c4 A* K) ` x[c]<-L
6 D# d3 m& h& s% M' ~6 r- P7 y x
Q1:写一个ell function compute log-likelihood od x.
# d( `' K( v  C$ H2 U* x6 |2 ^& FQ2用 R optim来maximize这个likelihood.

我现在写了一个ell的试子了，但是怎么也求不出最大值。
theta<-c(mu,sigma2)
7 B( {* u6 w6 |" ?$ P' J8 \ ell<-function(theta,x,L=-0.25){
stopifnot(is.vector(theta),length(theta)==2)
if(theta[2]<=0)
return(NA)
else(c<-(x==L))
. Z$ a9 H1 K1 ]6 i% V0 c7 usum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}
6 z5 q2 X* W: E5 k
+ L1 Y% f; c( f! X4 t! R我的求最大值试子是init<-c(mu=0,sigma2=1)

% j5 B9 L) @: w. g6 y% Mopt<-optim(par=init,fn=ell,method="L-BFGS-B",
/ Q' J# c% y7 _* ]. }& b$ j9 z  lower=c(-Inf,0.01),
  upper=c(Inf,Inf),hessian=TRUE,
  control=list(fnscale=-1))
/ Z8 C/ W: c$ h: `9 m谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧