ProxyConversion: Convert between Temperature in Degrees C and Proxy Units

In EarthSystemDiagnostics/sedproxy: Simulation of Sediment Archived Climate Proxy Records

Convert between Temperature in Degrees C and Proxy Units

Description

A wrapper function for accessing proxy - temperature conversion functions

Usage

ProxyConversion(
  temperature = NULL,
  proxy.value = NULL,
  calibration.type = "identity",
  slp.int.means = NULL,
  slp.int.vcov = NULL,
  calibration = switch(calibration.type, identity = NA, Uk37 = "Mueller global", MgCa =
    "Ten planktonic species_350-500"),
  point.or.sample = c("point", "sample"),
  n = 1
)

Arguments

temperature	Temperature in degrees C
proxy.value	Temperature in proxy units
calibration.type	Type of proxy, e.g. Uk37 or MgCa
slp.int.means	Optional user supplied vector of values for the slope and intercept of the calibration function. Overides the defaults.
slp.int.vcov	Optional user supplied variance covariance matrix calibration parameters. Overides the defaults.
calibration	The name of a specific calibration for which calibration parameters are provided by sedproxy. Currently applies only to calibration.type MgCa.
point.or.sample	Use the "best estimate" calibration parameters, or parameters sampled from the fitted calibration model
n	the number of replicate conversions to make in the case of sampled calibration parameters

Details

Valid entries for calibration are: "Ten planktonic species_350-500", "G. aequilateralis_350-500", "G. aequilateralis_500-1000", "G. conglobatus_350-500", "G. hirsuta_350-500", "G. inflata_350-500", "G. ruber pink_250-350", "G. ruber pink_350-500", "G. ruber white_250-350", "G. ruber white_350-500", "G. sacculifer with sac_350-500", "G. sacculifer without sac_350-500", "G. truncatulinoides_350-500", "G. truncatulinoides_500-1000", "N. dutertrei_350-500", "O. univesa_350-500", "P. obliquiloculata_350-500"

Value

a vector of temperatures or proxy values

Examples

# From temperature to UK'37
## With fixed calibration
ProxyConversion(temperature = c(10, 20), point.or.sample = "point",
                calibration.type = "Uk37")

## With random calibration, 5 replicates
ProxyConversion(temperature = c(1, 2), n = 5, point.or.sample = "sample",
                calibration.type = "Uk37")


## Back-transformation with same calibration
ProxyConversion(
  proxy.value = as.vector(
    ProxyConversion(
      temperature = c(21, 22),
      calibration.type = "Uk37",
      point.or.sample = "point"
    )
  ),
  point.or.sample = "point",
  calibration.type = "Uk37"
)

## Back-transformation with random calibration
ProxyConversion(
  proxy.value = as.vector(
    ProxyConversion(
      temperature = c(21, 22),
      calibration.type = "Uk37",
     point.or.sample = "point"
    )
  )
  ,
  n = 5,
  point.or.sample = "sample",
  calibration.type = "Uk37"
)

## Incompatible arguments
## Not run: 
ProxyConversion(temperature = 1, proxy.value = 1)

## End(Not run)

EarthSystemDiagnostics/sedproxy documentation built on Feb. 27, 2023, 10:11 a.m.