ci_standard {ciuupi}R Documentation

For given observed response vector yy, compute the standard 1α1 - \alpha confidence interval

Description

If σ\sigma is provided then compute the standard 1α1 - \alpha confidence interval for θ\theta. If σ\sigma is not provided then, as long as np30n-p \ge 30, replace σ\sigma by its estimate to compute an approximate 1α1 - \alpha confidence interval for θ\theta.

Usage

ci_standard(a, X, y, alpha, sig = NULL)

Arguments

a

The vector used to specify the parameter of interest θ=aβ\theta = a^{\top} \beta

X

The known n×pn \times p design matrix, with linearly independent columns

y

The nn-vector of observed responses

alpha

1α1 - \alpha is the coverage probability of the standard confidence interval

sig

Standard deviation of the random error. If a value is not specified then, provided that np30n-p \ge 30, sig is estimated from the data.

Details

Suppose that

y=Xβ+ε,y = X \beta + \varepsilon,

where yy is a random nn-vector of responses, XX is a known n×pn \times p matrix with linearly independent columns, β\beta is an unknown parameter pp-vector, and εN(0,σ2I)\varepsilon \sim N(0, \, \sigma^2 \, I), with σ2\sigma^2 assumed known. Suppose that the parameter of interest is θ=aβ\theta = a^{\top} \beta. The R function ci_standard computes the standard 1α1 - \alpha confidence interval for θ\theta.

The example below is described in Discussion 5.8 on p.3426 of Kabaila and Giri (2009). This example is obtained by extracting a 2×22 \times 2 factorial data set from the 232^3 factorial data set described in Table 7.5 of Box et al. (1963).

Value

If σ\sigma is provided then a data frame of the lower and upper endpoints of the standard 1α1 - \alpha confidence interval for θ\theta. If σ\sigma is not provided then, as long as np30n-p \ge 30, a data frame of the lower and upper endpoints of an approximation to this confidence interval.

References

Box, G.E.P., Connor, L.R., Cousins, W.R., Davies, O.L., Hinsworth, F.R., Sillitto, G.P. (1963) The Design and Analysis of Industrial Experiments, 2nd edition, reprinted. Oliver and Boyd, London.

Kabaila, P. and Giri, K. (2009) Confidence intervals in regression utilizing prior information. Journal of Statistical Planning and Inference, 139, 3419 - 3429.

Examples

y <- c(87.2, 88.4, 86.7, 89.2)
x1 <- c(-1, 1, -1, 1)
x2 <- c(-1, -1, 1, 1)
X <- cbind(rep(1, 4), x1, x2, x1*x2)
a <- c(0, 2, 0, -2)
ci_standard(a, X, y, 0.05, sig = 0.8)


[Package ciuupi version 1.2.3 Index]