s_ng | R Documentation |
Standard Normal Gammp prior calculates the posterior distribution for the fine mapping study. The number of individuals must be greater than the number of SNPs.
s_ng(y, data, medstar = 1, numb = 100, burnin = 1, every = 1)
y |
A vector of the pheontype. |
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. |
A list including:
alpha |
A vector of the posterior distribution of the intercept. |
beta |
A matrix of the posterior distribution of the effect sizes. |
sigmasq |
A vector of the posterior distribution of σ^2. |
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. |
Abulaziz Alenazi.
R implementation and documentation: Abulaziz Alenazi a.alenazi@nbu.edu.sa.
set.seed(1) data <- matrix( rnorm(500 * 30), ncol = 30) s_ng(y = rnorm(500), data = data)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.