powermarker的问题讨论
(2010-11-01 14:52:22)
标签:
杂谈 |
关于第一个问题的解答:生态类群间遗传距离是powermarker计算方法如下:
line | group | tt324 | tt24 | tt665 |
1 | 1 | 243/243 | 261/261 | 306/306 |
2 | 2 | 260/260 | 276/276 | 306/306 |
3 | 1 | 260/260 | 261/261 | 306/306 |
4 | 3 | 243/243 | 276/276 | 306/306 |
5 | 4 | 260/260 | 291/291 | 306/306 |
2.数据导入成功后,点击窗口顶部“Analysis”,在phylogene后拉框中选择compute frequence,在弹出的对话框中点击数据文件名(比如为:”*”)后,注意选择右侧相应的计算水平level-2,点击submit.
3. 点击窗口顶部“Analysis”,在Phylogene后拉框中选择“Frequence Based Distance”,在弹出的对话框中点击数据文件名(*)后,注意选择右侧距离计算方法,注意选择右侧相应的计算水平level-2,最后点击submit,此时distance文件夹里会出现“*.share allele”的类群间距离表。
4. 点击窗口顶部“Analysis”,在phylogene后拉框中选择“UPGMA/ NJ tree”,在弹出的对话框中点击数据文件名“*.share allele”后,注意选择右侧聚类模型,最后点击submit。
5. 点击窗口顶部“Data”,选择Batch export, 弹出的对话框中点击以NJ或upgma结尾的数据文件名,选择输出文件保存路径,再点击submit.
6.用Figtree软件打开输出树形图文件进行编辑,即可得到类群间聚类图,
文老师:
1、您的博士论文里的55页的对生态类群聚类时,生态类群间的 遗传距离是怎么算出来的呢,我用powermarker但是不知 道怎么算出来,我只能算出来生态类群见的Fst,能算AVOMA 。
3、对于TASSEL,我通过GLM模型算出来的结果有贡献率 ,但是MLM模型却找不到贡献率这一项,说明书看了再看,也还是 没有找到,如果不用MLM也就是没考虑kinship了。
4、我的群体是2个种,不是亚种,对于性状与标记的关联,作为 整个群体和分别作为2个群体来算,结果差异确实很大,从贡献率上 来开最明显,对于一个整体,对于相同的标记和性状,贡献率为5% ,分为2个群体来计算时,在一个群体中贡献率达到20%多(这2 个品种的性状差异比较大)