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.
First, you must prepare your raw data file as follow:
and read it into R with
library(radmixture) genotype <- read.table(file = '/path/to/file') # genotype <- read.csv(file = 'path/to/file')
tfrdpub()to transfer your raw data to a format understood by radmixture
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.
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