GWAS软件包:GAPIT3它来啦

news2024/11/14 7:10:32

GAPIT是一款非常老的而且非常流行的软件包,傻瓜式操作,一键出图出结果,一篮子的解决方案,是我最经常使用的GWAS分析软件包。

最近,GAPIT现在的版本是GAPIT3,速度比第二版有较大的提升:

更大的变化,终于有GAPIT这个软件包了,可以用library载入进去,而且安装方式可以用github安装,更符合R-style。

1. GAPIT3官网地址

官网地址:
https://zzlab.net/GAPIT/

github地址:

https://github.com/jiabowang/GAPIT

2. windows安装

推荐安装方式:

devtools::install_github("jiabowang/GAPIT3",force=TRUE)
library(GAPIT)

安装完成:

发现了一个官方文档的bug:

这里应该是librfary(GAPIT)而不是library(GAPIT3,否则会报错。

安装后的软件包:

3. Linux安装

安装代码不变:

在这里插入代码片

安装成功后的截图:

安装后的系统查看:

4. 分析测试数据

在官网上面下载示例数据:

下载的数据:
在这里插入图片描述

解压后的数据:


将路径设置为测试数据的文件夹,运行测试数据:

测试代码:

library(GAPIT)

#Tutorial 1: Basic Scenario of Compressed MLM by Zhang and et. al. (Nature Genetics, 2010) 
#----------------------------------------------------------------------------------------
#Step 1: Set data directory and import files
myY  <- read.table("mdp_traits.txt", head = TRUE)
myG <- read.delim("mdp_genotype_test.hmp.txt", head = FALSE)

#Step 2: Run GAPIT
myGAPIT <- GAPIT(
  Y=myY,
  G=myG,
  PCA.total=3
)

Windows系统测试GAPIT


分析结果:

在这里插入图片描述

Linux系统测试:

运行日志:

在这里插入图片描述

运行结果:

5. 安装和测试GAPIT3代码汇总

## 安装GAPIT代码
# install.packages("devtools") # 如果没有devtools,先运行本行代码安装
devtools::install_github("jiabowang/GAPIT3",force=TRUE)
library(GAPIT)

## 下载示例数据后,测试分析是否正常
library(GAPIT)

#Tutorial 1: Basic Scenario of Compressed MLM by Zhang and et. al. (Nature Genetics, 2010) 
#----------------------------------------------------------------------------------------
#Step 1: Set data directory and import files
myY  <- read.table("mdp_traits.txt", head = TRUE)
myG <- read.delim("mdp_genotype_test.hmp.txt", head = FALSE)

#Step 2: Run GAPIT
myGAPIT <- GAPIT(
  Y=myY,
  G=myG,
  PCA.total=3
)

6. 示例数据百度网盘下载

官网的示例数据,貌似不太方便下载,这里将示例数据和运行代码放到公众号后台,回复:gapit,获得下载链接。

在这里插入图片描述

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/1133971.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

快速解决安装ps打开找不到MSVCP140.dll问题,教你5个解决方法,

如果你在安装 Photoshop 时遇到找不到MSVCP140.dll的问题&#xff0c;MSVCP140.dll是Microsoft Visual C 2015 Redistributable的一个组件&#xff0c;它提供了许多常用的C函数库&#xff0c;用于支持各种软件的正常运行。当安装或运行某些软件时&#xff0c;如果系统中MSVCP14…

鱼眼图像去畸变python / c++

#鱼眼模型参考链接 本文假设去畸变后的图像与原图大小一样大。由于去畸变后的图像符合针孔投影模型&#xff0c;因此不同的去畸变焦距得到不同的视场大小&#xff0c;且物体的分辨率也不同。可以见上图&#xff0c;当焦距缩小为一半时&#xff0c;相同大小的图像&#xff08;横…

【办公常用软件分享】

在平时的工作生活中&#xff0c;经常会遇到各种各样的需求&#xff0c;没有合适的工具&#xff0c;不仅会降低效率&#xff0c;还会影响结果&#xff0c;有些工具的功能虽然能够满足&#xff0c;但是需要付费&#xff0c;偶尔用一次总显得不划算&#xff0c;所以今天就分享几个…

Modbus转Profinet网关连接三菱变频器博图快速配置案例

本案例将分享如何使用兴达易控的modbus转profinet网关&#xff08;XD-MDPN100&#xff09;来连接西门子1200系列plc&#xff0c;并实现三菱变频器的485通讯兼容转modbusTCP通信。通过在博图中进行配置&#xff0c;我们可以实现设备之间的连接和通信。 首先&#xff0c;我们需…

通信基础(三):多路复用技术

一、时分复用 时分复用造成线路资源的浪费: 使用时分复用系统传送计算机数据时&#xff0c;由于计算机数据的突发性质&#xff0c;用户对分配到的子信道的利用率一般是不高的。 二、 统计时分复用 STDM(Statistic TDM)

select distinct 语句详解

select distinct &#xff1a;当我们期望返回的数据不存在重复数据时&#xff08;每一行的数据都不一样&#xff09; 例如&#xff1a;表a select distinct a.* from a -- 可以查询出所有的信息 select distinct a.id from a -- 可以查询出id不同的信息&#xff0c;则还是全部…

【STM32】GPIO控制LED(HAL库版)

STM32最新固件库v3.5/Libraries/CMSIS/CM3/DeviceSupport/ST/STM32F10x/system_stm32f10x.c 林何/STM32F103C8 - 码云 - 开源中国 (gitee.com) STM32最新固件库v3.5/Libraries/STM32F10x_StdPeriph_Driver/src/stm32f10x_gpio.c 林何/STM32F103C8 - 码云 - 开源中国 (gitee.…

基于V/F控制的三相逆变器MATLAB仿真模型

微❤关注“电气仔推送”获得资料&#xff08;专享优惠&#xff09; 参考文献&#xff1a;张飞,刘亚,张玉杰.基于V/F控制的三相逆变器仿真模型的研究[J].自动化与仪器仪表,2015 关于V/F控制的论文非常多&#xff0c;随意下载&#xff01; 当分布式电源经过逆变器运行于孤岛模…

【Java 进阶篇】Java Servlet 入门指南

Java Servlet 是一种用于构建Web应用程序的Java技术&#xff0c;它允许您处理HTTP请求和生成HTTP响应。本篇博客将向您详细介绍Servlet的入门知识&#xff0c;无论您是初学者还是有一定经验的开发者&#xff0c;都能受益匪浅。 什么是 Servlet&#xff1f; Servlet 是 Java 技…

18日草稿

AI视野今日CS.CV 计算机视觉论文速览 Tue, 17 Oct 2023 Totally 158 papers &#x1f449;上期速览✈更多精彩请移步主页 Daily Computer Vision Papers HairCLIPv2: Unifying Hair Editing via Proxy Feature Blending Authors Tianyi Wei, Dongdong Chen, Wenbo Zhou, Jing …

rocksdb db.h

编译 https://github.com/facebook/rocksdb/blob/main/INSTALL.md 如果您计划在生产中运行 RocksDB&#xff0c;请不要使用默认 make 或 make all 进行编译。这将以调试模式编译 RocksDB&#xff0c;这比发布模式慢得多。 默认存在centos7 build和run rocksdb的脚本 https://g…

Linux 函数调用的用户态与内核态

在用户态中&#xff0c;程序的执行往往是一个函数调用另一个函数。函数调用都是通过栈来进行的。 在进程的内存空间里面&#xff0c;栈是一个从高地址到低地址&#xff0c;往下增长的结构&#xff0c;也就是上面是栈底&#xff0c;下面是栈顶&#xff0c;入栈和出栈的操作都是…

适用于 Windows 10 和 Windows 11 设备的笔记本电脑管理软件

便携式计算机管理软件使 IT 管理员能够简化企业中使用的便携式计算机的部署和管理&#xff0c;当今大多数员工使用Windows 笔记本电脑作为他们的主要工作机器&#xff0c;他们确实已成为几乎每个组织不可或缺的一部分。由于与台式机相比&#xff0c;笔记本电脑足够便携&#xf…

pycharm使用运行Docker容器的python解释器

根据官方的介绍&#xff1a;pycharm使用docker中的python解释器&#xff0c;都是基于镜像来做 也就是说docker中的python解释器不能使用现有的docker容器&#xff0c;而是必须基于镜像重新构建专属的python环境 如果我们每个项目都需要添加一个相同的package就会导致需要重新构…

php危险函数及rce漏洞

php代码执行语句 eval() eval()语句 eval() 会将符合PHP 语法规范字符串当作php 代码执行。 <meta charset"UTF-8"> <pre><?php$dd$_REQUEST[dd];eval($dd);?>可以执行php代码 也可以套一层system执行系统操作指令 assert()函数 assert() …

基于利用协议模拟工具解决工控CTF题

概述 对于参赛者而言&#xff0c;工控CTF题目往往感觉很头疼&#xff0c;不知道如何下手&#xff0c;闲来之时&#xff0c;从网上看到一道协议分析的题目&#xff0c;想着用模拟工具试下&#xff0c;发现有意向不到的效果&#xff0c;本文中的小工具为开源工具&#xff0c;读者…

RCE 远程代码执行漏洞分析

RCE 漏洞 1.漏洞描述 Remote Command/Code Execute 远程命令执行/远程代码执行漏洞 这种漏洞通常出现在应用程序或操作系统中&#xff0c;攻击者可以通过利用漏洞注入恶意代码&#xff0c;并在受攻击的系统上执行任意命令。 2.漏洞场景 PHP 代码执行PHP 代码注入OS 命令执…

Android 系统架构

首语 由于工作内容的转变&#xff0c;使得我向Android系统方向转变&#xff0c;对于一个Android系统工程师&#xff0c;了解Android整个系统架构是必然的。本篇是Android系统学习的开篇&#xff0c;Android系统庞大且复杂&#xff0c;但是能对Android的认识更深&#xff0c;更…

测开 ( 项目篇 )

正文 - 项目实践 - 文件压缩 1.项目启动&#xff0c;介入了解需求 项目背景&#xff1a;磁盘空间不够&#xff0c;需要一个压缩、解压缩软件目标&#xff1a;完成所有文件类型的压缩、解压缩。时间&#xff1a;2018-6-1到2018-6-7项目成员&#xff1a;项目经理&#xff1a;唐僧…

如何做一个无符号数识别程序

1.状态分析 我们可以把无符号数分为&#xff1a;整数&#xff0c;带小数&#xff0c;带指数部分三种形式。以此构建一个DFA。首先需识别输入是整数还是小数点&#xff0c;若是整数部分输入然后还要再循环识别一次是否有小数点&#xff0c;最后识别是否有指数部分&#xff0c;指…