Vous n'êtes pas identifié.
Vous êtes sur le forum du master ESA !
Le site du master ESA - description de la formation, notes de cours, contacts... vient de déménager !!!
Venez visiter notre nouveau site : www.master-esa.fr
Bonjour,
j'ai un problème d'autant plus énervant que je suis sure que la solution est évidente :
J'essai de tester l'adéquation de ma distribution à la loi lognormale et gamma que j'ai simulée. Pour ca je dois faire un test du chi2. J'ai trouvé à droite à gauche des commencements de solutions qui m'ont plus embrouillé qu'autre chose:
comme ca par exemple:
proc capability data=tmp1.tabzone ;
var mag;
by zone;
histogram / midpoints= 0.025 to 0.975 by 0.05
OUTHISTOGRAM = tmp1.sortie ;
run ;
data tmp1.chisqtest1 ;
set tmp1.sortie ;
if zone=1;
s =_n_ ;
retain cumul ;
if _n_=1 then cumul=(_OBSPCT_-5)**2/5 ;
else cumul=cumul+(_OBSPCT_-5)**2/5 ;
drop s ;
run ;
data tmp1.chisqtest2 ;
set tmp1.sortie ;
if zone=2;
s =_n_ ;
retain cumul ;
if _n_=1 then cumul=(_OBSPCT_-5)**2/5 ;
else cumul=cumul+(_OBSPCT_-5)**2/5 ;
drop s ;
run ;
proc print data=tmp1.chisqtest1 (obs=5) ;
run ;
alors c'est bien mais je ne vois pas quoi en faire... en plus ça m'étonne qu'il n'existe pas d'option pour faire ca (??!!).
Du coup si quelqu'un la déjà fait je veux bien un peu d'aide.
Merci!
Hors ligne
Vous découpez votre échantillon en classes (par exemple la variable classe prenant les valeurs 1, 2,..., et la variable eff contenant les effectifs de chacune des classes). Sous votre loi supposée vous évaluez les proportions théoriques (p1, p2,...) de chacune de ces classes puis vous faites quelque chose comme :
proc freq;
tables classe / testp=(p1,p2,...);
weight eff;
run;
Naturellement le choix des classes affecte les propriétés du test. En général, il est conseillé de les choisir de sorte que les p1, p2,.. soient égales.
Hors ligne