s_ng {NGBVS}

Standard NG prior

Description

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.

Usage

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

Arguments

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.

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.

sigmasq

A vector of the posterior distribution of \sigma^2.

psi

A matrix of the posterior distribution of \psi.

lambda

A vector of the posterior distribution of \lambda.

gammasq

A vector of the posterior distribution of \gamma^2.

Author(s)

Abulaziz Alenazi.

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

Examples

set.seed(1)
data <- matrix( rnorm(500 * 30), ncol = 30)
s_ng(y = rnorm(500), data = data)