Slice1: Slice1

In Zink-Antoine/TLpack: Analyse TL measurement using Bayesian approach, MCMC, Gibbs and slice samplers

Slice1

Description

MC Analysis TL (slice +gibbs) following gibbs1.R Hickey (2006) extrapolation y0=0 (gibbs1) additional scalar variable T (Slice sampler)

Usage

Slice1(Dose, df.T, df.y, mu_b, mu_0, var_b, var_0, var_y, n.iter)

Arguments

Dose	numeric (required) set of the irradiation doses
df.T	matrix (required) Luminescence data: the temperatures
df.y	matrix (required) Luminescence data; the luminescence signal
mu_b	numeric (required) Prior of the slope's average
mu_0	numeric (required) Prior of the true dose's average
var_b	numeric (required) Prior of the slope's variance
var_0	numeric (required) Prior of the true dose's variance
var_y	numeric (required) Prior of the intensity's variance
n.iter	numeric (required) the number of iteration

Value

an (n.iter x 3)-matrix

column	Type	Description
De	numeric	'True' Dose value
b	numeric	slope
T	numeric	Temperature

References

Gibbs sampler: Hickey, G. L. 2006. « The Linear Calibration Problem: A Bayesian Analysis ». PhD Thesis, PhD dissertation, University of Durham. 1–148. http://www.dur.ac.uk/g.l.hickey/dissertation.pdf.

chapter 3 and Appendix G.1

Slice sampler: Neal, R. 2003. « Slice Sampling ». Annals of Statistics 31 (3): 705‑67.

Examples

## Not run: 
##load data
if(dev.cur()!=1) dev.off()
data(multiTL, envir = environment())
attach(multiTL)
Slice1(Dose,df.T,df.y,mu_b,mu_0,var_b,var_0,var_y,n.iter)
detach(multiTL)

## End(Not run)

Zink-Antoine/TLpack documentation built on April 14, 2025, 1:58 p.m.