Computes the actual size of a intersection union test procedure that corresponds to a (1 – alpha) confidence set.

1 | ```
iutsize(p, n, alpha=0.1, sim=1e6)
``` |

`p` |
An integer giving the number of dimensions. |

`n` |
An integer giving the sample size. |

`alpha` |
A numeric value specifying the type I error level to be controlled. Default is 0.1. |

`sim` |
An integer giving the number of simulations to be carried out. Default is 1 million. |

A (1 – alpha) confidence set can be used to derive a two one-sided tests (TOST) procedure (Schuirmann 1987) whereby type I error rate control is ensured at level alpha due to the intersection union principle (Berger 1982). The actual test size, however, is often substantially lower than alpha i.e., the approach is conservative. It is well known for the one-dimensional case that the TOST corresponding to a (1 – alpha) confidence interval has size (1 – alpha/2). This function computes the achieved test size with dimension `p`

and `n`

according to the formula on p. 399 of Wang et al (1999).

A numeric value giving the actual size of the test.

Philip Pallmann (p.pallmann@lancaster.ac.uk)

Roger L. Berger (1982) Multiparameter hypothesis testing and acceptance sampling. Technometrics, 24(4), 295–300.

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.

Weizhen Wang, J. T. Gene Hwang, Anirban DasGupta (1999) Statistical tests for multivariate bioequivalence. Biometrika, 86(2), 395–402.

1 2 3 4 5 6 7 8 9 10 11 12 13 | ```
# For p=1 we get the well-known result that the 90% CI corresponds to the TOST at 5%:
#iutsize(p=1, n=20)
# With increasing dimension the test gets conservative:
#iutsize(p=2, n=20)
#iutsize(p=3, n=20)
# For p>1 the conservativeness also depends on sample size:
#iutsize(p=2, n=10)
#iutsize(p=2, n=1000)
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.