PostProcessChain          package:Geneland          R Documentation

_C_o_m_p_u_t_a_t_i_o_n _f_o_r  _m_a_p_s _o_f _p_o_s_t_e_r_i_o_r _p_r_o_b_a_b_i_l_i_t_y _o_f _p_o_p_u_l_a_t_i_o_n _m_e_m_b_e_r_s_h_i_p

_D_e_s_c_r_i_p_t_i_o_n:

     Computes posterior probabilities of population membership for each
     pixel of the spatial domain.

_U_s_a_g_e:

     PostProcessChain(coordinates,genotypes,allele.numbers,
     path.mcmc,nxdom, nydom,burnin)

_A_r_g_u_m_e_n_t_s:

coordinates: Spatial coordinates of individuals. A matrix with 2
          columns and one line per individual.

genotypes: Genotypes of individuals. A matrix with one line per
          individual and 2 columns per locus

allele.numbers: A vector of integer containing the number of possible
          allele for each locus

path.mcmc: Path to output files directory 

   nxdom: Number of pixel for discretization of the spatial domain in
          the horizontal direction

   nydom: Number of pixel for discretization of the spatial domain in
          the vertical direction

  burnin: Number of iterations of the chain to throw away. WARNING :
          this argument should be given the number of stored iterations
          (and not the number of computed iterations which differ if
          'burnin' !=1). If you have 'nit'=100000 and 'thinning'=100,
          then only 1000 iterations are stored. Then 'burnin'=10 will
          throw away 10 stored iterations, namely 100*10 computed
          iterations.

_V_a_l_u_e:

Posterior probability of population membership for each pixel:: They 
          are written in an ascii file called
          'proba.pop.membership.txt' (one column per population,
          'npopmax' values are computed for each pixel. Images in each
          column of 'proba.pop.membership.txt' are stored column-wise
          starting from the bottom left pixel. First line of
          'proba.pop.membership.txt'  = bottom left pixel , second line
          of 'proba.pop.membership.txt' = upward neighboor of the
          previous pixel, etc...)

          Another file called 'proba.pop.membership.perm.txt' tries to
          get rid of label switching issues by labelling the population
          according to a fixed constraint. (This has proved to be
          usefull with a small number of loci, (e.g. nloc=3), for well
          differentiated populations.)

Posterior probability of population membership for each
individual:: Th
          ey are written in a file named
          'proba.pop.membership.indiv.txt'. 

Label of modal population for individuals:: They are written in a file
          named 'modal.pop.indiv.txt'.  See the example section  of
          function 'mcmcFmodel' to see how they can be added in a plot.

_A_u_t_h_o_r(_s):

     Gilles Guillot

_R_e_f_e_r_e_n_c_e_s:

     A spatial statistical model for landscape genetics, Guillot,
     Estoup, Mortier, Cosson, Genetics, 2005

     Guillot, Mortier, Estoup, Geneland : A program for landscape
     genetics. Molecular Ecology  Notes, 2005.

_S_e_e _A_l_s_o:

     'PlotTessellation'

