dsquintana

# Write function

t.test2 <- function(m1,m2,s1,s2,n1,n2,m0=0,equal.variance=FALSE)
{
  if( equal.variance==FALSE ) 
  {
    se <- sqrt( (s1^2/n1) + (s2^2/n2) )
    # welch-satterthwaite df
    df <- ( (s1^2/n1 + s2^2/n2)^2 )/( (s1^2/n1)^2/(n1-1) + (s2^2/n2)^2/(n2-1) )
  } else

dsquintana

# Load TOSTER package

library(TOSTER)

# Perform equivalence test

TOSTtwo(m1=37.33,m2=39.65 ,sd1=13.40,sd2=15.91,n1=23,n2=32,
        low_eqbound_d=-0.26,high_eqbound_d=0.26,
        plot = TRUE)

dsquintana

# Load package

library(pwr)

# Perform power analysis

pwr.t2n.test(sig.level = 0.05, n1= 32,n2= 23, power = 0.8,
             alternative="two.sided")

dsquintana

### Contour-enhanced funnel plots using metafor ###

# Load metafor package

library("metafor")

# Load dataset

dat <- get(data(dat.molloy2014))