将vcf转换为hap文件，折叠基因型

Question

我用照片问一个问题。

• 我有一个vcf file。
• 还有一个R-code。
• 这就是我的输出：output

我希望我的输出文件仅包含数字/行/样本。

如何删除newdata ... C ...和“”

有什么建议吗？

注意：我是第一次在这里提出问题。我试图遵守规则。 我还在学习。

Answer 1

尝试以下示例：

library(data.table)

# example vcf
hap <- fread('
##fileformat=VCFv4.0
##fileDate=20090805
##source=myImputationProgramV3.1
##reference=1000GenomesPilot-NCBI36
##phasing=partial
##INFO=<ID=NS,Number=1,Type=Integer,Description="Number of Samples With Data">
##INFO=<ID=AN,Number=1,Type=Integer,Description="Total number of alleles in called genotypes">
##FILTER=<ID=q10,Description="Quality below 10">
##FILTER=<ID=s50,Description="Less than 50% of samples have data">
##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype">
#CHROM  POS ID  REF ALT QUAL    FILTER  INFO    FORMAT  NA00001 NA00002 NA00003
19  111 .   A   C   9.6 .   .   GT  0|0 0|0 0|1
19  112 .   A   G   10  .   .   GT  0|0 1|0 1|1
19  112 .   A   G   4   .   .   GT  1|0 1|0 1|1
')

data.table(gsub("|", "", do.call(paste0, hap[, -c(1:9)]), fixed = TRUE))
#        V1
# 1: 000001
# 2: 001011
# 3: 101011