24小时热门版块排行榜    

查看: 1763  |  回复: 4

inovermyhead

新虫 (小有名气)

[求助] 求助,第一次用MrBayes算进化树的问题。 已有2人参与

我用MrBayes运算进化树,MCMC抽样,已经运算了1000万代了(48小时),但是Average standard deviation of split frequencies 始终在0.02以上,从600万代以后就降不下了,请问如何能降下来? lset nst = 6 rates = invgamma ,samplefreq = 5000。
是否和比对后序列的剪切有关,我没有剪切序列。请问还能设定那些参数?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

growlywolf

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
inovermyhead(西门吹雪170代发): 金币+2, 鼓励回帖交流 2014-05-17 19:13:30
inovermyhead: 金币+5, ★★★很有帮助, 谢谢你,我剪切了一下序列,好点了,其他参数设置上有没有什么优化? 2014-05-19 08:45:12
序列差异较大的情况下做进化树是会出现这个问题。
如果是编码序列,可以考虑只用核心保守结构域。
另外,不要用clustal做序列比对了,准确性太低。推荐muscle
2楼2014-05-17 11:42:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

growlywolf

金虫 (小有名气)

【答案】应助回帖

★ ★
kx444555: 金币+2, 鼓励交流 2014-06-03 16:17:12
不要意思,我用这玩意儿做进化树还是七八年前的事情,具体设置早就忘记了
3楼2014-05-19 10:43:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lq102021

金虫 (小有名气)

【答案】应助回帖

★ ★
kx444555: 金币+2, 鼓励交流 2014-06-03 16:17:15
你使用核酸还是蛋白序列啊?你的参数时怎么设置的,我用的是蛋白序列,用modeltest检测最适模型为JTT+I+G,但是我不知道如何设置lset 和prset,你是怎么设的? 我估计你的问题和我的问题是一样的。还是求有经验的人解答一下吧
4楼2014-06-02 20:31:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

inovermyhead

新虫 (小有名气)

我用的核算序列,lset=6, 蛋白序列不会设啊.

[ 发自手机版 http://muchong.com/3g ]
5楼2014-06-02 20:33:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 inovermyhead 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 xvJS1BYDxvBS 2026-10-31 11/550 2026-11-02 00:21 by xzcv0Tehd8lQ
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-02 00:08 by xzcv0Tehd8lQ
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-02 00:08 by xzcv0Tehd8lQ
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 23:58 by xzcv0Tehd8lQ
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 23:55 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 15/750 2026-11-01 23:45 by xzcv0Tehd8lQ
[基金申请] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:11 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:11 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 22:05 by xzcv0Tehd8lQ
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 21:45 by xzcv0Tehd8lQ
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-01 21:19 by xzcv0Tehd8lQ
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 20:03 by xzcv0Tehd8lQ
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 16:52 by BCdvBAEpt3zm
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 16:45 by BCdvBAEpt3zm
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +9 zuesli 2026-10-27 15/750 2026-11-01 12:59 by Equinoxhua
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 12:40 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 03:59 by YXM1FDCSE8Bn
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
[考博] 申博求助 20+3 Katrinae 2026-10-26 4/200 2026-10-27 09:56 by 北京莱茵编辑
信息提示
请填处理意见