Populate INFO column with VariantAnnotation::writeVcf
Bewertung
- Schwierigkeit
- 3/5
- Geschätzter Aufwand
- 1-2 Tage
- Anfängerfreundlichkeit
- 35/100
- Issue-Typ
- Bug
- Klarheit
- Größtenteils klar
- Aktivitätsstatus
- Veraltet
- Tech-Stack
- r
- Bereich
- bioinformatics
Rechercherichtung
Beginnen Sie damit, das bereitgestellte R-Beispiel auszuführen, und untersuchen Sie das konstruierte VRanges-Objekt vor und nach VariantAnnotation::writeVcf. Verfolgen Sie den Einstiegspunkt von writeVcf, um zu sehen, wie INFO und zusätzliche Spalten behandelt werden. Als erledigt gilt die Aufgabe, wenn erklärt oder korrigiert wurde, warum die angeforderten RSID-Einträge in der erzeugten VCF fehlen.
Vom Indexierungsmodell aus dem Issue-Text verfasst.
Beschreibung
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
- Vorherrschende Sprache
- R
- Sterne
- 32
- Forks
- 21
- PR-Merge-Kennzahlen
- Keine gemergten PRs in 30 T.
Entwicklungsumgebung
Dieses Projekt bietet weder Dev-Container noch Dockerfile noch Beitragsleitfaden – die Einrichtung liegt bei Ihnen. Beginnen Sie mit der README; die allgemeinen Schritte stehen in unserem Leitfaden für den ersten Beitrag.
Erste Schritte
- Lesen Sie das ganze Issue und danach den Beitragsleitfaden des Projekts.
- Schreiben Sie ins Issue, dass Sie es übernehmen — das erspart doppelte Arbeit.
- Forken Sie das Repository und arbeiten Sie in einem Branch.
- Öffnen Sie einen Pull Request, der die Issue-Nummer nennt.
Mehr aus Bioconductor/VariantAnnotation
-
refactor the packageOffen
Schwierigkeit 5/5 Über eine Woche Anfängerfreundlichkeit 25/100
Bioconductor/VariantAnnotation#115 · 4 Kommentare ·
-
what is the "covr" condition of this package?Evtl. vergeben @vjcitn hat das vor 96 Tagen übernommen. Offen
Bioconductor/VariantAnnotation#114 · 1 Reaktion · 2 zugewiesene Personen ·
-
produce false positive and false negative PRsEvtl. vergeben @jmg421 hat das vor 86 Tagen übernommen. Offen
Schwierigkeit 5/5 Über eine Woche Anfängerfreundlichkeit 35/100
Bioconductor/VariantAnnotation#113 · 2 Kommentare ·
-
possible problem parsing INFO description fieldEvtl. vergeben @jmg421 hat das vor 118 Tagen übernommen. Offen
Schwierigkeit 3/5 1-2 Tage Anfängerfreundlichkeit 35/100
-
Schwierigkeit 5/5 Über eine Woche Anfängerfreundlichkeit 30/100
Bioconductor/VariantAnnotation#88 · 2 Kommentare ·
Alle Issues in Bioconductor/VariantAnnotation
Ähnliche Issues
-
Schwierigkeit 1/5 Unter einer Stunde Anfängerfreundlichkeit 78/100
-
Typos in schema fileOffen
Schwierigkeit 1/5 1-3 Stunden Anfängerfreundlichkeit 65/100
-
Schwierigkeit 2/5 1-3 Stunden Anfängerfreundlichkeit 69/100
Merck/metalite.ae#261 ·
-
Schwierigkeit 2/5 1-3 Stunden Anfängerfreundlichkeit 68/100
-
[R][Parquet] ParquetReaderProperties does not expose page_checksum_verification, so a corrupt page CRC reads silentlyEvtl. wieder frei Ein Pull Request dazu wurde ohne Merge geschlossen. OffenComponent: R Type: enhancement
Schwierigkeit 2/5 1-3 Stunden Anfängerfreundlichkeit 72/100
apache/arrow#53138 · 1 Kommentar ·
Maintainer antworten meist innerhalb von 2 Tagen