求高手指点写 ell,optim code

HongYu_柳

mu<-0
  @9 |$ E9 s! W! f* I3 { sigma2<-1
* L: Z  h: P0 K4 g) J3 [! k  n<-20
  L<--0.25
  ^) b4 y8 h. f9 C  y<-rnorm(n,mu,sqrt(sigma2))
7 }3 I( V/ U0 O; {  x<-y
6 h, p4 p* P; Y1 y; H$ f  Z/ J  c<-(y<=L)
x[c]<-L
x
" h5 Y; R6 m3 k5 h* a& sQ1:写一个ell function compute log-likelihood od x.
Q2用 R optim来maximize这个likelihood.

我现在写了一个ell的试子了，但是怎么也求不出最大值。
theta<-c(mu,sigma2)
! n0 |* b) h) O4 F4 | ell<-function(theta,x,L=-0.25){
; Z: Y, l4 D' g% d! U stopifnot(is.vector(theta),length(theta)==2)
if(theta[2]<=0)
7 v0 n7 R1 A/ t; r% k return(NA)
else(c<-(x==L))
/ o+ M. e& h& i; [5 W; {, L" gsum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}
: A8 W) V0 A5 G3 o0 @% I$ K
我的求最大值试子是init<-c(mu=0,sigma2=1)

+ r" m! p% B0 F- b  s2 r6 I7 `opt<-optim(par=init,fn=ell,method="L-BFGS-B",
  lower=c(-Inf,0.01),
  upper=c(Inf,Inf),hessian=TRUE,
& b' ?! N4 {( Q) Z( t  control=list(fnscale=-1))
2 P8 i1 |/ l; m& q谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧