Description Usage Arguments Details Value Author(s) References See Also Examples
BST is based on the test statistic of Goeman et al (2005) following a general empirical Bayes method to test on a large number of parameters in a logistic regression model. BST is closely related to the usual Score test, although it assumes an empirical Bayesian model with an independent prior on the genetic variant effects. The null distribution of the BST statistic is unknown and has to be estimated by permutation.
1 | BST(y, X, perm = 100)
|
y |
numeric vector with phenotype status: 0=controls, 1=cases. No missing data allowed |
X |
numeric matrix or data frame with genotype data coded as 0, 1, 2. Missing data is allowed |
perm |
positive integer indicating the number of permutations (100 by default) |
The BST statistic does not offer an asymptotic p-value. Permutations are required
There is no imputation for the missing data. Missing values are simply ignored in the computations.
An object of class "assoctest"
, basically a list with the following elements:
bst.stat |
bst statistic |
perm.pval |
permuted p-value |
args |
descriptive information with number of controls, cases, variants, and permutations |
name |
name of the statistic |
Gaston Sanchez
Goeman JJ, van de Geer SA, van Houwelingen HC (2006) Testing against a high dimensional alternative. Journal of the Royal Statistical Society, 68: 477-493
Chapman J, Whittaker J (2008) Analysis of Multiple SNPs in a Candidate Gene or Region. Genetic Epidemiology, 32: 560-566
Pan W (2009) Asymptotic Tests of Association with Multiple SNPs in Linkage Disequilibrium. Genetic Epidemiology, 33: 497-507
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 | ## Not run:
# number of cases
cases = 500
# number of controls
controls = 500
# total (cases + controls)
total = cases + controls
# phenotype vector
phenotype = c(rep(1,cases), rep(0,controls))
# genotype matrix with 10 variants (random data)
set.seed(123)
genotype = matrix(rbinom(total*10, 2, 0.05), nrow=total, ncol=10)
# apply BST with 500 permutations
mybst = BST(phenotype, genotype, perm=500)
mybst
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.