Arnhofen: Point pattern of mining pits from the Neolithic chert mine at...

In archdata: Example Datasets from Archaeological Research

Description

The list object contains two data tables of coordinates, one representing the centers of round mining pits, the other holding the vertices of the observed polygonal area. Data was generated by G. Roth in 2006 (Roth 2008). Spatial distance unit is meter. For converting the data to a point pattern see Examples.

Usage

data(Arnhofen)

Format

A list with two entries. The first is a data frame, points, with 216 observations of x and y coordinates. Each line represents the center of a round mining pit. The second is a data frame, window, with 100 observations and 3 variables, x, y, and vertex ID. Each line represents a vertex of the polygonal observation area.

points$x

(num) x coordinates of pit centers in m

points$y

(num) y coordinates of pit centers in m

window$x

(num) x coordinates of observation window vertices in m

window$y

(num) y coordinates of observation window vertices in m

window$id

(num) id for observation window vertices

Details

The coordinates in dataframe points represent the centres of 216 mining pits in the southeastern section of the 2001 excavation at the Neolithic chert mine of Abensberg-Arnhofen (Roth 2008). Direct dates for mining from the excavation place the site at 5300-4200 BC while use of mined material ends with the Bell Beaker Culture suggesting an end date for the mine of ca. 2200 BC. The regular pit pattern presented here dates to about 4200 BC, i.e. to the Münchshöfen Culture (4500-4000 BC). Arnhofen is the largest neolithic mine in Germany. The mining organization was analyzed by Roth (2008) using point pattern analysis (cf. Baddeley et al. 2016) which showed the neolithic mining to be conducted by farmers from surrounding villages (presumably on a seasonal basis).

The pit centers were located manually in a GIS using excavation maps from between 3 to 4 m below present surface. Mining pits were similar to vertical tubes with an average diameter of 1 m and a rounded horizontal section. A point therefore represents the center of such a vertical tube in the horizontal plane. A few of the pits reached a depth of nearly 8 m below surface. The vertices of the observation area polygon circumscribe a slightly smaller region than the excavated area. The list contains the additional attributes: reference for the data, short data description (site) and geographical coordinates (Lat/Lon) of the excavation.

Source

Roth, G. 2008. Geben und Nehmen, Eine Wirtschaftshistorische Studie zum Neolithischen Hornsteinbergbau von Abensberg-Arnhofen, Kr. Kelheim (Niederbayern) [in 4 volumes]. online PhD-thesis, University of Cologne 2008. http://kups.ub.uni-koeln.de/4176.

References

Baddeley, A., E. Rubak and R. Turner. 2016. Spatial Point Patterns: Methodology and Applications with R. CRC Press. Boca Raton.

Examples

# data and package spatstat by A. Baddeley et al. 2016 for point pattern analysis
# package spatstat is described and illustrated in Baddeley et al. (2016)
if (requireNamespace("spatstat", quietly = TRUE)) {
library(spatstat)
data(Arnhofen)
ap <- Arnhofen	      # to shorten the following code

# generate observation window object; note the polygonal outline.
arnwin <- owin(poly=ap$window[, 1:2])	

# generate point process pattern object from points and owin object
app <- ppp(ap$points$x, ap$points$y, arnwin) 
unitname(app) <- c("metre", "metres")   # optional, asign unitnames

# note that owin vertices traverse the polygon anticlockwise
plot(arnwin)
points(ap$window[, 1:2], pch=3, cex=.5)
text(ap$window[, 1], ap$window[, 2], ap$window[, 3], pos=3, cex=.7)

# visual inspection of the point process pattern 
plot(app)

# Computing the summary function "centered Besag's L" assuming
# homogeneous intensity. Centered Besag's L is just a conveniently 
# transformed Ripley's K. see references in ?Lest.

set.seed(1)	
Lcentrd <- envelope(app, Lest, nsim=49, nrank=1, global=TRUE, r=seq(0,7, 0.01),
     correction="translate", transform=expression(.-r))
# for the arguments see ?Kest and ?envelope. 

tm <- "Centered Besags's L for Arnhofen-Southeast"  # title

plot(Lcentrd, legendpos="bottomright", legendargs=list(bg="white"),
     main=tm, las=1)

# The deviations below envelopes suggest regular inter point distances
# at the 1 percent level - deviations above would have suggested clustering
# with r representing the radius of round clusters.

plot(Lcentrd, xlim=c(.5,2), legendpos="topright", legendargs=list(bg="white"),
     las=1, main=tm)
(inhibr <- Lcentrd$r[Lcentrd$obs<Lcentrd$lo])	

# significant inhibition between pits with an average diameter of 1 m pits
# were spaced at regular distances up to about 0.7 m apart:

max(inhibr) - 1

citation("spatstat") # don't forget to reference the method.
} else {
  cat("This example requires package spatstat.\n")
}

Example output

Loading required package: spatstat.data
Loading required package: nlme
Loading required package: rpart

spatstat 1.64-1       (nickname: ‘Help you I can, yes!’) 
For an introduction to spatstat, type ‘beginner’ 


Note: spatstat version 1.64-1 is out of date by more than 10 months; we recommend upgrading to the latest version.
Generating 49 simulations of CSR  ...
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.

Done.

To cite spatstat in publications use:

  Adrian Baddeley, Ege Rubak, Rolf Turner (2015). Spatial Point
  Patterns: Methodology and Applications with R. London: Chapman and
  Hall/CRC Press, 2015. URL
  http://www.crcpress.com/Spatial-Point-Patterns-Methodology-and-Applications-with-R/Baddeley-Rubak-Turner/9781482210200/

If you use hybrid models, please also cite:

  Adrian Baddeley, Rolf Turner, Jorge Mateu, Andrew Bevan (2013).
  Hybrids of Gibbs Point Process Models and Their Implementation.
  Journal of Statistical Software, 55(11), 1-43. URL
  http://www.jstatsoft.org/v55/i11/.

In survey articles, please cite the original paper on spatstat:

  Adrian Baddeley, Rolf Turner (2005). spatstat: An R Package for
  Analyzing Spatial Point Patterns. Journal of Statistical Software
  12(6), 1-42. URL http://www.jstatsoft.org/v12/i06/.

To see these entries in BibTeX format, use 'print(<citation>,
bibtex=TRUE)', 'toBibtex(.)', or set
'options(citation.bibtex.max=999)'.

archdata documentation built on Jan. 16, 2021, 5:35 p.m.