Geneland-package {Geneland}R Documentation

Simulation and inference of a spatial statistical model for landscape genetics.

Description

The main function mcmcFmodel takes geo-referenced individual multilocus genetic data and tries to detect population structure, i.e sub-populations making use of both spatial and genetic information.

Details

Package: Geneland
Type: Package
Version: 1.0.8
Date: 2006-12-01
License: GPL

The overall population is viewed as a co-existence of sub-populations at Hardy-Weiberg equilibrium. The sub-populations are supposed to be spatially organised through the so-called colored Poisson-Voronoi tessellation. Allele frequencies are assumed to be drawn from the F-model as described in (Falush 2003) although the particular case of independent Dirichlet allele frequencies, as described by Pritchard (2000), is also handled. Individuals within populations are assumed to be randomly located and Hardy-Weinberg equilibrium and linkage equilibrium are assumed.

The main purpose of the program is to perform Bayesian inference of all the parameters involved through Markov Chain Monte-Carlo simulation. This is achievied by the function mcmcFmodel. Function PostProcessChain read some output files of mcmcFmodel and computes some statistics suitable to print maps of inferred populations.

See Storage format section in mcmcFmodel help page.

The following functions are provided by the package:

simFmodel: simulation from the prior of the spatial F-model

simdata: Simulation of georeferenced genotypes under an IBD + barrier model

show.simdata: Graphical display of data simulated by simdata

mcmcFmodel: Full Bayesian Markov Chain Monte Carlo inference of parameters in the spatial F-model

PostProcessChain: Post-procesing of MCMC output for maps of posterior probability of populations subdomains

PlotTessellation: Graphical display of inferred sub-domains

The following functions are very basic and are only intended to be an aid for those not familiar with R. Most probably you may want to use directly the output files of mcmcFmodel and PostProcessChain to print your own figures.

PlotDrift: Graphical display of drift factors along MCMC run

PlotFreqA: Graphical display of allele frequencies in the ancestral population along MCMC run

PlotFreq: Graphical display of allele frequencies in the present time population along MCMC run

Plotnpop: Graphical display of number of populations along MCMC run

Plotntile: Graphical display of number of tiles along MCMC run

PosteriorMode: Computation and/or graphical display of mode in the posterior distribution of class membership at each pixel

Fstat: Computations of pairwise F statistics between known subpopulations of a dataset

Fstat.output: Computations of pairwise F statistics between inferred subpopulations

nullify: Alter genotpyes of a dataset by simulating null alleles.

FormatGenotypes: Internal function. Transform a file of genotypes into a format suitable for function mcmcFmodel

setplot: Internal function

Author(s)

Arnaud Estoup, Gilles Guillot, Filipe Santos

www.inapg.inra.fr/ens_rech/mathinfo/personnel/guillot/welcome.html

See also

www.inapg.inra.fr/ens_rech/mathinfo/personnel/guillot/Geneland.html

References

[1] G. Guillot, Estoup, A., Mortier, F. Cosson, J.F. A spatial statistical model for landscape genetics. Genetics, 170, 1261-1280, 2005.

[2] G. Guillot, Mortier, F., Estoup, A. Geneland : A program for landscape genetics. Molecular Ecology Notes, 5, 712-715, 2005.

[3] A. Coulon, Guillot G., Cosson J.-F., Angibault J.M.A., Aulagnier S. , Cargnelutti B., Galan M., Hewison A.J.M. Genetics structure is influenced by lansdcape features. Empirical evidence from a roe deer population. Molecular Ecology, 15 1669-1679, 2006.

[4] G. Guillot, F. Santos, A. Estoup. Inference in population genetics with Geneland: a sensitivity analysis to spatial sampling scheme, null alleles and isolation by distance. Submitted.

See Also

Examples






[Package Geneland version 2.0.0 Index]