equiv.boot {equivalence} | R Documentation |
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.
equiv.boot(x, y, alpha = 0.05, b0.ii = 0.25, b1.ii = 0.25, reps = 100)
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 |
In each case, if the two one-sided confidence interval is inside the region of similarity then the null hypothesis of dissimilarity is rejected.
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.}
Feedback from Mohammad Al-Ahmadi has been very useful for this function.
Andrew Robinson A.Robinson@ms.unimelb.edu.au
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
# Approximately reproduces the first row from Table 2 of Robinson et al. (2005) data(pref.4PG) equiv.boot(pref.4PG$volinc4PG, pref.4PG$stemvolinc)