psex: Probability of originating from distinct sexual events

In RClone: Partially Clonal Populations Analysis

Description

psex and psex_Fis compute the probability that repeated genotypes originate from distinct sexual events (i.e. being different genets and not ramets of the same MLG), with or without taking account of H-W equilibrium departures.

Usage

psex(data1, haploid = FALSE, vecpop = NULL, genet = FALSE, RR = FALSE, 
		MLGsim = FALSE,  nbrepeat = NULL, bar = FALSE)
psex_Fis(data1, vecpop = NULL, genet = FALSE, RR = FALSE, MLGsim = FALSE, 
		nbrepeat = NULL, bar = FALSE)

Arguments

data1	a Rclone table with one allele per column.
haploid	logical, option, haploid indicates the ploidy level of data1. Not edible for psex_Fis.
vecpop	vector, option, vecpop indicates the population name of each unit of data1, if data1 contains several populations. If data1 contains only one population, leave vecpop = NULL.
genet	option, if genet = TRUE, computes pgen on genet level.
RR	option, if RR = TRUE, computes pgen with Round-Robin method.
MLGsim	option, the method of psex calculation (see details).
nbrepeat	option, numeric, the population is simulated nbrepeat times, based on frequency values.
bar	option, if TRUE, a progression bar appears.

Details

We strongly recommand to use RR = TRUE option to compute allelic frequencies for clonal data. Otherwise, we let the options to work with frequencies at genet level (genet = TRUE) or ramet level (RR = FALSE and genet = FALSE).

if MLGsim = TRUE, psex are computed as probability for two units to be derived from distinct sexual reproductive event to be C(N,2) (Stenberg et al. 2003).

If MLGsim = FALSE, psex are computed with more conservative C(n,1) (Parks & Werth 1993) with n, "number of separated fragments with identical genotype to some previously encountered ramet".

The pvalue method calculation is largely inspired from MLGsim (Stenberg et al., 2003) and MLGsim2.0 (Ivens et al., 2012), with authors agreements.
For each repeat, a population is simulated with allelic frequencies.
If clones occurred, a simulated psex is computed and kept in memory.
At the end, a distribution of sim psex is constructed and p-value is computed as upper p-value (Monte Carlo).

psex and psex_Fis could be time consuming with a certain number of repeats.

Values must differ from MLGsim and MLGsim2.0 because of Round-Robin frequencies and Fis calculation (see freq_RR and Fis).

Value

For one population:

• if nbrepeat is not provided, a table with psex values,

• if nbrepeat is provided, a list of a table with psex values and p-values and a vector of sim psex.

If data1 is a multi-population table (vecpop != NULL), a list of either tables/tables and vectors for each population.

Warning

If sim_psex are less than 100, a warning message pops, as clones are not necessarily generated each simulation.

If no repeated genotype is generated during simulations, a warning message pops as well.

Author(s)

Creator/Author: Diane Bailleul <[e-mail]>
Author: Sophie Arnaud-Haond <[e-mail]>
Contributor: Solenn Stoeckel

The R implementation of RClone was written by Diane Bailleul.

The design was inspired by GenClone program described in Arnaud-Haond & Belkhir (2007).

References

Stenberg et al., 2003, MLGsim: a program for detecting clones using a simulation approach.

Arnaud-Haond et al., 2007, Standardizing methods to address clonality in population studies.

Ivens, A.B.F., van de Sanden, M. and Bakker, J. MLGsim 2.0: updated software for detecting clones from micro satellite data using a simulation approach.
In: The Evolutionary Ecology of Mutualism. PhD Thesis, 2012, University of Groningen. Pg 107-111
http://www.rug.nl/research/institute-evolutionary-life-sciences/tres/downloads for MLGsim 2.0.

See Also

Fis, freq_RR, pgen and pgen_Fis

Examples

data(posidonia)

psex(posidonia, RR = TRUE)
psex(posidonia, RR = TRUE, MLGsim = TRUE)
#psex(posidonia, RR = TRUE, nbrepeat = 1000, bar = TRUE)
##time consuming

Example output

   genet                 psex
1                            
2                            
3                            
4      3  0.00190284159898287
5      3 1.76851132496336e-06
6      3 1.06767920426143e-09
7      3  4.7070966236469e-13
8                            
9                            
10                           
11                           
12    11  0.00199659872286699
13                           
14                           
15    11 1.94726485008159e-06
16                           
17                           
18    17  0.00025371601828402
19    17 3.13892261776203e-08
20                           
21                           
22                           
23                           
24    23  0.00157820543710591
25                           
26    23 1.21615286466571e-06
27                           
28    20 0.000833166996686032
29                           
30    23  6.0875218010658e-10
31                           
32                           
33                           
34                           
35                           
36                           
37                           
38                           
39                           
40                           
   genet                 psex
1                            
2                            
3      3 1.61531120764607e-16
4      3 1.61531120764607e-16
5      3 1.61531120764607e-16
6      3 1.61531120764607e-16
7      3 1.61531120764607e-16
8                            
9                            
10                           
11    11 1.23363587525764e-09
12    11 1.23363587525764e-09
13                           
14                           
15    11 1.23363587525764e-09
16                           
17    17 2.52255742863493e-12
18    17 2.52255742863493e-12
19    17 2.52255742863493e-12
20    20 3.38688830725152e-07
21                           
22                           
23    23 2.22521655757152e-13
24    23 2.22521655757152e-13
25                           
26    23 2.22521655757152e-13
27                           
28    20 3.38688830725152e-07
29                           
30    23 2.22521655757152e-13
31                           
32                           
33                           
34                           
35                           
36                           
37                           
38                           
39                           
40                           

RClone documentation built on May 15, 2021, 5:07 p.m.