A function for TOST with all types of t-tests from summary statistics.
tsum_TOST( m1, sd1, n1, m2 = NULL, sd2 = NULL, n2 = NULL, r12 = NULL, hypothesis = c("EQU", "MET"), paired = FALSE, var.equal = FALSE, eqb, low_eqbound, high_eqbound, mu = 0, eqbound_type = c("raw", "SMD"), alpha = 0.05, bias_correction = TRUE, rm_correction = FALSE, glass = NULL, smd_ci = c("nct", "goulet", "t", "z") )
mean of group 1
standard deviation of group 1
sample size in group 1
mean of group 2
standard deviation of group 2
sample size in group 2
correlation of dependent variable between group 1 and group 2
'EQU' for equivalence (default), or 'MET' for minimal effects test, the alternative hypothesis.
a logical indicating whether you want a paired t-test.
a logical variable indicating whether to treat the two variances as being equal. If TRUE then the pooled variance is used to estimate the variance otherwise the Welch (or Satterthwaite) approximation to the degrees of freedom is used.
Equivalence bound. Can provide 1 value (negative value is taken as the lower bound) or 2 specific values that represent the upper and lower equivalence bounds.
lower equivalence bounds (deprecated).
upper equivalence bounds (deprecated).
a number indicating the true value of the mean for the two tailed test (or difference in means if you are performing a two sample test).
Type of equivalence bound. Can be set to "SMD" for standardized mean difference (i.e., Cohen's d) or "raw" for the mean difference. Default is "raw". Raw is strongly recommended as SMD bounds will produce biased results.
alpha level (default = 0.05)
Apply Hedges' correction for bias (default is TRUE).
Repeated measures correction to make standardized mean difference Cohen's d(rm). This only applies to repeated/paired samples. Default is FALSE.
A option to calculate Glass's delta as an alternative to Cohen's d type SMD. Default is NULL to not calculate Glass's delta, "glass1" will use the first group's SD as the denominator whereas "glass2" will use the 2nd group's SD.
Method for calculating SMD confidence intervals. Methods include Goulet, noncentral t (nct), central t (t), and normal method (z).
An S3 object of class
"TOSTt" is returned containing the following slots:
A table of class
"data.frame" containing two-tailed t-test and both one-tailed results.
A table of class
"data.frame" containing equivalence bound settings.
table of class
"data.frame" containing effect size estimates
String stating the hypothesis being tested
List containing the results of the standardized mean difference calculations (e.g., Cohen's d).Items include: d (estimate), dlow (lower CI bound), dhigh (upper CI bound), d_df (degrees of freedom for SMD), d_sigma (SE), d_lambda (non-centrality), J (bias correction), smd_label (type of SMD), d_denom (denominator calculation)
Alpha level set for the analysis.
Type of t-test.
List included text regarding the decisions for statistical inference.
For details on the calculations in this function see vignette("IntroTOSTt") & vignette("SMD_calcs").
# example code # One sample test tsum_TOST(m1 = 0.55, n1 = 18, sd1 = 4, eqb = 2) #> #> One-sample t-Test #> #> The equivalence test was non-significant, t(17) = -1.538, p = 7.12e-02 #> The null hypothesis test was non-significant, t(17) = 0.583, p = 5.67e-01 #> NHST: don't reject null significance hypothesis that the effect is equal to zero #> TOST: don't reject null equivalence hypothesis #> #> TOST Results #> t df p.value #> t-test 0.5834 17 0.567 #> TOST Lower 2.7047 17 0.008 #> TOST Upper -1.5380 17 0.071 #> #> Effect Sizes #> Estimate SE C.I. Conf. Level #> Raw 0.5500 0.9428 [-1.0901, 2.1901] 0.9 #> Hedges's g 0.1313 0.2367 [-0.2427, 0.5015] 0.9 #> Note: SMD confidence intervals are an approximation. See vignette("SMD_calcs").