The output file contains the snp data in the format expected by treemix -- see the treemix manual. The file needs to be gzipped before it will be recognised by treemix. Plotting functions can be obtained using source("C:/workspace/R/dartR/R/plotting_funcs.R").
gl2treemix(x, outfile = "treemix_input.gz", outpath = tempdir(),
v = 2)
-- name of the genlight object containing the SNP data [required]
-- file name of the output file (including extension).
-- path where to save the output file (set to tempdir by default)
-- verbosity: 0, silent or fatal errors; 1, begin and end; 2, progress log ; 3, progress and results summary; 5, full report [default 2]
Reference: Pickrell and Pritchard (2012). Inference of population splits and mixtures from genome-wide allele frequency data. PLoS Genetics https://doi.org/10.1371/journal.pgen.1002967
# NOT RUN {
gl2treemix(testset.gl)
# }
Run the code above in your browser using DataLab