palimpsestR: An R Package for the Identification and Probabilistic Decomposition of Archaeological Palimpsests

This paper presents palimpsestr, an R package and Shiny application dedicated to the identification and probabilistic decomposition of archaeological palimpsests - deposits in which material from multiple occupation phases is superimposed and partially intermixed. Such deposits represent one of the most persistent analytical challenges in field archaeology, and existing tools (Harris Matrix recording systems, k-means spatial clustering, GIS overlays) typically address only one evidence domain at a time. palimpsestr implements the Stratigraphic Entanglement Field (SEF) framework, which integrates four evidence domains - horizontal coordinates, vertical elevation, chronological range, and cultural class - into a single mixed-type finite mixture model fitted by Expectation-Maximisation: a diagonal-covariance Gaussian over the numeric features combined with a per-phase categorical (multinomial) distribution over the cultural class. The model is augmented with optional taphonomic weighting, a uniform noise/outlier component, and stratigraphic entanglement constraints derived from the Harris Matrix. Three interpretable diagnostics - the Stratigraphic Entanglement Index (SEI), Excavation Stratigraphic Energy (ESE), and Palimpsest Dissolution Index (PDI) - allow practitioners to assess deposit coherence, detect potentially redeposited finds, and evaluate the reliability of chronological attribution at the find, unit, and phase levels. The framework supports archaeologists in three principal interpretive tasks: assessing the reliability and coherence of stratigraphic units (distinguishing genuine palimpsests from recording errors) ; identifying intrusive finds with respect to both their typology and their stratigraphic position; and statistical estimation of type longevity (the duration of use of pottery and material categories). The package supports CSV/TSV/Excel/SQLite/PostgreSQL data import, GIS export to GeoPackage via the sf package, publication-quality plots via ggplot2 and interactive plots via plotly, and includes a built-in Shiny dashboard for non-programmatic use. Modelling the cultural class as a per-phase categorical distribution, rather than as one-hot Gaussian columns, prevents typology from splitting stratigraphic units across phases; the spatial and vertical similarity kernels are bounded and scale-invariant; the domain weights enter the likelihood and can be cross-validated; the noise component yields a genuine posterior probability of intrusion and shields the phase estimates from outliers; and optional treatments propagate per-find dating uncertainty and apply the stratigraphic constraint as a dynamically updated Neighborhood-EM field. The intrusion diagnostic also distinguishes the direction of the chronological mismatch (residual versus latent-feature), and a companion helper (recommendₛetup () ) inspects a dataset and reports when its recording resolution limits the achievable inference. We demonstrate the application through a case study at the multi-period Roman villa of Poggio Gramignano (Lugnano in Teverina, Italy). We also discuss the methodological assumptions of the SEF framework - in particular the implicit assumption of horizontal stratigraphy and the dependence of phase resolution on the spatial and chronological recording resolution of the data - and outline planned developments