The goal of segregatr is to provide segregation analysis for clinical variant classification.
You can install segregatr from CRAN as follows:
install.packages("segregatr")
Alternatively, obtain the latest development version from GitHub:
devtools::install_github("magnusdv/segregatr")
library(segregatr)
The family below shows four brothers, all affected with a rare dominant disease with 90% penetrance and phenocopy rate 1%. The parents have unknown affection status. All four brothers are shown to carry a candidate variant.
We will use segregatr to analyse the co-segregation of the variant and the disease in this pedigree. Specifically we want to compute the full-likelihood Bayes factor (FLB), quantifying the evidence that the variant is pathogenic.
To create the pedigree we use the nuclearPed()
function from the
pedtools package, which is automatically loaded together with
segregatr.
x = nuclearPed(4)
Then we run the FLB()
function, filling in the necessary data:
FLB(x, carriers = 3:6, affected = 3:6, unknown = 1:2,
freq = 0.0001, penetrances = c(0.01, 0.9, 0.9), proband = 3)
#> [1] 7.732161
The resulting FLB score is less than 8, which unfortunately only indicates suggestive evidence for pathogenicity.
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.