GetMarginalPosteriors: Get Marginal {*U*,*D*,*I*} Posteriors

In mturchin20/bmass: Bayesian Multivariate Analysis of Summary Statistics

Description

Get marginal posteriors for how much every individual phenotype belongs to categories {U,D,I} across each SNP

Usage

GetMarginalPosteriors(DataSources, ListSNPs, Models, LogFile)

Arguments

DataSources	A string indicating the variable names of the input datafiles and phenotypes.
ListSNPs	A list produced from running bmass containing the SNPs of interest to get marginal posteriors for.
Models	A matrix describing the models being explored (default output from running bmass).
LogFile	A matrix of string outputs for function logging purposes (default output from running bmass).

Value

A list containing three matrices of SNPs x Phenotypes marginal posteriors for each category {U,D,I}; this list is appended to the input ListSNPs as a new object, Marginals (the full returned object is a list containing the input ListSNPs and the input LogFile).

Examples

Phenotypes <- c("bmass_SimulatedData1", "bmass_SimulatedData2")
bmassOutput <- bmass(Phenotypes, bmass_SimulatedSigSNPs)
bmassOutput[c("PreviousSNPs", "LogFile")] <-
GetMarginalPosteriors(Phenotypes, bmassOutput$PreviousSNPs,
bmassOutput$Models, bmassOutput$LogFile)
bmassOutput$PreviousSNPs$Marginals

mturchin20/bmass documentation built on Sept. 21, 2020, 7:51 p.m.