SignesVitauxSAS

libname SOURCE  "/folders/myfolders/Miniprojet/source";    /* endroit où se trouvent tous les datasets sources et les formats */
options fmtsearch=( SOURCE  WORK);
options ls=256;

/* Signes vitaux*/
ODS GRAPHICS ON;

/* Création variable Tension*/
DATA SV;
    SET SOURCE.vital_signs;
    ATTRIB Tension
        label="Diagnostique tension"
        length=$25;
    IF SYS=. OR DIA=. THEN
        Tension="Non Calculé";
    ELSE IF SYS < 120 AND DIA < 80 THEN
        Tension="Tension optimale";
    ELSE IF SYS <= 129 AND DIA <= 84 THEN
        Tension="Tension normale";
    ELSE IF SYS <= 139 AND DIA <= 89 THEN
        Tension="Tension normale élevée";
    ELSE IF SYS <= 159 AND DIA <= 99 THEN
        Tension="Hypertension légère";
    ELSE IF SYS <= 179 AND DIA <= 109 THEN
        Tension="Hypertension modérée";
    ELSE IF SYS > 180 AND DIA > 110 THEN
        Tension="Hypertension sévère";
    ELSE
        Tension="Hypertension sévère";
RUN;

/*Fusion avec date de visite*/
PROC SQL;
    CREATE TABLE VS AS
    SELECT s.*, d.VISDT
    FROM SV as s LEFT JOIN SOURCE.date_of_visit as d
    ON s.USUBJID=d.USUBJID AND s.VISID=d.VISID;
QUIT;
/*Fusion entre VS et PE*/
/* A VERIFIER SI FONCTIONNE COMME IL FAUT*/
PROC SQL;
    CREATE TABLE VS_PE AS
    SELECT *
    FROM VS LEFT JOIN PE
    ON VS.USUBJID=PE.USUBJID;

/*Fusion avec numero de traitement*/
PROC SQL;
    CREATE TABLE VS_PE_Treatment AS
    SELECT v.*, t.TRTCD
    FROM VS as v LEFT JOIN SOURCE.treatment_assignment as t
    ON v.USUBJID=t.USUBJID;
QUIT;

PROC SORT DATA=VS_PE_Treatment OUT=VS_PE_Treatment_sorted;
    BY TRTCD;
RUN;

/*PROC FREQ : frequence differents parametres de signes vitaux*/
TITLE "Fréquence des parametres vitaux en fonction du groupe de traitement";
PROC FREQ DATA=VS_PE_Treatment_sorted NOPRINT;
    by TRTCD;
    tables SYS / out=FreqVSSys;
RUN;

PROC FREQ DATA=VS_PE_Treatment_sorted NOPRINT;
    by TRTCD;
    tables TEMP / out=FreqVSTemp;
RUN;

TITLE "FREQ SYS";
PROC PRINT DATA=FreqVSSys NOOBS LABEL;
RUN;

TITLE "FREQ TEMP";
PROC PRINT DATA=FreqVSTemp NOOBS LABEL;
RUN;

/*Test égalité des moyennes*/
TITLE "PROC ANOVA";
PROC ANOVA DATA=VS_PE_Treatment_sorted;
    CLASS TRTCD;
    MODEL SYS=TRTCD;
    MEANS TRTCD;
RUN;

PROC SQL;
    CREATE TABLE VS_PE_Mod as
    select TRTCD, avg(SYS) as moyenne, count(TRTCD) as count
    from VS_PE_Treatment_sorted
    GROUP BY TRTCD;
QUIT;

TITLE "PROC FREQ option X²";
PROC FREQ DATA=VS_PE_Mod;
    weight count;
    tables TRTCD*moyenne /nocum chisq;
RUN;

ODS GRAPHICS OFF;