BayesAB is an implementation of Bayes A and B for estimating direct SNP effects in
high dimensional data problems (p >> N).
BayesAB utilizes the C Function
cBaysABF for speed ..
BayesAB(ga, numiter = 5000, numMHIter = 10, Pi = .9, y)
A matrix of genotypes with a number of rows identical to the number of genotyped individuals and a number of columns identical to the number of SNPs. Values in the matrix are 0, 1, 2, & 5 for homozygous, heterozygous, other homozygous, & unknown genotypes, respectively.
Number of iterations
Number of Metropolis-Hastings iterations. Set to a value of 1 with Bayes A.
Proportion of SNP loci with 0 effect. Set to 0 to run Bayes A instead of B.
Trait phenotypes or conventional breeding values
This function runs Bayes A and B to estimate direct SNP effects based on a matrix of genotypes,
ga and a vector of adjusted phenotypes,
y, (Meuwissen et al., 2001; Genetics 157:1819-1829).
This is essential in high dimensional data problems with highly overparameterized models (p >> N).
To run Bayes A, set
Pi to 0 and
numMHIter to 1.
Other data management functions in
gdmp can be used to construct the integer matrix of genotypes,
ga, to use as input to
A list object with a vector of SNP estimates,
meanb, and a vector of genomic values for individuals,
aHat, are returned. It is also possible to extract the number of SNP loci in
Meuwissen et al. (2001). Prediction of Total Genetic Value Using Genome-Wide Dense Marker Maps. Genetics, 157, 1819–1829.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.