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)