login

An R Package for Analysis of Whole-Genome Association Studies

Human HeredityPublished 1 January 2007
David Clayton, H.T. Leung
Citations99
SJR quartileQ3
SJR score0.48
SNIP0.44

TL;DR

Data classes in which each genotype call is stored as a single byte are implemented to facilitate the analysis of whole genome association studies in the R language for statistical computing.

Abstract

OBJECTIVE: To provide data classes and methods to facilitate the analysis of whole genome association studies in the R language for statistical computing. METHODS: We have implemented data classes in which each genotype call is stored as a single byte. At this density, data for single chromosomes derived from large studies and new high-throughput gene chip platforms can be handled in memory. We use the object-oriented programming model introduced with version 4 of the S-plus package, usually termed 'S4 methods'. RESULTS: At the current state of development the package only supports population-based studies, although we would hope to provide support for family-based studies soon. Both quantitative and qualitative phenotypes may be analysed. Flexible association testing functions are provided which can carry out single SNP tests which control for potential confounding by quantitative and qualitative covariates. Tests involving several SNPs taken together as 'tags' are also supported. Efficient calculation of pair-wise linkage disequilibrium measures is implemented and data input functions include a function which can download data directly from the international HapMap project website.

Keywords

Biochemistry, Genetics and Molecular Biology