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
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.