如何用Plink2在命令行中运行GWAS？ - 腾讯云开发者社区 - 腾讯云

开发者社区

文档建议反馈控制台

文章/答案/技术大牛

发布

bioinfo05-GWAS学习

参考： PLINK 1.9 (cog-genomics.org)[1] 前言发现plink2 和plink 差别还是挺大的，没什么plink2 教程，还是用老版。...map map，与ped文件相伴随的文件，主要包含ped文件中SNP的位置信息。一般包含4列。分别是1....每行一个SNP，顺序与ped文件中的SNP相对应。因为纯文本格式占用大量储存空间，实际操作中尽量使用二进制格式，一组ped/map文件可转换成一组bed/bim/fam文件。...bed+bim+fam bed 不同于在基因组比对时，使用的记录位置信息的bed 文件，这里为二进制格式，存储基因型，可以想象成ped文件中除去前6列，剩下基因型数据组成的矩阵。...plink2 --bfile HapMap_3_r3_1 --missing 会生成两个文件： $ ls plink2* plink2.log plink2.smiss plink2.vmiss

7872 0

在python中运行命令行命令的四种方案

本文由腾讯云+社区自动同步，原文地址 https://stackoverflow.club/article/run_shell_command_in_python/ 简介毫无疑问，使用python运行命令行是最方便的将模型测试自动化的途径...方案一：os.system 仅仅在一个子终端运行系统命令，而不能获取命令执行后的返回信息如果在命令行下执行，结果直接打印出来。

36.9K2 0

您找到你想要的搜索结果了吗？

是的

没有找到

笔记 | GWAS 操作流程1：下载数据

❝相关的文章：https://onlinelibrary.wiley.com/doi/full/10.1002/mpr.1608 ❞ 教程中包括数据的过滤，SNP的过滤，样本的过滤，质控的标准等等，介绍的非常清楚...下载数据和代码首先，在linux环境下，新建一个文件夹，进入后运行下面命令： git clone https://github.com/MareesAT/GWA_tutorial.git 下载之后，...R：https://www.r-project.org/ plink：http://zzz.bwh.harvard.edu/plink/ https://www.cog-genomics.org/plink2...其中，1_Main_script_QC_GWAS.txt，2_Main_script_MDS.txt为运行代码。...，3_Main_script_association_GWAS.txt为运行代码。

3.1K3 1

snpQT-又一个人基因组SNP填充和GWAS流程

这个流程的目的是让你的SNP cute，为处理人类基因变异提供了帮助：基因组版本转换（b37->b38或者反过来）样本质控人群分层填充前质控本地填充填充后质控 GWAS 使用自动化的nextflow...流程，我们在Singularity容器或 Anaconda 环境中运行一系列版本的生物信息学软件，以提高可靠性和可重复性。...样本有二分类/数量性状如果你符合以上几点，看看文档吧：https://snpqt.readthedocs.io/en/latest/snpQT可能对你没用，如果你想：原始序列质控 call变异家系GWAS...非人基因组数据引用：好像并没有发表在好的杂志上，康奈尔大学团队做的。...GPL3 PLINK2 2.00a2.3 Chang CC, Chow CC, Tellier LCAM, Vattikuti S, Purcell SM, Lee JJ (2015) Second-generation

8502 0

plink2.0和plink1.9的忧伤笔记

3，plink帮助文档可以通过官网查询具体参数：https://www.cog-genomics.org/plink/2.0/ 也可以在命令行中调出帮助文档：比如直接键入plink2，出现基础参数...plink2 --help [flag name(s)...]..."plink2 --help | more" describes all functions....plink2 --ped yuanshi.ped --map yuanshi.map 或者写为： plink2 --pedmap yuanshi 默认输出文件： plink2.log plink2...比如读取plink1.9的二进制文件，输出bgen格式： plink2 --bfile a1 --export bgen-1.1 --out t1 日志： $ plink2 --bfile a1 --export

3K1 0

【直播】我的基因组57：最简陋的祖源分析

FastPop软件： https://sourceforge.net/projects/fastpop/files/ cd ~/biosoft# https://www.cog-genomics.org/plink2...这个脚本比较考验shell能力，而且运行非常慢，因为千人基因组计划的数据太多了！...尤其是本次的祖源分析系列和后面的GWAS解读系列，我相信公司在这一块会做的更好，毕竟，这个可是他们的饭碗！...人类的Y染色体除了在端粒上的拟常染色体区的少部分片段（只占有染色体长度约5%）能与相应的X染色体发生重组，其外都不能发生重组。这些片区是由原本X染色体与Y染色体同源的片段遗留下来的。...Y染色体中不能发生重组的其他部分被称为“NRY区”（non-recombining region，非重组区）。这个区域中的单核苷酸多态性被用于父系祖先的追溯。

2.8K9 0

统计遗传学：第八章，基因型数据质控

PLINK软件入门命令行我们在第7章介绍了PLINK[1.2]，它是用于QC和GWAS分析的最流行的开源自由软件程序之一。...在这些示例中，我们展示了在光标后运行命令。当您进行更高级的分析时，您将从脚本运行PLINK，就像在R或其他程序中一样。...注意，这些命令总是以两个破折号开头 plink用户指南： 1 PLINK是一个命令行程序，因此它需要在该环境中运行（即DOS窗口或文本中定义的Unix终端）。...在终端下运行脚本因为我们只运行简单的命令来让读者熟悉如何使用PLINK，所以我们通常只在命令行中键入所有命令。然而，一旦你们进行了更高级的分析，你们就会想要编写你们的脚本，然后在终端中运行它。...通常，IBS>0.1875的成对个体从分析中删除。基于线性混合模型（如BOLT-LMM）的最新方法不要求我们删除相关个体，因为它们在计算中考虑了相关性。

2.6K1 0

多基因风险评分（PRS）分析教程

GWAS 结果 QC Base data 和 target data 都应执行严格的质量控制步骤。低 MAF 和 INFO 的 SNPs 在执行下游分析之前通常会被去除。...> Height.gz 等位基因不匹配在 base data 和 target data 中不匹配的 SNP 可通过“链翻转”进行匹配，例如某个 SNP 在 base data 中为 A/C，target...因此，我们建议在 R 中执行转换，或者用 PRS 软件直接执行转换。...Clumping 尽管原则上所有常见的 SNP 都可以用于 PRS 分析中，但习惯上先将 GWAS 结果 clump，然后再计算风险评分。...计算 PRS 在 plink 中我们可使用 --score 和 --q-score-range 计算 PRS。

18.2K4 4

如何使用GCTA分析考虑kinship的二分类性状：--fastGWA-mlm-binary

1，适合性状二分类性状，比如疾病，患病编码为1，疾病编码为2，比如是否发芽，是否抗病等等，这些性状在遗传研究中非常常见，传统的GWAS软件，一般对这类二分类性状都是用Logistic回归，比如plink...软件的实现方法（plink分析二分类Logistic的GWAS模型，表型值编码以及OR值意义）（plink分析二分类性状的GWAS），这类分析，只能矫正群体结构（通过PCA放到协变量中），不能考虑个体间的亲缘关系...3，如何用GCTA的--fastGWA-mlm-binary分析二分类性状该方法发表在2021年的NG上面，权威可靠。下面是具体分析流程，亲测有效。 ## 1....sp_grm --make-bK-sparse 0.05：保留K_ij>0.05的个体对，其余的删除；输出：sp_grm.sparse.bin，稀疏矩阵 ## 3. fastGWA-GLMM二分类GWAS...GWAS 的最优选择。

1741 0

BGEN格式如何使用？有经验的家长已经给孩子收藏了。。。

基于这些数据的文本表示的传统数据格式（如IMPUTE输出的GEN格式或变量调用格式）有时不太适合这些数据量。事实上，对于简单的程序，解析这些格式所花费的时间可以支配程序执行时间。...例如，下图显示了在1号染色体上121668个SNP的18496个样本的数据集中，列出各种常见格式（Y轴）、文件大小（X轴）的变体识别数据（即基因组位置、ID字段和等位基因）所需的时间。...对于PLINK二进制（.bid）文件，标识数据存储在单独的文件（.bim文件）中，因此时间实际上为零。对于基于文本的格式，文件压缩的使用和读取性能之间存在显著的权衡。...BGEN以334Mb存储了22.5亿个基因型的整个数据集，每个基因型略多于一位，在该测试中耗时1.5秒。...plink2 --bgen t1.bgen 'ref-last' --sample t1.sample --export ped --out x1 • --bgen文件：指定t1.bgen，后面跟着参数

1.6K1 0

BOLT-LMM用户手册笔记

2.3 运行要运行 bolt 可执行文件，只需在 Linux 命令行（在 BOLT-LMM 安装目录中）调用....第二步只需要适度的额外计算，不需要额外的RAM，因为它只是对BOLT-LMM在模型拟合期间计算的残余表型执行实值剂量SNP的基因组扫描（如GRAMMAR-Gamma [8[24]]）。...请注意，与混合模型关联不同，线性回归容易受到总体分层的影响，因此在执行线性回归时，您可能希望将主成分（使用其他软件计算，例如，EIGENSOFT v6.0+ 中的 PLINK2 或 FastPCA [18...如果您希望在命令行上同时查看此输出的同时保存它，则可以使用 ....在两个单独的 BOLT-LMM 运行中分析常染色体和 chrX 变异（使用两次运行中的所有常染色体和 chrX 类型变异作为模型拟合的 PLINK 输入）。

3.4K4 1

搞孟德尔随机化热点的小伙伴数据分析能力有点弱啊

月）NC文章《Genome-wide association analysis of plasma lipidome identifies 495 genetic associations》里面的数据在GWAS...文章数据在GWAS catalog ，里面的索引号是 GCST90277238-GCST90277416 ，所以我们查询 GCST90277238试试看：查询 GCST90277238试试看里面就链接到了...只需要一个简单的循环即可，前面的文章数据在GWAS catalog ，里面的索引号是 GCST90277238-GCST90277416 ，一个简单的 Shell 脚本示例，使用了 wget 命令和一个...你可以将上述代码保存到一个文件（比如 download_files.sh），然后通过 chmod +x download_files.sh 赋予执行权限，并运行 ....在Python中，你可以使用requests库来下载文件。

1.8K1 0

GWAS分析中SNP解释百分比PVE | 第三篇，MLM模型中如何计算PVE？

GWAS分析中SNP解释百分比PVE | 第三篇，MLM模型中如何计算PVE？ #2021.12.24 1. R语言计算的PVE能否用于MLM模型？...GAPIT3.0增加了MLM模型计算PVE的方法：下面是GAPIT论坛中，作者的回复，所以，我们可以在GAPIT软件中使用MLM模型，进行GWAS分析，并得到每个SNP的PVE值。...这里，我们进行比较一下两种PVE的计算结果：「GLM模型计算PVE结果」 > # GAPIT 中的glm模型 > a1 = fread("04_glm_gapit/GAPIT.GLM.V3.GWAS.Results.csv...所以，在MLM模型的GWAS中，我们要选择MLM方法计算的PVE。问题来了，如果不用GAPIT软件，该如何手动计算PVE值呢？ 4....其它GWAS分析软件如何计算PVE 我们知道，其它GWAS软件中是没有PVE的结果的，比如： GEMMA GCTA中的fast-GWA 下一节介绍一下如何用R语言进行演示MLM的PVE计算方法。

2.2K1 0

基因组分析工具的瑞士军刀—BEDtools

它堪称是基因组分析工具中的瑞士军刀。其设计灵活，可以轻松地与其他命令行工具集成，如 awk、grep、sort 等，使得它成为基因组研究和数据分析中不可或缺的工具之一。...主要编程语言：C++ 旨在 UNIX、LINUX 上的“命令行”环境中运行，因此不支持Windows 平台 2发表文章题目：BEDTools: a flexible suite of utilities...这个命令在基因组注释、变异位点分析等方面非常有用。如何找到两个或多个基因组数据集（例如BED文件）中重叠的区域 intersect图解 “A intersect B”展示了A和B之间的交集区域。...通过使用额外的参数（如“-wa -wb”），可以展示两个数据集B1和B2中与A数据集重叠的所有区域 ## 输出A和B有交集的区域 bedtools intersect -a ....第二行“merge I (-d 10)”展示了在合并时允许10bp距离内的区间合并成一个连续区间的结果。

2.9K1 0

一文搞定单倍型Haploview分析（万字详解）

配置java环境 https://www.java.com/zh-CN/download/ 下载安装好之后，在终端运行java，出现帮助文档，说明配置成功。...5 windows报错解决方案注意，有时候会有问题，windows系统安装Haploview总是错误，解决方案是直接用jar文件，在终端cmd中运行：有个小伙伴说，Haploview软件在windows...进阶：LDblock+GWAS+Annotation+Locuszoom 可以通过-TopSite在GWAS图中显示最显著位点与其它位点的LD关系。...六、批量绘制单倍型图有小伙伴在星球询问GWAS分析的显著性位点，如何批量绘制LDblock，之前写过LDblock的安装和使用说明：（LDblock绘制连锁不平衡和单体型图）问题有3个核心点： 1...，然后运行就行了。

1.5K1 0

孟德尔分析：代谢疾病相关的GWAS数据库

继上周分享了血液中的蛋白组学相关网站后➡【孟德尔随机化】血液循环中的蛋白质组：常用网站一网打尽，今天我们继续扩充孟德尔随机化GWAS数据的来源吧~ 今天主要分享与糖尿病、血糖代谢、肥胖、高血压等代谢综合征相关的数据源...中的数据“为我所用”~ 写在后面月底菜编要开题，所以代码部分暂时搁置了【之前发过的全流程代码（一）不知道有没有同学尝试过，让我康康~】，如果在运行过程中有任何问题欢迎后台留言或者评论区提问，我已经提前踩过坑了...，或许可以为大家提供一丢丢帮助可以看到，学孟德尔随机化之前要学gwas数据分析全基因组关联研究（GWAS）是一种研究方法，用于在大量人群中寻找基因变异与疾病或某种表型特征之间的关联。...基因分型：然后，使用基因分型技术（如SNP芯片）对每个样本进行全基因组的基因分型。质控：在进行关联分析之前，需要进行数据质控，包括移除低质量的SNP、移除亲缘关系较近的样本、移除性别不符的样本等。...关联分析：然后，使用统计方法（如logistic回归或线性回归）对每个SNP进行关联分析，以确定其与疾病或表型的关联程度。

3.8K3 0

PRS多基因评分教程学习笔记(一）

从我读的几篇文章来看，多基因风险评分分为两个派别，一个是从GWAS中挑选显著差异的snp，进行评分，另一个则是倾向于使用尽可能多的位点，比如几万甚至更多。...下面，先来看下整体的步骤：从图中也可以看出，PRS分析需要Base数据（GWAS统计数据如P值，基因型-表型的SNP关系等）和Target数据。...如果在计算PRS时做出了错误的假设，则PRS在目标数据中的作用将指向错误的方向。...重复的SNP 如果在基础数据的生成中发生了错误，则基础数据文件中可能存在重复的SNP。大多数PRS软件不允许在基本数据输入中重复SNP，因此应将其删除。...样本亲属关联在本教程中，目标数据是模拟的，因此在基础数据和目标数据中必须没有紧密相关的个人。但是，用户应确保将基础数据和目标数据之间紧密相关的个人的可能性降到最低。

3K1 0

FUMA：基因关联的功能图谱和注释

请注意，所选参考面板中不存在的变异将不会包含在任何分析中。输入文件必要的列：输入文件「必须」包括「P 值」和 hg19 参考基因组上的「rsID」或「染色体 + 基因位置」。...染色体列可以是字符串，如 "chr1"，也可以是整数，如 1。当输入文件包含 X 染色体时，将编码为 23 染色体，但输入文件也可以包含 "X"。...如果输入文件有其他名称，可在指定输入文件时在相应的输入框中输入。需要注意的是，应避免使用名称如上但元素不同的列。...输入的 GWAS 统计摘要文件可以是 SNPs 的子集（例如，只有您的研究中感兴趣的 SNPs），但在这种情况下，MAGMA 结果不再相关。...运行结束以后，可以获取以下信息：可以根据leading SNP进行后续分析…… GENE2FUNC 比较简单~ 遇到问题如何解决 [Check-list for troubleshooting errors

3.4K1 3

如何用血常规发 Nature，临床常见指标的深度挖掘

在正常生理变化如怀孕和更年期后，甲状腺功能减退症或肝炎等慢性疾病的发展后，以及脾切除术后，平均设定点会发生变化（扩展数据图3）。...e，25,254名MGB患者的HGB设定点GWAS的曼哈顿图。 f，g，在使用HGB设定点和单次门诊值进行的GWAS中，设定点的P值（f）更频繁地显著，效应大小（g；每个SNP的beta系数）相似。...SNP遗传力估计是使用LDAK67 v.5.2中的sum-hers函数从GWAS汇总统计中得出的，并使用LDAK-Thin标签文件（如LDAK网站上提供的那样），排除任何解释超过1%表型变异的预测因子。...每个GWAS的重要结果和位点列在补充表3-5中，完整的GWAS结果可以通过GWAS目录访问（数据和代码可用性声明中给出了访问编号）。 MGB生物库队列的总结特征见补充表1。...基因组数据通过Windows Linux子系统中的bcftools 1.21进行分析，并使用plink和plink2。

5570 0

如何用 GCTA软件进行PCA分析并进行可视化

大家好，我是邓飞，今天介绍一下如何用GCTA计算PCA，并使用R语言进行可视化。...二、GCTA 做 PCA 的核心优势遗传针对性强：直接基于遗传标记计算相似度，主成分能精准反映个体间遗传差异，可有效校正 GWAS 中的种群分层偏倚。...计算效率高：优化了 GRM 计算算法，支持大样本（数万个体）和高密度 SNP（百万级）数据，运行速度优于传统 PCA 工具（如 EIGENSOFT）。...兼容性好：与 PLINK 数据格式无缝衔接，输出结果可直接导入 R、Python 等软件进行可视化，且能与 GCTA 的其他功能（如遗传力估算）联动。...")) + ylab(paste0("PC2 ",round(pcaal$por[2],4),"%")) PCA分析拓展 1，PCA分析，可以根据分组，绘制置信区间分组PCA 2，PCA分析中，

5571 0

点击加载更多

相关资讯

热门标签

活动推荐

运营活动

活动名称

广告关闭