marginalLikelihoodM2 | R Documentation |
Calculation of the marginal likelihood of Model 2, HET
x |
Dataframe output by tLOHDataImport |
a |
Alpha value |
b |
Beta value |
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.
The value returned from marginalLikelihoodM1 is numeric
Michelle Webb
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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.