Installation

This is an introduction to radmixture which could help you with estimating individual ancestries from large SNP genotype data.

This document introduces you how to use this package.

Install the latest development version from github with r devtools::install_github("wegene-llc/radmixture")
Install the latest version from CRAN r install.packages("radmixture")

First, you must prepare your raw data file as follow:

and read it into R with read.table or read.csv.

library(radmixture)
genotype <- read.table(file = '/path/to/file')
# genotype <- read.csv(file = 'path/to/file')

Calculate your ancestry components with public dataset

res <- tfrdpub(genotype, 13, global13.alleles, global13.13.F)

ances <- fFixQN(res$g, res$q, res$f, tol = 1e-4, method = 'BR', pubdata = "K13")

ances$q is your K13 result.

D.H. Alexander, J. Novembre, and K. Lange. Fast model-based estimation of ancestry in unrelated individuals. Genome Research, 19:1655–1664, 2009.
H. Zhou, D. H. Alexander, and K. Lange. A quasi-Newton method for accelerating the convergence of iterative optimization algorithms. Statistics and Computing, 2009.

We welcome contributions on radmixture. You can fork this repo and make your changes and submit a pull request. Report bug issues on issue page.

MIT + LICENSE