Parses and tidies DRAGEN variant-related metric outputs (variant caller, SV, CNV, ploidy, TMB, HRD, etc.).
Super classes
nemo::Tool -> tidydragen::DragenTool -> DragenVar
Public fields
flat_tidy_names(
logical(1))TRUE: fan-out sub-tables are named<tool>_<tidy_name>(parser token dropped). Needed for theploidystats/ploidyratiosplit.
Methods
Method new()
Create a new DragenVar object.
Usage
DragenVar$new(path = NULL, files_tbl = NULL)Arguments
path(
character(1))
Output directory of tool. Iffiles_tblis supplied, this is ignored.files_tbl(
tibble(n))
Tibble of files fromnemo::list_files_dir().
Method tidy_vc()
Tidy vc_metrics.csv. rg (sample id) kept raw; region moved
from metric names into a region column.
Method tidy_gvcf()
Tidy gvcf_metrics.csv (gVCF postfilter). Same shape as
vc_metrics; region moved into a region column.
Method tidy_ploidystats()
Tidy ploidy_estimation_metrics.csv into ploidystats (sample
scalars) and ploidyratio (long, one row per chromosome ratio).
Method parse_contamination()
Parse contamination.json (cttso cross-sample contamination).
Flat JSON; the per-SNP SNPsUsed array is dropped, NaN p-value -> NA.
Examples
cls <- DragenVar; tool <- "dragenvar"
indir <- system.file("extdata", tool, package = "tidydragen")
odir <- tempdir()
obj <- cls$new(indir)
obj$run(output_dir = odir, format = "parquet", input_id = "run1")
(lf <- list.files(odir, pattern = "dragenvar_.*parquet", full.names = FALSE))
#> [1] "sampleA_dragenvar_cnv.parquet"
#> [2] "sampleA_dragenvar_contamination.parquet"
#> [3] "sampleA_dragenvar_gvcf.parquet"
#> [4] "sampleA_dragenvar_hethom.parquet"
#> [5] "sampleA_dragenvar_hrd.parquet"
#> [6] "sampleA_dragenvar_microsat.parquet"
#> [7] "sampleA_dragenvar_nuctrans.parquet"
#> [8] "sampleA_dragenvar_ploidyratio.parquet"
#> [9] "sampleA_dragenvar_ploidystats.parquet"
#> [10] "sampleA_dragenvar_ploidyvcf.parquet"
#> [11] "sampleA_dragenvar_sv.parquet"
#> [12] "sampleA_dragenvar_tmb.parquet"
#> [13] "sampleA_dragenvar_vc.parquet"
#> [14] "sampleB_dragenvar_cnv.parquet"
#> [15] "sampleB_dragenvar_ploidyratio.parquet"
#> [16] "sampleB_dragenvar_ploidystats.parquet"
#> [17] "sampleB_dragenvar_vc.parquet"
