TekYonluAnova

library(tidyverse)

## Warning: package 'tidyverse' was built under R version 4.5.2

## Warning: package 'ggplot2' was built under R version 4.5.2

## Warning: package 'tidyr' was built under R version 4.5.2

## Warning: package 'readr' was built under R version 4.5.2

## Warning: package 'purrr' was built under R version 4.5.2

## Warning: package 'stringr' was built under R version 4.5.2

## Warning: package 'forcats' was built under R version 4.5.2

## Warning: package 'lubridate' was built under R version 4.5.2

## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.1.4     ✔ readr     2.1.5
## ✔ forcats   1.0.1     ✔ stringr   1.6.0
## ✔ ggplot2   4.0.0     ✔ tibble    3.3.0
## ✔ lubridate 1.9.4     ✔ tidyr     1.3.1
## ✔ purrr     1.2.0     
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors

library(car)

## Warning: package 'car' was built under R version 4.5.3

## Zorunlu paket yükleniyor: carData

## Warning: package 'carData' was built under R version 4.5.3

## 
## Attaching package: 'car'
## 
## The following object is masked from 'package:dplyr':
## 
##     recode
## 
## The following object is masked from 'package:purrr':
## 
##     some

library(rstatix)

## Warning: package 'rstatix' was built under R version 4.5.3

## 
## Attaching package: 'rstatix'
## 
## The following object is masked from 'package:stats':
## 
##     filter

set.seed(123)
veri_anova <- data.frame(
  Grup = as.factor(rep(c("Grup_A", "Grup_B", "Grup_C"), each = 20)),
  Skor = c(round(rnorm(20, 15, 3), 1), round(rnorm(20, 18, 3), 1), round(rnorm(20, 12, 3), 1))
)

veri_anova %>% group_by(Grup) %>% shapiro_test(Skor)

## # A tibble: 3 × 4
##   Grup   variable statistic     p
##   <fct>  <chr>        <dbl> <dbl>
## 1 Grup_A Skor         0.970 0.752
## 2 Grup_B Skor         0.955 0.444
## 3 Grup_C Skor         0.974 0.833

leveneTest(Skor ~ Grup, data = veri_anova)

## Levene's Test for Homogeneity of Variance (center = median)
##       Df F value Pr(>F)
## group  2  0.0342 0.9664
##       57

ModelAnova <- aov(Skor ~ Grup, data = veri_anova)
summary(ModelAnova)

##             Df Sum Sq Mean Sq F value   Pr(>F)    
## Grup         2  307.9  153.96   20.06 2.53e-07 ***
## Residuals   57  437.4    7.67                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

TukeyHSD(ModelAnova)

##   Tukey multiple comparisons of means
##     95% family-wise confidence level
## 
## Fit: aov(formula = Skor ~ Grup, data = veri_anova)
## 
## $Grup
##                 diff       lwr       upr     p adj
## Grup_B-Grup_A  2.425  0.317005  4.532995 0.0204278
## Grup_C-Grup_A -3.110 -5.217995 -1.002005 0.0022175
## Grup_C-Grup_B -5.535 -7.642995 -3.427005 0.0000001

ggplot(veri_anova, aes(x = Grup, y = Skor, fill = Grup)) +
  geom_boxplot() +
  theme_minimal()