使用admixture进行群体结构分析
软件下载网址:http://software.genetics.ucla.edu/admixture/download.html
#使用vcftools进行格式转换,plink.ped与plin.map文件
vcftools --vcf maftest.vcf --plink --out adsaixuan#过滤SNP文件
plink --noweb --file adsaixuan --hwe 0.0001 --make-bed --out QC#群体结构分析,设置K(如1-7)
for K in 1 2 3 4 5 6 7; do admixture_linux-1.3.0/admixture --cv QC.bed $K | tee log${K}.out; done#提取出CV值
grep -h CV log*.outCV error (K=1): 0.55248
CV error (K=2): 0.48190
CV error (K=3): 0.47835
CV error (K=4): 0.48236
CV error (K=5): 0.48985
#值越小越好
相关参考网站:
admixture 群体结构分析 - 简书
【一起学生信】群体结构图形—structure堆叠图_探索者v的博客-CSDN博客_structure群体结构
科学网—群体结构分析之Structure 篇完整版(By Raindy) - 高芳銮的博文
Admixture:一款快速分析群体遗传结构的软件_.ped
使用admixture进行群体结构分析相关推荐
- 群体进化,群体结构分析之STRUCTURE分析三款软件比较
STRUCTURE分析三款软件比较 三篇高引用文章 2005-STRUCTURE[1] 把选k值写的很清楚 2020/5/16 引用13119 这篇文章发表的时候二代测序还没兴起,ssr等的标记数量有 ...
- 基于SSR数据的群体结构分析和DeltaK值可视化
1.ssr做群体结构数据准备 1.1数据转换. 1.2原始数据是不符合structure的格式要求的,分享给大家我的python脚本做数据格式转换. $cat dat2structure.py #!/ ...
- 群体结构分析:用 phylip 构建进化树
用 phylip 构N-J树 在 linux 系统,可以用 conda 安装 phylip 软件 1. 将 SNP 文件转换为 phylip 格式 用 tassel 的格式转换功能将 plink 格式 ...
- GWAS全基因组关联分析流程(BWA+samtools+gatk+Plink+Admixture+Tassel)
我梳理了GWAS全基因组关联分析的整个流程,并提供了基本的命令,用到的软件包括BWA.samtools.gatk.Plink.Admixture.Tassel等,在此分享出来给大家提供参考. 一.BW ...
- 二代测序群体进化-GWAS分析及案例解析
群体进化-gwas分析 群体进化基础分析 PCA 分析原理 PCA(Principal Component Analysis),即主成分分析方法,是一种使用最广泛的数据降维算法.PCA的主要思想是将n ...
- 统计遗传学:第三章,群体遗传
3. 群体遗传 大家好,我是飞哥. 前几天推荐了这本书,可以领取pdf和配套数据代码.这里,我将各个章节介绍一下,总结也是学习的过程. 引文部分是原书的谷歌翻译,正文部分是我的理解. 第一部分基础,分 ...
- 【GWAS文献】基于GWAS与群体进化分析挖掘大豆相关基因
Resequencing 302 wild and cultivated accessions identifies genes related to domestication and improv ...
- 群体进化技术原理和研究思路
学习链接 https://www.bilibili.com/video/BV11X4y1G7ZM/?spm_id_from=333.788.recommend_more_video.10 不同生存条件 ...
- Nat Genet | 湖南农业大学群体基因组学研究揭示古老芥菜身世之谜、百变之源
2021年9月6日,国际著名学术期刊**Nature Genetics在线发表了湖南农业大学****刘忠松教授团队联合中国农业科学院油料作物研究所华玮研究员团队等国内外多家单位完成的题为"G ...
最新文章
- AI 初学者必须要了解的术语盘点
- BeanShell脚本接口之匿名类
- GT Transceiver的回环模式
- [architecture]-ARMV8的RAS Extension(Reliability、Availability、Serviceability)介绍
- dict是python语言的内置对象_Python内置了字典:dict的支持
- 自编码 Autoencoder
- linux mysql恢复数据_删库不跑路详解MySQL数据恢复
- Ps 初学者教程,如何使用色阶功能提高照片的对比度和亮度?
- iptables 简介
- 帷幕的帷是什么意思_美术生每天坚持练习|但越画越差的原因是什么?
- 基于单片机和C语言的毕业设计,毕业论文基于51单片机的C语言程序设计实训100例(1)(喜欢就下吧)...
- [转载]Altium Designer 破解,避免局域网冲突的ad9.ini
- 惠普笔记本安装Linux 7及NVIDIA Quandro P600显卡驱动
- Linux磁盘分区总结
- threejs-自定义着色器材质
- 知识图谱及其关键技术
- 淘宝店图片尺寸是多少?怎么修改淘宝图片的尺寸?
- IntelliJ IDEA如何整理代码格式
- 如何在 Amazon EKS 中部署 SR-IOV 设备插件
- 显卡超频稳定测试软件,显卡超频稳定性测试终极手段
热门文章
- 【Flutter web】内网网站如何发布?解决外网下canvaskit.js和字体无法加载问题
- 淘宝直播的“一哥”和“一姐”:当年有多卑微,今天就有多荣耀!
- 在2022年想创业就一定要注册公司吗?
- 使用docker搭建elk
- mysql 字段移动平均值_MYSQL简单移动平均值计算
- 从微信「拍一拍」,我想到了那些神奇的一行代码功能(文末福利)
- curl linux 测试接口_Linux测试端口的连通性的四种方法
- 在physx中创建射线ray,并获取碰撞点
- oracle 块大小,如何查看Windows上文件系统块的大小
- 一大波开发者福利来了,一份微软官方Github上发布的开源项目清单等你签收