Travailler avec les objets de sortie
R pour les utilisateurs et utilisatrices de SAS
Melinda Higgins, PhD
Research Professor/Senior Biostatistician Emory University
Sortie de summary
# Enregistrer le résumé et afficher davissmry
davissmry <- daviskeep %>%
select(weight, height, bmi) %>%
summary()
# Afficher l'objet de sortie summary
davissmry
weight height bmi
Min. : 39.0 Min. :148.0 Min. :15.82
1er Qu.: 55.0 1er Qu.:164.0 1er Qu.:20.22
Médiane: 63.0 Médiane:170.0 Médiane:21.80
Moy. : 65.3 Moy. :170.6 Moy. :22.26
3e Qu. : 73.5 3e Qu. :177.5 3e Qu. :23.94
Max. :119.0 Max. :197.0 Max. :36.73
Sortie de summary
# Classe de davissmry
class(davissmry)
# Vérifier que davissmry est aussi une matrice
is.matrix(davissmry)
# Afficher les dimensions de davissmry
dim(davissmry)
[1] "table"
[1] TRUE
[1] 6 3
Sélectionner des éléments de l'objet summary
# Afficher les colonnes 1 à 2 de la sortie
davissmry[,1:2]
weight height
Min. : 39.0 Min. :148.0
1er Qu.: 55.0 1er Qu.:164.0
Médiane: 63.0 Médiane:170.0
Moy. : 65.3 Moy. :170.6
3e Qu. : 73.5 3e Qu. :177.5
Max. :119.0 Max. :197.0
# Afficher les lignes 4 à 6 de la sortie
davissmry[4:6,]
weight height bmi
Moy. : 65.3 Moy. :170.6 Moy. :22.26
3e Qu. : 73.5 3e Qu. :177.5 3e Qu. :23.94
Max. :119.0 Max. :197.0 Max. :36.73
Sortie de dplyr summarise
# Enregistrer la sortie de summarise
davissmall <- daviskeep %>%
select(weight, height) %>%
summarise(across(everything(),
list(mean = ~mean(.x),
sd = ~sd(.x))))
# Afficher l'objet de sortie davissmall
davissmall
weight_mean weight_sd height_mean height_sd
1 65.29648 13.34346 170.5879 8.948848
# Afficher la structure de l'objet davissmall
str(davissmall)
'data.frame': 1 obs. of 4 variables:
$ weight_mean: num 65.3
$ weight_sd : num 13.3
$ height_mean: num 171
$ height_sd : num 8.95
Afficher des éléments de la sortie summarise
# Afficher la colonne 3 pour la taille moyenne
davissmall[, 3]
[1] 170.5879
# Afficher l'élément nommé height_mean
davissmall$height_mean
[1] 170.5879
Sortie de dplyr group_by
# Utiliser group_by pour obtenir les stats de bmi par sexe
davisbmisex <- daviskeep %>%
group_by(sex) %>%
select(bmi, sex) %>%
summarise(across(everything(),
list(mean = ~mean(.x),
sd = ~sd(.x))))
# Afficher l'objet de sortie
davisbmisex
# A tibble: 2 × 3
sex bmi_mean bmi_sd
<fct> <dbl> <dbl>
1 F 21.0 2.18
2 M 23.9 3.12
# Obtenir la structure de davisbmisex
str(davisbmisex)
tibble [2 × 3] (S3: tbl_df/tbl/data.frame)
$ sex : Factor w/ 2 levels "F","M": 1 2
$ bmi_mean: num [1:2] 21 23.9
$ bmi_sd : num [1:2] 2.18 3.12
1 La variable sex a été automatiquement codée comme de type Factor.
Afficher des éléments de la sortie group_by
# Afficher la colonne 2 de davisbmisex
davisbmisex[, 2]
# A tibble: 2 × 1
bmi_mean
<dbl>
1 21.0
2 23.9
# Afficher la ligne 1 de davisbmisex
davisbmisex[1,]
# A tibble: 1 × 3
sex bmi_mean bmi_sd
<fct> <dbl> <dbl>
1 F 21.0 2.18
Afficher des éléments de la sortie group_by
# Afficher l'élément sd de davisbmisex
davisbmisex$bmi_sd
[1] 2.176094 3.120285
# Afficher la ligne pour les hommes
davisbmisex %>%
filter(sex == "M")
# A tibble: 1 × 3
sex bmi_mean bmi_sd
<fct> <dbl> <dbl>
1 M 23.9 3.12
Sortie de psych describe
# Enregistrer la sortie de psych::describe() comme davispsych
davispsych <- daviskeep %>%
select(weight, height, bmi) %>%
psych::describe()
# Afficher l'objet de sortie davispsych
davispsych
vars n mean sd median trimmed mad min max range skew kurtosis se
weight 1 199 65.30 13.34 63.0 64.12 11.86 39.00 119.00 80.00 0.91 0.84 0.95
height 2 199 170.59 8.95 170.0 170.40 10.38 148.00 197.00 49.00 0.21 -0.38 0.63
bmi 3 199 22.26 3.01 21.8 22.08 2.55 15.82 36.73 20.91 0.91 1.91 0.21
Classe et structure de la sortie psych describe
str(davispsych)
Classes 'psych', 'describe' et 'data.frame':
3 obs. de 13 variables :
$ vars : int 1 2 3
$ n : num 199 199 199
$ mean : num 65.3 170.6 22.3
$ sd : num 13.34 8.95 3.01
$ median : num 63 170 21.8
$ trimmed : num 64.1 170.4 22.1
$ mad : num 11.86 10.38 2.55
$ min : num 39 148 15.8
$ max : num 119 197 36.7
$ range : num 80 49 20.9
$ skew : num 0.905 0.21 0.913
$ kurtosis: num 0.843 -0.379 1.914
$ se : num 0.946 0.634 0.213
Afficher des éléments de la sortie psych describe
# Afficher les moyennes tronquées
davispsych$trimmed
[1] 64.12422 170.40373 22.08452
# Afficher la ligne 2 pour la taille
davispsych %>%
slice(2)
vars n mean sd median trimmed mad min max range skew kurtosis se
1 2 199 170.59 8.95 170 170.4 10.38 148 197 49 0.21 -0.38 0.63
Afficher des éléments de la sortie psych describe
# Afficher la taille d'échantillon, la médiane, le min et le max
davispsych %>%
select(n, median, min, max)
n median min max
weight 199 63.0 39.00 119.00
height 199 170.0 148.00 197.00
bmi 199 21.8 15.82 36.73
Personnalisons les résultats pour la sortie statistique sur les ormeaux !
R pour les utilisateurs et utilisatrices de SAS
