Solution
1. Caractéristiques empiriques de l’échantillon
1.1. Moyenne
∑
xi : module pluviométrique ; n : nombre total d’observations
292.0
1.2. Variance et écart-type :
∑ (xi − x )
2
Variance s 2
= = 810.8
n
Ecart-type : s = s 2 = 90.0 mm
1.3. Coefficient de variation
s
CV = =0.31
x
2. Ajustement par la loi de Gauss
2.1. Calcul de la fréquence expérimentale de chaque valeur de l’échantillon : les
fréquences expérimentales au non dépassement sont calculées avec la formule
de Hazen qui s’écrit comme suit :
r − 0.5
F ( x) =
N
avec F(x) la fréquence au non dépassement ; r : le rang et N le nombre total d’observations
Numéro Module Fréquence Numéro Module Fréquence
(mm) expérimentale (mm) expérimentale
1 149.5 0.0082 32 288.5 0.5164
2 150.7 0.0246 33 291.0 0.5328
3 151.0 0.0410 34 294.1 0.5492
4 161.8 0.0574 35 294.3 0.5656
5 172.1 0.0738 36 308.0 0.5820
6 185.0 0.0902 37 308.6 0.5984
7 185.9 0.1066 38 310.1 0.6148
8 192.8 0.1200 39 319.1 0.6311
9 197.8 0.1393 40 321.5 0.6475
10 205.0 0.1557 41 326.9 0.6639
11 205.8 0.1721 42 329.6 0.6803
12 206.0 0.1885 43 333.5 0.6967
13 208.5 0.2049 44 342.0 0.7131
14 210.0 0.2213 45 357.2 0.7295
15 214.6 0.2377 46 360.2 0.7459
16 222.1 0.2541 47 364.5 0.7623
17 222.6 0.2705 48 366.5 0.7787
18 225.6 0.2869 49 367.4 0.7951
19 235.7 0.3033 50 380.9 0.8115
20 235.7 0.3197 51 382.0 0.8279
21 242.8 0.3361 52 394.4 0.8443
22 246.4 0.3525 53 399.5 0.8607
23 252.6 0.3689 54 402.1 0.8770
24 255.0 0.3852 55 408.8 0.8934
25 257.0 0.4016 56 413.4 0.9098
26 258.5 0.4180 57 421.2 0.9262
27 260.3 0.4344 58 421.8 0.9426
28 273.0 0.4508 59 428.7 0.9590
29 275.9 0.4672 60 500.8 0.9754
30 278.6 0.4836 61 556.4 0.9918
31 281.8 0.5000
2.2. Pointage sur papier Gauss des points expérimentaux : voir graphique.
2.3. Tracé de la droite de Henry : voir graphique.
3. Test d’adéquation de χ²
χ =∑
2
k (n
i − n pi )
2
1 n pi
N
n pi =
nk
k
∑n pi
P(x ) = 1
N
xk = x + s × u k
ni = effectif observé de la classe i ; npi = effectif théorique de cette même classe ; N = effectif total ; nk = nombre
de classes ; k = numéro de la classe. ; P(x) = fréquence au dépassement ; xk = limite supérieure de la classe
Classe Limites des ni npi ni-npi (ni-npi)² (ni-
classes npi)²/npi
1 ≤188.5 7 7.6250 -0.6250 0.3906 0.0512
2 188.5 -231.1 11 7.6250 3.3750 11.3906 1.4938
3 231.1 – 263.3 9 7.6250 1.3750 1.8906 0.2479
4 263.3 – 292.0 6 7.6250 -1.6250 2.6406 0.3463
5 292.0 – 320.7 6 7.6250 -1.6250 2.6406 0.3463
6 320.7 – 352.7 5 7.6250 -2.6250 6.8906 0.9037
7 352.7 – 395.6 8 7.6250 0.3750 0.1406 0.0184
8 ≥ 395.6 9 7.6250 1.3750 1.8906 0.2479
ν = k − p −1
χ cal
2
= 3.66
k = nombre de classes
χ 02.05 = 11.07 = χ 2 critique
p = nonbre de paramètres utilisés par la loi
Comme χ cal
2
< χ 02.05 on accepte l' hypothèse
4. Intervalles de confiance à 95% des paramètres empiriques
x = 292 mm
s = 90 mm
4.1. Taille de 61 ans (1906 – 1975) n> 30 pour la moyenne et n > 50 pour
l’écart-type
s
IC pour x = x ± u1−α
2 n
s
IC pour s = s ± u1−α
2 2n
α = 95% ⇒ u 1−α = 1.96 = u 0.025
2
90
x = 292.0 ± 1.96 = 292 ± 22.6
61
22.6
La moyenne est connue à ×100 ≈ 8% près
292
90
s = 90 ± 1.96 = 90 ± 16
61
16
L' écart - type est connu à × 100 ≈ 18% près
90
x = 293.5 mm
4.2. Taille de 23 ans (1906 – 1929)
s = 70.5 mm
4.2.1. Moyenne : n < 30
s
IC pour x = x ± t 1−α .
2 n
t : variable de Student
t à 95% = 2.074 pour ν = 22
70.5
IC pour x = 293.5 ± 2.074 = 293.5 ± 30.5
23
30.5
La moyenne est connue à ≈ 10% près
293.5
4.2.2. Ecart-type : n < 50
ns 2 ns 2
< s2 <
χ12 χ 22
[ ] [
p χ 2 < χ12 = p χ 2 < χ 22 = 1 − ] 2
α
α = 95% on a χ = χ 2
1
2
0.025 et χ = χ 0.975
2
2
2
r = N − 1 = 23 − 1 = 22
χ 0.025
2
= 36.78 et χ 0.975
2
= 10.98
23 × (70.5)² 23 × (70.5)²
< s² <
36.78 10.98
55.7 < s < 102.0
5. Estimation des quantiles
xp = x + s ×up
Module décennal humide
x0.90 = 292 + 90 ×1.28 = 407.2 mm
Module centenaire humide
x0.99 = 292 + 90 × 2.32 = 500.8 mm
Module décennal sec
x0.10 = 292 + 90 × (− 1.28) = 176.8 mm
Module centenaire sec
x0.10 = 292 + 90 × (− 2.32) = 83.2 mm
6. Intervalle de confiance des quantiles
s
IC pour x p = x p ± u 1−α u 2p + 2
2 2n
α = 95% ⇒ u1−α = u0.025 = 1.96
2
IC pour le module décennal humide
90
IC pour x0.90 = 407.2 ± 1.96 (1.28)² + 2 = 407.2 ± 30.5
2 × 61
30.5
Le module décennal humide est connu à ≈ 7%
407.2
IC pour le module centenaire humide
90
IC pour x0.99 = 500.8 ± 1.96 (2.32)² + 2 = 500.8 ± 43.4
2 × 61
43.4
Le module centenaire humide est connu à ≈ 9%
500.8
IC pour le module décennal sec
90
IC pour x0.10 = 176.8 ± 1.96 (1.28)² + 2 = 176.8 ± 30.5
2 × 61
30.5
Le module décennal sec est connu à ≈ 17%
176.8
IC pour le module centenaire sec
90
IC pour x0.01 = 83.2 ± 1.96 (2.32)² + 2 = 83.2 ± 43.4
2 × 61
43.4
Le module centenaire sec est connu à ≈ 52%
83.2