equiv.boot {equivalence}R Documentation

Regression-based TOST using bootstrap

Description

This function wraps the regression-based TOST equivalence test inside the bootstrap, extracts and reports the useful quantities, and reports the outcome of the test.

Usage

equiv.boot(x, y, alpha = 0.05, b0.ii = 0.25, b1.ii = 0.25, reps = 100)

Arguments

x the predictor variable
y the response variable
alpha the size of the test
b0.ii the half-length of the region of similarity for the intercept.
b1.ii the half-length of the region of similarity for the slope.
reps the number of bootstrap replicates required

Details

In each case, if the two one-sided confidence interval is inside the region of similarity then the null hypothesis of dissimilarity is rejected.

Value

A list of length 10, comprising

n The effective (non-missing) sample size
ci.b0 The intercept TOST confidence interval
rs.b0 The intercept region of similarity
q.b0 The proportions below, within, and above the intercept region of similarity
Test.b0 The outcome of the test of the null hypothesis of dissimilarity for the intercept (Reject/Not Reject)
ci.b1 The slope TOST confidence interval
rs.b1 The slope region of similarity
q.b1 The proportions below, within, and above the slope region of similarity
Test.b1 The outcome of the test of the null hypothesis of dissimilarity for the slope (Reject/Not Reject)
eff.alpha

{The corrected alpha for each of the two independent tests.}

Acknowledgements

Feedback from Mohammad Al-Ahmadi has been very useful for this function.

Author(s)

Andrew Robinson A.Robinson@ms.unimelb.edu.au

References

Robinson, A.P., R.A. Duursma, and J.D. Marshall. 2005. A regression-based equivalence test for model validation: shifting the burden of proof. Tree Physiology 25, 903-913

See Also

lm, boot, equiv.boot.lm

Examples


# Approximately reproduces the first row from Table 2 of Robinson et al. (2005)

data(pref.4PG)
equiv.boot(pref.4PG$volinc4PG, pref.4PG$stemvolinc)


[Package equivalence version 0.3.4 Index]