set.seed(57)
grp1=rnorm(50,0,1.5)
grp2=rnorm(50,0.8,3.2)
p = as.vector(BayesFactor::ttestBF(x=grp1,y=grp2,
posterior = TRUE, iterations = 3000,
rscale = "medium")[,4])
# flat reference function
res = fbst(posteriorDensityDraws = p, nullHypothesisValue = 0,
dimensionTheta = 3, dimensionNullset = 2)
# Return the e-value from an fbst object
res$eValue
Run the code above in your browser using DataLab