求高手指点写 ell,optim code

HongYu_柳

mu<-0
& z& {/ w- d  T: P- |- C sigma2<-1
& z9 J! i4 J' x$ b; D  n<-20
  L<--0.25
( l  g2 O. f/ s% c9 u( X  y<-rnorm(n,mu,sqrt(sigma2))
- c$ B9 u- X* h/ F  x<-y
+ l# b2 E2 @" q$ r/ }: D  c<-(y<=L)
x[c]<-L
x
  D, e* p" j$ `7 w6 Y, z  dQ1:写一个ell function compute log-likelihood od x.
$ Q) s& n1 X# r  G8 n" N  G% {6 zQ2用 R optim来maximize这个likelihood.
3 ?  r' w, P+ L- b" I# V
我现在写了一个ell的试子了，但是怎么也求不出最大值。
theta<-c(mu,sigma2)
5 a) J0 g8 D6 Y ell<-function(theta,x,L=-0.25){
6 t$ o) j  E0 o% A/ ^( d' | stopifnot(is.vector(theta),length(theta)==2)
4 [6 l" a( Z0 q8 ` if(theta[2]<=0)
return(NA)
9 j. z2 z. k( H5 L, k  x else(c<-(x==L))
- y  c4 B# I8 j% b) H6 ~sum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
( b. K& E* y: A4 e) y4 l4 m, r sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}

( r$ T8 H2 P* U5 k/ p我的求最大值试子是init<-c(mu=0,sigma2=1)
% G  L: Q5 l2 v  ^6 y9 b$ [: I
opt<-optim(par=init,fn=ell,method="L-BFGS-B",
  lower=c(-Inf,0.01),
  upper=c(Inf,Inf),hessian=TRUE,
6 D1 V# F: |+ V0 G* \+ r, w4 _  control=list(fnscale=-1))
谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧