PlotSections: Plot multiple corresponding age rangse for measured data

In MDSS: Modeling Human Dentin Serial Sectioning

Plot multiple corresponding age rangse for measured data

Description

PlotSections plots the corresponding age range for horizontal dentin serial sections with multiple class "cca" objects.

Usage

PlotSections(age, delta, obj,
  error = 0.1, col = "black", pch = 21, ...)

Arguments

age	List of user-assigned ages for each sections.
delta	List of vector observed delta values for each sections.
obj	List of class "cca" objects.
error	Vector of analytical error. Used for the height of "beans".
col	Vector of "col" values for plot.
pch	Vector of "pch" values for plot.
...	Additional arguments passed to plot.

Details

Function PlotSections plots corresponding age range with its measured delta values and equally assigned ages for a series of serial sections. Multiple trajectories can be plotted together into a single figure.

Value

A plot of multiple trajectories of dentin serial sections.

Author(s)

Takumi Tsutaya developed this model.

References

Tsutaya T 2020. Blurred time resolution of tooth dentin serial sections. American Journal of Physical Anthropology 173:748–759. DOI: 10.1002/ajpa.24113.

See Also

MDSS, mds, cca

Examples

## Example of the first molar of ST61 from the Unseiji temple.
unseiji_ST61_M1.mds <- mds(
  tooth = "M1",
  rt.y = unseiji_ST61_M1$rt.y,
  ecj = unseiji_ST61_M1$ecj,
  section.y = unseiji_ST61_M1$section.y,
  n.GL = 100)
unseiji_ST61_M1.cca <- cca(
  x = unseiji_ST61_M1.mds,
  t.start = unseiji_ST61_M1$t.start,
  t.end = unseiji_ST61_M1$t.end)
unseiji_ST61_C.mds <- mds(
  tooth = "C",
  rt.y = unseiji_ST61_C$rt.y,
  ecj = unseiji_ST61_C$ecj,
  section.y = unseiji_ST61_C$section.y,
  n.GL = 100)
unseiji_ST61_C.cca <- cca(
  x = unseiji_ST61_C.mds,
  t.start = unseiji_ST61_C$t.start,
  t.end = unseiji_ST61_C$t.end)

# Plot C and M1 modeled age ranges with stable nitrogen isotope ratios.
PlotSections(
  age = list(
    unseiji_ST61_C$section$age,
    unseiji_ST61_M1$section$age),
  delta = list(
    unseiji_ST61_C$section$d15N,
    unseiji_ST61_M1$section$d15N),
  obj = list(
    C = unseiji_ST61_C.cca,
    M1 = unseiji_ST61_M1.cca))

# Plot the same data with different colors.
PlotSections(
  age = list(
    unseiji_ST61_C$section$age,
    unseiji_ST61_M1$section$age),
  delta = list(
    unseiji_ST61_C$section$d15N,
    unseiji_ST61_M1$section$d15N),
  obj = list(
    C = unseiji_ST61_C.cca,
    M1 = unseiji_ST61_M1.cca),
  col = c("red", "blue"),
  pch = c(1, 2))

MDSS documentation built on Aug. 22, 2022, 1:05 a.m.