使用admixture进行群体结构分析「建议收藏」

使用admixture进行群体结构分析「建议收藏」使用admixture进行群体结构分析

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全家桶1年46,售后保障稳定

软件下载网址: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*.out

CV 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
#值越小越好

Jetbrains全家桶1年46,售后保障稳定

相关参考网站:

admixture 群体结构分析 – 简书

【一起学生信】群体结构图形—structure堆叠图_探索者v的博客-CSDN博客_structure群体结构

科学网—群体结构分析之Structure 篇完整版(By Raindy) – 高芳銮的博文

Admixture:一款快速分析群体遗传结构的软件_.ped

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/210766.html原文链接:https://javaforall.net

赞 (0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • JS 暂时性死区[通俗易懂]

    JS 暂时性死区[通俗易懂]JS暂时性死区ES6暂时性死区引用ES6暂时性死区只要块级作用域内存在let命令,它所声明的变量就“绑定”(binding)这个区域,不再受外部的影响。vartmp=123;if(true){tmp=’abc’;//ReferenceErrorlettmp;}上面代码中,存在全局变量tmp,但是块级作用域内let又声明了一个局部变量tmp,导致后…

    2022年6月30日
  • Redis的相关问题总结

    Redis的相关问题总结

    2021年11月7日
  • ios软件_iOS倒计时

    ios软件_iOS倒计时GitHub:https://github.com/remember17/CountdownTimer全局倒计时工具,可以维护任意多个倒计时当APP从后台进入前台的时候,倒计时不受影响,会根据时间线继续。使用方法: 把CountdownTimer.swift文件拖进项目 在CountdownTimer.swift的枚举CountDownKey中添加定时器 每个枚…

    2022年9月16日
  • Oracle修改表名「建议收藏」

    Oracle修改表名「建议收藏」rename原表名to新表名转载于:https://www.cnblogs.com/LeiYang5237/p/8549526.html

    2022年5月13日
  • ssm管理系统课题_p2实验室

    ssm管理系统课题_p2实验室开发目的方便高效地实验室设备统一管理,除了实现基本的增删改查,还提供借用、归还、购买和问题反馈功能,可实现对实验室设备的基本业务的处理本项目由本人负责开发完成,项目能保证正常运行,当然其中不免也会有缺漏或不完善的地方解决方案1.后端Java框架使用spring+springmvc+mybatisspring功能是实现参数参数注入,请求分发处理,对数据库操作进行事务控制,其中mybatis使用注解查询,整体上大部分使用xml配置,少部分使用注解2.前端使用HTML+javascript+css+j

    2022年10月13日
  • vim设置编码

    vim设置编码Vim设置编码让Vim能正确显示和保存各种编码文件内容是个头疼的事情,今天专注的把这个问题研究一下首先是Vim编码配置的三个选项:encoding(enc):此项目配置的是vim的默认显示编码,如果你打开的文件编码与该属性配置的编码不同,那么Vim会将编码自动转换到配置编码再显示fileencodings(fencs):此项目配置的是一个编

    2022年9月29日

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号