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.

Installation

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

Example

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.

References

Contributions

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.

License

MIT + LICENSE

TODO