Summarize the posterior model support by the number of SNPs contained in a model

Summarize the posterior model support by the number of SNPs contained in a model

pp.nsnp(results, plot = FALSE, expected = NULL, overdispersion = 1)

Arguments

results

object of class snpmod
plot

if TRUE, print a plot to current plotting device
expected

optional, if specified, this is the prior expectation of the number of causal SNPs underlying the trait in the studied region. Specifying will cause a line for the prior to be added to the plot
overdispersion

optional, default=1. Ignored unless expected is given. The overdispersion of the beta binomial distribution used for the prior relative to a binomial.

Value

a named list containing

  • "pp"a named vector summarizing the posterior for each count of SNPs

  • "plot"a ggplot object

See also

snpprior

Other plotting GUESSFM results: addlines, ggbed, ggchr, ggld, scalepos, signal.plot