Description Usage Arguments Value References See Also Examples

Calculating and selecting equivalent doses in a batch model according to the "standardised growth curve" (SGC) method suggested by Roberts and Duller (2004) or the "global standardised growth curve" (gSGC) method suggested by Li et al. (2015, 2016).

1 2 3 4 5 |

`obj_analyseBIN` |
list( |

`SGCpars` |
vector( |

`model` |
character( |

`origin` |
logical( |

`avgDev` |
numeric( |

`method` |
character(with default): method used for equivalent dose calculation, i.e., |

`SAR.Cycle` |
character(with default): SAR cycles used for SGC equivalent dose calculation. |

`errMethod` |
character(with default): method used for equivalent dose error assessment |

`Tn.above.3BG` |
logical(with default): logical value indicating if aliquot (grain) with Tn below 3 sigma BG should be rejected |

`TnBG.ratio.low` |
numeric(optional): lower limit on ratio of initial Tn signal to BG |

`rseTn.up` |
numeric(optional): upper limit on relative standard error of Tn in percent |

`FR.low` |
numeric(optional): lower limit on fast ratio of Tn |

`rseED.up` |
numeric(optional): upper limit on the relative standard error of equivalent dose in percent |

`use.se` |
logical(with default): logical value indicating if standard errors of values should be used during application of rejection criteria |

`outpdf` |
character(optional): if specified, results of SGC equivalent dose calculation will be written to a PDF file named |

`outfile` |
character(optional): if specified, SGC equivalent doses related quantities will be written to a CSV file named |

Return an invisible list that contains the following elements:

`scale.Ltx` |
scaled standardised natural-dose signals and associated standard errors used for SGC equivalent dose calculation. Note that standardised natural-dose signals will remain un-scaled if |

`sgcED` |
calculated SGC equivalent doses |

`ConfInt` |
68 percent (one sigma) and 95 percent (two sigma) confidence intervals of SGC equivalent doses |

`agID` |
aliquot (grain) ID of calculated (selected) SGC equivalent doses |

`summary.info` |
a summary of the SGC equivalent dose calculation |

Li B, Roberts RG, Jacobs Z, Li SH, 2015. Potential of establishing a "global standardised growth curve" (gSGC) for optical dating of quartz from sediments. Quaternary Geochronology, 27: 94-104.

Li B, Jacobs Z, Roberts RG, 2016. Investigation of the applicability of standardised growth curves for OSL dating of quartz from Haua Fteah cave, Libya. Quaternary Geochronology, 35: 1-15.

Roberts HM, Duller GAT, 2004. Standardised growth curves for optical dating of sediment using multiple-grain aliquots. Radiation Measurements, 38(2): 241-252.

fitGrowth; lsNORM; SARdata;

scaleSGCN; calED; calSARED

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 | ```
data(SARdata)
### (1) gSGC ED calculation:
### gSGCpars were obtained using function "lsNORM".
gSGCpars <- c(137.440874251, 0.007997863, 2.462035263, -0.321536177)
avg.error2 <- 0.1111623
res <- calSGCED(as_analyseBIN(SARdata), gSGCpars, method="gSGC",
model="gok", origin=FALSE, avgDev=avg.error2,
SAR.Cycle=c("N","R3"))
print(res$sgcED)
### (2) SGC ED calculation (not run):
### SGCpars were obtained using function "fitGrowth".
# SGCpars <- c(183.474322547, 0.007038048, 4.254287622, -0.337734151)
# avg.error <- 0.3156259
# calSGCED(as_analyseBIN(SARdata), SGCpars, method="SGC", model="gok",
# origin=FALSE, avgDev=avg.error, SAR.Cycle="N", outpdf="SGCED")
### (3) gSGC ED calculation and signal-related
### rejection criteria application (not run):
# data(BIN)
# res_pickBIN <-pickBINdata(BIN, LType="OSL")
# res_analyseBIN <- analyseBINdata(res_pickBIN, nfchn=4, nlchn=30)
# res_lsNORM <- lsNORM(res_analyseBIN$SARdata, model="gok", origin=FALSE)
# calSGCED(res_analyseBIN, SGCpars=res_lsNORM$LMpars2[,1],
# model="gok", origin=FALSE, avgDev=res_lsNORM$avg.error2,
# method="gSGC", SAR.Cycle=c("N","R3"), Tn.above.3BG=TRUE,
# TnBG.ratio.low=4, rseTn.up=10, outpdf="foo", outfile="foo")
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.