3D STATS by sublayer

In this document we compute several 3d stats with the package spatstat https://www.rdocumentation.org/packages/spatstat/versions/1.49-0 including the Pair Correlation Function of a Three-Dimensional Point Pattern https://www.rdocumentation.org/packages/spatstat/versions/1.49-0/topics/pcf3est function, the K-function of a Three-Dimensional Point Pattern https://www.rdocumentation.org/packages/spatstat/versions/1.49-0/topics/K3est, the Empty Space Function of a Three-Dimensional Point Pattern https://www.rdocumentation.org/packages/spatstat/versions/1.49-0/topics/F3est and the Nearest Neighbour Distance Distribution Function of a Three-Dimensional Point Pattern https://www.rdocumentation.org/packages/spatstat/versions/1.49-0/topics/G3est.

First we load the required libraries

require(spatstat)

We’re ready to load the files with coordinates

sh1<-read.csv("suelo_h1.csv",header=FALSE)
sh2<-read.csv("suelo_h2.csv",header=FALSE)
sh3<-read.csv("suelo_h3.csv",header=FALSE)
sh4<-read.csv("suelo_h4.csv",header=FALSE)
suh<-list(sh1,sh2,sh3,sh4)

We begin with the faunal remains and LAYER H1

suh<-list(sh1,sh2,sh3,sh4)
par(mfrow=c(2,2))
tit<-c("Faunal remains","Pottery","Lithics","Gasteropods")
for (i in 1:4){
fauh<-subset(suh[[1]],suh[[1]][,6]==i)
fauh3d<-pp3(fauh[,2],fauh[,3],fauh[,4],box3(c(999,1004),c(995,999),c(407,408)))
plot.pp3(fauh3d, main=tit[i])
}

suh<-list(sh1,sh2,sh3,sh4)
par(mfrow=c(2,2))
tit<-c("Faunal remains","Pottery","Lithics","Gasteropods")
for (i in 1:4){
  
fauh<-subset(suh[[1]],suh[[1]][,6]==i)
fauh3d<-pp3(fauh[,2],fauh[,3],fauh[,4],box3(c(999,1004),c(995,999),c(407,408)))
p3<-pcf3est(fauh3d)
plot(p3,main="3D Pair Correlation")
k3<-K3est(fauh3d)
plot(k3,main="K-function")
g3<-G3est(fauh3d)
plot(g3,main="Nearest Neighbour Distribution")
f3<-F3est(fauh3d)
plot(f3,main="Empty space")
print(toupper(tit[i]))
}

[1] "FAUNAL REMAINS"

[1] "POTTERY"

[1] "LITHICS"

[1] "GASTEROPODS"