VariantFiltering (version 1.4.3)

autosomalDominant: Autosomal dominant inheritance analysis

Description

This function identifies variants present in all the affected individual(s) discarding the ones that also occur in at least one of the unaffected subjects.

Usage

"autosomalDominant"(param, BPPARAM=bpparam("SerialParam"))

Arguments

param

A VariantFilteringParam object built from a multisample VCF file with at least one affected individual and zero or more unaffected ones, and from a PED file specyfing the family relationships among individuals as well as their gender and phenotype status (affected or unaffected).

BPPARAM

An object of class BiocParallelParam specifiying parameters related to the parallel execution of some of the tasks and calculations within this function. See function bpparam() from the BiocParallel package.

Value

An object of class VariantFilteringResults including functional annotations on variants segregating according to an autosomal dominant inheritance model.

Details

This function requires as an input a VariantFilteringParam class object built from an input multisample VCF file, along with a PED file.

References

Elurbe D.M., Mila, M., Castelo, R. The VariantFiltering package, in preparation.

See Also

autosomalRecessiveHomozygous autosomalRecessiveHeterozygous xLinked deNovo unrelatedIndividuals VariantFilteringResults

Examples

## Not run: 
# 
# CEUvcf <- file.path(system.file("extdata", package="VariantFiltering"), "CEUtrio.vcf.bgz")
# CEUped <- file.path(system.file("extdata", package="VariantFiltering"), "CEUtrio.ped")
# param <- VariantFilteringParam(vcfFilenames=CEUvcf,
#                      pedFilename=CEUped)
# aDo <- autosomalDominant(param)
# aDo
# ## End(Not run)