求高手指点写 ell,optim code

HongYu_柳

mu<-0
  e: P8 G, `0 s" z( y2 U sigma2<-1
  |' s  ^; Y4 Q. j: H2 \! }5 `  n<-20
4 t$ c4 H& Q0 z5 o1 M, U3 X7 }  L<--0.25
  y<-rnorm(n,mu,sqrt(sigma2))
  x<-y
# R$ r9 |1 U6 ]* E2 j  c<-(y<=L)
x[c]<-L
( I' T* w9 z8 @# _: W x
Q1:写一个ell function compute log-likelihood od x.
Q2用 R optim来maximize这个likelihood.
. j0 n2 ]; b5 c2 p& [9 h3 u
# H, u& `7 d3 g+ o  ]3 D6 g我现在写了一个ell的试子了，但是怎么也求不出最大值。
theta<-c(mu,sigma2)
4 C. x3 Z* K. ~- c# T4 P( C8 @ ell<-function(theta,x,L=-0.25){
- D  U: r$ k9 w! d5 B( T stopifnot(is.vector(theta),length(theta)==2)
' Q4 X8 }* k9 x0 ^! Y, s( k6 i" p if(theta[2]<=0)
: Y, p- s; [4 S; S0 g- t return(NA)
3 u/ p8 t9 q: |8 a3 l else(c<-(x==L))
" ]5 {' @  w2 D( |- ?sum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}
4 O9 ?) ?4 U* E' e
; @- q# x" K& @6 j7 \/ m' X我的求最大值试子是init<-c(mu=0,sigma2=1)

opt<-optim(par=init,fn=ell,method="L-BFGS-B",
- }" y& {) L& G! _2 |" O$ R9 P6 v  lower=c(-Inf,0.01),
, u! W5 e$ j* i* g  upper=c(Inf,Inf),hessian=TRUE,
; a) `4 t, H$ L! q" \1 b  control=list(fnscale=-1))
' y7 c4 _) O8 w( r谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧