MDSS (version 1.0-1)

PlotSections: 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, ...)

Value

A plot of multiple trajectories of dentin serial sections.

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.

Author

Takumi Tsutaya developed this model.

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.

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))