Hacktoberfest 2026:维护者为十月标记出来的 issue,仍然开放、适合新手。 浏览 Hacktoberfest issue

Populate INFO column with VariantAnnotation::writeVcf

未关闭
#61 4 条评论 0 个 reaction 已指派 0 人 在 GitHub 查看

还没有人认领这个 Issue。

评估

难度
3/5
预计耗时
1-2 天
新手友好度
35/100
Issue 类型
缺陷
描述清晰度
基本清楚
活跃度
停滞
技术栈
r

调研方向

首先运行提供的 R 示例,并检查 VariantAnnotation::writeVcf 前后构造的 VRanges 对象。跟踪 writeVcf 的入口点,以了解 INFO 和额外列是如何处理的。在解释或修正所请求的 RSID 条目为何未出现在生成的 VCF 中之后,即视为完成。

由索引模型根据 Issue 内容生成。

描述

Hi ,

I'm the creator of the MungeSumstats Bioconductor package in which we use VariantAnnotation::writeVcf as an output options for the formatted sumstats. I want to align the formatting of our sumstats to IEU OpenGWAS VCF, part of which means adding the RSIDs to the INFO column in the VCF. However, I can't seem to populate INFO with VariantAnnotation::writeVcf. Is there any way to do this?

An example:

#get example data 
sumstats_dt <- MungeSumstats::formatted_example()
#convert to VRanges
gr <- to_granges(sumstats_dt)
message("Converting summary statistics to VRanges.")
gr$dummy <- "GWAS"
vr <- VariantAnnotation::makeVRangesFromGRanges(gr,
    ref.field = "A1", 
    alt.field = "A2",
    keep.extra.columns = TRUE,
    sampleNames.field = "dummy"
)
#make SNP field into INFO
vr$INFO <- paste0("RSID=",vr$SNP)
vr$SNP <- NULL
#Output as VCF
VariantAnnotation::writeVcf(
            obj = vr,
            ### Must supply filename without compression suffix
            filename ="~/Downloads/test_vcf.vcf",
            index = FALSE
        )
#inspect
print(readLines("~/Downloads/test_vcf.vcf",100))
#INFO not populated
主要语言
R
星标
32
派生
21
PR 合并指标
30 天内没有已合并 PR

贡献指南

这个仓库没有索引到贡献指南

从这里开始

  1. 先读完整个 Issue,再读项目的贡献指南。
  2. 在 Issue 下留言说明你要接手 —— 这能避免两个人做同样的事。
  3. Fork 仓库,在一个分支上完成修改。
  4. 提交 Pull Request,并在描述里引用这个 Issue 编号。

Bioconductor/VariantAnnotation 的其他 Issue

查看 Bioconductor/VariantAnnotation 的全部 Issue

相似的 Issue

更多 R Issue

把新 issue 发到你的邮箱

精选适合新手参与的 GitHub issue 摘要。