1.vcf转ped/map:

plink --vcf spirits.recode.vcf --recode --out test

2.统计每个snp的MAF,并将maf < 0.05; miss > 0.2的SNP筛选出来并过滤掉。

plink --file 513lines_4429snps_110506 --freq --out test
plink --file 513lines_4429snps_110506 --maf 0.05 --geno 0.2 --out test --recode

3.针对508lines_traits.txt文件中的性状进行全基因组关联分析
用excel编辑获取fid,iid文件以及表型文件:
选择个体文件:
1 L1
2 L2
3 L3
4 L4
5 L5
6 L6
7 L7
8 L8
9 L9
10 L10

FID     IID     T1      T2      T3
1       L1      4.5514  4.0234  1.5055
2       L2      2.4888  2.2576  1.4878
3       L3      4.9849  3.413   1.1411
4       L4      7.9829  5.8398  1.31127
5       L5      1.7586  3.088   1.7824
6       L6      2.2273  4.0377  1.943
7       L7      8.1565  4.6783  0.6993
8       L8      2.8122  3.6162  1.33192
9       L9      2.3489  4.4685  1.9838

命令:

plink --bfile test --keep keep.txt --recode --out use_test
plink --allow-no-sex --file use_test --pheno phe.txt --assoc --adjust --all-pheno --out test

生成.qassoc.adjusted,.qassoc,分别是校正p值和原始p值

plink做SNP筛选和GWAS相关推荐

  1. 用plink做GWAS(PCA、关联分析)并用R绘图

    用plink做GWAS(PCA.关联分析)并用R绘图 GWAS 一.观察初始数据 二.质量控制 样本缺失率和位点缺失率过滤(产生.imiss和lmiss文件) 计算MAF 数据清理 三.检查亲缘关系 ...

  2. 使用Plink对SNP数据进行质量控制

    文章目录 前 中 对个体和SNP的检测和筛选 哈代平衡 次等位基因频率 性别检测 后 前 在做和基因相关的分析时,拿到基因型数据,首先要进行质量控制. 在使用Plink进行质量控制时,一般包含以下几步 ...

  3. 陈表达VBA笔记-Excel VBA 编程-在表格做自动筛选的功能

    Excel VBA 编程,在表格做自动筛选的功能 下面的例子是自动筛选85列的内容,筛选内容为处方 Private Sub CommandButton11_Click()     '筛选处方     ...

  4. Lasso回归系列一:用LASSO回归做特征筛选踩坑记

    用LASSO回归做特征筛选使用踩坑记 什么时候应该用LASSO回归进行特征筛选?使用时容易踩的坑你都避开了吗?Lasso回归的原理(请看我的另一篇博客: Lasso回归系列二:Lasso回归/岭回归的 ...

  5. Plink数据管理, 筛选个体和变异,合并不同的遗传文件,并附加一个表型

    0.2 数据管理 0.2.1 筛选个体和标记物(marker) 1)筛选个体 PLINK可用于筛选和排除某个个体,这需要一个含有个体IDs的文件.该文件必须是一个以空格/制表符分隔的文本文件,第一列是 ...

  6. javaweb和ajax使用查询出来的数据做下拉菜单_不会用Excel做数据筛选,老板叼的你没话说!...

    很多人会做漂亮的表格,自信满满地给老板汇报,老板突然冒出一句话:你给我看下公司员工的学历分布情况. 这个问题能瞬间让很多人傻眼,自然少不了老板的一顿叼.其实数据筛选并不难,只要认真学,10分钟轻松搞定 ...

  7. python运行excel高级筛选怎么用_懂点EXCEL就行!教你利用Python做数据筛选(上)...

    前言 Python的数据清洗功能有多厉害,相信不用我说大家都知道了,寥寥几行代码便可以把一份杂乱无章的表格给处理的干干净净.但是python也是不容易入门的,毕竟编程语言要理解和精通也是要花不少功夫的 ...

  8. python数据处理可以做什么菜_厌倦了EXCEL想玩点新花样?教你利用Python做数据筛选(下)...

    前言 回顾一下上次的案例模型,我们知道了只要写3行代码,就可以利用python对EXCEL数据做简单的筛选,说明了python是一门非常简单和实用的数据分析工具.前面学过的这几句简单的语言,不知道你们 ...

  9. python excel筛选统计_懂点EXCEL就行!教你利用Python做数据筛选(上)

    前言 Python的数据清洗功能有多厉害,相信不用我说大家都知道了,寥寥几行代码便可以把一份杂乱无章的表格给处理的干干净净.但是python也是不容易入门的,毕竟编程语言要理解和精通也是要花不少功夫的 ...

最新文章

  1. 《程序员代码面试指南》第三章 二叉树问题 二叉树节点间的最大距离问题
  2. 深度学习之优化详解:batch normalization
  3. c语言goto语句用法_C语言32个关键字9种控制语句34种运算符整理
  4. 数据结构—链表-单链表基本操作实现
  5. OpenShift 4 Hands-on Lab (10) 限制集群资源的使用量
  6. mysql then_mysql语法之case when then与列转行
  7. [leetcode] 1335. 工作计划的最低难度
  8. iOS 并发编程之 Operation Queues
  9. 从零开始学习Sencha Touch MVC应用之十三
  10. paip.python错误解决 0x64024e96 指令引用的 0x00000135 内存。该内存不能为 read。
  11. Mysql 8.0 | #08004Client does not support authentication protocol requested by server
  12. xposed android id修改,Xposed模块开发:修改设备IMEI
  13. 苹果妙控鼠标二代(Magic Mouse 2 )如何连接到 Window 10系统
  14. 100道经典AI面试题帮你搞定算法offer!
  15. 关键词搜索-关键词搜索引擎工具-关键词搜索排行榜
  16. 阿里云Landing Zone系列--场景说明
  17. JAVA 实现语音播报
  18. 任天堂服务器维护11.2,switch11.0.0版本更新内容汇总 NS11.0.0更新维护公告
  19. 阿里视频云web播放器常见问题汇总
  20. uniapp 搜索关键字高亮显示实现

热门文章

  1. 如何提升论文写作能力?
  2. unity拖拽UI生成Cube并跟随鼠标移动,放置点可放下
  3. 玩回合制手游《问道》心得
  4. 国产Google地球如何下载2021版谷歌高清卫星地图?
  5. 如何在谷歌浏览器中安装Dark Reader插件
  6. 超密集组网(UDN)
  7. win2008不能连接mysql_win2008 r2 安装sql server 2005/2008 无法连接服务器解决方法
  8. 思考:苟日新,日日新,又日新
  9. 一本图书的信息包括:图书编号(BookNo),书名(BookName),总册数(Total),借出数(Lend),价格(Value),利用结构体实现以下要求:
  10. AC-PEG-NH2,Acrylate-PEG-Amine,丙烯酸酯PEG氨基含有PEG间隔基