Marginal PP for models sharing information between traits, when trait correlation is zero

Marginal PP for models sharing information between traits, when trait correlation is zero

marginalpp0(STR, PP, mbeta, kappa, N, nsnps, NCORES)

Arguments

STR

list of models for traits 1, 2, ..., n, each given in the form of a character vector, with entries "snp1%snp2%snp3". The null model is given by "1" OR "0". It is assumed that all elements of ABF, PP and pr below follow this same order.

PP

list of posterir probabilities for the models in M

mbeta

list of joint beta estimates for each trait

kappa

single value or vector of values to consider for the sharing scale parameter. the value of kappa=1 must be included, and if not will be prepended.

N

number of individiduals with measurements for all traits

nsnps

number of snps in region

NCORES

number of cores for parallel computing; recommend NCORES=M, but if on Windows, use NCORES=1

Value

list of: - single.pp: list of pp for each model in STR[[i]] for trait i - shared.pp: list of pp for each model in STR[[i]] for trait i