Title: | Dynamic Range Boxes |
---|---|
Description: | Improves the concept of multivariate range boxes, which is highly susceptible for outliers and does not consider the distribution of the data. The package uses dynamic range boxes to overcome these problems. |
Authors: | Manuela Schreyer <manuelalarissa.schreyer@sbg.ac.at>, Wolfgang Trutschnig <Wolfgang.Trutschnig@sbg.ac.at>, Robert R. Junker <Robert.Junker@sbg.ac.at>, Jonas Kuppler <Jonas.Kuppler@uni-duesseldorf.de>, Arne Bathke <Arne.Bathke@sbg.ac.at>, Judith H. Parkinson <judith.parkinson@sbg.ac.at>, Raoul Kutil <Raoul_Hugo.Kutil@stud.sbg.ac.at> |
Maintainer: | Marco Tschimpke <marco.tschimpke@plus.ac.at> |
License: | GPL-2 |
Version: | 0.18 |
Built: | 2025-03-26 07:19:58 UTC |
Source: | CRAN |
The package DynRB improves the concept of multivariate range boxes, which is highly susceptible for outlines and does not consider the distribution of the data. The package uses dynamic range boxes to overcome these problems.
Package: | dynRB |
Type: | Package |
Version: | 0.16 |
Date: | 2021-05-11 |
Manuela Schreyer manuelalarissa.schreyer@sbg.ac.at,
Wolfgang Trutschnig Wolfgang.Trutschnig@sbg.ac.at,
Robert R. Junker Robert.Junker@sbg.ac.at (corresponding author),
Jonas Kuppler Jonas.Kuppler@uni-duesseldorf.de,
Arne Bathke Arne.Bathke@sbg.ac.at,
Judith H. Parkinson judith.parkinson@sbg.ac.at,
Raoul Kutil Raoul_Hugo.Kutil@stud.sbg.ac.at
Junker RR, Kuppler J, Bathke AC, Schreyer ML, Trutschnig W (2016) Dynamic range boxes - A robust non-parametric approach to quantify size and overlap of n-dimensional hypervolumes. Methods in Ecology and Evolution doi: 10.1111/2041-210X.12611
Judith H. Parkinson, Raoul Kutil, Jonas Kuppler, Robert R. Junker, Wolfgang Trutschnig, Arne C. Bathke: A Fast and Robust Way to Estimate Overlap of Niches and Draw Inference, International Journal of Biostatistics (2018)
# example function dynRB_VPa # for reliable results use steps = 201 data(finch2) r<-dynRB_VPa(finch2, steps = 101) r$result
# example function dynRB_VPa # for reliable results use steps = 201 data(finch2) r<-dynRB_VPa(finch2, steps = 101) r$result
Function returns pairwise overlaps for each dimension . Number of dynamic range boxes (
steps
) can be adjusted. Default: steps = 201
dynRB_Pn(A = A, steps = 201, correlogram = FALSE, row_col = c(2, 2))
dynRB_Pn(A = A, steps = 201, correlogram = FALSE, row_col = c(2, 2))
A |
Data frame, where the first column is a character vector containing the objects (e.g. species) and the other columns are numeric vectors (containing measurements). |
steps |
Number of range boxes. Default: |
correlogram |
If |
row_col |
Number of rows and columns of the figures (correlogram for each species). Default: |
Data frame containing the summarized overlaps for each pair of objects and dimension.
Manuela Schreyer manuelalarissa.schreyer@sbg.ac.at,
Wolfgang Trutschnig Wolfgang.Trutschnig@sbg.ac.at,
Robert R. Junker Robert.Junker@sbg.ac.at (corresponding author),
Jonas Kuppler Jonas.Kuppler@uni-duesseldorf.de,
Arne Bathke Arne.Bathke@sbg.ac.at
Junker RR, Kuppler J, Bathke AC, Schreyer ML, Trutschnig W (2016) Dynamic range boxes - A robust non-parametric approach to quantify size and overlap of n-dimensional hypervolumes. Methods in Ecology and Evolution doi: 10.1111/2041-210X.12611
# example function dynRB_Pn # for reliable results use steps = 201 data(finch2) r<-dynRB_Pn(finch2, steps = 101)
# example function dynRB_Pn # for reliable results use steps = 201 data(finch2) r<-dynRB_Pn(finch2, steps = 101)
Function returns Dynamic Range Box size of each dimension . Number of dynamic range boxes (
steps
) can be adjusted. Default: steps = 201
dynRB_Vn(A = A, steps = 201, correlogram = FALSE, row_col = c(2, 2))
dynRB_Vn(A = A, steps = 201, correlogram = FALSE, row_col = c(2, 2))
A |
Data frame, where the first column is a character vector and the other columns are numeric vectors. |
steps |
Number of range boxes. Default: |
correlogram |
If |
row_col |
Number of rows and columns of the figures (correlogram for each species). Default: |
Data frame containing the summarized niche length for each object and dimension.
Manuela Schreyer manuelalarissa.schreyer@sbg.ac.at,
Wolfgang Trutschnig Wolfgang.Trutschnig@sbg.ac.at,
Robert R. Junker Robert.Junker@sbg.ac.at (corresponding author),
Jonas Kuppler Jonas.Kuppler@uni-duesseldorf.de,
Arne Bathke Arne.Bathke@sbg.ac.at
Junker RR, Kuppler J, Bathke AC, Schreyer ML, Trutschnig W (2016) Dynamic range boxes - A robust non-parametric approach to quantify size and overlap of n-dimensional hypervolumes. Methods in Ecology and Evolution doi: 10.1111/2041-210X.12611
# example function dynRB_Vn # for reliable results use steps = 201 data(finch2) r<-dynRB_Vn(finch2, steps = 101)
# example function dynRB_Vn # for reliable results use steps = 201 data(finch2) r<-dynRB_Vn(finch2, steps = 101)
Function returns size and pairwise overlaps of niches or trait-spaces. Size or overlaps of dimensions can be aggregated by using either "product", "mean" or "geometric mean" as aggregation method. The results obtained by using the product are automatically printed. Number of dynamic range boxes (steps
) can be adjusted. Default: steps = 201
dynRB_VPa(A = A, steps = 201, correlogram = FALSE, row_col = c(2, 2), pca.corr = FALSE, var.thres = 0.9)
dynRB_VPa(A = A, steps = 201, correlogram = FALSE, row_col = c(2, 2), pca.corr = FALSE, var.thres = 0.9)
A |
Data frame, where the first column is a character vector and the other columns are numeric vectors. |
steps |
Number of range boxes. Default: |
correlogram |
If |
row_col |
Number of rows and columns of the figures (correlogram for each species). Default: |
pca.corr |
If |
var.thres |
Variance predicted by the PCA-axes, if |
Data frame containing the summarized niche overlap (and volume) for each pair of objects aggregated by all three possible choices (i.e. product, mean, geometric mean).
Manuela Schreyer manuelalarissa.schreyer@sbg.ac.at,
Wolfgang Trutschnig Wolfgang.Trutschnig@sbg.ac.at,
Robert R. Junker Robert.Junker@sbg.ac.at (corresponding author),
Jonas Kuppler Jonas.Kuppler@uni-duesseldorf.de,
Arne Bathke Arne.Bathke@sbg.ac.at
Junker RR, Kuppler J, Bathke AC, Schreyer ML, Trutschnig W (2016) Dynamic range boxes - A robust non-parametric approach to quantify size and overlap of n-dimensional hypervolumes. Methods in Ecology and Evolution doi: 10.1111/2041-210X.12611
# example function dynRB_VPa # for reliable results use steps = 201 data(finch2) r<-dynRB_VPa(finch2, steps = 101, correlogram = TRUE, row_col = c(1,1)) r$result
# example function dynRB_VPa # for reliable results use steps = 201 data(finch2) r<-dynRB_VPa(finch2, steps = 101, correlogram = TRUE, row_col = c(1,1)) r$result
To demonstrate the application of the functions for real world data, we used existing data sets on niches and trait-spaces and quantified their sizes and overlaps. The data set finch is a data set on morphological measurements of Darwin finches. The data set comprises quantitative measurements of nine traits characterizing five species of finches, each trait was measured at least in 10 individuals per species.
data("finch")
data("finch")
A data frame with 146 observations on the following 10 variables.
Species
a character vector of the Species Geospiza heliobates, Geospiza prosthemelas prosthemelas, Geospiza fuliginosa parvula, Geospiza fortis fortis and Geospiza fortis platyrhyncha
BodyL
a numeric vector
WingL
a numeric vector
TailL
a numeric vector
BeakW
a numeric vector
BeakH
a numeric vector
LBeakL
a numeric vector
UBeakL
a numeric vector
N.UBkL
a numeric vector
TarsusL
a numeric vector
Snodgrass R and Heller E (1904) Papers from the Hopkins-Stanford Galapagos Expedition, 1898-99. XVI. Birds. Proceedings of the Washington Academy of Sciences 5: 231-372.
data(finch) ## quick overview head(finch)
data(finch) ## quick overview head(finch)
To demonstrate the application of the functions for real world data, we used existing data sets on niches and trait-spaces and quantified their sizes and overlaps. The data set finch2 is a data set on morphological measurements of three Darwin finches. The data set comprises quantitative measurements of nine traits characterizing two species of finches, each trait was measured at least in 10 individuals per species.
data("finch2")
data("finch2")
A data frame with 103 observations on the following 10 variables.
Species
a character vector of the Species Geospiza fuliginosa parvula and Geospiza fortis fortis
BodyL
a numeric vector
WingL
a numeric vector
TailL
a numeric vector
BeakW
a numeric vector
BeakH
a numeric vector
LBeakL
a numeric vector
UBeakL
a numeric vector
N.UBkL
a numeric vector
TarsusL
a numeric vector
Snodgrass R and Heller E (1904) Papers from the Hopkins-Stanford Galapagos Expedition, 1898-99. XVI. Birds. Proceedings of the Washington Academy of Sciences 5: 231-372.
data(finch2) ## quick overview head(finch2)
data(finch2) ## quick overview head(finch2)
This functions can be used to show the graphics generated by the functions dynRB_Pn,dynRB_Vn and dynRB_VPa.
overview(r, row_col = c(3, 3))
overview(r, row_col = c(3, 3))
r |
Output of the function dynRB_Pn,dynRB_Vn or dynRB_VPa. |
row_col |
Number of rows and columns of the figures. Default: |
Manuela Schreyer manuelalarissa.schreyer@sbg.ac.at,
Wolfgang Trutschnig Wolfgang.Trutschnig@sbg.ac.at,
Robert R. Junker Robert.Junker@sbg.ac.at (corresponding author),
Jonas Kuppler Jonas.Kuppler@uni-duesseldorf.de,
Arne Bathke Arne.Bathke@sbg.ac.at
# example for the function dynRB_Pn # for reliable results use steps = 201 data(finch2) r<-dynRB_Pn(finch2, steps = 101) overview(r)
# example for the function dynRB_Pn # for reliable results use steps = 201 data(finch2) r<-dynRB_Pn(finch2, steps = 101) overview(r)
Function returns the asymmetric overlaps for each dimension, calculated by the method published by Parkinson et al. (2018) using ranks. Further two confidence intervals are returned for each estimate. The confidence level, as well as the repetitions for bootstrap can be adjusted.
ranks_OV(A = A, alpha = 0.05, reps4boot = 1000, digit = 3)
ranks_OV(A = A, alpha = 0.05, reps4boot = 1000, digit = 3)
A |
Data frame, where the first column contains two objects (e.g. species) and the other columns are numeric vectors (containing measurments). |
alpha |
The confidence level. Default: |
reps4boot |
Number of repetitions for the bootstrap. . Default: |
digit |
Number of digits after which the results are cut off. Default: |
Data Frame containing the two asymmetric overlaps for each dimension together with their confidence intervals. The last row contains the d-dimensional asymmetric overlaps.
Judith H. Parkinson judith.parkinson@sbg.ac.at,
Raoul Kutil Raoul_Hugo.Kutil@stud.sbg.ac.at,
Jonas Kuppler Jonas.Kuppler@uni-duesseldorf.de,
Robert R. Junker Robert.Junker@sbg.ac.at (corresponding author),
Wolfgang Trutschnig Wolfgang.Trutschnig@sbg.ac.at,
Arne Bathke Arne.Bathke@sbg.ac.at
Judith H. Parkinson, Raoul Kutil, Jonas Kuppler, Robert R. Junker, Wolfgang Trutschnig, Arne C. Bathke: A Fast and Robust Way to Estimate Overlap of Niches and Draw Inference, International Journal of Biostatistics (2018)
# example function ranks_OV data(finch2) head(finch2) ranks_OV(finch2[1:4], alpha = 0.05)
# example function ranks_OV data(finch2) head(finch2) ranks_OV(finch2[1:4], alpha = 0.05)