confint.pumping_test: confint.pumping_test

Description Usage Arguments Value Author(s) See Also Examples

View source: R/pumpingtest.R

Description

Calculates the confindence intervals of the hydraulic parameters, drawdown and its derivative by bootstrapping the residuals

Usage

1
2
3
4
## S3 method for class 'pumping_test'
confint(object, parm, level = 0.95, ...,
  slug = FALSE, method = "wald", d = 20, neval = 100,
  seed = 12345)

Arguments

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:

  • wald

  • bootstrap

  • jackniffe

.

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.

Value

In the case of method=wald the function returns a list with the following entries:

In the case of method=bootstrap the function returns a list with the following entries:

Author(s)

Oscar Garcia-Cabrejo khaors@gmail.com

See Also

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

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
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)

khaors/pumpingtest documentation built on Nov. 15, 2019, 8:10 p.m.