R code for equivalence test in one-way ANOVA
[TYPE INTO R THE DESIRED INPUTS RSQ, N, K, TDEL AND BETA USING VALUES IN FORM BELOW].
RSQ represents partial $$eta^text{2}$$ using one of Cohen's [:FAQ/effectSize:rules of thumb], n is a vector of group sizes, tdel represents the [:FAQ/mageq:criterion for t] and beta is the type II error.
rsq <- 0.006 n <- c(10,12,13,15) tdel <- 0.5 beta <- 0.05
[THEN COPY AND PASTE THE BELOW INTO R]
k <- length(n) ns <- sum(n) psi2 <- (rsq/(1-rsq))*k*((ns-k)/ns) cstats <- (k*(k-1)/ns)*qf(beta,k-1,ns-k,(ns/k)*tdel*tdel) ind <- 0 if (psi2 < cstats) ind = 1 print(ind)