yihui
不知道……也许只是想锻炼你吧……等研究比较超前时,你就可以自己变程序实现自己的算法了。
cran
[quote]引用第60楼summer2007于2007-06-21 22:44发表的“”:
向各位高手请教一个问题。现在用WinBUGS软件就能很好的实现Gibbs抽样,但我的导师不支持我用WinBUGS、SAS这些现成的软件,他认为搞研究最好用C++自己编程序。但自己编属实很困难,而且有这么好的软件为什么不用?我毕竟是个新手,或许这个问题看得不深,请各位高手多多指教。谢谢![/quote]
没错的,个人认为有几个原因吧。
1.作者在第一页就写了:Gibbs samplers can be dengerous.
2.自己编写MCMC可以更好的理解MCMC的内容。
3.BUGS的确是很简单,但是没有认真学过MCMC的人,很容易‘abuse’的。。。
4.自己编写,你可以加一些BUGS没有的功能在里面。
oliyiyi
[quote]引用第60楼summer2007于2007-06-21 22:44发表的“”:
向各位高手请教一个问题。现在用WinBUGS软件就能很好的实现Gibbs抽样,但我的导师不支持我用WinBUGS、SAS这些现成的软件,他认为搞研究最好用C++自己编程序。但自己编属实很困难,而且有这么好的软件为什么不用?我毕竟是个新手,或许这个问题看得不深,请各位高手多多指教。谢谢![/quote]
同意cran的观点,自己编程可能慢一点,但是对于你掌握和理解方法是很有好处的,如果我是导师我也会这样要求学生的,尤其是你研究领域的那一快统计方法,作为方法的研究者而不是应用者,你必须掌握黑匣子的内容,但不一定局限到c++!
summer2007
多谢各位的建议。我的研究方向是动物模型,我们从事生物领域研究的人,数学和计算机功底比较差。但听了各位高手的意见,我现在有信心去做了。多谢!
summer2007
[quote]引用第63楼oliyiyi于2007-06-22 07:52发表的“”:
同意cran的观点,自己编程可能慢一点,但是对于你掌握和理解方法是很有好处的,如果我是导师我也会这样要求学生的,尤其是你研究领域的那一快统计方法,作为方法的研究者而不是应用者,你必须掌握黑匣子的内容,但不一定局限到c++![/quote]
我也是才接触编程,导师建议我学C++或者fortran。你说可以不局限到C++,我主要是用它们实现Gibbs抽样,那还有别的更好的吗?可以给点建议吗?
redlou
现在做生物信息的人计算机很强的啊,会很多编程语言的。
oliyiyi
[quote]引用第65楼summer2007于2007-06-22 21:53发表的“”:
我也是才接触编程,导师建议我学C++或者fortran。你说可以不局限到C++,我主要是用它们实现Gibbs抽样,那还有别的更好的吗?可以给点建议吗?[/quote]
R, matlab, sas/iml, fortran, genstat,... 任何一门高级的语言都可以,不过最好和研究小组使用软件一致
summer2007
[quote]引用第67楼oliyiyi于2007-06-22 23:21发表的“”:
R, matlab, sas/iml, fortran, genstat,... 任何一门高级的语言都可以,不过最好和研究小组使用软件一致[/quote]
奥,多谢多谢!
cran
[quote]引用第64楼summer2007于2007-06-22 21:44发表的“”:
多谢各位的建议。我的研究方向是动物模型,我们从事生物领域研究的人,数学和计算机功底比较差。但听了各位高手的意见,我现在有信心去做了。多谢![/quote]
同行阿,我是作multispeices count modelling的。其实我只是接触了MCMC的理论和简单的MH 算法。我现在都是用OpenBUGS,不过还加一些R中间Coda的package,还需要做一些sensitivity analysis。不知道你是哪个方向的?
summer2007
[quote]引用第69楼cran于2007-06-23 17:03发表的“”:
同行阿,我是作multispeices count modelling的。其实我只是接触了MCMC的理论和简单的MH 算法。我现在都是用OpenBUGS,不过还加一些R中间Coda的package,还需要做一些sensitivity analysis。不知道你是哪个方向的?[/quote]
是吗?太意外了。你说的这个模型我还真不知道,是分类性状吗?我的研究方向也是多品种,至于具体用什么模型我还在思考、研究。我也是只接触了一些M-H算法、Gibbs抽样的理论。本来想用WinBUGS实现Gibbs抽样的,听了你的建议我打算学c++了。你试过自己编吗?我还没找到感觉呢。
cran
[quote]引用第70楼summer2007于2007-06-23 21:54发表的“”:
是吗?太意外了。你说的这个模型我还真不知道,是分类性状吗?我的研究方向也是多品种,至于具体用什么模型我还在思考、研究。我也是只接触了一些M-H算法、Gibbs抽样的理论。本来想用WinBUGS实现Gibbs抽样的,听了你的建议我打算学c++了。你试过自己编吗?我还没找到感觉呢。[/quote]
说真的,我的programming不好。其实我觉得这个还是看你的主要方向了。如果你想在MCMC理论上作点贡献的话,还是做自己的程序。我的目标比较简单,发展一个模型,然后忽悠生态学家而已。。。当然,一些posterior sample还是要认真分析的。
redlou
我在做生物遗传模型的,也想用MCMC来分析,不过发展模型也挺难的,
自己编程又差,入手很慢啊
summer2007
[quote]引用第72楼redlou于2007-06-24 21:20发表的“”:
我在做生物遗传模型的,也想用MCMC来分析,不过发展模型也挺难的,
自己编程又差,入手很慢啊[/quote]
同感啊!
longoR
才发现sas居然也有bayesian procedure了,只可惜07年才出来个preliminary的,晚了人家好几步。
summer2007
[quote]引用第74楼longoR++于2007-06-27 10:37发表的“”:
才发现sas居然也有bayesian procedure了,只可惜07年才出来个preliminary的,晚了人家好几步。[/quote]
好用吗?sas8.0有吗?
schyu
强烈支持!!!
ghlian
aukland 有好几位统计大牛哦
lz在那里出来 一定也很牛滴
我是金融数学的, 最近有用到mcmc才开始学, 好难懂