- 在线时间
- 6 小时
- 最后登录
- 2012-12-16
- 注册时间
- 2012-12-7
- 听众数
- 0
- 收听数
- 0
- 能力
- 0 分
- 体力
- 63 点
- 威望
- 0 点
- 阅读权限
- 20
- 积分
- 21
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 4
- 主题
- 1
- 精华
- 0
- 分享
- 0
- 好友
- 1
升级   16.84% 该用户从未签到 - 自我介绍
- 大学本科数学专业
 |
mu<-0
Y; [, T3 T8 v; x' ~ e sigma2<-1
; Q0 d9 [) d7 @7 V0 l3 X; a n<-20% S4 P& `; B, p6 o2 y7 }# M* S! @
L<--0.25
2 u4 G' x! ?( P; e y<-rnorm(n,mu,sqrt(sigma2))
7 U3 Z+ l O0 }: _ x<-y; y& n9 \+ G: r2 k
c<-(y<=L)- m, H7 M3 `2 }( I2 \
x[c]<-L
3 Z. E3 C$ i+ Q ?+ w x
; I# [: a% ]0 \' TQ1:写一个ell function compute log-likelihood od x.0 c1 ~4 w3 z4 Z; m! Y: W: Q
Q2用 R optim来maximize这个likelihood./ r# x% ]% @9 s0 K# j$ a
. g" X9 h% Y2 s9 ^) j我现在写了一个ell的试子了,但是怎么也求不出最大值。3 z" n, k. e* F- ^0 @" @' M; [
theta<-c(mu,sigma2)! _. W1 d. u8 d0 C! R
ell<-function(theta,x,L=-0.25){1 D, ?- O0 N+ E6 k' S) E+ V; L' |
stopifnot(is.vector(theta),length(theta)==2)# ]% a& p) m2 V+ o8 E
if(theta[2]<=0)4 d, k& J1 D0 R1 v U
return(NA)
* u2 r3 v0 i7 }8 p3 Q' @' s5 O! a7 y% V9 U else(c<-(x==L))# a9 X* i* U# @7 O: P+ w
sum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE))$ y, T9 Z/ l' O+ w3 W g! @/ u( c& ~
sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}
" H3 o8 \0 }6 b# c0 `7 r5 ~' q% g0 U0 D' k
我的求最大值试子是init<-c(mu=0,sigma2=1)
1 ` q4 N( X" C& f% ?$ j) J7 f( Q" f) J; D+ M) m: T
opt<-optim(par=init,fn=ell,method="L-BFGS-B",
4 {6 U5 Y: [- p6 R% M lower=c(-Inf,0.01),* X% [7 a* h$ ? W3 G6 s/ F
upper=c(Inf,Inf),hessian=TRUE,8 {; E7 e) B7 A" Z, O
control=list(fnscale=-1))
- t' ?% q) ]. Q. K谁能看看哪里错了。谢谢了。 |
zan
|