---
# System prepended metadata

title: Untitled

---

/** TP noté SAS du Mardi 17 novembre 2020 **/

/** Question 1 **/

LIBNAME resultat "P:\TP-note\Resultat";


/**Question 2**/

PROC IMPORT datafile="P:\TP-note\Csinges.xlsx"
    OUT=Csinges
    DBMS=XLSX ;
RUN;

DATA Csinges ;
	SET Csinges ;
	rename Longueur_bras__cm_=Longueur_bras ;
	rename Poids__kg_ = Poids ;
RUN ;


/** Question 3 **/
DATA resultat.Csinges ;
    set Csinges ;
RUN;

/**Question 4**/

PROC FORMAT library=resultat ;
    VALUE $espece
        "Espece 1" = "Douroucouli"
        "Espece 2" = "Tamarin";
RUN;

/*Nous allons changer le nom des variables pour des raisons pratiques*/


/**Question 5**/

DATA resultat.Csinges;
    set resultat.Csinges;
    format Espece $espece.;
RUN ;

/**Question 6**/

/**Etape DATA du calcul de la moyenne arithmétique**/
DATA resultat.Statpoid;
    set resultat.Csinges;
    retain cumul 0;
    cumul=sum(cumul,Poids);
    cumul1=cumul/100;
    drop cumul;
    rename cumul1=moyenne;
RUN;

/**Etape DATA du calcul de la variance corrigée**/
DATA resultat.Statpoid;
    set resultat.Statpoid;
    intermediaire1=Poids-moyenne;
    intermediaire2=intermediaire1*intermediaire1;
    retain cumul 0;
    cumul=sum(cumul,intermediaire2);
    cumul1=cumul/99;
    drop cumul intermediaire1 intermediaire2 ;
    rename cumul1=variance_corrigee;
RUN;

/**Vérifier si calcul fonctionne sur SAS**/


/**Etape DATA du calcul de l'écart-type**/
DATA resultat.Statpoid;
    set resultat.Statpoid;
    ecart-type = sqrt(variance);
RUN ;


PROC MEANS data=resultat.Csinges mean var std NOPRINT NOOBS ;
    var Poids ;
    OUTPUT OUT=resultat.Statpoid_means (drop=_TYPE_ _FREQ_) MEAN=moyenne  VAR=variance_corrigee STD=ecart-type ;
RUN ;

/**Commentaires de comparaison entre méthode DATA et méthode PROC MEANS à faire**/

/**Question 7**/

PROC UNIVARIATE data=resultat.Csinges ;
    qqplot Poids ;
RUN ;


/**Question 8**/



PROC MEANS data=resultat.Csinges Q1 median Q3 min max NOOBS NOPRINT;
    var Poids;
    output out=resultat.Statpoid_extremes_quartiles (drop=_TYPE_ _FREQ_) min=minimum max=maximum Q1=quartile1 median=mediane Q3=quartile3 ;
    CLASS espece;
RUN;

PROC PRINT data=resultat.Statpoid_extremes_quartiles;
RUN;


/**Question 9**/
PROC SGPLOT resultat.Csinges;
    vbox Poids / category=Nombre de blessures apparentes ;
RUN;




/**Question 10**/

PROC FREQ data=resultat.Csinges ;
    tables Poids*Nombre de blessures apparentes / nocol nofreq chisq ;
RUN ;

/*Comme la valeur de "prob" (0,0833) est supérieure au seuil de signification (0,05), on ne possède pas suffisamment de preuves pour conclure que les variables sont associées*/

/**Vérifier nom de la variable sur SAS pour les blessures**/


/**Question 11**/

PROC SGPLOT data=resultat.Csinges;
    scatter X=Poids Y=Nombre_de_blessures_apparentes ;
RUN ;


/**Question 12**/
DATA resultat.Correl;
    set resultat.Csinges;
    retain cumula 0;
    cumula=sum(cumula,Longueur_bras);
    cumulb=cumula/100;
    drop cumula;
    rename cumulb=moyennelb;
/*moyenne longueur bras*/
RUN;

DATA resultat.Correl;
    set resultat.Csinges;
    intermediaire3=Poids-moyenne;
    intermediaire4=intermediaire3*intermediaire3;
    retain cumulc 0;
    cumulc=sum(cumulc,intermediaire4);
    cumuld=cumulc/99;
    drop cumulc intermediaire3 intermediaire4 ;
    rename cumuld=variance_corrigeelb;
/*variance corrigée de longueur de bras*/
RUN;

DATA resultat.Correl;
    set resultat.Correl;
    ecart-typelb = sqrt(variance_corrigeelb);
/*ecart-type longueur bras*/
RUN ;

/*nous avons observé précédement que l'écart-type du poids était de 4,75, et nous venons de calculer l'écart-type de la longueur de bras*/    

DATA resultat.Correl;
    set resultat.Correl;
    cov=COV(Poids,longueur_bras); /*Pas de fonction pour calculer la covariance sur SAS donc nous n'avons pas réussi à la calculer manuellement*/
RUN;

/*Nous venons de calculer la covariance, et nous allons calculer le coefficient de Pearson*/
DATA resultat.Correl;
    set resultat.Correl;
    Pearson=cov/(ecart-type*ecart-typelb);
RUN;

PROC PRINT data=resultat.Correl;
RUN;


PROC CORR data=resultat.Csinges;
RUN;

/**Question 13**/

PROC REG data=resultat.Csinges;
    MODEL Poids=Longueur_bras ;
RUN ;

/** Question 14 **/

DATA resultat.SCR ;
    SET resultat.Csinges ;
    