我下载了一个VCF文件+索引(https://gwas.mrcieu.ac.uk/datasets/ebi-a-GCST90001942/),我试图在R上打开它(见下面的脚本)。但是,R只会加载而不会打开我的文件。我将VCF解压缩为CSV文件,R可以读取。然而,现在我需要获得汇总统计数据(beta,SE,EA,非EA,EA频率。我正在做以下工作
install.packages("tidyverse")
install.packages("remotes")
install.packages("gwasrapidd")
library(tidyverse)
library(remotes)
library(gwasrapidd)
remotes::install_github("mrcieu/ieugwasr")
remotes::install_github("MRCIEU/TwoSampleMR")
remotes::install_github("MRCIEU/gwasvcf")
remotes::install_github("MRCIEU/gwasglue")
library(ieugwasr)
library(TwoSampleMR)
library(gwasvcf)
library(gwasglue)
使用GWAS VCF文件
library(gwasvcf)
set_bcftools('/path/to/bcftools')
set_plink('/path/to/plink')
remotes::install_github('mrcieu/genetics.binaRies', force = TRUE)
set_plink()
set_bcftools()
suppressWarnings(suppressPackageStartupMessages({
library(gwasvcf)
library(VariantAnnotation)
library(dplyr)
library(magrittr)
}))
set_bcftools()
阅读在一切
要读取整个数据集,请使用readVcf函数。
vcffile <- "/Users/path/Desktop/ebi-a-GCST90001942.vcf.gz"
vcf_data <- readVcf(vcffile)
R永远加载并且不打开...
1条答案
按热度按时间bfrts1fy1#
我还没有尝试过
gwasvcf
包,但确实尝试过转换MungeSumstats
R包提供的摘要vcf
文件,它工作得很好。安装
MungeSumstats
包后,您可以简单地该软件包自动将
vcf
转换为标准的汇总统计文件。最终的输出是你所期望的
希望对你有帮助。