求高手指点写 ell,optim code

HongYu_柳

mu<-0
sigma2<-1
  n<-20
% k" `% h+ a" j$ |7 c) U& `  L<--0.25
' Z: }; D6 C: m7 `2 {/ J$ i- B  y<-rnorm(n,mu,sqrt(sigma2))
% G6 h2 _8 |. x, ]2 u9 N4 M  x<-y
/ u3 F9 ~9 T9 x  c<-(y<=L)
x[c]<-L
x
Q1:写一个ell function compute log-likelihood od x.
1 x! c/ I$ T7 vQ2用 R optim来maximize这个likelihood.
( V6 r; ]6 ^1 U0 N2 J
我现在写了一个ell的试子了，但是怎么也求不出最大值。
theta<-c(mu,sigma2)
ell<-function(theta,x,L=-0.25){
& P" q: [9 z5 a6 n# u8 g stopifnot(is.vector(theta),length(theta)==2)
4 Z- w* _, u0 \: ~  z if(theta[2]<=0)
return(NA)
+ h: x* D9 ?; \4 r5 \# e else(c<-(x==L))
, }8 j4 c* H! F+ S1 V, ]$ {sum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}
: M8 O7 a( I8 l7 F9 \) X
3 B! S0 `8 c8 `5 a0 T% U* f我的求最大值试子是init<-c(mu=0,sigma2=1)

opt<-optim(par=init,fn=ell,method="L-BFGS-B",
  lower=c(-Inf,0.01),
  upper=c(Inf,Inf),hessian=TRUE,
% E( Q5 b* Q2 v! b- g; D  H  control=list(fnscale=-1))
' A! Z% j& F; [0 Y8 v谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧