求救各位大佬,用之前在网上copy的代码解读vcf文件时出现了问题,原代码是用在GWAS 上下载的数据,数据格式直接就是vcf文件,但是芬兰数据库下载的文件没有后缀名,用VariantAnnotation包的readvcf函数读根本就读不出来,代码:

library(VariantAnnotation)

outcomeFile <- fread(input="finngen_R9_G6_MYASTHENIA.gz",sep="\t",header=T)

vcfRT <- readVcf(outcomeFile)

报错内容: 

函数‘readVcf’标签‘"data.table", "missing"’找不到继承方法

Logo

腾讯云面向开发者汇聚海量精品云计算使用和开发经验,营造开放的云计算技术生态圈。

更多推荐