krCI: Calculating confidence intervals for WTP using a parametric...

Description Usage Arguments Details Value References See Also Examples

View source: R/krCI.R

Description

This function calculates confidence intervals for WTP using the method of Krinsky and Robb (1986, 1990).

Usage

1
2
3
4
krCI(obj, nsim = 1000, CI = 0.95, individual = NULL)

## S3 method for class 'krCI'
print(x, ...)

Arguments

obj

an S3 class object "dbchoice" or "sbchoice".

nsim

the number of draws of the parameters.

CI

a percentile of the confidence intervals to be estimated.

individual

a data frame containing covariates that show an individual of which to estimate WTP.

x

an object of class "krCI".

...

optional arguments. Currently not in use.

Details

In the method of Krinsky and Robb (1986, 1990), a set of parameters is drawn nsim times from a multivariate normal distribution with a vector of the parameter estimates as a mean and the variance-covariance matrix of the parameter estimates. Then, various WTPs are computed for each draw of simulated parameters. As a result, we are able to build an empirical distribution of the WTPs concerned, and hence the confidence intervals. For each WTP, and when nsim = 1000, the lower and the upper bound of the 95% confidence interval (CI = 0.95) correspond to the 26th and the 975th sorted estimates, respectively.

Confidence intervals based on the bootstrap method are calculated by bootCI.

Hole (2007) conducted simulation experiments to compare the performance of the method of Krinsky and Robb (1986, 1990) with the bootstrap one.

A WTP of a specific individual (e.g., a representative respondent) can be estimated when assigning covariates to individual. See Example for details.

Value

The function krCI() returns an object of S3 class "krCI". An object of "krCI" is a list with the following components.

out

the output table with simulated confidence intervals as well as the four type of WTP estimates (mean, truncated mean, truncated mean with adjustment and median) from the ML estimation.

mWTP

a vector of simulated mean WTP. When |beta| < 1, this item is set to -999.

tr.mWTP

a vector of simulated mean WTP truncated at the maximum bid.

adj.tr.mWTP

a vector of simulated mean WTP truncated at the maximum bid with the adjustment.

medWTP

a vector of simulated median WTP.

When the parameter estimate on the bid does not satisfy the condition for the existence of the finite mean WTP (|beta|>1), the values of the lower and the upper bound of the confidence interval are coerced to set to -999.

The table contains the confidence intervals for the four types (mean, truncated mean, truncated mean with adjustment and median) of WTP estimate from the ML estimation. The adjustment for the truncated mean WTP is implemented by the method of Boyle et~al.(1988).

The generic function print() is available for the object of class "krCI" and displays the table of simulated confidence intervals.

References

Boyle KJ, Welsh MP, Bishop RC (1988). “Validation of Empirical Measures of Welfare Change: Comment.” Land Economics, 64(1), 94–98.

Hole AR (2007). “A Comparison of Approaches to Estimating Confidence Intervals for Willingness to Pay Measure.” Health Economics, 16, 827–840.

Krinsky I, Robb AL (1986). “On Approximating the Statistical Properties of Elasticities.” The Review of Economics and Statistics, 68, 715–719.

Krinsky I, Robb AL (1990). “On Approximating the Statistical Properties of Elasticities: A Correction.” The Review of Economics and Statistics, 72, 189–190.

See Also

bootCI, dbchoice, sbchoice

Examples

1
## See Examples in dbchoice and sbchoice.

Example output

Loading required package: MASS
Loading required package: interval
Loading required package: survival
Loading required package: perm
Loading required package: Icens
Loading required package: MLEcens
Loading required package: Formula

DCchoice documentation built on May 2, 2019, 4:44 p.m.