diff --git a/docs/articles/index.html b/docs/articles/index.html index 2466648a..d11c2a15 100644 --- a/docs/articles/index.html +++ b/docs/articles/index.html @@ -96,7 +96,7 @@
diff --git a/docs/articles/mlg.html b/docs/articles/mlg.html index 1904d990..7326239d 100644 --- a/docs/articles/mlg.html +++ b/docs/articles/mlg.html @@ -5,7 +5,7 @@ -Analysis of Multilocus Genotypes and Lineages in poppr 2.6.1 • poppr +Analysis of Multilocus Genotypes and Lineages in poppr 2.7.1 • poppr @@ -75,10 +75,10 @@
diff --git a/docs/articles/poppr_manual.html b/docs/articles/poppr_manual.html index 901578f7..df3f17ce 100644 --- a/docs/articles/poppr_manual.html +++ b/docs/articles/poppr_manual.html @@ -5,7 +5,7 @@ -Data import and manipulation in poppr version 2.6.1 • poppr +Data import and manipulation in poppr version 2.7.1 • poppr @@ -75,10 +75,10 @@
@@ -135,7 +135,7 @@

citation(package = "poppr")
## 
 ## To cite poppr in publications or presentations, please specify
-## poppr version 2.6.1 and with the following citation:
+## poppr version 2.7.1 and with the following citation:
 ## 
 ##   Kamvar ZN, Tabima JF, Grünwald NJ. (2014) Poppr: an R package
 ##   for genetic analysis of populations with clonal, partially
@@ -201,12 +201,12 @@ 

## Loading required package: adegenet
## Loading required package: ade4
## 
-##    /// adegenet 2.1.0 is loaded ////////////
+##    /// adegenet 2.1.1 is loaded ////////////
 ## 
 ##    > overview: '?adegenet'
 ##    > tutorials/doc/questions: 'adegenetWeb()' 
 ##    > bug reports/feature requests: adegenetIssues()
-
## This is poppr version 2.6.1. To get started, type package?poppr
+
## This is poppr version 2.7.1. To get started, type package?poppr
 ## OMP parallel support: available
x <- getfile()

A pop up window will appear like this1:

@@ -1655,13 +1655,13 @@

This allows you to produce publication quality graphs directly in R. Please see Hadley Wickham“s ggplot2 package for more details (Wickham 2009). Note that if you don”t like using ggplot2, you can access the data in the ggplot2 object and plot the data yourself:

head(p$data)
## # A tibble: 6 x 4
-##   Population MLG    count order 
-##   <fctr>     <chr>  <int> <fctr>
-## 1 Athena     MLG.20     9 1     
-## 2 Athena     MLG.66     5 2     
-## 3 Athena     MLG.14     3 3     
-## 4 Athena     MLG.35     3 4     
-## 5 Athena     MLG.13     2 5     
+##   Population MLG    count order
+##   <fct>      <chr>  <int> <fct>
+## 1 Athena     MLG.20     9 1    
+## 2 Athena     MLG.66     5 2    
+## 3 Athena     MLG.14     3 3    
+## 4 Athena     MLG.35     3 4    
+## 5 Athena     MLG.13     2 5    
 ## 6 Athena     MLG.16     2 6

diff --git a/docs/authors.html b/docs/authors.html index d3edfd49..22222878 100644 --- a/docs/authors.html +++ b/docs/authors.html @@ -170,6 +170,10 @@

Authors

Brian J. Knaus. Contributor.
0000-0003-1665-4343

+
  • +

    Patrick G. Meirmans. Contributor. +
    0000-0002-6395-8107

    +
  • Niklaus J. Grunwald. Thesis advisor.
    0000-0003-1656-7602

    diff --git a/docs/index.html b/docs/index.html index 419ff74e..a8401249 100644 --- a/docs/index.html +++ b/docs/index.html @@ -206,11 +206,11 @@

    Data import and manipulation -vignette("poppr_manual", "poppr") +vignette("poppr_manual", "poppr") Multilocus Genotype Analysis -vignette("mlg", "poppr") +vignette("mlg", "poppr") @@ -261,8 +261,9 @@

    Developers

    Dev status

      -
    • Build Status
    • -
    • Coverage Status
    • +
    • Build Status
    • +
    • AppVeyor build status
    • +
    • Coverage Status
    • CRAN version
  • diff --git a/docs/news/index.html b/docs/news/index.html index f4d486ed..af49df97 100644 --- a/docs/news/index.html +++ b/docs/news/index.html @@ -103,6 +103,51 @@

    Change log All releases

    +
    +

    +poppr 2.7.1

    +
    +

    +MISC

    +
      +
    • Missing documentation for poppr.amova() has been added.
    • +
    • Polysat is now listed in imports.
    • +
    +
    +
    +
    +

    +poppr 2.7.0

    +
    +

    +NEW FUNCTIONS

    +
      +
    • make_haplotypes() will split your data into pseudo-haplotypes for use in AMOVA-like analyses. This was a previously internal function, but has been promoted to a user-facing function in this version.

    • +
    • as.genambig() will convert genind/genclone objects to Polysat’s “genambig” class. Note that polysat must be installed for this to work.

    • +
    +
    +
    +

    +ALGORITHMIC CHANGE

    +
      +
    • AMOVA will now default to using euclidean distance. This affects all calculations made with within = FALSE or filter = TRUE without a user-supplied distance. This will not have affect those with haploid or diploid data using within = TRUE. The dissimilarity distance is equivalent to a squared euclidean distance for haploid genotypes, but not for any higher ploidy. Those using filter = TRUE without specifying a distance should use a euclidean threshold. This should not be an issue for those who simply want to group isolates with missing data, however as a zero distance is the same for euclidean and dissimilarity. Thanks goes to Patrick Meirmans for alerting me to this error.
    • +
    +
    +
    +

    +NEW FEATURES

    +
      +
    • AMOVA will now calculate within-individual variance for polyploid data.
    • +
    +
    +
    +

    +MISC

    +
      +
    • printing of AMOVA will now better handle any changes in methods from pegas or ade4.
    • +
    +
    +

    poppr 2.6.1

    @@ -117,16 +162,16 @@

    poppr 2.6.0

    -
    +

    -NEW FUNCTIONS

    +NEW FUNCTIONS

    -
    +

    -NEW FEATURES

    +NEW FEATURES
    • The function resample.ia() now can resample individuals weighted by their Psex value.
    • The minimum spanning networks will now scale nodes by area instead of radius. This gives a more accurate picture of the differences between MLGs. See https://github.com/grunwaldlab/poppr/issues/154 for details.
    • @@ -152,9 +197,9 @@

    • The minimum version of igraph has been set to 1.0.0.
    -
    +

    -MISC

    +MISC
    • The MSN is now plotted last in plot_poppr_msn() so additional legends can be added if necessary.
    @@ -163,9 +208,9 @@

    poppr 2.5.0

    -
    +

    -ALGORITHMIC CHANGE

    +ALGORITHMIC CHANGE

    • Identified in https://github.com/grunwaldlab/poppr/issues/139, Bruvo’s distance will now consider all possible combinations of ordered alleles in the calculation under the genome addition and loss models for missing data. This will affect those who have polyploid data that contain more than one missing allele at any genotype
    @@ -187,9 +232,9 @@

  • A bug in read.genalex() where removed samples would have incorrect strata labels was fixed. Thanks to Hernán Dario Capador-Barreto for identifying it. See https://github.com/grunwaldlab/poppr/issues/147.
  • -
    +

    -MISC

    +MISC
    • The internal plotting function for mlg.table now uses tidy evaluation for dplyr versions > 0.5.0
    • The package reshape2 was removed from imports and replaced with base functions (see https://github.com/grunwaldlab/poppr/issues/144 for details)
    • @@ -218,17 +263,17 @@

      poppr 2.4.0

      -
      +

      -NEW FUNCTIONS

      +NEW FUNCTIONS

      • jack.ia() will randomly jackknife your sample to a specified n (default is the number of MLG), and calculate the index of association over multiple iterations, giving a distribution of possible values at a given sample size.
    -
    +

    -NEW FEATURES

    +NEW FEATURES
    • The function mlg.table() gains new parameters, “color” and “background”. The “color” parameter will create a single barplot with colors representing populations while the “background” parameter will create a background plot showing the abundance of MLGs across populations within the facets.
    • The function win.ia() will now take into consideration chromosomal coordinates when constructing windows. It has additionally acquired a new parameter chromosome_buffer, which allows the user to specify whether or not the window should be limited to within chromosomes.
    • @@ -262,9 +307,9 @@

    • PROTECT statements were placed around allocation statements. For details, see https://github.com/grunwaldlab/poppr/issues/133.
    -
    +

    -MISC

    +MISC
    -
    +

    -NEW FEATURES

    +NEW FEATURES
    • There is now a plot parameter for the genotype curve to enable or suppress plotting.
    • Progress bars are now automatically suppressed when running non-interactively. to turn them on when running non-interactively, use options(poppr.debug = TRUE).
    • The progress bar for ia() and poppr() will now show estimated time. This is from dplyr’s progress_estimated().
    -
    +

    -MISC

    +MISC
    • The hist argument in the ia() is deprecated in favor of plot.
    • The x axis for the genotype_curve() plot is now numeric, allowing you to fit a smoothing function over the points without having to use the hack geom_smooth(aes(group = 1)). This is thanks to Kara Woo for pointing this out on twitter (https://twitter.com/kara_woo/status/783336540407685120).
    • @@ -325,18 +370,18 @@

    -
    +

    -NEW FEATURES

    +NEW FEATURES
    -
    +

    -MISC

    +MISC
    -
    +

    -NEW FEATURES

    +NEW FEATURES
    -
    +

    -MISC

    +MISC
    • Documentation for mlg.filter() was clarified and updated with more examples.
    • The vignette “Migration from poppr version 1” has been removed.
    • @@ -394,9 +439,9 @@

      poppr 2.1.1

      -
      +

      -NEW FEATURES

      +NEW FEATURES

    -
    +

    -MISC

    +MISC
    • mlg.filter() now defaults to using diss.dist() @@ -449,9 +494,9 @@

    • The internal code for the genotype_curve has been implemented in C for a 10x increase in speed.
    -
    +

    -NEW FEATURES

    +NEW FEATURES
    • poppr.msn, bruvo.msn, and plot_poppr_msn gain the ability to take character vectors for color palettes. See issue #55 (https://github.com/grunwaldlab/poppr/issues/55) for details.
    • @@ -475,9 +520,9 @@

      poppr.amova can now perform amova using the pegas implementation.

    -
    +

    -NEW FUNCTIONS

    +NEW FUNCTIONS
    • rrmlg will calculate round-robin multilocus genotypes for each locus.
    • @@ -496,9 +541,9 @@

    • because we’re through being cool.
    -
    +

    -MISC

    +MISC
    • Documentation for genclone and snpclone classes are more coherent.
    • Accessors added for internal MLG objects (for developers).
    • @@ -546,9 +591,9 @@

      read.genealex can now correctly import missing data for diploids (issue #42)

    -
    +

    -NEW FEATURES

    +NEW FEATURES
    • Startup message now tells you if poppr was compiled with OMP support.
    @@ -593,9 +638,9 @@

  • refreshing!
  • -
    +

    -NEW FEATURES

    +NEW FEATURES
    • The default plot for the index of association will now be a single histogram. The user has the option to visualize the standardized index of association (index = "rbarD", default) or the classic index of association (index = "Ia"). If the user uses the function ia with the argument valuereturn = TRUE, then the resulting object can be plotted with the plot function.
    • The function poppr will now plot all populations in a single faceted plot instead of one plot per population.
    • @@ -624,9 +669,9 @@

      genind2genalex gains the ability to selectively write different strata.

    -
    +

    -NEW FUNCTIONS

    +NEW FUNCTIONS
    • mlg.filter will contract multilocus genotypes given a genetic distance and threshold using one of three algorithms. It can report statistics such as the multilocus genotypes returned, the number of samples within each multilocus genotype, the thresholds at which multilocus genotypes were collapsed, and the genetic distance matrix that represents the new multilocus genotypes.
    • @@ -710,9 +755,9 @@

      poppr 1.1.3

      -
      +

      -NEW FEATURES

      +NEW FEATURES

      • new arguments to plot_poppr_msn to allow for easier manipulation of node sizes and of labeling
      • read.genalex can now take read text connections as input. Addresses issue #8
      • @@ -735,9 +780,9 @@

    -
    +

    -MISC

    +MISC
    • info_table will print a discrete scale as opposed to colorbar when type = “ploidy”
    • @@ -773,9 +818,9 @@

      poppr 1.1.0

      -
      +

      -NEW FEATURES

      +NEW FEATURES

      • Polyploids with ambiguous genotypes are now supported in poppr. See documentation for recode_polyploids for details.
      • Calculations of Bruvo’s distance now features correction for partial missing data utilizing genome addition and genome loss models as presented in Bruvo et al. 2004.
      • @@ -800,9 +845,9 @@

      • The genclone object is a new extension of the genind object from adegenet. This object contains slots containing population hierarchies and multilocus genotype definitions and will work with all analyses in adegenet and poppr.
    -
    +

    -NEW FUNCTIONS

    +NEW FUNCTIONS
    • [get,set,name,split,add]hierarchy - functions that will manipulate the hierarchy slot in a genclone object utilizing hierarchical formulae as arguments for simplification.
    • @@ -867,9 +912,9 @@

    • The function poppr will no longer return rounded results, but rather is printed with three significant digits.
    -
    +

    -MISC

    +MISC
    • Added unit tests.
    • The poppr user manual has been shortened to only include instructions on data manipulation.
    • @@ -905,9 +950,9 @@

    • Fixed bug for users who have downloaded ape version 3.1 or higher where bruvo.boot would throw an error.
    -
    +

    -MISC

    +MISC
    • Updated citation information.
    @@ -939,9 +984,9 @@

  • Changes to bruvo.boot allow for ever so slightly faster bootstrapping.
  • -
    +

    -MISC

    +MISC
    • Permutations for I_A and \bar{r}_d are now visualized as a progress bar as opposed to dots.
    @@ -959,9 +1004,9 @@

    (1994) normalization for NJ trees.

    -
    +

    -MISC

    +MISC
    • github repository for poppr has changed from github.com/poppr/poppr to github.com/grunwaldlab/poppr
    @@ -978,9 +1023,9 @@

  • Utilized rmultinom function to increase speed of bootstrap sampling methods for shufflepop and ia.
  • -
    +

    -NEW FEATURES

    +NEW FEATURES
    • Function informloci will remove phylogenetically uninformative loci.
    @@ -1020,9 +1065,9 @@

  • Expanded installation section to include installation instructions from github.
  • -
    +

    -MISC

    +MISC
    • internal permutation algorithm no longer lists permutations in reverse order
    @@ -1047,9 +1092,9 @@

  • Input values that are not multiples of the specified repeat length for Bruvo’s distance are now rounded (as opposed to being forced as integers).
  • -
    +

    -MISC

    +MISC
    • Vignette updated for aesthetics and to reflect algorithmic changes.
    @@ -1058,9 +1103,9 @@

    poppr 1.0.0

    -
    +

    -MISC

    +MISC

    • Poppr has been confirmed to work on Linux, Mac, and Windows systems with R 3.0.0.
    • Vignette poppr_manual now has cross-references to different sections.
    • @@ -1078,9 +1123,9 @@

      poppr 0.4.1

      -
      +

      -NEW FEATURES

      +NEW FEATURES

      • getfile has a new argument, “combine”, which will automatically add the path to the list of files, so they can be read without switching working directory.
      • @@ -1111,9 +1156,9 @@

        read.genalex can now correctly recognize regional formatting without an extra column.

    -
    +

    -NEW FEATURES

    +NEW FEATURES
    • read.genalex will now be able to take in a file that is formatted with both regional and geographic data.
    • @@ -1156,17 +1201,17 @@

      poppr 0.3

      -
      +

      -NEW FUNCTIONS

      +NEW FUNCTIONS

      • poppr.msn will draw a minimum spanning network for any distance matrix derived from your data set.
    -
    +

    -NEW FEATURES

    +NEW FEATURES
    • vignette now has sections describing poppr.msn, diss.dist, greycurve, and a section discussing how to export graphics.
    @@ -1184,9 +1229,9 @@

    poppr 0.2.2

    -
    +

    -NEW FEATURES

    +NEW FEATURES

    • index of association distributions will now feature a rug plot at the bottom as a better way to visualize the distribution of the index of association from the shuffled data sets.
    @@ -1195,9 +1240,9 @@

    poppr 0.2.1

    -
    +

    -NEW FUNCTIONS

    +NEW FUNCTIONS

    • diss.dist will produce a distance matrix based on discreet distances.
    • @@ -1205,9 +1250,9 @@

      greycurve will produce a grey scale adjusted to user-supplied parameters. This will be useful for future minimum spanning network functions.

    -
    +

    -NEW FEATURES

    +NEW FEATURES
    • bruvo.msn can now adjust the edge grey level to be weighted toward either closely or distantly weighted individuals.
    • @@ -1227,9 +1272,9 @@

      poppr 0.2

      -
      +

      -NEW FEATURES

      +NEW FEATURES

      • Added NEWS file and will now be incrementing version number (3/15/2013)
      @@ -1249,6 +1294,8 @@

      Contents

      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Warning: Infinite values detected.
      #> Running bootstraps: 100 / 100 -#> Calculating bootstrap values... done.
      +#> Calculating bootstrap values... done.
      set.seed(9999) # Generate a tree using custom distance bindist <- function(x) dist(tab(x), method = "binary") binnan <- aboot(nan9, dist = bindist)
      #> Running bootstraps: 100 / 100 -#> Calculating bootstrap values... done.
      +#> Calculating bootstrap values... done.
      # NOT RUN { # Distances from other packages. # diff --git a/docs/reference/bitwise.dist.html b/docs/reference/bitwise.dist.html index 679c58e8..eb16c5a2 100644 --- a/docs/reference/bitwise.dist.html +++ b/docs/reference/bitwise.dist.html @@ -196,7 +196,7 @@

      Examp #> @other: a list containing: ancestral.pops #>

    # Assess fraction of different alleles (finer measure, usually the most sensible) system.time(xd <- bitwise.dist(x))
    #> user system elapsed -#> 0.002 0.000 0.002
    xd
    #> 1 2 3 4 5 6 7 8 9 +#> 0.003 0.001 0.002
    xd
    #> 1 2 3 4 5 6 7 8 9 #> 2 0.2105 #> 3 0.2195 0.2180 #> 4 0.3950 0.3825 0.3855 @@ -208,7 +208,7 @@

    Examp #> 10 0.3980 0.3925 0.4005 0.2180 0.2035 0.2220 0.2185 0.2125 0.2140

    # Assess fraction of different loci (coarse measure) system.time(xdt <- bitwise.dist(x, differences_only = TRUE))
    #> user system elapsed -#> 0.001 0.000 0.002
    xdt
    #> 1 2 3 4 5 6 7 8 9 +#> 0.067 0.001 0.086
    xdt
    #> 1 2 3 4 5 6 7 8 9 #> 2 0.375 #> 3 0.394 0.385 #> 4 0.564 0.546 0.547 diff --git a/docs/reference/boot.ia.html b/docs/reference/boot.ia.html index 40caa3ff..7358c27b 100644 --- a/docs/reference/boot.ia.html +++ b/docs/reference/boot.ia.html @@ -161,7 +161,7 @@

    See a

    Examples

    data(Pinf) -boot.ia(Pinf, reps = 99)
    #> |============== | 27% ~0 s remaining |================================ | 60% ~0 s remaining |================================================= | 91% ~0 s remaining
    #> Ia rbarD +boot.ia(Pinf, reps = 99)
    #> |= | 2% ~3 s remaining |================= | 32% ~0 s remaining |==================================== | 68% ~0 s remaining |===================================================== | 99% ~0 s remaining
    #> Ia rbarD #> 1 0.6430049 0.07032330 #> 2 0.5705043 0.06235351 #> 3 0.5774512 0.06340112 diff --git a/docs/reference/bruvo.boot-1.png b/docs/reference/bruvo.boot-1.png index 2b1811d2..49c4506f 100644 Binary files a/docs/reference/bruvo.boot-1.png and b/docs/reference/bruvo.boot-1.png differ diff --git a/docs/reference/bruvo.boot.html b/docs/reference/bruvo.boot.html index ae492450..4862c683 100644 --- a/docs/reference/bruvo.boot.html +++ b/docs/reference/bruvo.boot.html @@ -225,7 +225,7 @@

    Examp #> (note: calculation of node labels can take a while even after the progress bar is full) #> #> Running bootstraps: 100 / 100 -#> Calculating bootstrap values... done.

    #> +#> Calculating bootstrap values... done.
    #> #> Phylogenetic tree with 10 tips and 9 internal nodes. #> #> Tip labels: diff --git a/docs/reference/bruvo.msn-1.png b/docs/reference/bruvo.msn-1.png index c9b43205..313c9bc8 100644 Binary files a/docs/reference/bruvo.msn-1.png and b/docs/reference/bruvo.msn-1.png differ diff --git a/docs/reference/bruvo.msn.html b/docs/reference/bruvo.msn.html index 935abd36..e10fb7c6 100644 --- a/docs/reference/bruvo.msn.html +++ b/docs/reference/bruvo.msn.html @@ -311,11 +311,11 @@

    Examp # View populations 8 and 9 with default colors. bruvo.msn(nancycats, replen = rep(2, 9), sublist=8:9, vertex.label="inds", - vertex.label.cex=0.7, vertex.label.dist=0.4)

    #> $graph -#> IGRAPH 7f82159 UNW- 19 18 -- + vertex.label.cex=0.7, vertex.label.dist=0.4)
    #> $graph +#> IGRAPH 9d04bb2 UNW- 19 18 -- #> + attr: name (v/c), size (v/n), shape (v/c), pie (v/x), pie.color #> | (v/x), label (v/c), weight (e/n), color (e/c), width (e/n) -#> + edges from 7f82159 (vertex names): +#> + edges from 9d04bb2 (vertex names): #> [1] N43 --N93 N92 --N112 N94 --N98 N95 --N96 N95 --N97 N98 --N99 #> [7] N98 --N100 N98 --N97 N98 --N111 N100--N108 N93 --N97 N104--N107 #> [13] N105--N109 N106--N107 N106--N109 N107--N108 N107--N112 N111--N113 diff --git a/docs/reference/coercion-methods.html b/docs/reference/coercion-methods.html index 307d6bce..83aca1f2 100644 --- a/docs/reference/coercion-methods.html +++ b/docs/reference/coercion-methods.html @@ -110,7 +110,9 @@

    Switch between genind and genclone objects.

    as.genclone(x, ..., mlg, mlgclass = TRUE) -genclone2genind(x)
    +genclone2genind(x) + +as.genambig(x)

    Arguments

    @@ -142,6 +144,8 @@

    Details

    genclone2genind will remove the mlg slot from the genclone object, creating a genind object.

    +

    as.genambig will convert a genind or genclone object to a polysat genambig +class.

    See also

    @@ -166,7 +170,9 @@

    Examp #> // Optional content #> @pop: population of each individual (group size range: 90-97) #> @other: a list containing: population_hierarchy -#>
    Aeut.gc <- as.genclone(Aeut) +#>
    +# Conversion to genclone -------------------------------------------------- +Aeut.gc <- as.genclone(Aeut) Aeut.gc
    #> #> This is a genclone object #> ------------------------- @@ -179,7 +185,9 @@

    Examp #> Population information: #> #> 0 strata. -#> 2 populations defined - Athena, Mt. Vernon

    Aeut.gi <- genclone2genind(Aeut.gc) +#> 2 populations defined - Athena, Mt. Vernon
    +# Conversion to genind ---------------------------------------------------- +Aeut.gi <- genclone2genind(Aeut.gc) Aeut.gi
    #> /// GENIND OBJECT ///////// #> #> // 187 individuals; 56 loci; 56 alleles; size: 65.8 Kb @@ -194,7 +202,22 @@

    Examp #> // Optional content #> @pop: population of each individual (group size range: 90-97) #> @other: a list containing: population_hierarchy -#>

    data(nancycats) +#>
    +# Conversion to polysat's "genambig" class -------------------------------- +if (require("polysat")) { + data(Pinf) + Pinf.gb <- as.genambig(Pinf) + summary(Pinf.gb) +}
    #> Loading required package: polysat
    #> Dataset with allele copy number ambiguity. +#> Insert dataset description here. +#> Number of missing genotypes: 10 +#> 86 samples, 11 loci. +#> 2 populations. +#> Ploidies: 2 3 NA +#> Length(s) of microsatellite repeats: NA
    +data(nancycats) + +# Conversion to bootgen for random sampling of loci ----------------------- nan.bg <- new("bootgen", nancycats[pop = 9]) nan.bg
    #> An object of class "bootgen" #> Slot "type": @@ -449,7 +472,9 @@

    Examp #> #> Slot "call": #> NULL -#>

    nan.gid <- bootgen2genind(nan.bg) +#>
    +# Conversion back to genind ----------------------------------------------- +nan.gid <- bootgen2genind(nan.bg) nan.gid
    #> /// GENIND OBJECT ///////// #> #> // 9 individuals; 9 loci; 108 alleles; size: 23.5 Kb @@ -464,7 +489,8 @@

    Examp #> @call: .local(.Object = .Object, tab = ..1) #> #> // Optional content -#> - empty -

    +#> - empty -
    +
    #> #> Confidence Intervals have been centered around observed statistic. -#> Please see ?diversity_ci for details.
    #> $obs +#> Please see ?diversity_ci for details.
    #> $obs #> Index #> Pop H G lambda E.5 #> South America 3.267944 23.29032 0.9570637 0.8825297 diff --git a/docs/reference/filter_stats-1.png b/docs/reference/filter_stats-1.png index 973419b2..47b7f3ef 100644 Binary files a/docs/reference/filter_stats-1.png and b/docs/reference/filter_stats-1.png differ diff --git a/docs/reference/filter_stats.html b/docs/reference/filter_stats.html index ea53a147..72d8ca17 100644 --- a/docs/reference/filter_stats.html +++ b/docs/reference/filter_stats.html @@ -201,7 +201,7 @@

    See a

    Examples

    data(Pinf) pinfreps <- fix_replen(Pinf, c(2, 2, 6, 2, 2, 2, 2, 2, 3, 3, 2)) -filter_stats(Pinf, distance = bruvo.dist, replen = pinfreps, plot = TRUE, threads = 1L)
    +filter_stats(Pinf, distance = bruvo.dist, replen = pinfreps, plot = TRUE, threads = 1L)
    #> Calculating genotypes for 1/8 loci. Completed iterations: 1% Calculating genotypes for 1/8 loci. Completed iterations: 2% Calculating genotypes for 1/8 loci. Completed iterations: 3% Calculating genotypes for 1/8 loci. Completed iterations: 4% Calculating genotypes for 1/8 loci. Completed iterations: 5% Calculating genotypes for 1/8 loci. Completed iterations: 6% Calculating genotypes for 1/8 loci. Completed iterations: 7% Calculating genotypes for 1/8 loci. Completed iterations: 8% Calculating genotypes for 1/8 loci. Completed iterations: 9% Calculating genotypes for 1/8 loci. Completed iterations: 10% Calculating genotypes for 1/8 loci. Completed iterations: 11% Calculating genotypes for 1/8 loci. Completed iterations: 12% Calculating genotypes for 1/8 loci. Completed iterations: 13% Calculating genotypes for 1/8 loci. Completed iterations: 14% Calculating genotypes for 1/8 loci. Completed iterations: 15% Calculating genotypes for 1/8 loci. Completed iterations: 16% Calculating genotypes for 1/8 loci. Completed iterations: 17% Calculating genotypes for 1/8 loci. Completed iterations: 18% Calculating genotypes for 1/8 loci. Completed iterations: 19% Calculating genotypes for 1/8 loci. Completed iterations: 20% Calculating genotypes for 1/8 loci. Completed iterations: 21% Calculating genotypes for 1/8 loci. Completed iterations: 22% Calculating genotypes for 1/8 loci. Completed iterations: 23% Calculating genotypes for 1/8 loci. Completed iterations: 24% Calculating genotypes for 1/8 loci. Completed iterations: 25% Calculating genotypes for 1/8 loci. Completed iterations: 26% Calculating genotypes for 1/8 loci. Completed iterations: 27% Calculating genotypes for 1/8 loci. Completed iterations: 28% Calculating genotypes for 1/8 loci. Completed iterations: 29% Calculating genotypes for 1/8 loci. Completed iterations: 30% Calculating genotypes for 1/8 loci. Completed iterations: 31% Calculating genotypes for 1/8 loci. Completed iterations: 32% Calculating genotypes for 1/8 loci. Completed iterations: 33% Calculating genotypes for 1/8 loci. Completed iterations: 34% Calculating genotypes for 1/8 loci. Completed iterations: 35% Calculating genotypes for 1/8 loci. Completed iterations: 36% Calculating genotypes for 1/8 loci. Completed iterations: 37% Calculating genotypes for 1/8 loci. Completed iterations: 38% Calculating genotypes for 1/8 loci. Completed iterations: 39% Calculating genotypes for 1/8 loci. Completed iterations: 40% Calculating genotypes for 1/8 loci. Completed iterations: 41% Calculating genotypes for 1/8 loci. Completed iterations: 42% Calculating genotypes for 1/8 loci. Completed iterations: 43% Calculating genotypes for 1/8 loci. Completed iterations: 44% Calculating genotypes for 1/8 loci. Completed iterations: 45% Calculating genotypes for 1/8 loci. Completed iterations: 46% Calculating genotypes for 1/8 loci. Completed iterations: 47% Calculating genotypes for 1/8 loci. Completed iterations: 48% Calculating genotypes for 1/8 loci. Completed iterations: 49% Calculating genotypes for 1/8 loci. Completed iterations: 50% Calculating genotypes for 1/8 loci. Completed iterations: 51% Calculating genotypes for 1/8 loci. Completed iterations: 52% Calculating genotypes for 1/8 loci. Completed iterations: 53% Calculating genotypes for 1/8 loci. Completed iterations: 54% Calculating genotypes for 1/8 loci. Completed iterations: 55% Calculating genotypes for 1/8 loci. Completed iterations: 56% Calculating genotypes for 1/8 loci. Completed iterations: 57% Calculating genotypes for 1/8 loci. Completed iterations: 58% Calculating genotypes for 1/8 loci. Completed iterations: 59% Calculating genotypes for 1/8 loci. Completed iterations: 60% Calculating genotypes for 1/8 loci. Completed iterations: 61% Calculating genotypes for 1/8 loci. Completed iterations: 62% Calculating genotypes for 1/8 loci. Completed iterations: 63% Calculating genotypes for 1/8 loci. Completed iterations: 64% Calculating genotypes for 1/8 loci. Completed iterations: 65% Calculating genotypes for 1/8 loci. Completed iterations: 66% Calculating genotypes for 1/8 loci. Completed iterations: 67% Calculating genotypes for 1/8 loci. Completed iterations: 68% Calculating genotypes for 1/8 loci. Completed iterations: 69% Calculating genotypes for 1/8 loci. Completed iterations: 70% Calculating genotypes for 1/8 loci. Completed iterations: 71% Calculating genotypes for 1/8 loci. Completed iterations: 72% Calculating genotypes for 1/8 loci. Completed iterations: 73% Calculating genotypes for 1/8 loci. Completed iterations: 74% Calculating genotypes for 1/8 loci. Completed iterations: 75% Calculating genotypes for 1/8 loci. Completed iterations: 76% Calculating genotypes for 1/8 loci. Completed iterations: 77% Calculating genotypes for 1/8 loci. Completed iterations: 78% Calculating genotypes for 1/8 loci. Completed iterations: 79% Calculating genotypes for 1/8 loci. Completed iterations: 80% Calculating genotypes for 1/8 loci. Completed iterations: 81% Calculating genotypes for 1/8 loci. Completed iterations: 82% Calculating genotypes for 1/8 loci. Completed iterations: 83% Calculating genotypes for 1/8 loci. Completed iterations: 84% Calculating genotypes for 1/8 loci. Completed iterations: 85% Calculating genotypes for 1/8 loci. Completed iterations: 86% Calculating genotypes for 1/8 loci. Completed iterations: 87% Calculating genotypes for 1/8 loci. Completed iterations: 88% Calculating genotypes for 1/8 loci. Completed iterations: 89% Calculating genotypes for 1/8 loci. Completed iterations: 90% Calculating genotypes for 1/8 loci. Completed iterations: 91% Calculating genotypes for 1/8 loci. Completed iterations: 92% Calculating genotypes for 1/8 loci. Completed iterations: 93% Calculating genotypes for 1/8 loci. Completed iterations: 94% Calculating genotypes for 1/8 loci. Completed iterations: 95% Calculating genotypes for 1/8 loci. Completed iterations: 96% Calculating genotypes for 1/8 loci. Completed iterations: 97% Calculating genotypes for 1/8 loci. Completed iterations: 98% Calculating genotypes for 1/8 loci. Completed iterations: 99% Calculating genotypes for 1/8 loci. Completed iterations: 100% Calculating genotypes for 2/8 loci. Completed iterations: 1% Calculating genotypes for 2/8 loci. Completed iterations: 2% Calculating genotypes for 2/8 loci. Completed iterations: 3% Calculating genotypes for 2/8 loci. Completed iterations: 4% Calculating genotypes for 2/8 loci. Completed iterations: 5% Calculating genotypes for 2/8 loci. Completed iterations: 6% Calculating genotypes for 2/8 loci. Completed iterations: 7% Calculating genotypes for 2/8 loci. Completed iterations: 8% Calculating genotypes for 2/8 loci. Completed iterations: 9% Calculating genotypes for 2/8 loci. Completed iterations: 10% Calculating genotypes for 2/8 loci. Completed iterations: 11% Calculating genotypes for 2/8 loci. Completed iterations: 12% Calculating genotypes for 2/8 loci. Completed iterations: 13% Calculating genotypes for 2/8 loci. Completed iterations: 14% Calculating genotypes for 2/8 loci. Completed iterations: 15% Calculating genotypes for 2/8 loci. Completed iterations: 16% Calculating genotypes for 2/8 loci. Completed iterations: 17% Calculating genotypes for 2/8 loci. Completed iterations: 18% Calculating genotypes for 2/8 loci. Completed iterations: 19% Calculating genotypes for 2/8 loci. Completed iterations: 20% Calculating genotypes for 2/8 loci. Completed iterations: 21% Calculating genotypes for 2/8 loci. Completed iterations: 22% Calculating genotypes for 2/8 loci. Completed iterations: 23% Calculating genotypes for 2/8 loci. Completed iterations: 24% Calculating genotypes for 2/8 loci. Completed iterations: 25% Calculating genotypes for 2/8 loci. Completed iterations: 26% Calculating genotypes for 2/8 loci. Completed iterations: 27% Calculating genotypes for 2/8 loci. Completed iterations: 28% Calculating genotypes for 2/8 loci. Completed iterations: 29% Calculating genotypes for 2/8 loci. Completed iterations: 30% Calculating genotypes for 2/8 loci. Completed iterations: 31% Calculating genotypes for 2/8 loci. Completed iterations: 32% Calculating genotypes for 2/8 loci. Completed iterations: 33% Calculating genotypes for 2/8 loci. Completed iterations: 34% Calculating genotypes for 2/8 loci. Completed iterations: 35% Calculating genotypes for 2/8 loci. Completed iterations: 36% Calculating genotypes for 2/8 loci. Completed iterations: 37% Calculating genotypes for 2/8 loci. Completed iterations: 38% Calculating genotypes for 2/8 loci. Completed iterations: 39% Calculating genotypes for 2/8 loci. Completed iterations: 40% Calculating genotypes for 2/8 loci. Completed iterations: 41% Calculating genotypes for 2/8 loci. Completed iterations: 42% Calculating genotypes for 2/8 loci. Completed iterations: 43% Calculating genotypes for 2/8 loci. Completed iterations: 44% Calculating genotypes for 2/8 loci. Completed iterations: 45% Calculating genotypes for 2/8 loci. Completed iterations: 46% Calculating genotypes for 2/8 loci. Completed iterations: 47% Calculating genotypes for 2/8 loci. Completed iterations: 48% Calculating genotypes for 2/8 loci. Completed iterations: 49% Calculating genotypes for 2/8 loci. Completed iterations: 50% Calculating genotypes for 2/8 loci. Completed iterations: 51% Calculating genotypes for 2/8 loci. Completed iterations: 52% Calculating genotypes for 2/8 loci. Completed iterations: 53% Calculating genotypes for 2/8 loci. Completed iterations: 54% Calculating genotypes for 2/8 loci. Completed iterations: 55% Calculating genotypes for 2/8 loci. Completed iterations: 56% Calculating genotypes for 2/8 loci. Completed iterations: 57% Calculating genotypes for 2/8 loci. Completed iterations: 58% Calculating genotypes for 2/8 loci. Completed iterations: 59% Calculating genotypes for 2/8 loci. Completed iterations: 60% Calculating genotypes for 2/8 loci. Completed iterations: 61% Calculating genotypes for 2/8 loci. Completed iterations: 62% Calculating genotypes for 2/8 loci. Completed iterations: 63% Calculating genotypes for 2/8 loci. Completed iterations: 64% Calculating genotypes for 2/8 loci. Completed iterations: 65% Calculating genotypes for 2/8 loci. Completed iterations: 66% Calculating genotypes for 2/8 loci. Completed iterations: 67% Calculating genotypes for 2/8 loci. Completed iterations: 68% Calculating genotypes for 2/8 loci. Completed iterations: 69% Calculating genotypes for 2/8 loci. Completed iterations: 70% Calculating genotypes for 2/8 loci. Completed iterations: 71% Calculating genotypes for 2/8 loci. Completed iterations: 72% Calculating genotypes for 2/8 loci. Completed iterations: 73% Calculating genotypes for 2/8 loci. Completed iterations: 74% Calculating genotypes for 2/8 loci. Completed iterations: 75% Calculating genotypes for 2/8 loci. Completed iterations: 76% Calculating genotypes for 2/8 loci. Completed iterations: 77% Calculating genotypes for 2/8 loci. Completed iterations: 78% Calculating genotypes for 2/8 loci. Completed iterations: 79% Calculating genotypes for 2/8 loci. Completed iterations: 80% Calculating genotypes for 2/8 loci. Completed iterations: 81% Calculating genotypes for 2/8 loci. Completed iterations: 82% Calculating genotypes for 2/8 loci. Completed iterations: 83% Calculating genotypes for 2/8 loci. Completed iterations: 84% Calculating genotypes for 2/8 loci. Completed iterations: 85% Calculating genotypes for 2/8 loci. Completed iterations: 86% Calculating genotypes for 2/8 loci. Completed iterations: 87% Calculating genotypes for 2/8 loci. Completed iterations: 88% Calculating genotypes for 2/8 loci. Completed iterations: 89% Calculating genotypes for 2/8 loci. Completed iterations: 90% Calculating genotypes for 2/8 loci. Completed iterations: 91% Calculating genotypes for 2/8 loci. Completed iterations: 92% Calculating genotypes for 2/8 loci. Completed iterations: 93% Calculating genotypes for 2/8 loci. Completed iterations: 94% Calculating genotypes for 2/8 loci. Completed iterations: 95% Calculating genotypes for 2/8 loci. Completed iterations: 96% Calculating genotypes for 2/8 loci. Completed iterations: 97% Calculating genotypes for 2/8 loci. Completed iterations: 98% Calculating genotypes for 2/8 loci. Completed iterations: 99% Calculating genotypes for 2/8 loci. Completed iterations: 100% Calculating genotypes for 3/8 loci. Completed iterations: 1% Calculating genotypes for 3/8 loci. Completed iterations: 2% Calculating genotypes for 3/8 loci. Completed iterations: 3% Calculating genotypes for 3/8 loci. Completed iterations: 4% Calculating genotypes for 3/8 loci. Completed iterations: 5% Calculating genotypes for 3/8 loci. Completed iterations: 6% Calculating genotypes for 3/8 loci. Completed iterations: 7% Calculating genotypes for 3/8 loci. Completed iterations: 8% Calculating genotypes for 3/8 loci. Completed iterations: 9% Calculating genotypes for 3/8 loci. Completed iterations: 10% Calculating genotypes for 3/8 loci. Completed iterations: 11% Calculating genotypes for 3/8 loci. Completed iterations: 12% Calculating genotypes for 3/8 loci. Completed iterations: 13% Calculating genotypes for 3/8 loci. Completed iterations: 14% Calculating genotypes for 3/8 loci. Completed iterations: 15% Calculating genotypes for 3/8 loci. Completed iterations: 16% Calculating genotypes for 3/8 loci. Completed iterations: 17% Calculating genotypes for 3/8 loci. Completed iterations: 18% Calculating genotypes for 3/8 loci. Completed iterations: 19% Calculating genotypes for 3/8 loci. Completed iterations: 20% Calculating genotypes for 3/8 loci. Completed iterations: 21% Calculating genotypes for 3/8 loci. Completed iterations: 22% Calculating genotypes for 3/8 loci. Completed iterations: 23% Calculating genotypes for 3/8 loci. Completed iterations: 24% Calculating genotypes for 3/8 loci. Completed iterations: 25% Calculating genotypes for 3/8 loci. Completed iterations: 26% Calculating genotypes for 3/8 loci. Completed iterations: 27% Calculating genotypes for 3/8 loci. Completed iterations: 28% Calculating genotypes for 3/8 loci. Completed iterations: 29% Calculating genotypes for 3/8 loci. Completed iterations: 30% Calculating genotypes for 3/8 loci. Completed iterations: 31% Calculating genotypes for 3/8 loci. Completed iterations: 32% Calculating genotypes for 3/8 loci. Completed iterations: 33% Calculating genotypes for 3/8 loci. Completed iterations: 34% Calculating genotypes for 3/8 loci. Completed iterations: 35% Calculating genotypes for 3/8 loci. Completed iterations: 36% Calculating genotypes for 3/8 loci. Completed iterations: 37% Calculating genotypes for 3/8 loci. Completed iterations: 38% Calculating genotypes for 3/8 loci. Completed iterations: 39% Calculating genotypes for 3/8 loci. Completed iterations: 40% Calculating genotypes for 3/8 loci. Completed iterations: 41% Calculating genotypes for 3/8 loci. Completed iterations: 42% Calculating genotypes for 3/8 loci. Completed iterations: 43% Calculating genotypes for 3/8 loci. Completed iterations: 44% Calculating genotypes for 3/8 loci. Completed iterations: 45% Calculating genotypes for 3/8 loci. Completed iterations: 46% Calculating genotypes for 3/8 loci. Completed iterations: 47% Calculating genotypes for 3/8 loci. Completed iterations: 48% Calculating genotypes for 3/8 loci. Completed iterations: 49% Calculating genotypes for 3/8 loci. Completed iterations: 50% Calculating genotypes for 3/8 loci. Completed iterations: 51% Calculating genotypes for 3/8 loci. Completed iterations: 52% Calculating genotypes for 3/8 loci. Completed iterations: 53% Calculating genotypes for 3/8 loci. Completed iterations: 54% Calculating genotypes for 3/8 loci. Completed iterations: 55% Calculating genotypes for 3/8 loci. Completed iterations: 56% Calculating genotypes for 3/8 loci. Completed iterations: 57% Calculating genotypes for 3/8 loci. Completed iterations: 58% Calculating genotypes for 3/8 loci. Completed iterations: 59% Calculating genotypes for 3/8 loci. Completed iterations: 60% Calculating genotypes for 3/8 loci. Completed iterations: 61% Calculating genotypes for 3/8 loci. Completed iterations: 62% Calculating genotypes for 3/8 loci. Completed iterations: 63% Calculating genotypes for 3/8 loci. Completed iterations: 64% Calculating genotypes for 3/8 loci. Completed iterations: 65% Calculating genotypes for 3/8 loci. Completed iterations: 66% Calculating genotypes for 3/8 loci. Completed iterations: 67% Calculating genotypes for 3/8 loci. Completed iterations: 68% Calculating genotypes for 3/8 loci. Completed iterations: 69% Calculating genotypes for 3/8 loci. Completed iterations: 70% Calculating genotypes for 3/8 loci. Completed iterations: 71% Calculating genotypes for 3/8 loci. Completed iterations: 72% Calculating genotypes for 3/8 loci. Completed iterations: 73% Calculating genotypes for 3/8 loci. Completed iterations: 74% Calculating genotypes for 3/8 loci. Completed iterations: 75% Calculating genotypes for 3/8 loci. Completed iterations: 76% Calculating genotypes for 3/8 loci. Completed iterations: 77% Calculating genotypes for 3/8 loci. Completed iterations: 78% Calculating genotypes for 3/8 loci. Completed iterations: 79% Calculating genotypes for 3/8 loci. Completed iterations: 80% Calculating genotypes for 3/8 loci. Completed iterations: 81% Calculating genotypes for 3/8 loci. Completed iterations: 82% Calculating genotypes for 3/8 loci. Completed iterations: 83% Calculating genotypes for 3/8 loci. Completed iterations: 84% Calculating genotypes for 3/8 loci. Completed iterations: 85% Calculating genotypes for 3/8 loci. Completed iterations: 86% Calculating genotypes for 3/8 loci. Completed iterations: 87% Calculating genotypes for 3/8 loci. Completed iterations: 88% Calculating genotypes for 3/8 loci. Completed iterations: 89% Calculating genotypes for 3/8 loci. Completed iterations: 90% Calculating genotypes for 3/8 loci. Completed iterations: 91% Calculating genotypes for 3/8 loci. Completed iterations: 92% Calculating genotypes for 3/8 loci. Completed iterations: 93% Calculating genotypes for 3/8 loci. Completed iterations: 94% Calculating genotypes for 3/8 loci. Completed iterations: 95% Calculating genotypes for 3/8 loci. Completed iterations: 96% Calculating genotypes for 3/8 loci. Completed iterations: 97% Calculating genotypes for 3/8 loci. Completed iterations: 98% Calculating genotypes for 3/8 loci. Completed iterations: 99% Calculating genotypes for 3/8 loci. Completed iterations: 100% Calculating genotypes for 4/8 loci. Completed iterations: 1% Calculating genotypes for 4/8 loci. Completed iterations: 2% Calculating genotypes for 4/8 loci. Completed iterations: 3% Calculating genotypes for 4/8 loci. Completed iterations: 4% Calculating genotypes for 4/8 loci. Completed iterations: 5% Calculating genotypes for 4/8 loci. Completed iterations: 6% Calculating genotypes for 4/8 loci. Completed iterations: 7% Calculating genotypes for 4/8 loci. Completed iterations: 8% Calculating genotypes for 4/8 loci. Completed iterations: 9% Calculating genotypes for 4/8 loci. Completed iterations: 10% Calculating genotypes for 4/8 loci. Completed iterations: 11% Calculating genotypes for 4/8 loci. Completed iterations: 12% Calculating genotypes for 4/8 loci. Completed iterations: 13% Calculating genotypes for 4/8 loci. Completed iterations: 14% Calculating genotypes for 4/8 loci. Completed iterations: 15% Calculating genotypes for 4/8 loci. Completed iterations: 16% Calculating genotypes for 4/8 loci. Completed iterations: 17% Calculating genotypes for 4/8 loci. Completed iterations: 18% Calculating genotypes for 4/8 loci. Completed iterations: 19% Calculating genotypes for 4/8 loci. Completed iterations: 20% Calculating genotypes for 4/8 loci. Completed iterations: 21% Calculating genotypes for 4/8 loci. Completed iterations: 22% Calculating genotypes for 4/8 loci. Completed iterations: 23% Calculating genotypes for 4/8 loci. Completed iterations: 24% Calculating genotypes for 4/8 loci. Completed iterations: 25% Calculating genotypes for 4/8 loci. Completed iterations: 26% Calculating genotypes for 4/8 loci. Completed iterations: 27% Calculating genotypes for 4/8 loci. Completed iterations: 28% Calculating genotypes for 4/8 loci. Completed iterations: 29% Calculating genotypes for 4/8 loci. Completed iterations: 30% Calculating genotypes for 4/8 loci. Completed iterations: 31% Calculating genotypes for 4/8 loci. Completed iterations: 32% Calculating genotypes for 4/8 loci. Completed iterations: 33% Calculating genotypes for 4/8 loci. Completed iterations: 34% Calculating genotypes for 4/8 loci. Completed iterations: 35% Calculating genotypes for 4/8 loci. Completed iterations: 36% Calculating genotypes for 4/8 loci. Completed iterations: 37% Calculating genotypes for 4/8 loci. Completed iterations: 38% Calculating genotypes for 4/8 loci. Completed iterations: 39% Calculating genotypes for 4/8 loci. Completed iterations: 40% Calculating genotypes for 4/8 loci. Completed iterations: 41% Calculating genotypes for 4/8 loci. Completed iterations: 42% Calculating genotypes for 4/8 loci. Completed iterations: 43% Calculating genotypes for 4/8 loci. Completed iterations: 44% Calculating genotypes for 4/8 loci. Completed iterations: 45% Calculating genotypes for 4/8 loci. Completed iterations: 46% Calculating genotypes for 4/8 loci. Completed iterations: 47% Calculating genotypes for 4/8 loci. Completed iterations: 48% Calculating genotypes for 4/8 loci. Completed iterations: 49% Calculating genotypes for 4/8 loci. Completed iterations: 50% Calculating genotypes for 4/8 loci. Completed iterations: 51% Calculating genotypes for 4/8 loci. Completed iterations: 52% Calculating genotypes for 4/8 loci. Completed iterations: 53% Calculating genotypes for 4/8 loci. Completed iterations: 54% Calculating genotypes for 4/8 loci. Completed iterations: 55% Calculating genotypes for 4/8 loci. Completed iterations: 56% Calculating genotypes for 4/8 loci. Completed iterations: 57% Calculating genotypes for 4/8 loci. Completed iterations: 58% Calculating genotypes for 4/8 loci. Completed iterations: 59% Calculating genotypes for 4/8 loci. Completed iterations: 60% Calculating genotypes for 4/8 loci. Completed iterations: 61% Calculating genotypes for 4/8 loci. Completed iterations: 62% Calculating genotypes for 4/8 loci. Completed iterations: 63% Calculating genotypes for 4/8 loci. Completed iterations: 64% Calculating genotypes for 4/8 loci. Completed iterations: 65% Calculating genotypes for 4/8 loci. Completed iterations: 66% Calculating genotypes for 4/8 loci. Completed iterations: 67% Calculating genotypes for 4/8 loci. Completed iterations: 68% Calculating genotypes for 4/8 loci. Completed iterations: 69% Calculating genotypes for 4/8 loci. Completed iterations: 70% Calculating genotypes for 4/8 loci. Completed iterations: 71% Calculating genotypes for 4/8 loci. Completed iterations: 72% Calculating genotypes for 4/8 loci. Completed iterations: 73% Calculating genotypes for 4/8 loci. Completed iterations: 74% Calculating genotypes for 4/8 loci. Completed iterations: 75% Calculating genotypes for 4/8 loci. Completed iterations: 76% Calculating genotypes for 4/8 loci. Completed iterations: 77% Calculating genotypes for 4/8 loci. Completed iterations: 78% Calculating genotypes for 4/8 loci. Completed iterations: 79% Calculating genotypes for 4/8 loci. Completed iterations: 80% Calculating genotypes for 4/8 loci. Completed iterations: 81% Calculating genotypes for 4/8 loci. Completed iterations: 82% Calculating genotypes for 4/8 loci. Completed iterations: 83% Calculating genotypes for 4/8 loci. Completed iterations: 84% Calculating genotypes for 4/8 loci. Completed iterations: 85% Calculating genotypes for 4/8 loci. Completed iterations: 86% Calculating genotypes for 4/8 loci. Completed iterations: 87% Calculating genotypes for 4/8 loci. Completed iterations: 88% Calculating genotypes for 4/8 loci. Completed iterations: 89% Calculating genotypes for 4/8 loci. Completed iterations: 90% Calculating genotypes for 4/8 loci. Completed iterations: 91% Calculating genotypes for 4/8 loci. Completed iterations: 92% Calculating genotypes for 4/8 loci. Completed iterations: 93% Calculating genotypes for 4/8 loci. Completed iterations: 94% Calculating genotypes for 4/8 loci. Completed iterations: 95% Calculating genotypes for 4/8 loci. Completed iterations: 96% Calculating genotypes for 4/8 loci. Completed iterations: 97% Calculating genotypes for 4/8 loci. Completed iterations: 98% Calculating genotypes for 4/8 loci. Completed iterations: 99% Calculating genotypes for 4/8 loci. Completed iterations: 100% Calculating genotypes for 5/8 loci. Completed iterations: 1% Calculating genotypes for 5/8 loci. Completed iterations: 2% Calculating genotypes for 5/8 loci. Completed iterations: 3% Calculating genotypes for 5/8 loci. Completed iterations: 4% Calculating genotypes for 5/8 loci. Completed iterations: 5% Calculating genotypes for 5/8 loci. Completed iterations: 6% Calculating genotypes for 5/8 loci. Completed iterations: 7% Calculating genotypes for 5/8 loci. Completed iterations: 8% Calculating genotypes for 5/8 loci. Completed iterations: 9% Calculating genotypes for 5/8 loci. Completed iterations: 10% Calculating genotypes for 5/8 loci. Completed iterations: 11% Calculating genotypes for 5/8 loci. Completed iterations: 12% Calculating genotypes for 5/8 loci. Completed iterations: 13% Calculating genotypes for 5/8 loci. Completed iterations: 14% Calculating genotypes for 5/8 loci. Completed iterations: 15% Calculating genotypes for 5/8 loci. Completed iterations: 16% Calculating genotypes for 5/8 loci. Completed iterations: 17% Calculating genotypes for 5/8 loci. Completed iterations: 18% Calculating genotypes for 5/8 loci. Completed iterations: 19% Calculating genotypes for 5/8 loci. Completed iterations: 20% Calculating genotypes for 5/8 loci. Completed iterations: 21% Calculating genotypes for 5/8 loci. Completed iterations: 22% Calculating genotypes for 5/8 loci. Completed iterations: 23% Calculating genotypes for 5/8 loci. Completed iterations: 24% Calculating genotypes for 5/8 loci. Completed iterations: 25% Calculating genotypes for 5/8 loci. Completed iterations: 26% Calculating genotypes for 5/8 loci. Completed iterations: 27% Calculating genotypes for 5/8 loci. Completed iterations: 28% Calculating genotypes for 5/8 loci. Completed iterations: 29% Calculating genotypes for 5/8 loci. Completed iterations: 30% Calculating genotypes for 5/8 loci. Completed iterations: 31% Calculating genotypes for 5/8 loci. Completed iterations: 32% Calculating genotypes for 5/8 loci. Completed iterations: 33% Calculating genotypes for 5/8 loci. Completed iterations: 34% Calculating genotypes for 5/8 loci. Completed iterations: 35% Calculating genotypes for 5/8 loci. Completed iterations: 36% Calculating genotypes for 5/8 loci. Completed iterations: 37% Calculating genotypes for 5/8 loci. Completed iterations: 38% Calculating genotypes for 5/8 loci. Completed iterations: 39% Calculating genotypes for 5/8 loci. Completed iterations: 40% Calculating genotypes for 5/8 loci. Completed iterations: 41% Calculating genotypes for 5/8 loci. Completed iterations: 42% Calculating genotypes for 5/8 loci. Completed iterations: 43% Calculating genotypes for 5/8 loci. Completed iterations: 44% Calculating genotypes for 5/8 loci. Completed iterations: 45% Calculating genotypes for 5/8 loci. Completed iterations: 46% Calculating genotypes for 5/8 loci. Completed iterations: 47% Calculating genotypes for 5/8 loci. Completed iterations: 48% Calculating genotypes for 5/8 loci. Completed iterations: 49% Calculating genotypes for 5/8 loci. Completed iterations: 50% Calculating genotypes for 5/8 loci. Completed iterations: 51% Calculating genotypes for 5/8 loci. Completed iterations: 52% Calculating genotypes for 5/8 loci. Completed iterations: 53% Calculating genotypes for 5/8 loci. Completed iterations: 54% Calculating genotypes for 5/8 loci. Completed iterations: 55% Calculating genotypes for 5/8 loci. Completed iterations: 56% Calculating genotypes for 5/8 loci. Completed iterations: 57% Calculating genotypes for 5/8 loci. Completed iterations: 58% Calculating genotypes for 5/8 loci. Completed iterations: 59% Calculating genotypes for 5/8 loci. Completed iterations: 60% Calculating genotypes for 5/8 loci. Completed iterations: 61% Calculating genotypes for 5/8 loci. Completed iterations: 62% Calculating genotypes for 5/8 loci. Completed iterations: 63% Calculating genotypes for 5/8 loci. Completed iterations: 64% Calculating genotypes for 5/8 loci. Completed iterations: 65% Calculating genotypes for 5/8 loci. Completed iterations: 66% Calculating genotypes for 5/8 loci. Completed iterations: 67% Calculating genotypes for 5/8 loci. Completed iterations: 68% Calculating genotypes for 5/8 loci. Completed iterations: 69% Calculating genotypes for 5/8 loci. Completed iterations: 70% Calculating genotypes for 5/8 loci. Completed iterations: 71% Calculating genotypes for 5/8 loci. Completed iterations: 72% Calculating genotypes for 5/8 loci. Completed iterations: 73% Calculating genotypes for 5/8 loci. Completed iterations: 74% Calculating genotypes for 5/8 loci. Completed iterations: 75% Calculating genotypes for 5/8 loci. Completed iterations: 76% Calculating genotypes for 5/8 loci. Completed iterations: 77% Calculating genotypes for 5/8 loci. Completed iterations: 78% Calculating genotypes for 5/8 loci. Completed iterations: 79% Calculating genotypes for 5/8 loci. Completed iterations: 80% Calculating genotypes for 5/8 loci. Completed iterations: 81% Calculating genotypes for 5/8 loci. Completed iterations: 82% Calculating genotypes for 5/8 loci. Completed iterations: 83% Calculating genotypes for 5/8 loci. Completed iterations: 84% Calculating genotypes for 5/8 loci. Completed iterations: 85% Calculating genotypes for 5/8 loci. Completed iterations: 86% Calculating genotypes for 5/8 loci. Completed iterations: 87% Calculating genotypes for 5/8 loci. Completed iterations: 88% Calculating genotypes for 5/8 loci. Completed iterations: 89% Calculating genotypes for 5/8 loci. Completed iterations: 90% Calculating genotypes for 5/8 loci. Completed iterations: 91% Calculating genotypes for 5/8 loci. Completed iterations: 92% Calculating genotypes for 5/8 loci. Completed iterations: 93% Calculating genotypes for 5/8 loci. Completed iterations: 94% Calculating genotypes for 5/8 loci. Completed iterations: 95% Calculating genotypes for 5/8 loci. Completed iterations: 96% Calculating genotypes for 5/8 loci. Completed iterations: 97% Calculating genotypes for 5/8 loci. Completed iterations: 98% Calculating genotypes for 5/8 loci. Completed iterations: 99% Calculating genotypes for 5/8 loci. Completed iterations: 100% Calculating genotypes for 6/8 loci. Completed iterations: 1% Calculating genotypes for 6/8 loci. Completed iterations: 2% Calculating genotypes for 6/8 loci. Completed iterations: 3% Calculating genotypes for 6/8 loci. Completed iterations: 4% Calculating genotypes for 6/8 loci. Completed iterations: 5% Calculating genotypes for 6/8 loci. Completed iterations: 6% Calculating genotypes for 6/8 loci. Completed iterations: 7% Calculating genotypes for 6/8 loci. Completed iterations: 8% Calculating genotypes for 6/8 loci. Completed iterations: 9% Calculating genotypes for 6/8 loci. Completed iterations: 10% Calculating genotypes for 6/8 loci. Completed iterations: 11% Calculating genotypes for 6/8 loci. Completed iterations: 12% Calculating genotypes for 6/8 loci. Completed iterations: 13% Calculating genotypes for 6/8 loci. Completed iterations: 14% Calculating genotypes for 6/8 loci. Completed iterations: 15% Calculating genotypes for 6/8 loci. Completed iterations: 16% Calculating genotypes for 6/8 loci. Completed iterations: 17% Calculating genotypes for 6/8 loci. Completed iterations: 18% Calculating genotypes for 6/8 loci. Completed iterations: 19% Calculating genotypes for 6/8 loci. Completed iterations: 20% Calculating genotypes for 6/8 loci. Completed iterations: 21% Calculating genotypes for 6/8 loci. Completed iterations: 22% Calculating genotypes for 6/8 loci. Completed iterations: 23% Calculating genotypes for 6/8 loci. Completed iterations: 24% Calculating genotypes for 6/8 loci. Completed iterations: 25% Calculating genotypes for 6/8 loci. Completed iterations: 26% Calculating genotypes for 6/8 loci. Completed iterations: 27% Calculating genotypes for 6/8 loci. Completed iterations: 28% Calculating genotypes for 6/8 loci. Completed iterations: 29% Calculating genotypes for 6/8 loci. Completed iterations: 30% Calculating genotypes for 6/8 loci. Completed iterations: 31% Calculating genotypes for 6/8 loci. Completed iterations: 32% Calculating genotypes for 6/8 loci. Completed iterations: 33% Calculating genotypes for 6/8 loci. Completed iterations: 34% Calculating genotypes for 6/8 loci. Completed iterations: 35% Calculating genotypes for 6/8 loci. Completed iterations: 36% Calculating genotypes for 6/8 loci. Completed iterations: 37% Calculating genotypes for 6/8 loci. Completed iterations: 38% Calculating genotypes for 6/8 loci. Completed iterations: 39% Calculating genotypes for 6/8 loci. Completed iterations: 40% Calculating genotypes for 6/8 loci. Completed iterations: 41% Calculating genotypes for 6/8 loci. Completed iterations: 42% Calculating genotypes for 6/8 loci. Completed iterations: 43% Calculating genotypes for 6/8 loci. Completed iterations: 44% Calculating genotypes for 6/8 loci. Completed iterations: 45% Calculating genotypes for 6/8 loci. Completed iterations: 46% Calculating genotypes for 6/8 loci. Completed iterations: 47% Calculating genotypes for 6/8 loci. Completed iterations: 48% Calculating genotypes for 6/8 loci. Completed iterations: 49% Calculating genotypes for 6/8 loci. Completed iterations: 50% Calculating genotypes for 6/8 loci. Completed iterations: 51% Calculating genotypes for 6/8 loci. Completed iterations: 52% Calculating genotypes for 6/8 loci. Completed iterations: 53% Calculating genotypes for 6/8 loci. Completed iterations: 54% Calculating genotypes for 6/8 loci. Completed iterations: 55% Calculating genotypes for 6/8 loci. Completed iterations: 56% Calculating genotypes for 6/8 loci. Completed iterations: 57% Calculating genotypes for 6/8 loci. Completed iterations: 58% Calculating genotypes for 6/8 loci. Completed iterations: 59% Calculating genotypes for 6/8 loci. Completed iterations: 60% Calculating genotypes for 6/8 loci. Completed iterations: 61% Calculating genotypes for 6/8 loci. Completed iterations: 62% Calculating genotypes for 6/8 loci. Completed iterations: 63% Calculating genotypes for 6/8 loci. Completed iterations: 64% Calculating genotypes for 6/8 loci. Completed iterations: 65% Calculating genotypes for 6/8 loci. Completed iterations: 66% Calculating genotypes for 6/8 loci. Completed iterations: 67% Calculating genotypes for 6/8 loci. Completed iterations: 68% Calculating genotypes for 6/8 loci. Completed iterations: 69% Calculating genotypes for 6/8 loci. Completed iterations: 70% Calculating genotypes for 6/8 loci. Completed iterations: 71% Calculating genotypes for 6/8 loci. Completed iterations: 72% Calculating genotypes for 6/8 loci. Completed iterations: 73% Calculating genotypes for 6/8 loci. Completed iterations: 74% Calculating genotypes for 6/8 loci. Completed iterations: 75% Calculating genotypes for 6/8 loci. Completed iterations: 76% Calculating genotypes for 6/8 loci. Completed iterations: 77% Calculating genotypes for 6/8 loci. Completed iterations: 78% Calculating genotypes for 6/8 loci. Completed iterations: 79% Calculating genotypes for 6/8 loci. Completed iterations: 80% Calculating genotypes for 6/8 loci. Completed iterations: 81% Calculating genotypes for 6/8 loci. Completed iterations: 82% Calculating genotypes for 6/8 loci. Completed iterations: 83% Calculating genotypes for 6/8 loci. Completed iterations: 84% Calculating genotypes for 6/8 loci. Completed iterations: 85% Calculating genotypes for 6/8 loci. Completed iterations: 86% Calculating genotypes for 6/8 loci. Completed iterations: 87% Calculating genotypes for 6/8 loci. Completed iterations: 88% Calculating genotypes for 6/8 loci. Completed iterations: 89% Calculating genotypes for 6/8 loci. Completed iterations: 90% Calculating genotypes for 6/8 loci. Completed iterations: 91% Calculating genotypes for 6/8 loci. Completed iterations: 92% Calculating genotypes for 6/8 loci. Completed iterations: 93% Calculating genotypes for 6/8 loci. Completed iterations: 94% Calculating genotypes for 6/8 loci. Completed iterations: 95% Calculating genotypes for 6/8 loci. Completed iterations: 96% Calculating genotypes for 6/8 loci. Completed iterations: 97% Calculating genotypes for 6/8 loci. Completed iterations: 98% Calculating genotypes for 6/8 loci. Completed iterations: 99% Calculating genotypes for 6/8 loci. Completed iterations: 100% Calculating genotypes for 7/8 loci. Completed iterations: 1% Calculating genotypes for 7/8 loci. Completed iterations: 2% Calculating genotypes for 7/8 loci. Completed iterations: 3% Calculating genotypes for 7/8 loci. Completed iterations: 4% Calculating genotypes for 7/8 loci. Completed iterations: 5% Calculating genotypes for 7/8 loci. Completed iterations: 6% Calculating genotypes for 7/8 loci. Completed iterations: 7% Calculating genotypes for 7/8 loci. Completed iterations: 8% Calculating genotypes for 7/8 loci. Completed iterations: 9% Calculating genotypes for 7/8 loci. Completed iterations: 10% Calculating genotypes for 7/8 loci. Completed iterations: 11% Calculating genotypes for 7/8 loci. Completed iterations: 12% Calculating genotypes for 7/8 loci. Completed iterations: 13% Calculating genotypes for 7/8 loci. Completed iterations: 14% Calculating genotypes for 7/8 loci. Completed iterations: 15% Calculating genotypes for 7/8 loci. Completed iterations: 16% Calculating genotypes for 7/8 loci. Completed iterations: 17% Calculating genotypes for 7/8 loci. Completed iterations: 18% Calculating genotypes for 7/8 loci. Completed iterations: 19% Calculating genotypes for 7/8 loci. Completed iterations: 20% Calculating genotypes for 7/8 loci. Completed iterations: 21% Calculating genotypes for 7/8 loci. Completed iterations: 22% Calculating genotypes for 7/8 loci. Completed iterations: 23% Calculating genotypes for 7/8 loci. Completed iterations: 24% Calculating genotypes for 7/8 loci. Completed iterations: 25% Calculating genotypes for 7/8 loci. Completed iterations: 26% Calculating genotypes for 7/8 loci. Completed iterations: 27% Calculating genotypes for 7/8 loci. Completed iterations: 28% Calculating genotypes for 7/8 loci. Completed iterations: 29% Calculating genotypes for 7/8 loci. Completed iterations: 30% Calculating genotypes for 7/8 loci. Completed iterations: 31% Calculating genotypes for 7/8 loci. Completed iterations: 32% Calculating genotypes for 7/8 loci. Completed iterations: 33% Calculating genotypes for 7/8 loci. Completed iterations: 34% Calculating genotypes for 7/8 loci. Completed iterations: 35% Calculating genotypes for 7/8 loci. Completed iterations: 36% Calculating genotypes for 7/8 loci. Completed iterations: 37% Calculating genotypes for 7/8 loci. Completed iterations: 38% Calculating genotypes for 7/8 loci. Completed iterations: 39% Calculating genotypes for 7/8 loci. Completed iterations: 40% Calculating genotypes for 7/8 loci. Completed iterations: 41% Calculating genotypes for 7/8 loci. Completed iterations: 42% Calculating genotypes for 7/8 loci. Completed iterations: 43% Calculating genotypes for 7/8 loci. Completed iterations: 44% Calculating genotypes for 7/8 loci. Completed iterations: 45% Calculating genotypes for 7/8 loci. Completed iterations: 46% Calculating genotypes for 7/8 loci. Completed iterations: 47% Calculating genotypes for 7/8 loci. Completed iterations: 48% Calculating genotypes for 7/8 loci. Completed iterations: 49% Calculating genotypes for 7/8 loci. Completed iterations: 50% Calculating genotypes for 7/8 loci. Completed iterations: 51% Calculating genotypes for 7/8 loci. Completed iterations: 52% Calculating genotypes for 7/8 loci. Completed iterations: 53% Calculating genotypes for 7/8 loci. Completed iterations: 54% Calculating genotypes for 7/8 loci. Completed iterations: 55% Calculating genotypes for 7/8 loci. Completed iterations: 56% Calculating genotypes for 7/8 loci. Completed iterations: 57% Calculating genotypes for 7/8 loci. Completed iterations: 58% Calculating genotypes for 7/8 loci. Completed iterations: 59% Calculating genotypes for 7/8 loci. Completed iterations: 60% Calculating genotypes for 7/8 loci. Completed iterations: 61% Calculating genotypes for 7/8 loci. Completed iterations: 62% Calculating genotypes for 7/8 loci. Completed iterations: 63% Calculating genotypes for 7/8 loci. Completed iterations: 64% Calculating genotypes for 7/8 loci. Completed iterations: 65% Calculating genotypes for 7/8 loci. Completed iterations: 66% Calculating genotypes for 7/8 loci. Completed iterations: 67% Calculating genotypes for 7/8 loci. Completed iterations: 68% Calculating genotypes for 7/8 loci. Completed iterations: 69% Calculating genotypes for 7/8 loci. Completed iterations: 70% Calculating genotypes for 7/8 loci. Completed iterations: 71% Calculating genotypes for 7/8 loci. Completed iterations: 72% Calculating genotypes for 7/8 loci. Completed iterations: 73% Calculating genotypes for 7/8 loci. Completed iterations: 74% Calculating genotypes for 7/8 loci. Completed iterations: 75% Calculating genotypes for 7/8 loci. Completed iterations: 76% Calculating genotypes for 7/8 loci. Completed iterations: 77% Calculating genotypes for 7/8 loci. Completed iterations: 78% Calculating genotypes for 7/8 loci. Completed iterations: 79% Calculating genotypes for 7/8 loci. Completed iterations: 80% Calculating genotypes for 7/8 loci. Completed iterations: 81% Calculating genotypes for 7/8 loci. Completed iterations: 82% Calculating genotypes for 7/8 loci. Completed iterations: 83% Calculating genotypes for 7/8 loci. Completed iterations: 84% Calculating genotypes for 7/8 loci. Completed iterations: 85% Calculating genotypes for 7/8 loci. Completed iterations: 86% Calculating genotypes for 7/8 loci. Completed iterations: 87% Calculating genotypes for 7/8 loci. Completed iterations: 88% Calculating genotypes for 7/8 loci. Completed iterations: 89% Calculating genotypes for 7/8 loci. Completed iterations: 90% Calculating genotypes for 7/8 loci. Completed iterations: 91% Calculating genotypes for 7/8 loci. Completed iterations: 92% Calculating genotypes for 7/8 loci. Completed iterations: 93% Calculating genotypes for 7/8 loci. Completed iterations: 94% Calculating genotypes for 7/8 loci. Completed iterations: 95% Calculating genotypes for 7/8 loci. Completed iterations: 96% Calculating genotypes for 7/8 loci. Completed iterations: 97% Calculating genotypes for 7/8 loci. Completed iterations: 98% Calculating genotypes for 7/8 loci. Completed iterations: 99% Calculating genotypes for 7/8 loci. Completed iterations: 100% Calculating genotypes for 8/8 loci. Completed iterations: 1% Calculating genotypes for 8/8 loci. Completed iterations: 2% Calculating genotypes for 8/8 loci. Completed iterations: 3% Calculating genotypes for 8/8 loci. Completed iterations: 4% Calculating genotypes for 8/8 loci. Completed iterations: 5% Calculating genotypes for 8/8 loci. Completed iterations: 6% Calculating genotypes for 8/8 loci. Completed iterations: 7% Calculating genotypes for 8/8 loci. Completed iterations: 8% Calculating genotypes for 8/8 loci. Completed iterations: 9% Calculating genotypes for 8/8 loci. Completed iterations: 10% Calculating genotypes for 8/8 loci. Completed iterations: 11% Calculating genotypes for 8/8 loci. Completed iterations: 12% Calculating genotypes for 8/8 loci. Completed iterations: 13% Calculating genotypes for 8/8 loci. Completed iterations: 14% Calculating genotypes for 8/8 loci. Completed iterations: 15% Calculating genotypes for 8/8 loci. Completed iterations: 16% Calculating genotypes for 8/8 loci. Completed iterations: 17% Calculating genotypes for 8/8 loci. Completed iterations: 18% Calculating genotypes for 8/8 loci. Completed iterations: 19% Calculating genotypes for 8/8 loci. Completed iterations: 20% Calculating genotypes for 8/8 loci. Completed iterations: 21% Calculating genotypes for 8/8 loci. Completed iterations: 22% Calculating genotypes for 8/8 loci. Completed iterations: 23% Calculating genotypes for 8/8 loci. Completed iterations: 24% Calculating genotypes for 8/8 loci. Completed iterations: 25% Calculating genotypes for 8/8 loci. Completed iterations: 26% Calculating genotypes for 8/8 loci. Completed iterations: 27% Calculating genotypes for 8/8 loci. Completed iterations: 28% Calculating genotypes for 8/8 loci. Completed iterations: 29% Calculating genotypes for 8/8 loci. Completed iterations: 30% Calculating genotypes for 8/8 loci. Completed iterations: 31% Calculating genotypes for 8/8 loci. Completed iterations: 32% Calculating genotypes for 8/8 loci. Completed iterations: 33% Calculating genotypes for 8/8 loci. Completed iterations: 34% Calculating genotypes for 8/8 loci. Completed iterations: 35% Calculating genotypes for 8/8 loci. Completed iterations: 36% Calculating genotypes for 8/8 loci. Completed iterations: 37% Calculating genotypes for 8/8 loci. Completed iterations: 38% Calculating genotypes for 8/8 loci. Completed iterations: 39% Calculating genotypes for 8/8 loci. Completed iterations: 40% Calculating genotypes for 8/8 loci. Completed iterations: 41% Calculating genotypes for 8/8 loci. Completed iterations: 42% Calculating genotypes for 8/8 loci. Completed iterations: 43% Calculating genotypes for 8/8 loci. Completed iterations: 44% Calculating genotypes for 8/8 loci. Completed iterations: 45% Calculating genotypes for 8/8 loci. Completed iterations: 46% Calculating genotypes for 8/8 loci. Completed iterations: 47% Calculating genotypes for 8/8 loci. Completed iterations: 48% Calculating genotypes for 8/8 loci. Completed iterations: 49% Calculating genotypes for 8/8 loci. Completed iterations: 50% Calculating genotypes for 8/8 loci. Completed iterations: 51% Calculating genotypes for 8/8 loci. Completed iterations: 52% Calculating genotypes for 8/8 loci. Completed iterations: 53% Calculating genotypes for 8/8 loci. Completed iterations: 54% Calculating genotypes for 8/8 loci. Completed iterations: 55% Calculating genotypes for 8/8 loci. Completed iterations: 56% Calculating genotypes for 8/8 loci. Completed iterations: 57% Calculating genotypes for 8/8 loci. Completed iterations: 58% Calculating genotypes for 8/8 loci. Completed iterations: 59% Calculating genotypes for 8/8 loci. Completed iterations: 60% Calculating genotypes for 8/8 loci. Completed iterations: 61% Calculating genotypes for 8/8 loci. Completed iterations: 62% Calculating genotypes for 8/8 loci. Completed iterations: 63% Calculating genotypes for 8/8 loci. Completed iterations: 64% Calculating genotypes for 8/8 loci. Completed iterations: 65% Calculating genotypes for 8/8 loci. Completed iterations: 66% Calculating genotypes for 8/8 loci. Completed iterations: 67% Calculating genotypes for 8/8 loci. Completed iterations: 68% Calculating genotypes for 8/8 loci. Completed iterations: 69% Calculating genotypes for 8/8 loci. Completed iterations: 70% Calculating genotypes for 8/8 loci. Completed iterations: 71% Calculating genotypes for 8/8 loci. Completed iterations: 72% Calculating genotypes for 8/8 loci. Completed iterations: 73% Calculating genotypes for 8/8 loci. Completed iterations: 74% Calculating genotypes for 8/8 loci. Completed iterations: 75% Calculating genotypes for 8/8 loci. Completed iterations: 76% Calculating genotypes for 8/8 loci. Completed iterations: 77% Calculating genotypes for 8/8 loci. Completed iterations: 78% Calculating genotypes for 8/8 loci. Completed iterations: 79% Calculating genotypes for 8/8 loci. Completed iterations: 80% Calculating genotypes for 8/8 loci. Completed iterations: 81% Calculating genotypes for 8/8 loci. Completed iterations: 82% Calculating genotypes for 8/8 loci. Completed iterations: 83% Calculating genotypes for 8/8 loci. Completed iterations: 84% Calculating genotypes for 8/8 loci. Completed iterations: 85% Calculating genotypes for 8/8 loci. Completed iterations: 86% Calculating genotypes for 8/8 loci. Completed iterations: 87% Calculating genotypes for 8/8 loci. Completed iterations: 88% Calculating genotypes for 8/8 loci. Completed iterations: 89% Calculating genotypes for 8/8 loci. Completed iterations: 90% Calculating genotypes for 8/8 loci. Completed iterations: 91% Calculating genotypes for 8/8 loci. Completed iterations: 92% Calculating genotypes for 8/8 loci. Completed iterations: 93% Calculating genotypes for 8/8 loci. Completed iterations: 94% Calculating genotypes for 8/8 loci. Completed iterations: 95% Calculating genotypes for 8/8 loci. Completed iterations: 96% Calculating genotypes for 8/8 loci. Completed iterations: 97% Calculating genotypes for 8/8 loci. Completed iterations: 98% Calculating genotypes for 8/8 loci. Completed iterations: 99% Calculating genotypes for 8/8 loci. Completed iterations: 100%
    # NOT RUN { # Marker Type Comparison -------------------------------------------------- # With AFLP data, it is often necessary to include more markers for resolution data(Aeut) diff --git a/docs/reference/greycurve.html b/docs/reference/greycurve.html index 93238278..bae38905 100644 --- a/docs/reference/greycurve.html +++ b/docs/reference/greycurve.html @@ -153,7 +153,7 @@

    Value

    Examples

    # Normal grey curve with an adjustment of 3, an upper limit of 0.8, and # weighted towards smaller values. -greycurve()
    # NOT RUN { +greycurve()
    # NOT RUN { # 1:1 relationship grey curve. greycurve(gadj=1, glim=1:0) diff --git a/docs/reference/ia.html b/docs/reference/ia.html index 70c8d9c3..430dd537 100644 --- a/docs/reference/ia.html +++ b/docs/reference/ia.html @@ -343,10 +343,10 @@

    Examp #> 0.17207262 0.02178965
    # Pairwise over all loci: data(partial_clone) -res <- pair.ia(partial_clone)
    #> | | | 0% | |== | 2% | |=== | 4% | |===== | 7% | |====== | 9% | |======== | 11% | |========= | 13% | |=========== | 16% | |============ | 18% | |============== | 20% | |================ | 22% | |================= | 24% | |=================== | 27% | |==================== | 29% | |====================== | 31% | |======================= | 33% | |========================= | 36% | |========================== | 38% | |============================ | 40% | |============================== | 42% | |=============================== | 44% | |================================= | 47% | |================================== | 49% | |==================================== | 51% | |===================================== | 53% | |======================================= | 56% | |======================================== | 58% | |========================================== | 60% | |============================================ | 62% | |============================================= | 64% | |=============================================== | 67% | |================================================ | 69% | |================================================== | 71% | |=================================================== | 73% | |===================================================== | 76% | |====================================================== | 78% | |======================================================== | 80% | |========================================================== | 82% | |=========================================================== | 84% | |============================================================= | 87% | |============================================================== | 89% | |================================================================ | 91% | |================================================================= | 93% | |=================================================================== | 96% | |==================================================================== | 98% | |======================================================================| 100%
    plot(res, low = "black", high = "green", index = "Ia")
    +res <- pair.ia(partial_clone)
    #> | | | 0% | |== | 2% | |=== | 4% | |===== | 7% | |====== | 9% | |======== | 11% | |========= | 13% | |=========== | 16% | |============ | 18% | |============== | 20% | |================ | 22% | |================= | 24% | |=================== | 27% | |==================== | 29% | |====================== | 31% | |======================= | 33% | |========================= | 36% | |========================== | 38% | |============================ | 40% | |============================== | 42% | |=============================== | 44% | |================================= | 47% | |================================== | 49% | |==================================== | 51% | |===================================== | 53% | |======================================= | 56% | |======================================== | 58% | |========================================== | 60% | |============================================ | 62% | |============================================= | 64% | |=============================================== | 67% | |================================================ | 69% | |================================================== | 71% | |=================================================== | 73% | |===================================================== | 76% | |====================================================== | 78% | |======================================================== | 80% | |========================================================== | 82% | |=========================================================== | 84% | |============================================================= | 87% | |============================================================== | 89% | |================================================================ | 91% | |================================================================= | 93% | |=================================================================== | 96% | |==================================================================== | 98% | |======================================================================| 100%
    plot(res, low = "black", high = "green", index = "Ia")
    # Resampling data(Pinf) -resample.ia(Pinf, reps = 99)
    #> |=============== | 29% ~0 s remaining |=============================== | 59% ~0 s remaining |============================================== | 86% ~0 s remaining
    #> Ia rbarD +resample.ia(Pinf, reps = 99)
    #> |================ | 31% ~0 s remaining |=============================== | 59% ~0 s remaining |=================================================== | 95% ~0 s remaining
    #> Ia rbarD #> 1 0.6191763 0.06783205 #> 2 0.6798527 0.07473909 #> 3 0.5200972 0.05695292 diff --git a/docs/reference/index.html b/docs/reference/index.html index 0c2b5cb8..1018756e 100644 --- a/docs/reference/index.html +++ b/docs/reference/index.html @@ -100,7 +100,7 @@ @@ -134,7 +134,7 @@

    bootgen2genind as.genclone genclone2genind

    +

    bootgen2genind as.genclone genclone2genind as.genambig

    @@ -168,7 +168,7 @@

    bootgen2genind as.genclone genclone2genind

    +

    bootgen2genind as.genclone genclone2genind as.genambig

    @@ -215,6 +215,12 @@

    recode_polyploids

    + + + +

    Switch between genind and genclone objects.

    Switch between genind and genclone objects.

    Recode polyploid microsatellite data for use in frequency based statistics.

    +

    make_haplotypes

    +

    Split samples from a genind object into pseudo-haplotypes

    diff --git a/docs/reference/info_table.html b/docs/reference/info_table.html index 5cae1068..fce503b9 100644 --- a/docs/reference/info_table.html +++ b/docs/reference/info_table.html @@ -204,8 +204,8 @@

    Details

    Examples

    data(nancycats) -nancy.miss <- info_table(nancycats, plot = TRUE, type = "missing")
    data(Pinf) -Pinf.ploid <- info_table(Pinf, plot = TRUE, type = "ploidy")
    +nancy.miss <- info_table(nancycats, plot = TRUE, type = "missing")
    data(Pinf) +Pinf.ploid <- info_table(Pinf, plot = TRUE, type = "ploidy")
    # Get a table -atab <- mlg.table(Aeut, color = TRUE)
    atab
    #> MLG.1 MLG.2 MLG.3 MLG.4 MLG.5 MLG.6 MLG.7 MLG.8 MLG.9 MLG.10 MLG.11 +atab <- mlg.table(Aeut, color = TRUE)
    atab
    #> MLG.1 MLG.2 MLG.3 MLG.4 MLG.5 MLG.6 MLG.7 MLG.8 MLG.9 MLG.10 MLG.11 #> Athena 1 0 0 0 0 0 1 1 1 1 1 #> Mt. Vernon 0 2 1 1 1 1 0 0 0 0 0 #> MLG.12 MLG.13 MLG.14 MLG.15 MLG.16 MLG.17 MLG.18 MLG.19 MLG.20 diff --git a/docs/reference/plot_poppr_msn-1.png b/docs/reference/plot_poppr_msn-1.png index dc91f7e5..40dd3272 100644 Binary files a/docs/reference/plot_poppr_msn-1.png and b/docs/reference/plot_poppr_msn-1.png differ diff --git a/docs/reference/plot_poppr_msn.html b/docs/reference/plot_poppr_msn.html index 51810cca..fabfa64a 100644 --- a/docs/reference/plot_poppr_msn.html +++ b/docs/reference/plot_poppr_msn.html @@ -320,9 +320,17 @@

    Examp amsn <- poppr.msn(Aeut, adist, showplot = FALSE) # Default -library("igraph") # To get all the layouts. -set.seed(500) -plot_poppr_msn(Aeut, amsn, gadj = 15)

    +library("igraph") # To get all the layouts.
    #> +#> Attaching package: ‘igraph’
    #> The following object is masked from ‘package:poppr’: +#> +#> diversity
    #> The following object is masked from ‘package:testthat’: +#> +#> compare
    #> The following objects are masked from ‘package:stats’: +#> +#> decompose, spectrum
    #> The following object is masked from ‘package:base’: +#> +#> union
    set.seed(500) +plot_poppr_msn(Aeut, amsn, gadj = 15)
    # NOT RUN { # Different layouts (from igraph) can be used by supplying the function name. set.seed(500) diff --git a/docs/reference/poppr.amova.html b/docs/reference/poppr.amova.html index 2398f431..7e54c212 100644 --- a/docs/reference/poppr.amova.html +++ b/docs/reference/poppr.amova.html @@ -104,73 +104,79 @@

    Perform Analysis of Molecular Variance (AMOVA) on genind or genclone objects

    This function simplifies the process necessary for performing AMOVA in R. It gives user the choice of utilizing either the ade4 or the pegas -implementation of AMOVA. See amova (ade4) and -amova (pegas) for details on the specific -implementation.

    +implementation of AMOVA. See ade4::amova() (ade4) and pegas::amova() +(pegas) for details on the specific implementation.

    poppr.amova(x, hier = NULL, clonecorrect = FALSE, within = TRUE,
    -  dist = NULL, squared = TRUE, correction = "quasieuclid", sep = "_",
    -  filter = FALSE, threshold = 0, algorithm = "farthest_neighbor",
    -  missing = "loci", cutoff = 0.05, quiet = FALSE, method = c("ade4",
    -  "pegas"), nperm = 0)
    + dist = NULL, squared = TRUE, freq = TRUE, correction = "quasieuclid", + sep = "_", filter = FALSE, threshold = 0, + algorithm = "farthest_neighbor", missing = "loci", cutoff = 0.05, + quiet = FALSE, method = c("ade4", "pegas"), nperm = 0)

    Arguments

    - + - + - + - + +set to NULL (default), the raw pairwise distances will be calculated via +dist().

    + + + + - + - + - + @@ -193,23 +199,22 @@

    Ar

    +options given in the function missingno(). Default is "loci".

    +removed/modified. See missingno() for details.

    +corrections will be printed to the screen. If TRUE, no messages will be +printed.

    - @@ -222,79 +227,88 @@

    Ar

    Value

    -

    a list of class amova from the ade4 package. See - amova for details.

    +

    a list of class amova from the ade4 or pegas package. See +ade4::amova() or pegas::amova() for details.

    Details

    The poppr implementation of AMOVA is a very detailed wrapper for the - ade4 implementation. The output is an amova class list - that contains the results in the first four elements. The inputs are - contained in the last three elements. The inputs required for the ade4 - implementation are:

      +ade4 implementation. The output is an ade4::amova() class list that +contains the results in the first four elements. The inputs are contained +in the last three elements. The inputs required for the ade4 implementation +are:

      1. a distance matrix on all unique genotypes (haplotypes)

      2. a data frame defining the hierarchy of the distance matrix

      3. a genotype (haplotype) frequency table.

      4. -

      All of this data can be constructed from a genind - object, but can be daunting for a novice R user. This function - automates the entire process. Since there are many variables regarding - genetic data, some points need to be highlighted:

      -

      On Hierarchies:

      -The hierarchy is defined by different - population strata that separate your data hierarchically. These strata are - defined in the strata slot of genind and - genclone

      objects. They are useful for defining the - population factor for your data. See the function strata for - details on how to properly define these strata.

      -

      On Within Individual Variance:

      +
    +

    All of this data can be constructed from a genind object, +but can be daunting for a novice R user. This function automates the +entire process. Since there are many variables regarding genetic data, +some points need to be highlighted:

    On Hierarchies:

    +The hierarchy is defined by different +population strata that separate your data hierarchically. These strata are +defined in the strata slot of genind and +genclone objects. They are useful for defining the +population factor for your data. See the function strata() for details on +how to properly define these strata. +

    On Within Individual Variance:

    Heterozygosities within - diploid genotypes are sources of variation from within individuals and can - be quantified in AMOVA. When within = TRUE, poppr will split diploid - genotypes into haplotypes and use those to calculate within-individual - variance. No estimation of phase is made. This acts much like the default - settings for AMOVA in the Arlequin software package. Within individual - variance will not be calculated for haploid individuals or dominant - markers. -

    On Euclidean Distances:

    - AMOVA, as defined by - Excoffier et al., utilizes an absolute genetic distance measured in the - number of differences between two samples across all loci. With the ade4 - implementation of AMOVA (utilized by poppr), distances must be Euclidean - (due to the nature of the calculations). Unfortunately, many genetic - distance measures are not always euclidean and must be corrected for before - being analyzed. Poppr automates this with three methods implemented in - ade4, quasieuclid, lingoes, and - cailliez. The correction of these distances should not - adversely affect the outcome of the analysis. -

    On Filtering:

    +genotypes are sources of variation from within individuals and can be +quantified in AMOVA. When within = TRUE, poppr will split genotypes into +haplotypes with the function make_haplotypes() and use those to calculate +within-individual variance. No estimation of phase is made. This acts much +like the default settings for AMOVA in the Arlequin software package. +Within individual variance will not be calculated for haploid individuals +or dominant markers as the haplotypes cannot be split further. Setting +within = FALSE uses the euclidean distance of the allele frequencies +within each individual +

    On Euclidean Distances:

    + With the ade4 implementation of AMOVA +(utilized by poppr), distances must be Euclidean (due to the nature of the +calculations). Unfortunately, many genetic distance measures are not always +euclidean and must be corrected for before being analyzed. Poppr automates +this with three methods implemented in ade4, quasieuclid(), lingoes(), +and cailliez(). The correction of these distances should not adversely +affect the outcome of the analysis. +

    On Filtering:

    Filtering multilocus genotypes is performed by - mlg.filter. This can necessarily only be done AMOVA tests - that do not account for within-individual variance. The distance matrix used - to calculate the amova is derived from using mlg.filter with - the option stats = "distance", which reports the distance between - multilocus genotype clusters. One useful way to utilize this feature is to - correct for genotypes that have equivalent distance due to missing data. - (See example below.) -

    On Methods:

    - Both ade4 and pegas have - implementations of AMOVA, both of which are appropriately called "amova". - The ade4 version is faster, but there have been questions raised as to the - validity of the code utilized. The pegas version is slower, but careful - measures have been implemented as to the accuracy of the method. It must be - noted that there appears to be a bug regarding permuting analyses where - within individual variance is accounted for (within = TRUE) in the - pegas implementation. If you want to perform permutation analyses on the - pegas implementation, you must set within = FALSE. In addition, - while clone correction is implemented for both methods, filtering is only - implemented for the ade4 version. - -

    Note

    - -

    The ade4 function randtest.amova contains a slight - bug as of version 1.7.4 which causes the wrong alternative hypothesis to be - applied on every 4th heirarchical level. Luckily, there is a way to fix it - by re-converting the results with the function - as.krandtest. See examples for details.

    +mlg.filter(). This can necessarily only be done AMOVA tests that do not +account for within-individual variance. The distance matrix used to +calculate the amova is derived from using mlg.filter() with the option +stats = "distance", which reports the distance between multilocus +genotype clusters. One useful way to utilize this feature is to correct for +genotypes that have equivalent distance due to missing data. (See example +below.) +

    On Methods:

    + Both ade4 and pegas have +implementations of AMOVA, both of which are appropriately called "amova". +The ade4 version is faster, but there have been questions raised as to the +validity of the code utilized. The pegas version is slower, but careful +measures have been implemented as to the accuracy of the method. It must be +noted that there appears to be a bug regarding permuting analyses where +within individual variance is accounted for (within = TRUE) in the pegas +implementation. If you want to perform permutation analyses on the pegas +implementation, you must set within = FALSE. In addition, while clone +correction is implemented for both methods, filtering is only implemented +for the ade4 version. +

    On Polyploids:

    + As of poppr version 2.7.0, this +function is able to calculate phi statistics for within-individual variance +for polyploid data with full dosage information. When a data set does +not contain full dosage information for all samples, then the resulting +pseudo-haplotypes will contain missing data, which would result in an +incorrect estimate of variance. + Instead, the AMOVA will be performed on the distance matrix derived from +allele counts or allele frequencies, depending on the freq option. This +has been shown to be robust to estimates with mixed ploidy (Ronfort et al. +1998; Meirmans and Liu 2018). If you wish to brute-force your way to +estimating AMOVA using missing values, you can split your haplotypes with +the make_haplotypes() function. + One strategy for addressing ambiguous dosage in your polyploid data set +would be to convert your data to polysat's genambig class with the +as.genambig(), estimate allele frequencies with polysat::deSilvaFreq(), +and use these frequencies to randomly sample alleles to fill in the +ambiguous alleles.

    References

    @@ -302,12 +316,17 @@

    R molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics, 131, 479-491.

    +

    Ronfort, J., Jenczewski, E., Bataillon, T., and Rousset, F. (1998). Analysis +of population structure in autotetraploid species. Genetics, 150, +921–930.

    +

    Meirmans, P., Liu, S. (2018) Analysis of Molecular Variance (AMOVA) for +Autopolyploids Submitted.

    See also

    -

    amova (ade4) amova (pegas) - clonecorrect diss.dist missingno - is.euclid strata

    +

    ade4::amova(), pegas::amova(), clonecorrect(), diss.dist(), +missingno(), ade4::is.euclid(), strata(), make_haplotypes(), +as.genambig()

    Examples

    @@ -350,7 +369,7 @@

    Examp #> Phi-samples-Pop 0.2803128 #> Phi-Pop-total 0.7000679 #>
    amova.test <- randtest(amova.result) # Test for significance -plot(amova.test)
    amova.test
    #> class: krandtest lightkrandtest +plot(amova.test)
    amova.test
    #> class: krandtest lightkrandtest #> Monte-Carlo tests #> Call: randtest.amova(xtest = amova.result) #> @@ -383,13 +402,6 @@

    Examp poppr.amova(monpop, ~Symptom/Year) # gets a warning of zero distances poppr.amova(monpop, ~Symptom/Year, filter = TRUE, threshold = 0.1) # no warning -# Correcting incorrect alternate hypotheses with >2 heirarchical levels -# -mon.amova <- poppr.amova(monpop, ~Symptom/Year/Tree) -mon.test <- randtest(mon.amova) -mon.test # Note alter is less, greater, greater, less -alt <- c("less", "greater", "greater", "greater") # extend this to the number of levels -with(mon.test, as.krandtest(sim, obs, alter = alt, call = call, names = names)) # }

    @@ -402,8 +414,6 @@

    Contents

  • Details
  • -
  • Note
  • -
  • References
  • See also
  • diff --git a/docs/reference/poppr.msn-1.png b/docs/reference/poppr.msn-1.png index fe285902..13536646 100644 Binary files a/docs/reference/poppr.msn-1.png and b/docs/reference/poppr.msn-1.png differ diff --git a/docs/reference/poppr.msn.html b/docs/reference/poppr.msn.html index 627e4484..8d79555a 100644 --- a/docs/reference/poppr.msn.html +++ b/docs/reference/poppr.msn.html @@ -288,7 +288,7 @@

    Examp A.dist <- diss.dist(Aeut) # Graph it. -A.msn <- poppr.msn(Aeut, A.dist, gadj = 15, vertex.label = NA)
    +A.msn <- poppr.msn(Aeut, A.dist, gadj = 15, vertex.label = NA)
    # Find the sizes of the nodes (number of individuals per MLL): igraph::vertex_attr(A.msn$graph, "size")^2
    #> [1] 2 1 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 9 1 1 5 1 1 #> [26] 2 1 1 1 1 2 3 1 1 2 1 1 1 2 1 1 1 1 1 1 2 1 2 1 1 diff --git a/docs/reference/psex.html b/docs/reference/psex.html index 6f98664f..ac6a6beb 100644 --- a/docs/reference/psex.html +++ b/docs/reference/psex.html @@ -244,7 +244,7 @@

    Examp # With multiple encounters Pram_psex <- psex(Pram, by_pop = FALSE, method = "multiple") -plot(Pram_psex, log = "y", col = ifelse(Pram_psex > 0.05, "red", "blue"))

    abline(h = 0.05, lty = 2)
    title("Probability of multiple encounters")
    # NOT RUN { +plot(Pram_psex, log = "y", col = ifelse(Pram_psex > 0.05, "red", "blue"))
    abline(h = 0.05, lty = 2)
    title("Probability of multiple encounters")
    # NOT RUN { # For a single encounter (default) Pram_psex <- psex(Pram, by_pop = FALSE) plot(Pram_psex, log = "y", col = ifelse(Pram_psex > 0.05, "red", "blue")) diff --git a/docs/reference/samp.ia.html b/docs/reference/samp.ia.html index 9888e076..feae2556 100644 --- a/docs/reference/samp.ia.html +++ b/docs/reference/samp.ia.html @@ -186,12 +186,12 @@

    Examp n.snp.struc = 5e2, ploidy = 2, parallel = FALSE) position(x) <- sort(sample(1e4, 1e3)) -res <- samp.ia(x)
    #> | | | 0% | |= | 1% | |= | 2% | |== | 3% | |=== | 4% | |==== | 5% | |==== | 6% | |===== | 7% | |====== | 8% | |====== | 9% | |======= | 10% | |======== | 11% | |======== | 12% | |========= | 13% | |========== | 14% | |========== | 15% | |=========== | 16% | |============ | 17% | |============= | 18% | |============= | 19% | |============== | 20% | |=============== | 21% | |=============== | 22% | |================ | 23% | |================= | 24% | |================== | 25% | |================== | 26% | |=================== | 27% | |==================== | 28% | |==================== | 29% | |===================== | 30% | |====================== | 31% | |====================== | 32% | |======================= | 33% | |======================== | 34% | |======================== | 35% | |========================= | 36% | |========================== | 37% | |=========================== | 38% | |=========================== | 39% | |============================ | 40% | |============================= | 41% | |============================= | 42% | |============================== | 43% | |=============================== | 44% | |================================ | 45% | |================================ | 46% | |================================= | 47% | |================================== | 48% | |================================== | 49% | |=================================== | 50% | |==================================== | 51% | |==================================== | 52% | |===================================== | 53% | |====================================== | 54% | |====================================== | 55% | |======================================= | 56% | |======================================== | 57% | |========================================= | 58% | |========================================= | 59% | |========================================== | 60% | |=========================================== | 61% | |=========================================== | 62% | |============================================ | 63% | |============================================= | 64% | |============================================== | 65% | |============================================== | 66% | |=============================================== | 67% | |================================================ | 68% | |================================================ | 69% | |================================================= | 70% | |================================================== | 71% | |================================================== | 72% | |=================================================== | 73% | |==================================================== | 74% | |==================================================== | 75% | |===================================================== | 76% | |====================================================== | 77% | |======================================================= | 78% | |======================================================= | 79% | |======================================================== | 80% | |========================================================= | 81% | |========================================================= | 82% | |========================================================== | 83% | |=========================================================== | 84% | |============================================================ | 85% | |============================================================ | 86% | |============================================================= | 87% | |============================================================== | 88% | |============================================================== | 89% | |=============================================================== | 90% | |================================================================ | 91% | |================================================================ | 92% | |================================================================= | 93% | |================================================================== | 94% | |================================================================== | 95% | |=================================================================== | 96% | |==================================================================== | 97% | |===================================================================== | 98% | |===================================================================== | 99% | |======================================================================| 100%
    hist(res, breaks = "fd")
    +res <- samp.ia(x)
    #> | | | 0% | |= | 1% | |= | 2% | |== | 3% | |=== | 4% | |==== | 5% | |==== | 6% | |===== | 7% | |====== | 8% | |====== | 9% | |======= | 10% | |======== | 11% | |======== | 12% | |========= | 13% | |========== | 14% | |========== | 15% | |=========== | 16% | |============ | 17% | |============= | 18% | |============= | 19% | |============== | 20% | |=============== | 21% | |=============== | 22% | |================ | 23% | |================= | 24% | |================== | 25% | |================== | 26% | |=================== | 27% | |==================== | 28% | |==================== | 29% | |===================== | 30% | |====================== | 31% | |====================== | 32% | |======================= | 33% | |======================== | 34% | |======================== | 35% | |========================= | 36% | |========================== | 37% | |=========================== | 38% | |=========================== | 39% | |============================ | 40% | |============================= | 41% | |============================= | 42% | |============================== | 43% | |=============================== | 44% | |================================ | 45% | |================================ | 46% | |================================= | 47% | |================================== | 48% | |================================== | 49% | |=================================== | 50% | |==================================== | 51% | |==================================== | 52% | |===================================== | 53% | |====================================== | 54% | |====================================== | 55% | |======================================= | 56% | |======================================== | 57% | |========================================= | 58% | |========================================= | 59% | |========================================== | 60% | |=========================================== | 61% | |=========================================== | 62% | |============================================ | 63% | |============================================= | 64% | |============================================== | 65% | |============================================== | 66% | |=============================================== | 67% | |================================================ | 68% | |================================================ | 69% | |================================================= | 70% | |================================================== | 71% | |================================================== | 72% | |=================================================== | 73% | |==================================================== | 74% | |==================================================== | 75% | |===================================================== | 76% | |====================================================== | 77% | |======================================================= | 78% | |======================================================= | 79% | |======================================================== | 80% | |========================================================= | 81% | |========================================================= | 82% | |========================================================== | 83% | |=========================================================== | 84% | |============================================================ | 85% | |============================================================ | 86% | |============================================================= | 87% | |============================================================== | 88% | |============================================================== | 89% | |=============================================================== | 90% | |================================================================ | 91% | |================================================================ | 92% | |================================================================= | 93% | |================================================================== | 94% | |================================================================== | 95% | |=================================================================== | 96% | |==================================================================== | 97% | |===================================================================== | 98% | |===================================================================== | 99% | |======================================================================| 100%
    hist(res, breaks = "fd")
    # with unstructured snps assuming 1e4 positions set.seed(999) x <- glSim(n.ind = 10, n.snp.nonstruc = 1e3, ploidy = 2) position(x) <- sort(sample(1e4, 1e3)) -res <- samp.ia(x)
    #> | | | 0% | |= | 1% | |= | 2% | |== | 3% | |=== | 4% | |==== | 5% | |==== | 6% | |===== | 7% | |====== | 8% | |====== | 9% | |======= | 10% | |======== | 11% | |======== | 12% | |========= | 13% | |========== | 14% | |========== | 15% | |=========== | 16% | |============ | 17% | |============= | 18% | |============= | 19% | |============== | 20% | |=============== | 21% | |=============== | 22% | |================ | 23% | |================= | 24% | |================== | 25% | |================== | 26% | |=================== | 27% | |==================== | 28% | |==================== | 29% | |===================== | 30% | |====================== | 31% | |====================== | 32% | |======================= | 33% | |======================== | 34% | |======================== | 35% | |========================= | 36% | |========================== | 37% | |=========================== | 38% | |=========================== | 39% | |============================ | 40% | |============================= | 41% | |============================= | 42% | |============================== | 43% | |=============================== | 44% | |================================ | 45% | |================================ | 46% | |================================= | 47% | |================================== | 48% | |================================== | 49% | |=================================== | 50% | |==================================== | 51% | |==================================== | 52% | |===================================== | 53% | |====================================== | 54% | |====================================== | 55% | |======================================= | 56% | |======================================== | 57% | |========================================= | 58% | |========================================= | 59% | |========================================== | 60% | |=========================================== | 61% | |=========================================== | 62% | |============================================ | 63% | |============================================= | 64% | |============================================== | 65% | |============================================== | 66% | |=============================================== | 67% | |================================================ | 68% | |================================================ | 69% | |================================================= | 70% | |================================================== | 71% | |================================================== | 72% | |=================================================== | 73% | |==================================================== | 74% | |==================================================== | 75% | |===================================================== | 76% | |====================================================== | 77% | |======================================================= | 78% | |======================================================= | 79% | |======================================================== | 80% | |========================================================= | 81% | |========================================================= | 82% | |========================================================== | 83% | |=========================================================== | 84% | |============================================================ | 85% | |============================================================ | 86% | |============================================================= | 87% | |============================================================== | 88% | |============================================================== | 89% | |=============================================================== | 90% | |================================================================ | 91% | |================================================================ | 92% | |================================================================= | 93% | |================================================================== | 94% | |================================================================== | 95% | |=================================================================== | 96% | |==================================================================== | 97% | |===================================================================== | 98% | |===================================================================== | 99% | |======================================================================| 100%
    hist(res, breaks = "fd")
    +res <- samp.ia(x)
    #> | | | 0% | |= | 1% | |= | 2% | |== | 3% | |=== | 4% | |==== | 5% | |==== | 6% | |===== | 7% | |====== | 8% | |====== | 9% | |======= | 10% | |======== | 11% | |======== | 12% | |========= | 13% | |========== | 14% | |========== | 15% | |=========== | 16% | |============ | 17% | |============= | 18% | |============= | 19% | |============== | 20% | |=============== | 21% | |=============== | 22% | |================ | 23% | |================= | 24% | |================== | 25% | |================== | 26% | |=================== | 27% | |==================== | 28% | |==================== | 29% | |===================== | 30% | |====================== | 31% | |====================== | 32% | |======================= | 33% | |======================== | 34% | |======================== | 35% | |========================= | 36% | |========================== | 37% | |=========================== | 38% | |=========================== | 39% | |============================ | 40% | |============================= | 41% | |============================= | 42% | |============================== | 43% | |=============================== | 44% | |================================ | 45% | |================================ | 46% | |================================= | 47% | |================================== | 48% | |================================== | 49% | |=================================== | 50% | |==================================== | 51% | |==================================== | 52% | |===================================== | 53% | |====================================== | 54% | |====================================== | 55% | |======================================= | 56% | |======================================== | 57% | |========================================= | 58% | |========================================= | 59% | |========================================== | 60% | |=========================================== | 61% | |=========================================== | 62% | |============================================ | 63% | |============================================= | 64% | |============================================== | 65% | |============================================== | 66% | |=============================================== | 67% | |================================================ | 68% | |================================================ | 69% | |================================================= | 70% | |================================================== | 71% | |================================================== | 72% | |=================================================== | 73% | |==================================================== | 74% | |==================================================== | 75% | |===================================================== | 76% | |====================================================== | 77% | |======================================================= | 78% | |======================================================= | 79% | |======================================================== | 80% | |========================================================= | 81% | |========================================================= | 82% | |========================================================== | 83% | |=========================================================== | 84% | |============================================================ | 85% | |============================================================ | 86% | |============================================================= | 87% | |============================================================== | 88% | |============================================================== | 89% | |=============================================================== | 90% | |================================================================ | 91% | |================================================================ | 92% | |================================================================= | 93% | |================================================================== | 94% | |================================================================== | 95% | |=================================================================== | 96% | |==================================================================== | 97% | |===================================================================== | 98% | |===================================================================== | 99% | |======================================================================| 100%
    hist(res, breaks = "fd")
    #> | | | 0% | |== | 3% | |==== | 6% | |====== | 9% | |======== | 12% | |========== | 15% | |============ | 18% | |============== | 21% | |================ | 24% | |=================== | 26% | |===================== | 29% | |======================= | 32% | |========================= | 35% | |=========================== | 38% | |============================= | 41% | |=============================== | 44% | |================================= | 47% | |=================================== | 50% | |===================================== | 53% | |======================================= | 56% | |========================================= | 59% | |=========================================== | 62% | |============================================= | 65% | |=============================================== | 68% | |================================================= | 71% | |=================================================== | 74% | |====================================================== | 76% | |======================================================== | 79% | |========================================================== | 82% | |============================================================ | 85% | |============================================================== | 88% | |================================================================ | 91% | |================================================================== | 94% | |==================================================================== | 97% | |======================================================================| 100%
    plot(res, type = "l")
    +res <- win.ia(x, window = 300L) # Calculate for windows of size 300
    #> | | | 0% | |== | 3% | |==== | 6% | |====== | 9% | |======== | 12% | |========== | 15% | |============ | 18% | |============== | 21% | |================ | 24% | |=================== | 26% | |===================== | 29% | |======================= | 32% | |========================= | 35% | |=========================== | 38% | |============================= | 41% | |=============================== | 44% | |================================= | 47% | |=================================== | 50% | |===================================== | 53% | |======================================= | 56% | |========================================= | 59% | |=========================================== | 62% | |============================================= | 65% | |=============================================== | 68% | |================================================= | 71% | |=================================================== | 74% | |====================================================== | 76% | |======================================================== | 79% | |========================================================== | 82% | |============================================================ | 85% | |============================================================== | 88% | |================================================================ | 91% | |================================================================== | 94% | |==================================================================== | 97% | |======================================================================| 100%
    plot(res, type = "l")
    # NOT RUN { # unstructured snps set.seed(999)

    x

    a genind or genclone -object

    a genind or genclone object

    hier

    a hierarchical formula that defines your population -hierarchy. (e.g.: ~Population/Subpopulation). See Details below.

    a hierarchical formula that defines your population +hierarchy. (e.g.: ~Population/Subpopulation). See Details below.

    clonecorrect

    logical if TRUE, the data set will be clone -corrected with respect to the lowest level of the hierarchy. The default is -set to FALSE. See clonecorrect for details.

    logical if TRUE, the data set will be clone corrected +with respect to the lowest level of the hierarchy. The default is set to +FALSE. See clonecorrect() for details.

    within

    logical. When this is set to TRUE (Default), -variance within individuals are calculated as well. If this is set to -FALSE, The lowest level of the hierarchy will be the sample level. -See Details below.

    logical. When this is set to TRUE (Default), variance +within individuals are calculated as well. If this is set to FALSE, The +lowest level of the hierarchy will be the sample level. See Details below.

    dist

    an optional distance matrix calculated on your data. If this is -set to NULL (default), the raw pairwise distances will be calculated -via diss.dist.

    squared

    if a distance matrix is supplied, this indicates whether or not it represents squared distances.

    freq

    logical. If within = FALSE, the parameter rho is calculated +(Ronfort et al. 1998; Meirmans and Liu 2018). By setting freq = TRUE, +(default) allele counts will be converted to frequencies before the +distance is calculated, otherwise, the distance will be calculated on +allele counts, which can bias results in mixed-ploidy data sets. Note that +this option has no effect for haploid or presence/absence data sets.

    correction

    a character defining the correction method for -non-euclidean distances. Options are quasieuclid -(Default), lingoes, and cailliez. -See Details below.

    a character defining the correction method for +non-euclidean distances. Options are ade4::quasieuclid() (Default), +ade4::lingoes(), and ade4::cailliez(). See Details below.

    sep

    Deprecated. As of poppr version 2, this argument serves no purpose.

    Deprecated. As of poppr version 2, this argument serves no +purpose.

    filter

    logical When set to TRUE, mlg.filter will be run -to determine genotypes from the distance matrix. It defaults to -FALSE. You can set the parameters with algorithm and -threshold arguments. Note that this will not be performed when -within = TRUE. Note that the threshold should be the number of -allowable substitutions if you don't supply a distance matrix.

    logical When set to TRUE, mlg.filter will be run to +determine genotypes from the distance matrix. It defaults to FALSE. You +can set the parameters with algorithm and threshold arguments. Note +that this will not be performed when within = TRUE. Note that the +threshold should be the number of allowable substitutions if you don't +supply a distance matrix.

    threshold
    missing

    specify method of correcting for missing data utilizing -options given in the function missingno. Default is -"loci".

    cutoff

    specify the level at which missing data should be -removed/modified. See missingno for details.

    quiet

    logical If FALSE (Default), messages regarding any -corrections will be printed to the screen. If TRUE, no messages will -be printed.

    method

    Which method for calculating AMOVA should be used? Choices +

    Which method for calculating AMOVA should be used? Choices refer to package implementations: "ade4" (default) or "pegas". See details for differences.