CoPheScan with fixed priors |
|
---|---|
per.snp.priors |
|
adjust_priors |
|
hypothesis.priors |
|
Bayesian cophescan analysis using Approximate Bayes Factors |
|
cophe.single.lbf |
|
run |
|
cophe.susie.lbf |
|
combine.bf |
|
Run cophescan on multiple traits at once |
|
print the summary of results from cophescan single or susie |
|
Simplifying the output obtained from |
|
logsum |
|
CoPheScan with hierarchical priors |
|
Run the hierarchical Metropolis Hastings model to infer priors |
|
Average of priors: pnk, pak and pck |
|
Average of priors: pnk, pak and pck from list (memory intensive) |
|
Average of posterior probabilities: Hn, Ha and Hc |
|
Average of posterior probabilities: Hn, Ha and Hc from list (memory intensive) |
|
Calculation of the posterior prob of Hn, Ha and Hc |
|
Extract beta and p-values of queried variant |
|
sample alpha |
|
sample beta |
|
sample gamma |
|
dnorm for alpha |
|
dgamma for beta |
|
dgamma for gamma |
|
Log likelihood calculation |
|
Log posterior calculation |
|
Log sum |
|
Run the hierarchical mcmc model to infer priors |
|
Conversion of parameters alpha, beta and gamma to pnk, pak and pck |
|
List of priors: pn, pa and pc over all iterations |
|
List of posterior probabilities: Hn, Ha and Hc over all iterations |
|
Calculate log priors |
|
Initiate parameters alpha, beta and gamma |
|
Proposal distribution |
|
Target distribution |
|
Predict hypothesis |
|
Predict cophescan hypothesis for tested associations |
|
Estimate the Hc.cutoff for the required FDR |
|
Visualization |
|
Plot region Manhattan for a trait highlighting the queried variant |
|
cophe_plots showing the Ha and Hc of all traits and labelled above the specified threshold |
|
Heatmap of multi-trait cophescan results |
|
Prepare data for plotting |
|
Test data |
|
Simulated multi-trait data |
|
Package |
|
The 'cophescan' package. |