Skip to contents

TOST Independent Samples T-Test for jamovi. This function is not meant to be utilized in R. See t_TOST function.

Usage

dataTOSTtwo(
  data,
  deps,
  group,
  var_equal = FALSE,
  hypothesis = "EQU",
  low_eqbound = -0.5,
  high_eqbound = 0.5,
  eqbound_type = "raw",
  alpha = 0.05,
  desc = FALSE,
  plots = FALSE,
  descplots = FALSE,
  low_eqbound_d = -999999999,
  high_eqbound_d = -999999999,
  smd_type = "g"
)

Arguments

data

the data as a data frame

deps

a vector of strings naming dependent variables in data

group

a string naming the grouping variable in data; must have two levels

var_equal

TRUE or FALSE (default), assume equal variances

hypothesis

'EQU' for equivalence (default), or 'MET' for minimal effects test, the alternative hypothesis.

low_eqbound

a number (default: -0.5) the lower equivalence/MET bounds

high_eqbound

a number (default: 0.5) the upper equivalence/MET bounds

eqbound_type

'SMD' (default) or 'raw'; whether the bounds are specified in Cohen's d or raw units respectively

alpha

alpha level (default = 0.05)

desc

TRUE or FALSE (default), provide descriptive statistics

plots

TRUE or FALSE (default), provide effect size plots

descplots

TRUE or FALSE (default), provide plots

low_eqbound_d

deprecated

high_eqbound_d

deprecated

smd_type

'd' (default) or 'g'; whether the calculated effect size is biased (d) or bias-corrected (g).

Value

A results object containing:

results$texta html
results$tosta table
results$eqba table
results$effsizea table
results$desca table
results$plotsan array of images
results$descplotsan array of images

Tables can be converted to data frames with asDF or as.data.frame. For example:

results$tost$asDF

as.data.frame(results$tost)

References

Berger, R. L., & Hsu, J. C. (1996). Bioequivalence Trials, Intersection-Union Tests and Equivalence Confidence Sets. Statistical Science, 11(4), 283-302.

Gruman, J. A., Cribbie, R. A., & Arpin-Cribbie, C. A. (2007). The effects of heteroscedasticity on tests of equivalence. Journal of Modern Applied Statistical Methods, 6(1), 133-140, formula for Welch's t-test on page 135

Examples

library(TOSTER)

## Load iris dataset, remove one of the three groups so two are left

data<-iris[which(iris$Species!="versicolor"),]

## TOST procedure on the raw data

dataTOSTtwo(data, deps="Sepal.Width", group="Species", var_equal = TRUE, low_eqbound = -0.5,
            high_eqbound = 0.5, alpha = 0.05, desc = TRUE, plots = TRUE)
#> 
#>  TOST INDEPENDENT SAMPLES T-TEST
#> 
#>  Hypothesis Tested: Equivalence
#> 
#>  Null Hypothesis: -0.5 &ge; (Mean1 - Mean2) or (Mean1 - Mean2) &ge; 0.5
#>  Alternative: -0.5 < (Mean1 - Mean2) < 0.5
#>  &#9989 NHST: reject null significance hypothesis that the effect is
#>  equal to zero
#>  &#10060 TOST: don't reject null equivalence hypothesis
#> 
#>  Note: SMD confidence intervals are an approximation. See our <a
#>  href="https://aaroncaldwell.us/TOSTERpkg/articles/SMD_calcs.html">documentation.
#> 
#>  TOST Results                                                
#>  ─────────────────────────────────────────────────────────── 
#>                                 t          df       p        
#>  ─────────────────────────────────────────────────────────── 
#>    Sepal.Width    t-test          6.450    98.00    < .001   
#>                   TOST Lower      13.55    98.00    < .001   
#>                   TOST Upper    -0.6536    98.00      .257   
#>  ─────────────────────────────────────────────────────────── 
#> 
#> 
#>  Equivalence Bounds                                 
#>  ────────────────────────────────────────────────── 
#>                                 Low        High     
#>  ────────────────────────────────────────────────── 
#>    Sepal.Width    Hedges's g     -1.421     1.421   
#>                   Raw           -0.5000    0.5000   
#>  ────────────────────────────────────────────────── 
#> 
#> 
#>  Effect Sizes                                                  
#>  ───────────────────────────────────────────────────────────── 
#>                                 Estimate    Lower     Upper    
#>  ───────────────────────────────────────────────────────────── 
#>    Sepal.Width    Hedges's g       1.280    0.9356     1.666   
#>                   Raw             0.4540    0.3371    0.5709   
#>  ───────────────────────────────────────────────────────────── 
#>    Note. Denominator set to the pooled SD
#> 
#> 
#>  Descriptives                                              
#>  ───────────────────────────────────────────────────────── 
#>               N     Mean     Median    SD        SE        
#>  ───────────────────────────────────────────────────────── 
#>    Group 1    50    3.428     3.400    0.3791    0.05361   
#>    Group 2    50    2.974     3.000    0.3225    0.04561   
#>  ───────────────────────────────────────────────────────── 
#>