confints: Confidence intervals around multivariate normal means

In PhilipPallmann/simbe: Simultaneous Assessment of Bioequivalence

Description

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

Usage

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 ([e-mail])

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

# 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.