[PDF] Procdure de tlchargement du logiciel R

View - Download Procdure de tlchargement du logiciel R

Previous PDF

Next PDF

R-Commander : Notions du chapitre 3

Estimation et tests d'hypothèses : problèmes à un échantillon

1) Présentation du jeu de données...........................

.....................2

2) Estimation et test sur une moyenne ȝ.........................................3

a) Test bilatéral et intervalle de confiance.........................3 b) Tests unilatéraux....................................... ............4 c) Quantiles et probabilités de la loi de Student...................5 d) Vérification de la normalité........................................ 7

Graphique quantile-quantile

Test de Shapiro et Wilk

3) Estimation d'une variance ı

2 a) Estimation ponctuelle....................................... .......8 b) Quantiles de la loi du khi-deux....................................8

4) Estimation et test sur une proportion p.....................................10

a) Test bilatéral et intervalle de confiance.........................10 b) Tests unilatéraux....................................... .............12 c) Quantiles et probabilités de la loi normale......................12 1

1) Présentation du jeu de données

Considérons le jeu de données

Oiseaux.xls utilisé dans les exercices du chapitre 3. Il est constitué de diverses mesures sur un échantillon de 49 oiseaux : - la survie ou non à un traitement - la longueur de l'oiseau - l'étendue des ailes - la longueur de la tête - la longueur de l'humérus - la longueur du sternum.

En voici un aperçu.

À l'aide des statistiques descriptives (voir chapitre 1), on peut calculer la moyenne et l'écart type des variables continues : On pourrait éventuellement construire des graphiques pour illustrer la dispersion de chaque variable. Il est important de faire une bonne exploration des données avant de se lancer dans l'analyse ou l' inférence sur les paramètres. 2

2) Estimation et test sur une moyenne ȝ

Les tests d'hypothèses et les intervalles de confiance s'obtiennent à partir des mêmes fenêtres de commandes. Dans le cas d'une moyenne ȝ, le seul test disponible est le test

utilisant la loi de Student (appelé t-test univarié), pour lequel on suppose que les données

sont issues d'une loi normale dont on ignore la variance ı 2 (cette dernière est estimée par la variance échantillonnale S 2 a) Test bilatéral et intervalle de confiance L'intervalle de confiance et le test bilatéral pour l'étendue moyenne ȝ s'obtiennent par la même série de commandes. sélectionner la variable d'intérêt ; spécifier l'hypothèse alternative Moyenne de la population = mu0 (on devrait lire Moyenne de la population mu0 pour bien représenter H 1 préciser la valeur de ȝ 0 dans la case Hypothèse nulle (ceci influence le test seulement, et non l'intervalle de confiance) ; spécifier la valeur de 1-Į désirée dans la case Niveau de confiance (ainsi la valeur du seuil Į en sera automatiquement déduite) ; cliquer sur OK. Testons par exemple si l'étendue moyenne diffère de 240 mm, au seuil de 1%. 3 Les résultats apparaissent dans la fenêtre de sortie : L'intervalle de confiance à 99% est [239.38 ; 243,27]. Pour réaliser le test bilatéral, la valeur observée de la statistique T obs doit être positionnée par rapport à la loi de Student à 48 degrés de liberté.

8323.149/0678.52403265.241

0 nsxT obs

Puisque le test est bilatéral et que T

obs est positif, le p-value a été obtenu en calculant la quantité suivante :

07312.003656.02)8323.1(2)(2

481
tPTtP obsn

L'hypothèse H

0 : ȝ = 240 n'est pas rejetée au seuil de Į = 0.01, car le p-value est supérieur à 0.01. (H 0 ne serait pas rejetée à 5% non plus, mais serait rejetée

à 10%.)

b) Tests unilatéraux Les commandes à entrer sont presque les mêmes que pour les tests bilatéraux, à

l'exception de l'hypothèse alternative, qui doit cette fois être spécifiée selon le test :

Moyenne de la population < mu0 pour un test unilatéral à gauche (H 1 0 Moyenne de la population > mu0 pour un test unilatéral à droite (H 1 0 On montre ci-dessous un exemple pour tester si l'étendue moyenne est supérieure à

240, au seuil de 1%. (H

0 : ȝ=240 versus H 1 : ȝ > 240, avec Į = 0.01). Remarquons que le seuil observé est la moitié du seuil pour le test bilatéral : . )8323.1( 48
tP Notons que l'intervalle de confiance calculé ici n'est pas l'intervalle dont nous avons discuté. Nous n'en tiendrons pas compte. 4 c) Quantiles et probabilités de la loi de Student On peut aussi calculer séparément les quantités formant l'intervalle de confiance : nstx n2/,1 La moyenne, l'écart-type et la taille d'échantillon s'obtiennent en calculant les statistiques descriptives. Il ne manque que le quantile d'ordre 1-Į/2 de la loi de

Student à n-1 degrés de liberté, t

n-1,Į/2 On montre ci-dessous que, pour l'intervalle de confiance à 99%, t

48, 0.005

= 2.6822. 5 De même, on peut calculer la probabilité associée au p-value d'un test, une fois que la statistique T obs est calculée. Voici comment obtenir (le p- value est le double de cette probabilité pour le test bilatéral). )8323.1( 48
tP 6 d) Vérification de la normalité Le test de Student a comme postulat la normalité des observations. Pour s'assurer de la validité des conclusions, il est bon de vérifier si nos données sont normales. Nous avons vu 3 techniques pour poser un diagnostic : l'histogramme, le graphique quantile-quantile et le test de Shapiro et Wilk.

Graphique quantile-quantile

Un tel graphique peut servir à vérifier l'adéquation des données à n'importe quel modèle statistique. Celui qui nous intéresse est la loi normale. Il suffit donc de cocher la distribution normale et de juger si le graphique forme bien une droite ou non. 7

Test de Shapiro et Wilk

On effectue le test par la série de commandes Statistiques - Résumés - Test de normalité Shapiro-Wilk. On rejettera la normalité des données si le p-value est inférieur au seuil Į que l'on se fixe pour ce test (qui peut être différent du Į choisi pour le test de Student). Ici, on accepterait le modèle normal pour nos données pour tous les choix raisonnables de Į.

3) Estimation d'une variance ı

2 L'intervalle de confiance et les tests sur une variance ne sont pas disponibles dans R- Commander. Il faut donc calculer les quantités séparément. a) Estimation ponctuelle

La variance de la population, ı

2 , est estimée par la variance échantillonnale S 2 . On obtient cette valeur via les statistiques descriptives. Dans le cas de la variable Étendue, l'écart-type s = 5.0678, donc l'estimation de la variance est s 2 = 25.683. b) Quantiles de la loi du khi-deux

L'intervalle de confiance pour ı

2 de niveau 1- Į contient les quantiles d'ordre Į/2 et

1-Į /2 de la loi du khi-deux à n-1 degrés de liberté.

2

2/1,12

2 2/,12 )1(,)1( nn SnSn 8 On montre ci-dessous comment obtenir ces quantiles pour n=49 et Į = 0.05 :

On a calculé :

75.30
2

975.0,482

2/1,1 n et 02.69 2

025.0,482

2/,1 n

L'intervalle de confiance à 95% serait alors :

2

2/1,12

2 2/,12 nn SnSn 9

4) Estimation et test sur une proportion p

Les tests d'hypothèses et les intervalles de confiance s'obtiennent à partir des mêmes fenêtres de commandes. Dans le cas d'une proportion, le test dont nous avons discuté est le test approximatif utilisant le théorème limite central, donc valide pour les grandes tailles d'échantillon. a) Test bilatéral et intervalle de confiance On s'intéresse à la proportion de décès après un traitement dans la population : p. sélectionner la variable binaire d'intérêt (SURVIE) ; spécifier l'hypothèse alternative Proportion de la population = p0 (on devrait lire Proportion de la population p0 pour bien représenter H 1 préciser la valeur de p 0 dans la case Hypothèse nulle (ceci influence le test seulement, et non sur l'intervalle de confiance) ; spécifier la valeur de 1-Į désirée dans la case Niveau de confiance (ainsi la valeur du seuil Į sera automatiquement déduite) ; préciser que l'on souhaite le test avec approximation normale ; cliquer sur OK. 10 Les résultats apparaissent dans la fenêtre de sortie : On précise d'abord le nombre de oui et de non dans l'échantillon. La proportion échantillonnale de non (i.e. de décès) est .5714.049/28ˆp L'intervalle de confiance à 95% (en jaune) est [0.433 ; 0.700]. Pour le test bilatéral (en vert), la valeur observée de la statistique Z obs est présentée au carré, identifiée X-squared en référence à la loi du khi-deux. 2 000

1667.0408.0

49)6.01(6.06.05714.0

)1(ˆ obsobs nppppZ

Puisque le test est bilatéral et que Z

obs est négatif, le p-value a été obtenu en calculant la quantité suivante :

683.03415.02)408.0(2)(2ZPZZP

obs

L'hypothèse H

0 : p = 0.6 n'est pas rejetée au seuil de Į = 0.05, car le p-value est supérieur à 0.05. 11 b) Tests unilatéraux Les commandes à entrer sont presque les mêmes que pour le test bilatéral, à

l'exception de l'hypothèse alternative, qui doit cette fois être spécifiée selon le test :

Proportion de la population < p0 pour un test unilatéral à gauche (H 1 : p < p 0 Proportion de la population > p0 pour un test unilatéral à droite (H 1 : p > p 0 On montre ci-dessous un exemple pour tester si la proportion de décès est inférieure à

0.6, au seuil de 5%. (H

0 : p=0.6 versus H 1 :p < 0.6 , avec Į = 0.05). Remarquons que le seuil observé est la moitié du seuil pour le test bilatéral : = 0.3415. )408.0(ZP Notons que l'intervalle de confiance calculé ici n'est pas l'intervalle dont nous avons discuté. Nous n'en tiendrons pas compte. c) Quantiles et probabilités de la loi normale

Comme pour les autres tests, les quantités

requises pour conduire un test d'hypothèses ou pour construire un interva lle de confiance peuvent être calculées séparément. Pour obtenir les quantiles ou les probabilités associées à la loi normale, voir le document sur le chapitre 2. 12quotesdbs_dbs35.pdfusesText_40

[PDF] intervalle de confiance d'une moyenne excel

[PDF] unité commerciale définition

[PDF] climat définition cycle 3

[PDF] definition de meteorologie

[PDF] unité commerciale physique et virtuelle complémentaire

[PDF] definition meteo

[PDF] dispense cap petite enfance

[PDF] deaes

[PDF] formule variance

[PDF] problème du second degré seconde

[PDF] bpjeps

[PDF] moyenne nationale bac francais 2017

[PDF] moyenne nationale math bac s

[PDF] moyenne nationale bac philo 2015

[PDF] moyenne nationale bac physique 2016