Write a gsi_sim file from a data frame (wide or long/tidy). Used internally in assigner and radiator and might be of interest for users.
write_gsi_sim(
data,
pop.levels = NULL,
pop.labels = NULL,
strata = NULL,
filename = "gsi_sim.unname.txt"
)A tidy genomic data set in the working directory tidy formats.
How to get a tidy data frame ?
Look for radiator tidy_genomic_data.
(option, string) This refers to the levels in a factor. In this
case, the id of the pop.
Use this argument to have the pop ordered your way instead of the default
alphabetical or numerical order. e.g. pop.levels = c("QUE", "ONT", "ALB")
instead of the default pop.levels = c("ALB", "ONT", "QUE").
Default: pop.levels = NULL. If you find this too complicated, there is also the
strata argument that can do the same thing, see below.
(optional, string) Use this argument to rename/relabel
your pop or combine your pop. e.g. To combine "QUE" and "ONT"
into a new pop called "NEW":
(1) First, define the levels for your pop with pop.levels argument:
pop.levels = c("QUE", "ONT", "ALB").
(2) then, use pop.labels argument:
pop.levels = c("NEW", "NEW", "ALB").#'
To rename "QUE" to "TAS":
pop.labels = c("TAS", "ONT", "ALB").
Default: pop.labels = NULL. If you find this too complicated, there is also the
strata argument that can do the same thing, see below.
(optional) A tab delimited file with 2 columns with header:
INDIVIDUALS and STRATA.
Default: strata = NULL. Use this argument to rename or change
the populations id with the new STRATA column.
The STRATA column can be any hierarchical grouping.
The name of the file written to the working directory.
A gsi_sim input file is saved to the working directory.
Anderson, Eric C., Robin S. Waples, and Steven T. Kalinowski. (2008) An improved method for predicting the accuracy of genetic stock identification. Canadian Journal of Fisheries and Aquatic Sciences 65, 7:1475-1486.
Anderson, E. C. (2010) Assessing the power of informative subsets of loci for population assignment: standard methods are upwardly biased. Molecular ecology resources 10, 4:701-710.
/hrefhttps://github.com/eriqande/gsi_simgsi_sim and /hrefhttps://github.com/eriqande/rubiasrubias: genetic stock identification (GSI) in the tidyverse.