plotCIsigmaP {soilphysics}

Percentile Confidence Intervals for Simulated Preconsolidation Stress

Description

Build and plot percentile confidence intervals for preconsolidation stress simulated from simSigmaP.

Usage

plotCIsigmaP(msim, conf.level = 0.95, shade.col = "orange", 
	ordered = TRUE, xlim = NULL, xlab = expression(sigma[P]), 
	las = 1, mar = c(4.5, 6.5, 2, 1), ...)

Arguments

msim

an object of class "simSigmaP".

conf.level

the confidence level for the intervals.

shade.col

a character or number indicating the color of the shaded area delimiting each CI. See colors.

ordered

logical; should the intervals be displayed according to the value of the simulated mean?

xlim

optional; a numeric vector of length two containing the limits of the x-axis.

xlab

optional; a character indicating the x-axis label.

las

optional; see par.

mar

optional; see par.

...

further graphical parameters; see par.

Value

A numeric matrix containing the simulated mean, coefficient of variation, lower and upper CI limits and the name of the method used to calculate the preconsolidation stress.

Author(s)

Anderson Rodrigo da Silva <[e-mail]>

See Also

simSigmaP, sigmaP

Examples

# input data: stress and void ratio
pres <- c(1, 12.5, 25, 50, 100, 200, 400, 800, 1600)
VR <- c(1.43, 1.41, 1.40, 1.39, 1.35, 1.31, 1.25, 1.18, 1.12)

# simulation (may take a few seconds)
simres <- simSigmaP(VR, pres, nsim = 30)
head(simres)

# percentile confidence intervals
ci <- plotCIsigmaP(simres, conf.level = 0.95, 
	shade.col = "blue", ordered = TRUE)
print(ci)

# End (Not run)