asym_s_ng {NGBVS}

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 \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(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)