confints: Confidence intervals around multivariate normal means

Description Usage Arguments Details Value Note Author(s) References See Also Examples

Description

Computes boundaries of (simultaneous) confidence intervals around multivariate normal means using different methods.

Usage

1
confints(dat, method, alpha=0.1, steps=100)

Arguments

dat

A matrix or data.frame with independent units in rows and multivariate outcomes in columns.

method

A character string specifying the method to be used. See details for available methods.

alpha

A numeric value giving the type I error level to be controlled. Default is 0.1.

steps

An integer setting the initial number of steps for the search algorithm.

Details

x

Available methods are: tost for xxxxxxxxxxxxxx; expanded for xxxxxxxxxxxx; fixseq for xxxxxxxxxxxxxxxxxxx.

Value

x

Note

x

Author(s)

Philip Pallmann ([email protected])

References

Eve Bofinger (1992) Expanded confidence intervals, one-sided tests, and equivalence testing. Journal of Biopharmaceutical Statistics, 2(2), 181–188.

Jason C. Hsu, J. T. Gene Hwang, Hung-Kung Liu, Stephen J. Ruberg (1994) Confidence intervals associated with tests for bioequivalence. Biometrika, 81(1), 103–114.

Willi Maurer, Ludwig A. Hothorn, Walter Lehmacher (1995) Multiple comparisons in drug clinical trials and preclinical assays: a priori ordered hypotheses. In: Joachim Vollmar (editor), Biometrie in der Chemisch-Pharmazeutischen Industrie, vol. 6, pp. 3–18. Fischer-Verlag, Stuttgart, Germany.

Philip Pallmann & Thomas Jaki (2016) Simultaneous confidence regions and intervals for multivariate bioequivalence. Submitted to Statistics in Medicine.

Donald J. Schuirmann (1987) A comparison of the two one-sided tests procedure and the power approach for assessing the equivalence of average bioavailability. Journal of Pharmacokinetics and Biopharmaceutics, 15(6), 657–680.

See Also

x

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
# Generate trivariate normal data

mydata <- mvtnorm::rmvnorm(n=50, mean=rep(0.1, 3), sigma=diag(3) * 0.05)
colnames(mydata) <- LETTERS[1:3]

# Compute 90% confidence intervals

confints(dat=mydata, method="tost", alpha=0.1)
confints(dat=mydata, method="expanded", alpha=0.1)
#confints(dat=mydata, method="fixseq", alpha=0.1)

PhilipPallmann/simbe documentation built on May 8, 2019, 1:34 a.m.