GWAMA:全球荟萃分析工具
下载和安装
从官网(https://genomics.ut.ee/en/tools/gwama-download)上下载或:
wget http:/www.geenivaramu.ee/tools/sample.zip
unzip gwama.zip (解压)
make (安装)
运行命令:
../GWAMA --filelist gwama.in --output gwama_xp_wz_jp --quantitative
其中,
gwama.in为输入文件,内容如下
每一行对应一个人群的gwas分析结果,需要的文件格式如下:
MARKERNAME
代表SNP位点的名称,EA
代表effect allele, NEA
代表other allele,EAF代表EA的频率, 对于连续型的性状,要求输入BETA
和SE
两列,对于二分类的性状,要求输入OR
, OR_95L
, OR_95U
3列,代表odd ratio以及95%置信区间。
以上是最基础的数据,如果有额外的数据也可以加入,比如STRAND
代表链的方向,IMPUTED
代表该位点是否是填充得到的等
gwama_xp_wz_jp为输出文件的前缀,输出结果保存在gwama_xp_wz_jp.out文件中,内容如下
--quantitative 表示使用的表型数据是连续变量。
其它需要请看官方文档:https://genomics.ut.ee/en/tools/gwama-tutorial
为什么需要做meta分析?
群体分层是GWAS研究中一个比较常见的假阳性来源。
具体来说,GWAS的核心在于SNP与表型的关联,具体来说,我们观测到每一个genome位点会有许多不同的SNP,如果某个SNP总是与某种表型/疾病同时出现,那么我们可以推测这个SNP极有可能会导致这个表型/疾病。如果存在群体分层,那么发现的SNP可能是群体的一个特征而不是与表型相关,也就是假阳性位点。
当群体出现分层时,常规手段时将分层的群体独立分析,最后做meta分析。
如何判断群体分层?
用plink计算PCA,画PC1和PC2在不同群体的散点图,看是否将人群分开,如果能观察到几个group,就是出现群体分层。
什么是meta分析?
GWAS meta analysis应该翻译成整合分析,略区别于医学的荟萃分析。GWAS中的meta analysis,是将不同的GWAS整合在一起,减少群体分层的影响,提升了power。(还是不知道原理)
一般用的软件有:plink、GWAMA、METAL
关于Meta的在不同场景下的意思
Meta的字面意思:Meta means about the thing itself. It's seeing the thing from a higher perspective instead of from within the thing, like being self-aware
Metadata is "data that provides information about other data"
在计算机里,meta就是描述数据的数据。
META是一种思想概念,一种抽象思维,用来描述数据的数据,比如有一张学生表,记录着学生的基本信息,我们通过表可以获取学生信息(数据),但是有时候也要得到表本身的信息数据(比如表结构信息:字段名称,字段数据类型,长度等信息),对于这种基础信息的描述,就会使用META的概念,使用META元数据来描述表本身。放到HTML中也是一样的,HTML用来描述网页信息,但是HTML自己也有一些信息(比如网页标题,网页描述,搜索关键字),这些信息也就称之为HTML META信息,并且HTML也定义了专门的META标签。
在医学里,meta就是荟萃分析。
Meta分析是对具有相同目的且相互独立的多个研究结果进行系统的综合评价和定量分析的一种研究方法。它不仅仅是搜集很多以往的研究成果,同时进行全面系统的质量评价,另外还要对符合条件的研究进行定量分析。
在微生物里,宏基因组学(Metagenomics),又称元基因组学。是以特定生境中的整个微生物群落作为研究对象,无需分离培养,直接提取环境样本DNA进行测序,研究环境微生物的群落结构、物种分类、系统进化、基因功能及代谢网络等,已广泛应用于微生物领域。
Reference:
https://genomics.ut.ee/en/tools/gwama
https://www.cnblogs.com/chenwenyan/p/10912521.html
https://www.cnblogs.com/leezx/p/13369846.html
推荐阅读
-
35 岁实现财务*,腾讯程序员手握2300万提前退休?-1000万房产、1000万腾讯股票、加上300万的现金,一共2300万的财产。有网友算了一笔账,假设1000万的房产用于自住,剩下1300万资产按照平均税后20-50万不等进行计算,大约花上26-60年左右的时间才能赚到这笔钱。也就是说,普通人可能奋斗一辈子,才能赚到这笔钱。在很多人还在为中年危机而惶惶不可终日的时候,有的人的35岁,就已经安全着陆,试问哪个打工人不羡慕?但问题是有这样财富积累必然有像样的实力做靠山。没有人可以不劳而获。 看到这里,肯定有人说,那么对于普通人来说,卷可能真就成了唯一的出路。但是卷也有轻松的卷,“偷懒”的卷法,对于程序员而言,刨除掉一时无法改掉的开会传统占用的大部分时间,如何把有限的时间和精力放在真正重要的架构设计、需求设计上,而不是重复的造*,编码、改bug、手动测试。因此在科技改变生活的今天,学会使用AI工具成为程序员们的必备技能。 以全栈式全自动的软件开发工具飞算SoFlu软件机器人为例,作为全球首款面向微服务架构设计和最佳实践的软件机器人,SoFlu软件机器人改变了原来手工编码的作业模式,通过可视化拖拽方式以及参数配置就能实现等同于编写复杂代码的业务逻辑,在设计业务逻辑时就完成了微服务应用开发,做到“业务即图,图即代码”,实现“软件开发,十倍提效”。 而作为飞算SoFlu软件机器人的一个重要组成部分,FuncGPT(慧函数)支持所有类型函数创建。通过自然语言描述Java函数需求,实时生成高质量、高可读性的Java函数代码。生成代码可直接复制到IDEA,或一键导入Java全自动开发工具函数库。 FuncGPT(慧函数)具有五大核心能力: ● 自然语言:通过自然语言即可生成函数,降低软件开发门槛。 ● 秒级函数开发:全自动秒级开发函数,秒级完成,效率千倍提升。 ● 一站式函数开发应用:函数生成即可直接应用,效率、安全有保障。 ● 质量可靠:生成的函数代码符合业界规范、具备良好的可读性,符合最优解。 ● 灵活生成:可以根据开发人员的具体需求生成并支持调整和修改。 以“将文件流转换为File对象”这个需求为例,以下是FuncGPT(慧函数)秒级生成的代码: / 类名:FileConverter // 函数名:convertStreamToFile // 函数功能:将文件流转换为File对象 // POM依赖包:无 import java.io.File; import java.io.FileOutputStream; import java.io.IOException; import java.io.InputStream; public class FileConverter { /** * 将文件流转换为File对象 * @param inputStream 文件流 * @param filePath 文件路径 * @param fileName 文件名 * @return 转换后的File对象 * @throws IOException 如果发生I/O异常 */ public static File convertStreamToFile(InputStream inputStream, String filePath, String fileName) throws IOException { File file = new File(filePath + File.separator + fileName); // 创建File对象 try (FileOutputStream outputStream = new FileOutputStream(file)) { // 创建文件输出流 byte buffer = new byte[1024]; int bytesRead; while ((bytesRead = inputStream.read(buffer)) != -1) { // 从文件流读取数据并写入文件 outputStream.write(buffer, 0, bytesRead); } } return file; // 返回转换后的File对象 } } // 函数示例 // 将文件流转换为File对象示例 // 入参:inputStream,文件流 // 入参:filePath,文件路径 // 入参:fileName,文件名 // 出参:file,转换后的File对象 // 调用示例: // InputStream inputStream = new FileInputStream("example.txt"); // String filePath = "C:\\Users\\User\\Documents"; // String fileName = "example.txt"; // File file = FileConverter.convertStreamToFile(inputStream, filePath, fileName); // System.out.println(file.getAbsolutePath); // 输出结果:例如,将文件流转换为File对象后,文件的绝对路径为:C:\Users\User\Documents\example.txt // 则输出结果为:C:\Users\User\Documents\example.txt 通过分析,不难发现以上代码:
-
GWAMA:全球荟萃分析工具
-
GWAMA:强大的基因组关联研究荟萃分析工具