求高手指点写 ell,optim code

HongYu_柳

mu<-0
sigma2<-1
  e: _3 p  T$ H; I0 [1 E  n<-20
5 B' L1 ]. }' N4 {2 m' b  L<--0.25
' ]1 G( g7 f' u2 G( B5 B. c  y<-rnorm(n,mu,sqrt(sigma2))
- B! ~% ?& S8 G3 D$ _0 I$ E  x<-y
2 Y: h( @! o$ W  c<-(y<=L)
/ h! E# y7 `9 f# V) o5 o x[c]<-L
x
Q1:写一个ell function compute log-likelihood od x.
Q2用 R optim来maximize这个likelihood.
: m& X- l; q$ d6 {
5 }4 f5 ?' ^5 @9 S9 p我现在写了一个ell的试子了，但是怎么也求不出最大值。
theta<-c(mu,sigma2)
6 a- m' L- N4 v4 F. v4 k3 x ell<-function(theta,x,L=-0.25){
stopifnot(is.vector(theta),length(theta)==2)
6 c8 d2 x2 N2 p2 ~ if(theta[2]<=0)
" z& l& j7 h6 D0 m4 B return(NA)
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))}
# V+ u/ M% U3 a4 Q5 s! ?0 @
1 v1 L1 r$ L+ _! V* ^. F6 a4 d' @我的求最大值试子是init<-c(mu=0,sigma2=1)

$ N5 n1 F" ~. d0 O7 s) p% Dopt<-optim(par=init,fn=ell,method="L-BFGS-B",
% X- T+ G- _8 L6 i- M2 }  lower=c(-Inf,0.01),
  upper=c(Inf,Inf),hessian=TRUE,
% U0 s, E+ p, Y  control=list(fnscale=-1))
谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧