求高手指点写 ell,optim code

HongYu_柳

mu<-0
sigma2<-1
  n<-20
+ `# G; Z* y! C. t  X. r) `  L<--0.25
  }) V& s0 i7 K# [1 a8 P+ c  y<-rnorm(n,mu,sqrt(sigma2))
  x<-y
  c<-(y<=L)
x[c]<-L
8 B  Q2 ?, I2 N0 S3 y9 I x
Q1:写一个ell function compute log-likelihood od x.
Q2用 R optim来maximize这个likelihood.

我现在写了一个ell的试子了，但是怎么也求不出最大值。
- m# `& `% @$ G& |7 Qtheta<-c(mu,sigma2)
4 l& E8 Y0 e. [3 n4 D) U ell<-function(theta,x,L=-0.25){
stopifnot(is.vector(theta),length(theta)==2)
if(theta[2]<=0)
) |1 R5 ~( `# t( I: R 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))}

# b/ {; C: U- f! X) |, k我的求最大值试子是init<-c(mu=0,sigma2=1)
3 x& E5 t' p4 Y9 I( }! T
- `  X$ |" |+ U* c) c* Aopt<-optim(par=init,fn=ell,method="L-BFGS-B",
6 N7 V0 Q% q# F  K  w  lower=c(-Inf,0.01),
  upper=c(Inf,Inf),hessian=TRUE,
  control=list(fnscale=-1))
. t) ?% S: [/ f, `6 i2 [谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧