Home My Page Projects Code Snippets Project Openings

Tag Cloud

Android Annuaire C C++ Data Mining ESUP Esup-Commons Experiments Git JSF Java Linux MATLAB MTGAUE Markov Decision Processes Matlab Neuroscience ODF Octave Open Educational Resources OpenMP PHP Perl Primer Python Python3 R Satellite Image Time Series Scilab Script Simulation Simulink TEI TEI. WIMS XML XSLT acquisition affiche planning algorithm alignment analyse texte annuaire applications astronomie astrophysics astrophysique atmosphere auth bash bibliotheque bibliotheques bioinformatics bioinformatique c c++ calcul scientifique calendrier climate cmake cnrs comparative genomics computer science controle crb data data analysis data mining data processing database django docker docker-compose e-education editeur editeurs editor environment esup exercices interactifs file formations fortran framework genetics gestion git hydrologie ical ics image processing informatique interface inventaire isotopic pattern java javascript jsonrpc k-mers kernel lamp library linguistique logging mapping mass spectrometry mathematics matlab mission model modelisation monitoring mudslide pedagogie peptide performance permutation method php pl portlet postgresql production protein proteomics python python3 qooxdoo qt5 quality quality alignment quantification read recommender systems satellite optical images science script scripting scripts security sequencing signal processing simulation spectrum statistical tests statistics statistiques toolbox unitary events user friendly visualisation web wordpress workflow zone critique
1 project in result set.
Random Forest ABC for IM model (micros) - Script in R to perform an Approximate Bayesian Computation (ABC) analysis under the Isolation with Migration model from microsatellite data. Demographic history was inferred from microsatellite data using an approximate Bayesian computation (ABC) approach via random forests (Marin et al. 2016; Pudlo et al. 2016). In this approach, data is simulated from the demographic model with parameter values taken from prior probability distributions and data is transformed into summary statistics. Random forests are used to learn about the parameters from the simulated summary statistics. The resulting random forests can then be used to estimate the posterior probability distributions of parameters. In this analysis, a model of two populations (eastern and western clusters) is evaluated. Each population is characterized by a parameter θ (θW=4NWμ and θE=4NEμ, where NW is the effective population size of the western population, NE is the effective population size of the eastern population, and μ is the mutation rate). Western population is founded by individuals from the eastern population at time T=t/4NW (time tF measured in number of generations). Two concurrent models are evaluated regarding the presence or absence of gene flow between the two populations. In the case of presence of gene flow there is an additional parameter, the scaled migration rate M=4NWm. Microsatellites are assumed to mutate following a generalized stepwise mutation model (GSM), in which the number of repeat units gained or lost in each mutation is taken from a geometric distribution with parameter PGSM. Data under this model is generated by simulation using coalescent simulator ms (Hudson 2002) with a custom script (see below) to transform its output into microsatellite data. Each simulated data set was summarized by statistics used in population genetics to characterize microsatellite genetic diversity and population differentiation, known to be informative about demography. Data input: file in STRUCTURE format (.str) with alleles coded in number of repeats (microsatellite data), two rows per individual, second column indicating the population, missing data coded as -9
Tags: population genetics, R
This project has not yet categorized itself in the Trove Software Map

Activity Percentile: 0
Activity Ranking: 0
Registered: 2017-01-24 09:53