genclone
computes main genetic/genotypic diversity/richness indices.
GenClone(data1, haploid = FALSE, coords = NULL, vecpop = NULL, listMLL = NULL,
nbrepeat = NULL, bar = FALSE)
a Rclone
table with one allele per column.
logical, option, haploid
indicates the ploidy level of data1
.
a table with coordinates of every units in data1
.
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
.
numeric, option, if pvalue = TRUE
, nbrepeat
is the number of
resampling to enable pvalues computation.
option, a custom list of MLL.
option, if TRUE
, displays a progression bar.
GenClone
returns a table with:
N, the number of units in data1,
Lineage, MLG or MLL,
nb_L, the number of MLG/MLL,
nb_all, the mean number of alleles,
SE, the standard error of nb_all
,
Fis, on ramets if diploid data
pval_2sides, the two-sided p-value of Fis
if nbrepeat
,
Fis_WR, on genets if diploid data
pval_2sides, the two-sided p-value of Fis_WR
if nbrepeat
,
R, the clonal diversity index (Dorken & Eckert 2001; Ellstrand & Roose 1987),
Pareto_index, the index of Pareto
Sp_Loiselle, Sp
index computed on ramets with Loiselle kinship results
used to quantify Spatial Genetic Structure (Vekemans and Hardy, 2004)
pval_2sides, the two-sided p-value of Sp_Loiselle
if nbrepeat
,
Sp_Ritland, Sp
index computed on ramets with Ritland kinship results
used to quantify SGS
pval_2sides, the two-sided p-value of Sp_Ritland
if nbrepeat
,
Sp_L_WR, Sp
index computed on genets with Loiselles kinship results
used to quantify SGS
pval_2sides, the two-sided p-value of Sp_L_WR
if nbrepeat
,
Sp_R_WR, , Sp
index computed on genets with Ritland kinship results
used to quantify SGS
pval_2sides, the two-sided p-value of Sp_R_WR
if nbrepeat
,
H'', the Shannon-Wiener index estimator (Pielou 1966),
J', the Pielou evenness index(Pielou 1975),
D', the Simpson complement unbiased (Pielou 1969 ; Gini 1912 ; Peet 1974),
V, the Simpson complement index (Hurlbert 1971 ; Fager 1972),
Hill, the reciprocal of Simpson index unbiased (Hurlbert 1971 ; Hill 1973).
GenClone
returns results of several functions of RClone
:
a summary of MLG_tab
, Fis
on ramets and genets with pvalues
(resample the population nbrepeat
times, with simulated sexual events),
B_Pareto
from Pareto_index
, Sp
from
autocorrelation
and indexes from clonal_index
.
If no coordinate at all are available, let coords = NULL
as it or
create a table with always the same number (i.e. "999", "-1", etc.).
If coordinates are available for some populations only, for the population
with missing coordinates: replace all the coordinates by the same number,
as "999".
GenClone
cannot handle mix situation with missing coordinates
only for some units of the population.
Arnaud-Haond et al., 2007, Standardizing methods to address clonality in population studies.
# NOT RUN {
data(posidonia)
data(coord_posidonia)
#GenClone(posidonia) #without coordinates
#GenClone(posidonia, coords = coord_posidonia) #with coordinates
#GenClone(posidonia, coords = coord_posidonia, nbrepeat = 1000)
##time consuming
# }
Run the code above in your browser using DataLab