BagimsizOrneklem

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_bagimsiz <- data.frame(
  Grup = as.factor(rep(c("Grup_A", "Grup_B"), each = 25)),
  Skor = c(round(rnorm(25, 18, 4), 1), round(rnorm(25, 14, 4), 1))
)

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

## # A tibble: 2 × 4
##   Grup   variable statistic     p
##   <fct>  <chr>        <dbl> <dbl>
## 1 Grup_A Skor         0.967 0.574
## 2 Grup_B Skor         0.977 0.811

leveneTest(Skor ~ Grup, data = veri_bagimsiz)

## Levene's Test for Homogeneity of Variance (center = median)
##       Df F value Pr(>F)
## group  1   0.001 0.9753
##       48

t.test(Skor ~ Grup, data = veri_bagimsiz, var.equal = TRUE)

## 
##  Two Sample t-test
## 
## data:  Skor by Grup
## t = 3.2658, df = 48, p-value = 0.002018
## alternative hypothesis: true difference in means between group Grup_A and group Grup_B is not equal to 0
## 95 percent confidence interval:
##  1.325206 5.570794
## sample estimates:
## mean in group Grup_A mean in group Grup_B 
##               17.860               14.412

veri_bagimsiz %>% cohens_d(Skor ~ Grup, var.equal = TRUE)

## # A tibble: 1 × 7
##   .y.   group1 group2 effsize    n1    n2 magnitude
## * <chr> <chr>  <chr>    <dbl> <int> <int> <ord>    
## 1 Skor  Grup_A Grup_B   0.924    25    25 large

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