Populate INFO column with VariantAnnotation::writeVcf
还没有人认领这个 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
贡献指南
这个仓库没有索引到贡献指南
从这里开始
- 先读完整个 Issue,再读项目的贡献指南。
- 在 Issue 下留言说明你要接手 —— 这能避免两个人做同样的事。
- Fork 仓库,在一个分支上完成修改。
- 提交 Pull Request,并在描述里引用这个 Issue 编号。
Bioconductor/VariantAnnotation 的其他 Issue
-
难度 5/5 一周以上 新手友好度 25/100
Bioconductor/VariantAnnotation#115 · 4 条评论 ·
-
Bioconductor/VariantAnnotation#114 · 1 个 reaction · 已指派 2 人 ·
-
难度 5/5 一周以上 新手友好度 35/100
Bioconductor/VariantAnnotation#113 · 2 条评论 ·
-
难度 3/5 1-2 天 新手友好度 35/100
-
难度 5/5 一周以上 新手友好度 30/100
Bioconductor/VariantAnnotation#88 · 2 条评论 ·
查看 Bioconductor/VariantAnnotation 的全部 Issue
相似的 Issue
-
难度 2/5 1-3 小时 新手友好度 75/100
briandconnelly/airnow#9 ·
-
难度 2/5 1-3 小时 新手友好度 75/100
OHDSI/CohortConstructor#774 ·
-
pre-review R TeX Track: 5 (DSAIS)
难度 1/5 1 小时以内 新手友好度 60/100
openjournals/joss-reviews#11330 · 7 条评论 ·
-
难度 2/5 1-3 小时 新手友好度 75/100
-
难度 2/5 1-3 小时 新手友好度 75/100