Description Usage Arguments Details Value Author(s) See Also Examples
batchChisqTest
calculates Chi-square values for batches from
2-by-2 tables of SNPs, comparing each batch with the other batches.
batchFisherTest
calculates Fisher's exact test values.
1 2 3 4 5 6 7 8 9 | batchChisqTest(genoData, batchVar, snp.include = NULL,
chrom.include = 1:22, sex.include = c("M", "F"),
scan.exclude = NULL, return.by.snp = FALSE,
correct = TRUE, verbose = TRUE)
batchFisherTest(genoData, batchVar, snp.include = NULL,
chrom.include = 1:22, sex.include = c("M", "F"),
scan.exclude = NULL, return.by.snp = FALSE,
conf.int = FALSE, verbose = TRUE)
|
genoData |
|
batchVar |
A character string indicating which annotation variable should be used as the batch. |
snp.include |
A vector containing the IDs of SNPs to include. |
chrom.include |
Integer vector with codes for chromosomes to
include. Ignored if |
sex.include |
Character vector with sex to include. Default is
c("M", "F"). If sex chromosomes are present in |
scan.exclude |
A vector containing the IDs of scans to be excluded. |
return.by.snp |
Logical value to indicate whether snp-by-batch matrices should be returned. |
conf.int |
Logical value to indicate if a confidence interval should be computed. |
correct |
Logical value to specify whether to apply the Yates continuity correction. |
verbose |
Logical value specifying whether to show progress information. |
Because of potential batch effects due to sample processing and genotype calling, batches are an important experimental design factor.
batchChisqTest
calculates the Chi square values from 2-by-2
table for each SNP, comparing each batch with the other batches.
batchFisherTest
calculates Fisher's Exact Test from 2-by-2
table for each SNP, comparing each batch with the other batches.
For each SNP and each batch, batch effect is evaluated by a 2-by-2 table:
# of A alleles, and # of B alleles in the batch, versus
# of A alleles, and # of B alleles in the other batches.
Monomorphic SNPs are set to NA
for all batches.
The default behavior is to combine allele frequencies from males and
females and return results for autosomes only. If results for sex
chromosomes (X or Y) are desired, use chrom.include
with values
23 and/or 25 and sex.include
="M" or "F".
If there are only two batches, the calculation is only performed once and the values for each batch will be identical.
batchChisqTest
returns a list with the following elements:
mean.chisq |
a vector of mean chi-squared values for each batch. |
lambda |
a vector of genomic inflation factor computed as |
chisq |
a matrix of chi-squared values with SNPs as rows and
batches as columns. Only returned if |
batchFisherTest
returns a list with the following elements:
mean.or |
a vector of mean odds-ratio values for each
batch. |
lambda |
a vector of genomic inflation factor computed as
|
Each of the following is a matrix with SNPs as rows and batches as
columns, and is only returned if return.by.snp=TRUE
:
pval |
P value |
oddsratio |
Odds ratio |
confint.low |
Low value of the confidence interval for the odds
ratio. Only returned if |
confint.high |
High value of the confidence interval for the odds
ratio. Only returned if |
batchChisqTest
and batchFisherTest
both also return the
following if return.by.snp=TRUE
:
allele.counts |
matrix with total number of A and B alleles over all batches. |
min.exp.freq |
matrix of minimum expected allele frequency with SNPs as rows and batches as columns. |
Xiuwen Zheng, Stephanie Gogarten
GenotypeData
, chisq.test
,
fisher.test
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 | library(GWASdata)
file <- system.file("extdata", "illumina_geno.gds", package="GWASdata")
gds <- GdsGenotypeReader(file)
data(illuminaScanADF)
genoData <- GenotypeData(gds, scanAnnot=illuminaScanADF)
# autosomes only, sexes combined (default)
res.chisq <- batchChisqTest(genoData, batchVar="plate")
res.chisq$mean.chisq
res.chisq$lambda
# X chromosome for females
res.chisq <- batchChisqTest(genoData, batchVar="status",
chrom.include=23, sex.include="F", return.by.snp=TRUE)
head(res.chisq$chisq)
# Fisher exact test of "status" on X chromosome for females
res.fisher <- batchFisherTest(genoData, batchVar="status",
chrom.include=23, sex.include="F", return.by.snp=TRUE)
qqPlot(res.fisher$pval)
close(genoData)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.