Quantifying missing data of all samples in VCF
Details
This function generates a dataframe of absolute and relative missing data per sample, and a plot of relative missing data % per sample.
Examples
f <- system.file("extdata", "vaillantii_discosnp_sub.vcf.gz",
package = "VCFArrow")
vcf <- read_vcf(f) |> set_vcf_groups(dirname(f))
#> ℹ VCF is being read in chunks of 50000 variants
#> ✔ VCF successfully read into a VCFArrow object
assess_vcf_missing_data(vcf, res_path = tempdir(),
species = "Phyllomedusa vaillantii",
project = "vaillantii")