求高手指点写 ell,optim code

HongYu_柳

mu<-0
3 Q5 j2 p7 R' e+ z4 U sigma2<-1
  n<-20
% S: T1 t3 M7 n( G  L<--0.25
8 i2 f' E, C: ~3 Z  y<-rnorm(n,mu,sqrt(sigma2))
7 T& \7 _% e' y* W5 L, a6 \  F4 b  x<-y
7 s: C* ~' `' q6 n  r' ^  c<-(y<=L)
x[c]<-L
x
Q1:写一个ell function compute log-likelihood od x.
Q2用 R optim来maximize这个likelihood.

. P# f+ r" i- H9 B7 q我现在写了一个ell的试子了，但是怎么也求不出最大值。
theta<-c(mu,sigma2)
- X5 m6 R  N8 \) W ell<-function(theta,x,L=-0.25){
* [% O) [: y# s' L7 G+ z stopifnot(is.vector(theta),length(theta)==2)
( d$ I8 X) z3 m  o if(theta[2]<=0)
return(NA)
else(c<-(x==L))
7 ^9 s/ i7 o" Lsum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}

2 K/ T  K+ T& w2 _: Z我的求最大值试子是init<-c(mu=0,sigma2=1)

: [& t& O# N8 Q, f9 zopt<-optim(par=init,fn=ell,method="L-BFGS-B",
  lower=c(-Inf,0.01),
; |; x6 n( D3 H% v0 k* c/ V/ f+ ?$ ^  upper=c(Inf,Inf),hessian=TRUE,
& O9 ~  M/ B2 c7 i0 ~* x1 j  control=list(fnscale=-1))
谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧