cats is a fork of the R package CATS developed by Anders Albrechtsen <aalbrechtsen@bio.ku.dk> for power estimation in two-stage genome-wide association designs. See the vignette for code examples.
# install.packages("devtools")
devtools::install_github("nanxstats/cats")The method for power analysis implemented in CATS is based on Skol et al. 2006, but more generalized so that the ratio between cases and controls may vary between stages. The allele frequencies, disease prevalence, and relative risk can also vary.
The changes made to the original R package CATS are only for fixing the package build and check issues (under R 3.6.0 as of May 2019), for example:
- Used unified lowercase names for package and functions
- Fixed C code build errors and warnings (now passes clang)
- Fixed other package formatting issues that blocks building and checking
- Used roxygen2 to automate namespace and documentation generation
- Organized the existing code examples as a vignette
This package is forked and modified from these resources. They are also available as archives under bin/:
- Skol AD, Scott LJ, Abecasis GR, Boehnke M: Joint analysis is more efficient than replication-based analysis for two-stage genome-wide association studies. Nat Genet 38: 209-213, 2006.