Français

/*étape1: description des variables*/

proc univariate

data=DMSTAT.Ue2score1;

var DatIDM Agei IMCi EtatSante ScoreHygienoM24;

run;

proc freq

data=DMSTAT.Ue2score1;

Tables SEXE NivEtud AntHypCho AntHypArt Tabac ActPhy ExpoOmeg;

run;

/*étape2: vérification des hypothèses*/

/* a-Indépendance: pas de raison ici de penser que les observations ne sont pas indépendantes*/

/* b-Normalité*/

Proc univariate data=DMSTAT.UE2SCORE1 normal;

var ScoreHygienoM24;

run;

/* pour avoir histogramme: proc univariate data=DMSTAT.UE2SCORE1;

var ScoreHygienoM24;

Histogramm Score HygienoM24 / normal;*/

/*run; */

/*c-Linéarité*/

data DMSTAT.Ue2score1;

set DMSTAT.Ue2score1;

/*âge*/

if Agei<55 then Agei_classe=1;

if Agei>=55 and Agei<65 then Agei_classe=2;

if Agei>=65 and Agei<75 then Agei_classe=3;

if Agei>=75 then Agei_classe=4;

if Agei=. then Agei_classe=.;

if Agei_classe=1 then Agei_55=1;

if Agei_classe not in (. 1) then Agei_55=0;

if Agei_classe=2 then Agei_5565=1;

if Agei_classe not in (. 2) then Agei_5565=0;

if Agei_classe=3 then Agei_6575=1;

if Agei_classe not in (. 3)then Agei_6575=0;

if Agei_classe=4 then Agei_75=1;

if Agei_classe not in (. 4)then Agei_75=0;

run;

proc freq data=DMSTAT.UE2score2;

tables Agei_classe Agei_55 Agei_5565 Agei_6575 Agei_75;

run;

/*IMCi*/

data DMSTAT.Ue2score1;

set DMSTAT.Ue2score1;

if IMCi<25 then IMCi_classe=1;

if IMCi>=25 and IMCi<27 then IMCi_classe=2;

if IMCi>=27 and IMCi<30 then IMCi_classe=3;

if IMCi>=30 then IMCi_classe=4;

if IMCi=. then IMCi_classe=.;

if IMCi_classe=1 then IMCi_25=1;

if IMCi_classe not in (. 1)then IMCi_25=0;

if IMCi_classe=2 then IMCi_2527=1;

if IMCi_classe not in (. 2) then IMCi_2527=0;

if IMCi_classe=3 then IMCi_2730=1;

if IMCi_classe not in(. 3)then IMCi_2730=0;

if IMCi_classe=4 then IMCi_30=1;

if IMCi_classe not in (. 4)then IMCi_30=0;

run;

proc freq data=DMSTAT.UE2score2;

tables IMCi_classe IMCi_25 IMCi_2527 IMCi_2730 IMCi_30;

run;

/*EtatSante*/

data DMSTAT.Ue2score1;

set DMSTAT.Ue2score1;

if EtatSante<6 then EtatSante_classe=1;

if EtatSante>=6 and EtatSante<8 then EtatSante_classe=2;

if EtatSante>=8 and EtatSante<9 then EtatSante_classe=3;

if EtatSante>=9 then EtatSante_classe=4;

if EtatSante=. then EtatSante_classe=.;

if EtatSante_classe=1 then EtatSante_6=1;

if EtatSante_classe not in (. 1) then EtatSante_6=0;

if EtatSante_classe=2 then EtatSante_68=1;

if EtatSante_classe not in (. 2) then EtatSante_68=0;

if EtatSante_classe=3 then EtatSante_89=1;

if EtatSante_classe not in(. 3)then EtatSante_89=0;

if EtatSante_classe=4 then EtatSante_9=1;

if EtatSante_classe not in (. 4) then EtatSante_9=0;

run;

proc freq data=DMSTAT.UE2score2;

tables EtatSante_classe EtatSante_6 EtatSante_68 EtatSante_89 EtatSante_9;

run;

/*vérification*/

proc freq data=DMSTAT.UE2score1;

tables Agei_classe Agei_55 Agei_5565 Agei_6575 Agei_75

IMCi_classe IMCi_25 IMCi_2527 IMCi_2730 IMCi_30

EtatSante_classe EtatSante_6 EtatSante_68 EtatSante_89 EtatSante_9;

run;

/*contrôle de la linéarité*/

/* âge*/

Proc reg data=DMSTAT.Ue2score2;

Model ScoreHygienoM24=Agei;

plot ScoreHygienoM24*Agei;

run;

Proc reg data=DMSTAT.Ue2score1;

Model ScoreHygienoM24= Agei_55 Agei_5565 Agei_6575 Agei_75;

ScoreHygienoM24: test Agei_55 Agei_5565 Agei_6575 Agei_75;

run;

Proc reg data=DMSTAT.Ue2score1;

Model ScoreHygienoM24= Agei_5565 Agei_6575 Agei_75; /*Réf: Agei_55*/

ScoreHygienoM24: test Agei_5565 Agei_6575 Agei_75;

run;

/*===>la linéarité est vérifié*/

/*===>on conserve Agei en continu*/

/*IMCi*/

Proc reg data=DMSTAT.Ue2score2;

Model ScoreHygienoM24=IMCi;

plot ScoreHygienoM24*IMCi;

run;

Proc reg data=DMSTAT.Ue2score1;

Model ScoreHygienoM24=IMCi_25 IMCi_2527 IMCi_2730 IMCi_30;

ScoreHygienoM24:test IMCi_25 IMCi_2527 IMCi_2730 IMCi_30;

run;

/*linéarité non verifié on prend une classe de référence: IMCi_25*/

Proc reg data=DMSTAT.Ue2score2;

Model ScoreHygienoM24=IMCi_2527 IMCi_2730 IMCi_30; /* Réf:IMCi_25*/

ScoreHygienoM24: test IMCi_2527 IMCi_2730 IMCi_30;

Run;

Proc reg data=DMSTAT.Ue2score1;

Model ScoreHygienoM24=IMCi_25 IMCI_2527 IMCi_27 IMCi_30;

ScoreHygienoM24: test IMCi_25 IMCi_2527 IMCi_27 IMCi_30;

run;

/*linéarité non

/*===>linéarité verifié,regroupement des classes IMCi_2527 et IMCi_2730 ensemble?*/

/*EtatSante*/

Proc reg data=DMSTAT.Ue2score2;

Model ScorehygienoM24=EtatSante;

plot ScoreHygienioM24*EtatSante;

run;

Proc reg data=DMSTAT.Ue2score2;

Model ScoreHygienoM24=EtatSante_6 EtatSante_68 EtatSante_89 EtatSante_9;

ScoreHygienoM24: test EtatSante_6 EtatSante_68, EtatSante_89, EtatSante_9;

run;

/*linéarité non verifié on prend une classe de référence: EtatSante_6*/

Proc reg data=DMSTAT.Ue2score2;

Model ScoreHygienoM24= EtatSante_68