# to find out what's in this package, run the following
?archdata

# or go to
# https://cran.rstudio.com/web/packages/archdata/archdata.pdf

library("archdata")

data("RBGlass1")
RBGlass1


RBGlass1.pca <- prcomp(RBGlass1[, -1], scale.=TRUE)

biplot(RBGlass1.pca, xlabs=abbreviate(RBGlass1$Site, 1), cex=.75)

PitHouses

data(PitHouses)
# Crosstabulation of Hearths with Size
PitHouses.tbl <- xtabs(~Hearths+Size, PitHouses)
PitHouses.tbl
barplot(PitHouses.tbl, ylab="Frequency", main="Arctic Norway Pithouses", beside=TRUE,legend.text=TRUE, args.legend=list(title="Hearths"))


data(Fibulae)
Fibulae

t(sapply(Fibulae[, 3:16], quantile, na.rm=TRUE))
plot(density(Fibulae$Length, bw="SJ"), main="Kernel Density Plot of Length")

data(EndScrapers)
xtabs(Freq~Site+Curvature, EndScrapers)
xtabs(Freq~Curvature+Sides+Site, EndScrapers)

data(Olorgesailie.maj)
Olorgesailie.maj
# Compute percentages over the localities
Olor.pct <- prop.table(as.matrix(Olorgesailie.maj), 1)*100
boxplot(Olor.pct)

data(Olorgesailie.sub)

Olorgesailie.sub

chisq.test(Olorgesailie.sub[,3:18], simulate.p.value=TRUE)

BSR_XRF <- read.csv(file="./data-import/xrd-majors-bsr-brickstamps.csv", header=TRUE, sep=",")

BSR_XRF

BSR.pca <- prcomp(BSR_XRF[, -1], scale.=TRUE)
biplot(BSR.pca, xlabs=abbreviate(BSR_XRF$Sample, 1), cex=.75)