A collection of classes that represent archaeological data. This package provides a set of S4 classes that represent different special types of matrix (absolute/relative frequency, presence/absence data, co-occurrence matrix, etc.) upon which package developers can build subclasses. It also provides a set of generic methods (mutators and coercion mechanisms) and functions (e.g. predicates). In addition, a few classes of general interest (e.g. that represent stratigraphic relationships) are implemented.
You can install the released version of arkhe from CRAN with:
Or install the development version from GitHub with:
# install.packages("remotes") remotes::install_github("nfrerebeau/arkhe")
# Load the package library(arkhe)
arkhe provides a set of S4 classes that represent different special types of matrix.
CountMatrixrepresents absolute frequency data,
AbundanceMatrixrepresents relative frequency data,
OccurrenceMatrixrepresents a co-occurrence matrix,
SimilarityMatrixrepresents a (dis)similarity matrix,
IncidenceMatrixrepresents presence/absence data,
StratigraphicMatrixrepresents stratigraphic relationships.
It assumes that you keep your data tidy: each variable (type/taxa) must be saved in its own column and each observation (assemblage/sample) must be saved in its own row.
These new classes are of simple use, on the same way as the base
# Define a count data matrix # (data will be rounded to zero decimal places, then coerced with as.integer) quanti <- CountMatrix(data = sample(0:10, 100, TRUE), nrow = 10, ncol = 10) # Define a logical matrix # (data will be coerced with as.logical) quali <- IncidenceMatrix(data = sample(0:1, 100, TRUE), nrow = 10, ncol = 10)
arkhe uses coercing mechanisms (with validation methods) for data type conversions:
## Create a count matrix A0 <- matrix(data = sample(0:10, 100, TRUE), nrow = 10, ncol = 10) ## Coerce to absolute frequencies A1 <- as_count(A0) ## Coerce to relative frequencies B <- as_abundance(A1) ## Row sums are internally stored before coercing to a frequency matrix ## (use get_totals() to get these values) ## This allows to restore the source data A2 <- as_count(B) all(A1 == A2) #>  TRUE ## Coerce to presence/absence C <- as_incidence(A1) ## Coerce to a co-occurrence matrix D <- as_occurrence(A1)
Many familiar methods and group generic functions are available for all
*Matrix classes (such as
Logic…). In addition, all functions that call
as.data.frame first on their main argument should work (e. g.
Please note that all
*Matrix classes extend the R base
matrix, but the S3 part of the object does not store the data. Values are stored in a specific slot (allowing type checking).
X <- CountMatrix(data = sample(0:10, 25, TRUE), nrow = 5, ncol = 5) ## Get the S3 object S3Part(X, strictS3 = TRUE) #> col1 col2 col3 col4 col5 #> row1 1 6 11 16 21 #> row2 2 7 12 17 22 #> row3 3 8 13 18 23 #> row4 4 9 14 19 24 #> row5 5 10 15 20 25 ## Coerce to an S3 matrix as.matrix(X) #> col1 col2 col3 col4 col5 #> row1 5 3 1 0 5 #> row2 9 2 5 0 7 #> row3 4 7 0 3 10 #> row4 9 10 2 3 3 #> row5 5 9 3 9 9
Please note that the arkhe project is released with a Contributor Code of Conduct. By contributing to this project, you agree to abide by its terms.