Datasets for teaching quantitative approaches and modeling in archaeology and paleontology. This package provides several types of data related to broad topics (cultural evolution, radiocarbon dating, paleoenvironments, etc.), which can be used to illustrate statistical methods in the classroom (multivariate data analysis, compositional data analysis, diversity measurement, etc.).

This package is for teaching purposes only: some datasets may be outdated.

When using **folio**, you must cite the relevant
publications for each dataset (see the documentation).

You can install the released version of **folio** from
CRAN with:

`install.packages("folio")`

And the development version from GitHub with:

```
# install.packages("remotes")
::install_github("tesselle/folio") remotes
```

```
## Arnold and Libby's Curve of Knowns
## Replicate fig. 1 from Arnold and Libby (1949)
data("arnold1949")
<- lm(activity_found ~ age_expected, data = arnold1949)
fit summary(fit)
#>
#> Call:
#> lm(formula = activity_found ~ age_expected, data = arnold1949)
#>
#> Residuals:
#> 1 2 3 4 5 6
#> 0.31341 -0.27618 -0.04573 -0.19344 0.09780 0.10414
#>
#> Coefficients:
#> Estimate Std. Error t value Pr(>|t|)
#> (Intercept) 12.2665215 0.2903726 42.24 1.88e-06 ***
#> age_expected -0.0011588 0.0000938 -12.35 0.000247 ***
#> ---
#> Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
#>
#> Residual standard error: 0.2421 on 4 degrees of freedom
#> Multiple R-squared: 0.9745, Adjusted R-squared: 0.9681
#> F-statistic: 152.6 on 1 and 4 DF, p-value: 0.0002467
plot(
x = arnold1949$age_expected,
y = arnold1949$activity_found,
type = "p",
pch = 16,
xlab = "Historical age (years)",
ylab = "Specific activity (C/m/g)",
main = "Samples of known ages"
)segments(
x0 = arnold1949$age_expected,
y0 = arnold1949$activity_found - arnold1949$activity_found_error,
x1 = arnold1949$age_expected,
y1 = arnold1949$activity_found + arnold1949$activity_found_error
)segments(
x0 = arnold1949$age_expected - arnold1949$age_expected_error,
y0 = arnold1949$activity_found,
x1 = arnold1949$age_expected + arnold1949$age_expected_error,
y1 = arnold1949$activity_found
)abline(fit, col = "red")
```

Please note that the **folio** project is released with
a Contributor Code of
Conduct. By contributing to this project, you agree to abide by its
terms.