asym_s_ng: Standard NG prior

In NGBVS: Bayesian Variable Selection for SNP Data using Normal-Gamma

Standard NG prior

Description

Standard Normal Gammp prior calculates the posterior distribution for the fine mapping cases-controls study. The number of case-controls must be greater than the number of SNPs.

Usage

asym_s_ng(y, data,  medstar = 1, numb = 100, burnin = 1, every = 1)

Arguments

y	A vector of the pheontype, where takes 0s and 1s.
data	An N \times p finemap data, where N and p denote the samples and number of SNPs respectively.
medstar	The value of M.
numb	Number of samples for each SNP.
burnin	The amount of burn-in for the MCMC sample.
every	The amount of thining for the MCMC sample.

Value

A list including:

alpha	A vector of the posterior distribution of the intercept.
beta	A matrix of the posterior distribution of the effect sizes.
psi	A matrix of the posterior distribution of ψ.
lambda	A vector of the posterior distribution of λ.
gammasq	A vector of the posterior distribution of γ^2.

Author(s)

Abulaziz Alenazi.

R implementation and documentation: Abulaziz Alenazi [e-mail].

Examples

set.seed( 1 )
data <- matrix(sample( c( 0, 1, 2 ), 500 * 30, replace = TRUE,
prob = c( 0.35, 0.35, 0.3)), ncol = 30)
asym_s_ng(y = rbinom(500, 1, 0.5), data = data)

NGBVS documentation built on Sept. 16, 2022, 5:06 p.m.