admixture软件_使用ADMIXTURE进行群体结构分析
大家好,今天给大家推荐一款群体结构分析软件——AMIXTURE,以前大家经常用STRUCTURE软件进行群体结构分析,但因STRUCTURE运行速度极慢,而被不少研究者所诟病。ADMIXTURE使用与STRUCTURE相同的统计模型,但使用快速数值优化算法计算估计值,因此大大提升其运行速度,逐渐成为群体结构分析的主流软件。
接下来,给大家分享一下如何使用这款软件进行分析。
第一步:必须是下载该软件了
wget http://dalexander.github.io/admixture/binaries/admixture_linux-1.3.0.tar.gz
tar -zvxf admixture_linux-1.3.0.tar.gz
cd admixture_linux-1.3.0
./admixture --version
显示以上信息,证明软件安装成功,接下来就是准备输入文件了。
第二步:准备输入文件
ADMIXTURE输入文件为二进制的PLINK文件(.bed),同时将对应的bim和fam文件放在同一个文件夹下,对于我们得到到非bed格式文件,使用Plink1.9软件进行转化,举例,我们现在有example.gen文件,使用以下命令,即可得到example.bed, example.bim, example.fam文件
% plink1.9 --file example.gen --make-bed --out example
第三步:使用ADMIXTURE进行计算
% for K in 1 2 3 4 5; do admixture --cv example.bed $K | tee log${K}.out; done
第四步:选择交叉熵最小值(best K)就是最佳的分群
% grep -h CV log*.out
CV error (K=1): 0.94522
CV error (K=2): 0.60260
CV error (K=3): 0.58336
CV error (K=4): 0.59707
CV error (K=5): 0.59939
该组最佳分群数,k=3
第五步:使用R做图
输出结果为example.3.P和example.3.Q两个文件,我们使用example.3.Q做图
>tbl=read.table("example.3.Q")
>barplot(t(as.matrix(tbl)), col=rainbow(3),
xlab="", ylab="Ancestry",border = NA)
你学会了吗?
欢迎大家关注公众号,并提出宝贵意见。
参考:
http://dalexander.github.io/admixture/admixture-manual.pdf
admixture软件_使用ADMIXTURE进行群体结构分析相关推荐
- 使用admixture进行群体结构分析
软件下载网址:http://software.genetics.ucla.edu/admixture/download.html #使用vcftools进行格式转换,plink.ped与plin.ma ...
- 群体遗传学---admixture软件快速群体分群
群体遗传学中测的很多个个体,得到了最终的SNP vcf文件,需要将其分成群体,看那几个物种聚在一起,一般使用的软件就是STRUCTURE,但是STREUTURE运行速度极慢,后面frappe软件提升了 ...
- plink质控及转换文件、admixture软件学习记录
plink质控及转换文件.admixture软件学习记录 目的:自定义ped文件和map文件,用plink质控及转换文件,使用admixture做群体结构学分析. ped文件和map文件根据格式用pa ...
- 心理测评软件php mysql_心理测量软件_心理测评系统
心理测量软件_心理测评系统 时间:2020-08-17 16:28:24 浏览:次 你也许会说,心理健康与否,只不过是个人的小事,和别人.社会没有什么牵连.你说错了.比方说,如果一个人手脚不灵便 ...
- animiz动画制作软件_学习动画制作需要掌握的常用软件合集
影视动画越来越深入我们的工作和生活,其中动画的发展也深深吸引着我们的眼球. 在当下创作一部动画,大多已放下传统纸上动画的形式,开始尝试使用电脑绘图和3D方式来制作.所以,就让我们来看看常用的动画制作的 ...
- 录制电脑屏幕的软件_电脑屏幕录制操作方法,太实用了!
电脑屏幕录制操作应该怎么做?本期将和大家分享三种电脑屏幕录制的操作方法,适用于win7.win10和mac三种电脑系统,包括windows7系统自带的步骤记录器和windows10操作系统自带的屏幕录 ...
- 上网时间监控软件_电脑监控系统软件有哪些作用?主要功能有哪些?
随着IT技术和互联网的发展,企业开始利用网络办公来降低管理成本.信息化办公给企业带来了高效快捷的同时,机密泄露事件也层出不穷,网络安全问题变得日益严重.正是在这种背景下,以员工计算机上网行为监管为主的 ...
- 工程计算软件_同望BIM工程量计算软件—土石方
天工造价关注我们,给您想要的工程类资讯 复制以下链接,我要参与 https://www.tgcost.com/member/activity/20190705/index.html 推荐阅读 ...
- 软件_手把手教vscode配置c++,python开发环境
原创:软件_手把手教vscode配置c++,python开发环境 之前主用Python作为项目开发语言,将项目迁移到arm边缘盒子上后发现arm的cpu不给力,软件速度低于预期,所以计划将部分程序改为 ...
- 软件_搭建rtsp视频推送环境
原创博客地址:软件_搭建rtsp视频推送环境 live555编译安装启动 编译 1 2 3 4 5 wget http://www.live555.com/liveMedia/public/live5 ...
最新文章
- Python运算符与编码
- html和jsp的差异
- 报错,null [java.lang.IndexOutOfBoundsException,Index: 5, Size: 5]
- application实现网页计数_手把手教你利用爬虫爬网页(Python代码)
- Python解决模块无法引用问题之pymysql模块导入
- Apache虚拟机和访问控制
- 关于ES6 Class语法相关总结
- python中高阶函数与装饰器(3)
- c#【 网络日志解析工具】开发篇-【媒体预览功能】--base64,webp响应解码为图片-》可做复制base64进行图片查看
- mysql sql 分析工具下载_DB Query Analyzer下载
- android刷机教程 华为,华为Mate20X 刷机教程 华为Mate20X 强刷升级教程
- Android Studio生成APP方法及其所在位置
- 《软件工具》分享2款好用的时序图工具
- 素数之魂——黎曼和他的伟大猜想
- 李密(582~619)
- 全面提高你的搜索技巧
- java 各种编码间转换
- 数字IC后端实现40天速成篇(中)
- 如何给自己的公司做一个网站
- 查询电话联系人(包括sim卡)
热门文章
- 【CSON原创】CSS的障眼法:利用border实现图片的翻转
- jQuery jqModal弹出层
- 独热向量编码(one-hot encoding)原理详解与实现
- 10、斐波那契数列,跳台阶问题(Python)
- 【JavaEE】CH10 JDBC
- java生成16位唯一性的订单号
- exawear能运行java_VirSCAN.org-多引擎在线病毒扫描网 v1.02,当前支持 47 款杀毒引擎...
- python下载完是什么样的图片_[宅男福利]用Python下载页面中所有的图片
- 常平计算机培训班,常平大朗CNC编程培训速成班,一个月学会UG编程
- Home Assistant系列 -- 设置界面语言与地理位置