estimate_range {ArchaeoPhases}  R Documentation 
Calculates the ranges of summary statistics from the output of two or more runs of the MCMC algorithm. Results are given in calendar years for statistics that estimate them.
estimate_range(
mcmc,
position,
app = "bcal",
estimates = c("mean", "q1", "median", "q3", "ci.inf", "ci.sup"),
quiet = "partial",
bin_width = 1,
decimal = ".",
separator = ","
)
mcmc 
A vector of path names to the MCMC files. 
position 
Numeric vector containing the positions of the columns corresponding to the MCMC chains of interest, or a vector of column names. 
app 
Name of the application that created the MCMC files,
one of 
estimates 
Numeric vector containing the positions of the columns
corresponding to the statistics of interest returned by the

quiet 
One of 
bin_width 
If 
decimal 
If 
separator 
If 
This function is useful for estimating the sensitivity of calibration results to different model parameters.
A list with the following components:
A matrix of estimate ranges.
The mean of the ranges in range_table
.
The standard deviation of the ranges in range_table
.
The minimum of the ranges in range_table
.
The median of the ranges in range_table
.
The maximum value of the ranges in range_table
.
Thomas S. Dye, tsd@tsdye.online
## Not run:
## Generate 0's
res < estimate_range(mcmc = c("http://tsdye.online/AP/ox.csv",
"http://tsdye.online/AP/ox.csv"), position = c(1, 2),
app = "oxcal", quiet = "yes")
sum(res$range_table)
## End(Not run)