Population Assignment using Genetic, Non-Genetic or Integrated Data in a Machine-learning Framework
This R package helps perform population assignment and infer population structure using a machine-learning framework. It employs supervised machine-learning methods to evaluate the discriminatory power of your data collected from source populations, and is able to analyze large genetic, non-genetic, or integrated (genetic plus non-genetic) data sets. This framework is designed for solving the upward bias issue discussed in previous studies. Main features are listed as follows.
- Use principle component analysis (PCA) for dimensionality reduction (or data transformation)
- Use Monte-Carlo cross-validation to estimate mean and variance of assignment accuracy
- Use K-fold cross-validation to estimate membership probability
- Allow to resample various sizes of training datasets (proportions or fixed numbers of individuals and proportions of loci)
- Allow to choose from various proportions of training loci either randomly or based on locus Fst values
- Provide several machine-learning classification algorithms, including LDA, SVM, naive Bayes, decision tree, and random forest, to build tunable predictive models.
- Output results in publication-quality plots that can be modified using ggplot2 functions
You can install the released version from CRAN or the up-to-date version from this Github respository.
-
To install from CRAN
- Simply enter
install.packages("assignPOP")
in your R console
- Simply enter
-
To install from Github
- step 1. Install devtools package by entering
install.packages("devtools")
- step 2. Import the library,
library(devtools)
- step 3. Then enter
install_github("alexkychen/assignPOP")
- step 1. Install devtools package by entering
Note: When you install the package from Github, you may need to install additional packages before the assignPOP can be successfully installed. Follow the hints that R provided and then re-run install_github("alexkychen/assignPOP")
.
Please visit our tutorial website for more infomration
Changes in ver. 1.1.9 (2020.3.16)
- Fix input non-genetic data (x1) error in assign.X
History
Changes in ver. 1.1.8 (2020.2.28)
- update following functions to work with R 4.0.0
- accuracy.MC, accuracy.kfold, assign.matrix, compile.data, membership.plot
- add stringsAsFactor=T to read.table and read.csv
- temporarily turn off testthat due to its current failure to pass test in Debian system
Changes in ver. 1.1.7 (2019.8.26)
- add broken-stick method for principal component selection in assign.MC, assign.kfold, and assign.X functions
- update accuracy.MC, accuracy.kfold, assign.matrix to handle missing levels of predicted population in test results
- update assign. and accuracy. functions to handle numeric population names
Changes in ver. 1.1.6 (2019.6.8)
- fix multiprocess issue in assign.kfold function
Changes in ver. 1.1.5 (2018.3.23)
- Update assign.MC & assign.kfold to detect pop size and train.inds/k.fold setting
- Update accuracy.MC & assign.matrix to handle test individuals not from every pop
- Slightly modify levels method in accuracy.kfold
- fix bugs in accuracy.plot for K-fold results
- fix membership.plot title positioning and set text size to default
Changes in ver. 1.1.4 (2018.3.8)
- Fix missing assign.matrix function
Changes in ver. 1.1.3 (2017.6.15)
- Add unit tests (using package testthat)
Changes in ver. 1.1.2 (2017.5.13)
- Change function name read.genpop to read.Genepop; Add function read.Structure.
- Update read.genpop function, now can read haploid data
Chen, K. Y., Marschall, E. A., Sovic, M. G., Fries, A. C., Gibbs, H. L., & Ludsin, S. A. (2018). assign POP: An R package for population assignment using genetic, non-genetic, or integrated data in a machine-learning framework. Methods in Ecology and Evolution. 9(2)439-446. https://doi.org/10.1111/2041-210X.12897
Previous packages can be found and downloaded at the releases page