Test d'adequation : Récupération p-value pour macro

[SAStableau]

Bonjour,

Je me contente pour l'instant de faire des tests d'adéquation aux lois continues avec proc univariate. Afin de réaliser une macro, j'aimerais récupérer la p-value de chaque statistique de test (Kolmogorov, Cramer et Anderson). J'ai regardé toutes les sorties possibles de output de SAS et je n'ai rien trouvé.

Si quelqu'un a une idée, ça serait génial.

Merci d'avance

[fafabzh6]

Bonjour,
Pourrais-tu nous donner le code permettant de calculer tes statistiques de tests d'adéquation.

[SAStableau]

Voilà avec ce code, on a les tests d'adéquation :

Code :

1
2
3
4
5
6
7
proc univariate DATA=avoine;
var yield;
histogram / normal (l=1 color=red)
			/*weibull (l=2 color=blue)
			gamma (l=3 color=green)
			exponential (l=4 color=yellow)*/;
run;

S'il faut aussi des données, pas de problème, je peux les passer. Mais bon, avec n'importe quelle variable continue ça devrait marcher.

Merci d'avance

[fafabzh6]

Bonjour,
Si tu pouvais nous envoyer tes données cela pourrait aider.
Fafabzh

[SAStableau]

Voici mes données :

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
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
 
DATA avoine;
input Block Variety $ nitro yield;
cards;
1 Victory 0 111
1 Victory 0.2 130
1 Victory 0.4 157
1 Victory 0.6 174
1 GoldenRain 0 117
1 GoldenRain 0.2 114
1 GoldenRain 0.4 161
1 GoldenRain 0.6 141
1 Marvellous 0 105
1 Marvellous 0.2 140
1 Marvellous 0.4 118
1 Marvellous 0.6 156
2 Victory 0	61
2 Victory 0.2 91
2 Victory 0.4 97
2 Victory 0.6 100
2 GoldenRain 0 70
2 GoldenRain 0.2 108
2 GoldenRain 0.4 126
2 GoldenRain 0.6 149
2 Marvellous 0 96
2 Marvellous 0.2 124
2 Marvellous 0.4 121
2 Marvellous 0.6 144
3 Victory 0	68
3 Victory 0.2 64
3 Victory 0.4 112
3 Victory 0.6 86
3 GoldenRain 0 60
3 GoldenRain 0.2 102
3 GoldenRain 0.4 89
3 GoldenRain 0.6 96
3 Marvellous 0 89
3 Marvellous 0.2 129
3 Marvellous 0.4 132
3 Marvellous 0.6 124
4 Victory 0	74
4 Victory 0.2 89
4 Victory 0.4 81
4 Victory 0.6 122
4 GoldenRain 0 64
4 GoldenRain 0.2 103
4 GoldenRain 0.4 132
4 GoldenRain 0.6 133
4 Marvellous 0 70
4 Marvellous 0.2 89
4 Marvellous 0.4 104
4 Marvellous 0.6 117
5 Victory 0	62
5 Victory 0.2 90
5 Victory 0.4 100
5 Victory 0.6 116
5 GoldenRain 0 80
5 GoldenRain 0.2 82
5 GoldenRain 0.4 94
5 GoldenRain 0.6 126
5 Marvellous 0 63
5 Marvellous 0.2 70
5 Marvellous 0.4 109
5 Marvellous 0.6 99
6 Victory 0	53
6 Victory 0.2 74
6 Victory 0.4 118
6 Victory 0.6 113
6 GoldenRain 0 89
6 GoldenRain 0.2 82
6 GoldenRain 0.4 86
6 GoldenRain 0.6 104
6 Marvellous 0 97
6 Marvellous 0.2 99
6 Marvellous 0.4 119
6 Marvellous 0.6 121
;
run;

Encore merci de chercher.

[fafabzh6]

Bonjour,
Tu devrais utiliser l'option qqplot pour tester tes distributions.

[SAStableau]

Qqplot me donne une interprétation graphique et j'ai besoin d'une p-value pour automatiser une macro.
J'aimerais que ma macro prenne en paramètre ma variable puis compare la p-value du test à 0.05.
Il faut donc que je récupère la p-value. De ce fait, qqplot ne me donne aucun résultat théorique.

Je commence à douter de l'accessibilité de cette p-value puisqu'elle n'est pas programmée pour sortir par Output.

J'ai même réfléchis à une sortie ODS dans un classeur Excel et récupérer le numéro de la case pour atteindre la valeur de la p-value mais ce n'est pas possible.

Encore merci de ton aide. Si tu as la moindre idée, je suis preneuse !

[fafabzh6]

Bonjour,

Je pense moi aussi que SAS ne délivre pas de p-value pour les tests d'adéquation ... après je ne connais bien évidemment pas toutes les proc SAS.
En tout cas pour la proc univariate ... il fait juste des tests graphiques de type qqplot et une superposition de loi sur les données.

Dans tous les cas tiens nous au courant.

[SAStableau]

Pas de problème. En tout cas merci pour toute ton aide.

Je crois que je vais devoir passer sous R qui donne tous ces résultats rapidement.

Bonne journée

[fafabzh6]

Bonjour,
C'est aussi l'éventualité que je voyais.
Tu peux poster tes questions sur ce forum :
http://www.developpez.net/forums/f27/autres-langages/autres-langages/
Tu auras très certainement de très bon retour ... comme j'ai pu en avoir.

[green_fr]

Envoyé par SAStableau

Pas de problème. En tout cas merci pour toute ton aide.

Je crois que je vais devoir passer sous R qui donne tous ces résultats rapidement.

Bonne journée

Ou sous MatLab. C'est ce qu'on a fait, nous.

[fafabzh6]

Dans ce cas là, le très bon et fourni forum Matlab :
http://www.developpez.net/forums/f148/environnements-developpement/matlab/

[fafabzh6]

Un dernière solution consisterai à programmer via des macro programmes les p-values ....