Description Usage Arguments Value Author(s) References See Also Examples

`tdtGxG`

and `colGxG`

perform the genotypic TDT for the interaction of two SNPs or of each pair of columns
of a genotype matrix, respectively.

`fastGxG`

provides a fast implementation for the genotypic TDT for two-way interactions when considering the
simplest conditional logistic regression model only containing one parameter for the interaction effect. It thus leads
to the same results as `colGxG`

with `test = "screen"`

. In `fastGxGrec`

,
an analytic solution to the genotypic TDT based on the simplest model for testing a recessive x recessive model is
implemented, which is even faster than `fastGxG`

with `model = "recessive"`

. In future versions of this
package, `fastGxG`

and `fastGxGrec`

will be joint with `colGxG`

.

The genotypic TDT for testing two-way interactions makes use of the 16 possible genotypes that can be obtained from combining the parents' genotypes of the two considered SNPs. Thus, for each family, genotypes for one case (i.e. the affected offspring) and 15 pseudo-controls are used.

1 2 3 4 5 6 7 8 9 10 11 | ```
tdtGxG(snp1, snp2, test = c("epistatic", "lrt", "full", "screen"),
model = c("additive", "dominant", "recessive"))
colGxG(mat.snp, test = c("epistatic", "lrt", "full", "screen"), genes = NULL,
maf = FALSE, model = c("additive", "dominant", "recessive"))
fastGxG(mat.snp, model = c("additive", "dominant", "recessive"),
genes = NULL, interval = c(-10, 10), tol = 10^-8, maxiter = 1000,
size = 20)
fastGxGrec(mat.snp, genes = NULL, size = 20)
``` |

`snp1, snp2` |
numeric vectors of length |

`mat.snp` |
a numeric matrix in which each column represents a SNP. Each of the SNPs must have
the same structure as |

`test` |
character string naming the GxG test that should be performed. If |

`genes` |
a character vector containing the names of the genes to which the SNPs belong. If specified, only the two-way interactions
between SNPs from different genes are tested. If |

`maf` |
logical indicating whether the minor allele frequency (computed by considering the genotypes of only the parents) should be added to the output. |

`model` |
type of model that should be considered. Abbreviations are allowed. Thus, e.g., |

`interval` |
the end-points of the interval to be searched for the root. For details, see |

`tol` |
the desired accuracy/convergence tolerance. For details, see |

`maxiter` |
the maximum number of iterations. For details, see |

`size` |
the number of interactions considered simultaneously when computing the parameter estimates. |

Depending on `test`

, the output contains statistics and p-values either of a likelihood ratio test (`test = "epistatic"`

or
`test = "lrt"`

) or the Wald statistics and the corresponding p-values for the interaction term in the conditional logistic
regression model (`test = "full"`

or `test = "screen"`

). If `maf = TRUE`

, a vector `maf`

containing the minor allele frequencies of each SNP and a matrix `mat.maf`

with two columns containing the SNP-wise minor allele
frequencies for each tested pair of SNPs are added to the output of `colGxG`

.

Holger Schwender, holger.schwender@udo.edu

Cordell, H. J. (2002). Epistasis: What it Means, what it Doesn't mean, and Statistical Methods to Detect it in Humans. Human Molecular Genetics, 11, 2463-2468.

Schwender, H., Taub, M.A., Beaty, T.H., Marazita, M.L., and Ruczinski, I. (2011).
Rapid Testing of SNPs and Gene-Environment Interactions in Case-Parent Trio Data Based on
Exact Analytic Parameter Estimation. *Biometrics*, 68, 766-773.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 | ```
# Load the simulated data.
data(trio.data)
# The interaction between a particular pair of SNPs
# (e.g., the ones in the first and second column of
# mat.test) can be tested by
gxg.out <- tdtGxG(mat.test[,1], mat.test[,2])
# All pairs of SNPs in mat.test can be tested by
gxg.out2 <- colGxG(mat.test)
# By default, Cordell's likelihood ratio test for
# epistatistic interactions is used. This is the
# most sophisticated, but also most time-consuming
# test. If another test, e.g., the one considering
# a conditional logistic regression model only
# containing a term for the interaction, should
# be used, then this can be done by
gxg.out3 <- colGxG(mat.test, test = "screen")
# In this case, different modes of inheritance can
# be considered (by default, the additive mode is
# considered). If a dominant model (for both SNPs)
# should be tested, this can be done by
gxg.out4 <- colGxG(mat.test, test = "screen", model ="dom")
# If just a subset of all pairs of SNPs should be
# tested, e.g., only pairs of SNPs belonging to different
# genes, then this can be done by first specifying a
# vector specifying which SNP belongs to which genes.
# If we, e.g., assume that the first two SNPs in mat.test
# belong to gene G1 and the other four SNPs to G2, then
# this vector can be specified by
genes <- paste("G", rep(1:2, c(2,4)), sep="")
# and only the pairs of SNPs in which the two SNPs belong
# to different genes can be tested with Cordell's
# likelihood ratio test by
gxg.out5 <- colGxG(mat.test, genes = genes)
``` |

trio documentation built on Nov. 8, 2020, 7:41 p.m.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.