Calcul boxplot avec poids
Bonjour,
La fonction wtd.boxplot ne donne pas les mêmes résultats que la fonction wtd.quantile.
Il semblerait que pour tracer le boxplot, la fonction wtd multiplie par 100 la colonne poids.
J’ai essayé avec la fonction syvboxplot. Sauf que là j’ai un tout autre résultat.
Je ne comprends pas, quelqu’un pourrait m’éclairer à ce sujet ?
Aussi, si je décide de conserver ces boxplot, sont-elles justes ?
Une dernière question, peut-on construire une boxplot manuellement (avec des lignes de codes tout de même toujours sur R).
Code:
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
| > bd
vitesse poids poids2
1 1 0.12 12
2 2 0.36 36
3 3 0.25 25
4 4 0.32 32
5 5 0.85 85
> wtd.quantile(bd$vitesse, weights = bd$poids)
0% 25% 50% 75% 100%
4.000 4.225 4.450 4.675 4.900
> wtd.quantile(bd$vitesse, weights = bd$poids2)
0% 25% 50% 75% 100%
1.00 2.25 4.00 5.00 5.00
> wtd.boxplot(bd$vitesse, weights = bd$poids)
> wtd.boxplot(bd$vitesse, weights = bd$poids2)
> dw <- svydesign(ids = ~1, data = bd, weights = ~bd$poids)
> svyboxplot(vitesse ~ 1, dw, main = "Boxplot simple", ylab = "Âge")
> svyquantile(~vitesse, dw, quantile = c(0,0.25, 0.5, 0.75,1), ci = TRUE)
$quantiles
0 0.25 0.5 0.75 1
vitesse 1 1.986111 3.6875 4.441176 5
$CIs
, , vitesse
0 0.25 0.5 0.75 1
(lower 1.000000 1.199615 1 1.389876 5
upper) 1.453163 3.087933 5 5.000000 5 |
