Description Usage Arguments Details Value Note Author(s) References See Also Examples
The Genomic Distance-Based Regression has been developed by Wessel et al (2006). This approach captures genotype information across multiple loci through a similarity measure between any two individuals. GDBR is unique in its regression analysis relating variation in the measure of genomic similarity to variation in their trait values. Note that this approach is computationally expensive.
1 |
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. |
distance |
character string indicating the type of distance to be used. Possible options are "IBS" or "wIBS" ( |
weights |
optional numeric vector with weights for the genetic variants ( |
perm |
positive integer indicating the number of permutations ( |
The argument distance
is used to specify the similarity distance. "IBS"
indicates Identity-By-Share, "wIBS"
indicates weighted IBS.
An object of class "assoctest"
, basically a list with the following elements:
gdbr.stat |
gdbr statistic |
perm.pval |
permuted p-value |
args |
descriptive information with number of controls, cases, variants, permutations, and selected distance |
name |
name of the statistic |
This method is computationally expensive
Gaston Sanchez
Wessel J, Schork NJ (2006) Generalized Genomic Distance-Based Regression Methodology for Multilocus Association Analysis. The American Journal of Human Genetics, 79: 792-806
Schaid DJ (2010) Genomic Similarity and Kernel Methods I: Advancements by Building on Mathematical and Statistical Foundations. The American Journal of Human Heredity, 70: 109-131
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 | ## Not run:
# number of cases
cases = 250
# number of controls
controls = 250
# 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 GDBR with 50 permutations
# (it takes some time to run the permutations!)
mygdbr = GDBR(phenotype, genotype, perm=50)
mygdbr
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.