#One-way ANOVA

jmv::anovaOneW(
  formula = anxiety ~ groupcat,
  data = spider,
  welchs = FALSE,
  fishers = TRUE,
  desc = TRUE,
  descPlot = TRUE,
  eqv = TRUE)

## 
##  ONE-WAY ANOVA
## 
##  One-Way ANOVA (Fisher's)                           
##  -------------------------------------------------- 
##               F           df1    df2    p           
##  -------------------------------------------------- 
##    anxiety    2.826923      1     22    0.1068392   
##  -------------------------------------------------- 
## 
## 
##  Group Descriptives                                                 
##  ------------------------------------------------------------------ 
##               groupcat    N     Mean        SD           SE         
##  ------------------------------------------------------------------ 
##    anxiety    spider      12    47.00000    11.028888    3.183766   
##               toupee      12    40.00000     9.293204    2.682717   
##  ------------------------------------------------------------------ 
## 
## 
##  ASSUMPTION CHECKS
## 
##  Homogeneity of Variances Test (Levene's)            
##  --------------------------------------------------- 
##               F            df1    df2    p           
##  --------------------------------------------------- 
##    anxiety    0.7815619      1     22    0.3862236   
##  ---------------------------------------------------

library(ggplot2)
library(Hmisc)

## Loading required package: lattice

## Loading required package: survival

## Loading required package: Formula

## 
## Attaching package: 'Hmisc'

## The following objects are masked from 'package:base':
## 
##     format.pval, units

ggplot(spider, aes(y=anxiety, x=groupcat))+ 
  geom_jitter(position=position_jitter(w=0.15, h=1)) +
  stat_summary(fun="mean", geom="point", size=2.5, color="red") +
  stat_summary(fun.data = mean_cl_normal, geom = "errorbar", 
               fun.args = list(mult = 1.96), color="red", width=.4) +
  labs(y = "Anxiety", x = "Phobia group") +
  theme_bw()