Calculates the confindence intervals of the hydraulic parameters, drawdown and its derivative by bootstrapping the residuals
# S3 method for pumping_test confint(object, parm, level = 0.95, ..., slug = FALSE, method = "wald", d = 20, neval = 100, seed = 12345)
object | An estimated pumping_test object |
---|---|
parm | Name of parameters to be analyzed |
level | A numeric vector with the probabilities used to calculate the confidence intervals. |
... | additional argument(s) for methods. |
slug | Logical flag to indicate a slug test. |
method | A character string with the method to estimate the confidence intervals. The methods currently supported are:
. |
d | Number of points used in the calculation of the derivative. |
neval | Number of bootstrap realizations used in the estimation of confidence intervals. |
seed | Random seed. |
In the case of method=wald the function returns a list with the following entries:
hydraulic.parameters.names: A character vector with the names of the hydraulic parameters
hydraulic.parameters.ci: Confidence intervals of the hydraulic parameters
In the case of method=bootstrap the function returns a list with the following entries:
fit.parameters: A matrix with the values of the parameters of the nonlinear regression model
total.drawdown: A matrix with the values of the drawdown
total_dsdlogt: A matrix with the values of the derivative of drawdown
sci: A matrix with the confidence intervals of drawdown
dsci: A matrix with the confidence intervals of the derivative of drawdown
hydraulic.parameters: A matrix with the hydraulic parameter estimated from the boostrapp realizations
parameter.names: A vector with the name of the fitting parameters
hydraulic.parameters.names: A character vector with the names of the hydraulic parameters
hydraulic.parameters.ci: Confidence intervals of the hydraulic parameters
Other base functions: additional.parameters<-
,
confint_bootstrap
,
confint_jackniffe
,
confint_wald
, estimated<-
,
evaluate
, fit.optimization
,
fit.parameters<-
,
fit.sampling
, fit
,
hydraulic.parameter.names<-
,
hydraulic.parameters<-
,
model.parameters
, model<-
,
plot.pumping_test
,
plot_model_diagnostic
,
plot_sample_influence
,
plot_uncert
,
print.pumping_test
,
pumping_test
, simulate
,
summary.pumping_test
data(theis) ptest.theis <- pumping_test("Well1", Q = 1.388e-2, r = 250, t= theis$t, s = theis$s) ptest.theis.fit <- fit(ptest.theis, "theis") hydraulic.parameters(ptest.theis) <- ptest.theis.fit$hydraulic_parameters fit.parameters(ptest.theis) <- ptest.theis.fit$parameters model(ptest.theis) <- "theis" estimated(ptest.theis) <- TRUE ptest.theis.ci1 <- confint(ptest.theis, level = 0.975) ptest.theis.ci2 <- confint(ptest.theis, level = c(0.025, 0.975), method = 'bootstrap', d = 20, neval = 100)