Description Usage Arguments Details Value Note Author(s) References See Also Examples
Computes boundaries of (simultaneous) confidence intervals around multivariate normal means using different methods.
1 | confints(dat, method, alpha=0.1, steps=100)
|
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 |
alpha |
A numeric value giving the type I error level to be controlled. Default is |
steps |
An integer setting the initial number of steps for the search algorithm. |
x
Available method
s are: tost
for xxxxxxxxxxxxxx; expanded
for xxxxxxxxxxxx; fixseq
for xxxxxxxxxxxxxxxxxxx.
x
x
Philip Pallmann (p.pallmann@lancaster.ac.uk)
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.
x
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)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.