marginalLikelihoodM2: Marginal M2 Calculation

View source: R/functions.R

marginalLikelihoodM2R Documentation

Marginal M2 Calculation

Description

Calculation of the marginal likelihood of Model 2, HET

Arguments

x

Dataframe output by tLOHDataImport

a

Alpha value

b

Beta value

Details

The reference and total counts should come from a .csv output by the spatial LOH pre-processing pipeline. The recommended values for both Alpha2 and Beta2 is 500.

Value

The value returned from marginalLikelihoodM1 is numeric

Author(s)

Michelle Webb

Examples

test <- data.frame(REF=c(10,2,3,4,5,10),TOTAL=c(20,20,20,20,20,20))
apply(test, MARGIN = 1, FUN = marginalLikelihoodM1, a = 500, b = 500)

USCDTG/tLOH documentation built on Oct. 23, 2022, 8:05 p.m.