Structure是由斯坦福大学Pritchard实验室开发的一款群体结构分析软件。structure可通过贝叶斯聚类方法对每个样本的来源进行判断,进而反应群体的遗传结构。真正跑时burnin设为100000,跑100000代,k设为1-20,每轮运行10次。
序列文件用xmfa2struct转为structure可识别的格式,SSR用Genalex转为structure可识别的格式。
操作步骤打开structure选择新建任务输入任务名,保存路径与输入文件输入个体数,倍性,位点数,missing data的比重(一般都是-9)输入除了数据外多余的行的内容输入除了数据外多余的列代表的内容,点击finish确认信息,点击proceed
点击parameter set下的New 输入burnin多少代和burnin后再跑多少代,可以是5000与50000,其他模型之类的按照默认,点击OK
输入配置名,点确定
点击start a job 设定k为1到10,number of iterations为20次
点击开始跑,等待跑完 找到结果的目录,找到results文件夹,压缩为zip文件 进入网站 http://taylor0.biology.ucla.edu/structureHarvester/# 选则results.zip,点击上传 等待跑完 查看最近的K值是多少 返回structure,点击具有最近k值的结果文件 点击bar plot下的show 结果出来了,图像可以以多种形式显示,并可保存。 可以标上居群ID 还可以根据Q值来划分 还可以分成多层来显示 此外,还有FST的图
直方图
三角图 树图等 另外,也可用结果中的文件利用distruct软件进行美化处理
好了,以上就是Structure教程。如果文章对你有所帮助,请转发给你身边需要的人噢!
没有服务器想快速查看高通量比对数据?Tablet可以帮你
FASTQ格式详解
如何用Pathview画出高大上的基因与代谢通路热图?
使用PAML计算选择压力
测序数据到手后,如何高效的进行质量统计?
2021-03-15
用DOGMA对叶绿体或线粒体基因组进行注释
扫码关注我们Get更多科研小工具- 本文固定链接: https://maimengkong.com/kyjc/1638.html
- 转载请注明: : 萌小白 2023年10月8日 于 卖萌控的博客 发表
- 百度已收录