大家好,今天给大家推荐一款群体结构分析软件——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进行群体结构分析相关推荐

  1. 使用admixture进行群体结构分析

    软件下载网址:http://software.genetics.ucla.edu/admixture/download.html #使用vcftools进行格式转换,plink.ped与plin.ma ...

  2. 群体遗传学---admixture软件快速群体分群

    群体遗传学中测的很多个个体,得到了最终的SNP vcf文件,需要将其分成群体,看那几个物种聚在一起,一般使用的软件就是STRUCTURE,但是STREUTURE运行速度极慢,后面frappe软件提升了 ...

  3. plink质控及转换文件、admixture软件学习记录

    plink质控及转换文件.admixture软件学习记录 目的:自定义ped文件和map文件,用plink质控及转换文件,使用admixture做群体结构学分析. ped文件和map文件根据格式用pa ...

  4. 心理测评软件php mysql_心理测量软件_心理测评系统

    心理测量软件_心理测评系统 时间:2020-08-17 16:28:24    浏览:次 你也许会说,心理健康与否,只不过是个人的小事,和别人.社会没有什么牵连.你说错了.比方说,如果一个人手脚不灵便 ...

  5. animiz动画制作软件_学习动画制作需要掌握的常用软件合集

    影视动画越来越深入我们的工作和生活,其中动画的发展也深深吸引着我们的眼球. 在当下创作一部动画,大多已放下传统纸上动画的形式,开始尝试使用电脑绘图和3D方式来制作.所以,就让我们来看看常用的动画制作的 ...

  6. 录制电脑屏幕的软件_电脑屏幕录制操作方法,太实用了!

    电脑屏幕录制操作应该怎么做?本期将和大家分享三种电脑屏幕录制的操作方法,适用于win7.win10和mac三种电脑系统,包括windows7系统自带的步骤记录器和windows10操作系统自带的屏幕录 ...

  7. 上网时间监控软件_电脑监控系统软件有哪些作用?主要功能有哪些?

    随着IT技术和互联网的发展,企业开始利用网络办公来降低管理成本.信息化办公给企业带来了高效快捷的同时,机密泄露事件也层出不穷,网络安全问题变得日益严重.正是在这种背景下,以员工计算机上网行为监管为主的 ...

  8. 工程计算软件_同望BIM工程量计算软件—土石方

         天工造价关注我们,给您想要的工程类资讯 复制以下链接,我要参与 https://www.tgcost.com/member/activity/20190705/index.html 推荐阅读 ...

  9. 软件_手把手教vscode配置c++,python开发环境

    原创:软件_手把手教vscode配置c++,python开发环境 之前主用Python作为项目开发语言,将项目迁移到arm边缘盒子上后发现arm的cpu不给力,软件速度低于预期,所以计划将部分程序改为 ...

  10. 软件_搭建rtsp视频推送环境

    原创博客地址:软件_搭建rtsp视频推送环境 live555编译安装启动 编译 1 2 3 4 5 wget http://www.live555.com/liveMedia/public/live5 ...

最新文章

  1. Python运算符与编码
  2. html和jsp的差异
  3. 报错,null [java.lang.IndexOutOfBoundsException,Index: 5, Size: 5]
  4. application实现网页计数_手把手教你利用爬虫爬网页(Python代码)
  5. Python解决模块无法引用问题之pymysql模块导入
  6. Apache虚拟机和访问控制
  7. 关于ES6 Class语法相关总结
  8. python中高阶函数与装饰器(3)
  9. c#【 网络日志解析工具】开发篇-【媒体预览功能】--base64,webp响应解码为图片-》可做复制base64进行图片查看
  10. mysql sql 分析工具下载_DB Query Analyzer下载
  11. android刷机教程 华为,华为Mate20X 刷机教程 华为Mate20X 强刷升级教程
  12. Android Studio生成APP方法及其所在位置
  13. 《软件工具》分享2款好用的时序图工具
  14. 素数之魂——黎曼和他的伟大猜想
  15. 李密(582~619)
  16. 全面提高你的搜索技巧
  17. java 各种编码间转换
  18. 数字IC后端实现40天速成篇(中)
  19. 如何给自己的公司做一个网站
  20. 查询电话联系人(包括sim卡)

热门文章

  1. 【CSON原创】CSS的障眼法:利用border实现图片的翻转
  2. jQuery jqModal弹出层
  3. 独热向量编码(one-hot encoding)原理详解与实现
  4. 10、斐波那契数列,跳台阶问题(Python)
  5. 【JavaEE】CH10 JDBC
  6. java生成16位唯一性的订单号
  7. exawear能运行java_VirSCAN.org-多引擎在线病毒扫描网 v1.02,当前支持 47 款杀毒引擎...
  8. python下载完是什么样的图片_[宅男福利]用Python下载页面中所有的图片
  9. 常平计算机培训班,常平大朗CNC编程培训速成班,一个月学会UG编程
  10. Home Assistant系列 -- 设置界面语言与地理位置