Le test T peut être défini comme un test statistique utilisé pour comparer deux moyennes. Ce guide fournit de multiples tutoriels décrivant les différents types de tests T, notamment:
- tests-t pour échantillon unique. Comparer la moyenne d’un échantillon à une moyenne standard connue.
- t-test pour échantillons indépendants : Test t de Student et test t de Welch. Comparer deux groupes indépendants
- test-t pour échantillons appariés. Comparer deux échantillons reliés entre eux.
Vous apprendrez la formule du test t et comment:
- Calculez les différents t-tests dans R
- Vérifier les hypothèses du test t
- Calculez et rapportez la taille de l’effet du test t en utilisant le d de Cohen.
Sommaire:
Livre Apparenté
Pratique des Statistiques dans R II - Comparaison de Groupes: Variables NumériquesPrérequis
Assurez-vous d’avoir installé les paquets R suivants:
tidyverse
pour la manipulation et la visualisation des donnéesggpubr
pour créer facilement des graphiques prêts à la publicationrstatix
contient des fonctions R facilitant les analyses statistiques.datarium
: contient les jeux de données requis pour ce chapitre.
Commencez par charger les packages requis suivants:
library(tidyverse)
library(ggpubr)
library(rstatix)
Exemples de codes R
Comparaison de deux groupes indépendants:
# Préparation des données
data("genderweight", package = "datarium")
head(genderweight, 3)
## # A tibble: 3 x 3
## id group weight
## <fct> <fct> <dbl>
## 1 1 F 61.6
## 2 2 F 64.6
## 3 3 F 66.2
# Test statistique
stat.test <- genderweight %>%
t_test(weight ~ group) %>%
add_significance()
stat.test
## # A tibble: 1 x 9
## .y. group1 group2 n1 n2 statistic df p p.signif
## <chr> <chr> <chr> <int> <int> <dbl> <dbl> <dbl> <chr>
## 1 weight F M 20 20 -20.8 26.9 4.30e-18 ****
# Visualisation :
# Créer un box-plot
bxp <- ggboxplot(
genderweight, x = "group", y = "weight",
ylab = "Weight", xlab = "Groups", add = "jitter"
)
# Ajouter la p-value et les niveaux de significativité
stat.test <- stat.test %>% add_xy_position(x = "group")
bxp +
stat_pvalue_manual(stat.test, tip.length = 0) +
labs(subtitle = get_test_label(stat.test, detailed = TRUE))
Version: English
No Comments