求高手指点写 ell,optim code

HongYu_柳

mu<-0
6 r5 E" R/ u; z- L sigma2<-1
) s& W3 l$ Y: l# t/ Z  n<-20
  L<--0.25
  y<-rnorm(n,mu,sqrt(sigma2))
  x<-y
# X7 t6 _0 [' V  c<-(y<=L)
; H: b  X7 t: |& X. j x[c]<-L
x
Q1:写一个ell function compute log-likelihood od x.
Q2用 R optim来maximize这个likelihood.
" ~2 y8 a+ B4 e  V, T" q+ l
% }/ n( N6 F$ V2 c; t我现在写了一个ell的试子了，但是怎么也求不出最大值。
- W8 a. D3 h, w! ^. \9 j0 ntheta<-c(mu,sigma2)
* {. \$ z& z1 |" z ell<-function(theta,x,L=-0.25){
stopifnot(is.vector(theta),length(theta)==2)
if(theta[2]<=0)
( @+ j1 N6 C  l% O5 Q return(NA)
4 w4 r! M4 L4 _" z( [ else(c<-(x==L))
' W5 i6 c1 l5 j- Lsum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))
sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}
+ Q8 r- O% ^/ A2 D! y  D/ P
' C; |1 A& H. Z: e2 a5 L( H) m我的求最大值试子是init<-c(mu=0,sigma2=1)

opt<-optim(par=init,fn=ell,method="L-BFGS-B",
  lower=c(-Inf,0.01),
: q9 g- @! ]7 f2 g  upper=c(Inf,Inf),hessian=TRUE,
  control=list(fnscale=-1))
谁能看看哪里错了。谢谢了。

zxtttxxqm

这个的话你去问版主吧