aion

aion est un package R qui fournit un système de classes et de méthodes permettant de représenter et de manipuler des séries chronologiques et des intervalles de temps archéologiques. Les dates sont représentées sous forme de rata die et peuvent être converties en n’importe quel calendrier défini par Reingold et Dershowitz (2018).

Une fois qu’une série chronologique est créée avec aion, n’importe quel calendrier peut être utilisé pour afficher ou visualiser les données (le calendrier grégorien est utilisé par défaut ; voir vignette("aion")).

aion ne fournit pas d’outils pour la modélisation chronologique. Il offre plutôt une API simple qui peut être utilisée par d’autres packages spécialisés (voir kairos ou ananke).

Pour citer aion dans une publication :

Frerebeau, N. (2024). aion : An R Package to Represent Archaeological Time Series. Journal of Open Source Software, 9(96). DOI : 10.21105/joss.06210.

Ce package fait partie du projet tesselle.

Installation

Vous pouvez installer la dernière version depuis le CRAN avec :

install.packages("aion")

Et la version de développement depuis Codeberg avec :

# install.packages("remotes")
remotes::install_git("https://codeberg.org/tesselle/aion")

Utilisation

## Installer les packages supplémentaires
## (si nécessaire)
# install.packages("folio")

## Charger le package
library(aion)

Données de comptage comme séries temporelles :

## Charger les données de comptage
## (données issues de Husi 2022)
data("loire", package = "folio")

## Conserver les variables dont la somme est supérieure à 600
keep <- c("01f", "01k", "01L", "08e", "08t", "09b", "15i", "15q")

## Calculer les moyennes
mid <- rowMeans(loire[, c("lower", "upper")])

## Création des séries temporelles
X <- series(
  object = loire[, keep],
  time = mid,
  calendar = AD()
)

## Visualiser (calendrier par défaut)
plot(
  x = X, 
  type = "h" # barres verticales
)

Visualiser des intervalles de temps :

## Création des intervalles
Y <- intervals(
  start = loire$lower,
  end = loire$upper,
  names = rownames(loire),
  calendar = AD()
)

## Visualiser (calendrier par défaut)
plot(x = Y)