QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3044|回复: 1
打印 上一主题 下一主题

求高手指点写 ell,optim code

[复制链接]
字体大小: 正常 放大

1

主题

0

听众

21

积分

升级  16.84%

该用户从未签到

自我介绍
大学本科数学专业
跳转到指定楼层
1#
发表于 2012-12-10 04:20 |只看该作者 |正序浏览
|招呼Ta 关注Ta
mu<-0% f8 B0 Z& [4 S  R* W8 `
sigma2<-1
9 W; d# e5 X, x8 Y  n<-20# B; O* {' ?" S- N% E4 _
  L<--0.25% h& _0 Y) \% W( s1 S- I' ~% X
  y<-rnorm(n,mu,sqrt(sigma2))
# }! S. w& B4 E1 }  x<-y2 \% ?0 V; H4 }( a9 O
  c<-(y<=L)
; P  I+ O4 O# j4 [- {* n# |# m7 F x[c]<-L
: J2 m1 B/ h3 n0 ?7 e" ] x. X# C5 \" g2 i+ e  B7 q
Q1:写一个ell function compute log-likelihood od x.
  v! |+ G2 w9 NQ2用 R optim来maximize这个likelihood.- l2 V' l- i! J
* a) Y& X8 c% u9 c1 ?6 m% M
我现在写了一个ell的试子了,但是怎么也求不出最大值。# _! x: p7 R4 U; k: f5 G9 X- Q3 \
theta<-c(mu,sigma2)" y' j6 ^+ N. i( n2 K$ H
ell<-function(theta,x,L=-0.25){
7 _# B/ e" J+ }+ q! f7 f$ B, l stopifnot(is.vector(theta),length(theta)==2)
) o2 S6 e2 o( T. o7 R if(theta[2]<=0)
4 A5 u8 f2 M  \8 h return(NA)7 }4 O7 g% Q1 r7 M& q2 t
else(c<-(x==L))
% S+ ^9 O( v5 V+ M0 \, T/ `sum(dnorm(x[c],theta[1],sqrt(theta[2]),log=TRUE)). ~- d; k* A" e/ d: H
sum2(dnorm(x[!c],theta[1],sqrt(theta[2]),lop=true))}& S" I6 m- t, _! N

  r' `3 h8 _" l/ B& ~; L. N7 r我的求最大值试子是init<-c(mu=0,sigma2=1)
6 Y% n; h5 R- e. d: S0 _* Y$ w3 X8 C) |! z  l/ c
opt<-optim(par=init,fn=ell,method="L-BFGS-B",, |* H. }6 I4 S  f3 Z
  lower=c(-Inf,0.01),
6 P5 @/ G8 _3 k$ _) u# {$ [: N  upper=c(Inf,Inf),hessian=TRUE,
, K  ?; k, R: x& f  control=list(fnscale=-1))0 J) O/ P) ]2 |% \- x4 r- e
谁能看看哪里错了。谢谢了。
zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
zxtttxxqm 实名认证       

9

主题

11

听众

461

积分

升级  53.67%

  • TA的每日心情

    2014-10-6 23:12
  • 签到天数: 142 天

    [LV.7]常住居民III

    社区QQ达人 新人进步奖

    群组数学建模培训课堂1

    群组学术交流B

    群组学术交流C

    群组学术交流D

    群组学术交流A

    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-8-26 07:51 , Processed in 0.373913 second(s), 63 queries .

    回顶部