Skip to contents

Assign Priors to Major/Minor Allele Combinations Based on Biological Difficulty

Source: R/Modeling.R
AssignPriors.Rd

Assigns prior probabilities to major/minor allele combinations using an exponential decay model that incorporates biological plausibility.

Usage

AssignPriors(combinations, lambda)

Arguments

combinations

Data frame. Output from GenerateCombinations(), with columns major, minor, and CN.

lambda

Numeric. Decay rate parameter for the exponential prior.

Value

A data frame with the input columns plus:

Bio\_diff

Assigned biological difficulty score

prior

Normalized prior probability for each combination

Details

The function assigns a biological difficulty score to each combination. For unlisted combinations, the score is set according to rules based on CN, major, and minor. The prior is then calculated as \(\exp(-\lambda \times \text{Bio\_diff})\) and normalized to sum to 1.

Examples

combos <- GenerateCombinations(3)
AssignPriors(combos, lambda = 0.5)
#> # A tibble: 2 × 5
#>   major minor    CN Bio_diff prior
#>   <dbl> <dbl> <int>    <dbl> <dbl>
#> 1     3     0     3        4 0.269
#> 2     2     1     3        2 0.731