应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 生物医药区 » 微生物 » 综合其他 » 请各位来看棵树
52/1返回列表
查看: 2686  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 MicEPI ，点击这里进行查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

秋电工

新虫 (小有名气)

[求助] 请各位来看棵树

近日学习用mega4.0建系统进化树，大概会操作了，但仍然有很多疑问。坛子里关于mega的帖子看了遍，迷茫发问的朋友多，解答的大牛人也有，我的问题和很多求助的虫友相似，现系统提出，渴望知道的大侠尽早看到此贴，尽早给若干小虫指点下迷津，不甚感激！
一，附件中的小树是用mega4.0，Bootsrap Value设置1000次建成。具体步骤如下：
a、打开mega软件—Alignment 下拉—Alignment exposure/clustal—create a new alignment --然后根据实际情况(DNA or 蛋白序列)选择Yes
b、Edit—Insert sequence from file（.fastal的格式）；
c、Algin by clustal W---Data—Export Alignment—Mega fomat；
d、打开上一步保存的文件—Distance—compute pairwise—compute；
e、Phylogeny—Bootstrap Test of phylogeny—Neighbor-Joining—将Boostrap（500 replicates；seed=64238修改为1000 replicates）--compute---获得进化树—设好格式---image –save as TIF file
问题来了，在坛子里看见各位大侠给出的步骤，要用Clustal X软件进行序列比对之后才能建树，请问，mega4.0中是不是嵌合了clustal么，为什么有大侠说还要下载一个clustal软件呢？上面的建树步骤C是不是算经过了clustal呢？
二，附件中小树节点处数值的含义，请各位解释下为什么HX4和AB546196.1（bootstrap值100）与HX6和AJ308316.1（bootstrap值77）两类群聚一起后又97呢?是否意味着HX4和AB546196.1亲缘关系很近，甚至是同一个东西，HX6和AJ308316.1关系较近，这两类群又聚一起，而且bootstrap值97，大于77，很费解呀很费解。
三，树下的标尺有什么用？去算枝长吗？枝长又代表什么，进化距离吗？如果枝长是进化距离，那bootstrap值又代表什么？
四，在坛子里逛的时候，看到有大侠说“遗传标度法和步长值法是用来表示两种不同的结果时用的 步长值常用在鉴定菌种的进化树中 而遗传距离法常用在表示不同菌株间进化关系的进化树中 ”。什么是遗传标度法，什么是步长值法？本人建树目的是想得到菌株间的进化关系
五，看文献过程中，发现有部分文献的进化树要么只有标尺，要么只有bootstrap值，而且也没有对数值，标尺的含义做出解释，究竟什么样的树才算规范？
昨晚码了一个求助贴，满怀希望发的时候，系统维护，伤啊伤！今天抖擞一下，再次情绪饱满的发一个，敬候各位的佳音哈！
回复此楼

» 猜你喜欢
1楼 2011-04-18 09:54:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xp198766

铁杆木虫 (著名写手)

小木虫职业打酱油滴~~！
引用回帖:
1283233楼: Originally posted by 黑藻先生 at 2012-11-15 09:57:48
有关5
我觉得貌似BOOTSTRAP值这个可信度是不是有待商榷
我做的时候跟老板商量的据说是75以下的都容易变
这个变是只插入或删除别的序列后容易摆动
所以具体这个数值是多少我也没查文献

看很多人说50 我觉得对 ...

我觉得截齐了做树的话，避免了不同长度的序列可能对做树造成的误差；但是针对一些蛋白基因，有些本身就很长，有些本身短，这种情况，我现在也不确定是不是也要截成一样长，有时，很多序列都很长，就1－2条比较短，这时候我也不知道该怎么办了，可以确定的是，截成一样长做树，做出来的树应该是不会有疑问，但是可能对序列的完整度有影响，得到的结果可能会有变化。
bootstrap值低于50确定是不可信的，但高于50只能说针对当前数据，在bootstrap方法检验下，是高于一半的支持率的，这个值当然越高越好，你选择不同数据，不同模型，bootstrap值一般都会有变化的，bootstrap值为100%也只能说是两个序列聚在一起的概率非常大，是个推测值，也不能绝对说他们就是一类的，我个人觉得就是一个统计分析结果，低于50%一般是不可信，高于50%可信度就高一点而已……
赞一下 回复此楼
13楼2012-11-15 12:05:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 19 个回答

xp198766

铁杆木虫 (著名写手)

小木虫职业打酱油滴~~！

【答案】应助回帖

★ ★ ★
秋电工(金币+6): 很感谢，解答了我大部分问题，但是方法中有些具体的操作我还是不会，希望能更详尽 2011-04-18 11:08:26
scelab(金币+3, EPI+1): good 2011-04-18 11:54:10
1、clustalx 也是为了alignment，如果在MEGA里面alignment了，就不用clustalx了；

2、bootstrap值一般是将你的序列保留一部分，把剩下部分随机打乱，拼成不同的序列，组成1000个你的alignment文件，做树，显示的77表明，1000次做树的过程有有77%次，也就是770次得到HX6和AJ308316.1聚在一起这个结果；

3、下面的标尺是枝长，也是进化距离，bootstrap值是可信度；

4、这个我水平有限，也比较迷惑，个人觉得你现在做成的树适合的进化研究；

5、一般是要有bootstrap值的，这个表示可信度，要是值低于50,一般别人是不认同的，有的文章没有标尺，可能是只想得到树的分枝情况，没有想要计算各个物种的距离，只要得到的一个树的拓扑结果吧，但是bootstrap值一般是要的。

对了，你的alignment，排列好了之后，一般要将首尾序列截成一样齐……

只是个人意见，可能有偏差，有错误欢迎指出。
赞一下(8人) 回复此楼
3楼2011-04-18 10:52:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

秋电工

新虫 (小有名气)
引用回帖:
Originally posted by xp198766 at 2011-04-18 10:52:28:
1、clustalx 也是为了alignment，如果在MEGA里面alignment了，就不用clustalx了；

2、bootstrap值一般是将你的序列保留一部分，把剩下部分随机打乱，拼成不同的序列，组成1000个你的alignment文件，做树，显示 ...

“对了，你的alignment，排列好了之后，一般要将首尾序列截成一样齐……
”弱弱问问大侠，这个具体怎么操作？mega4.0里嵌合的clustal没有截序列的功能么？
赞一下 回复此楼
4楼2011-04-18 11:13:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xp198766

铁杆木虫 (著名写手)

小木虫职业打酱油滴~~！
引用回帖:
Originally posted by 秋电工 at 2011-04-18 11:13:05:
“对了，你的alignment，排列好了之后，一般要将首尾序列截成一样齐……
”弱弱问问大侠，这个具体怎么操作？mega4.0里嵌合的clustal没有截序列的功能么？

MEGA里面的不知道有没有，没有怎么用MEGA，好像clustalx里面也没有，不过bioedit里面有……可以先在BIOEDIT里面截，再导入到MEGA里面，再转换成MEGA格式，做树
赞一下 回复此楼
5楼2011-04-18 11:26:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 19 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭