libname doc "C:/doc2/ologit";


/*https://www.researchgate.net/post/Reference_level_in_a_multiple_regression*/


proc freq data=Doc.data3;

table faculte;

run;

/*Paris 6, 15%*/



proc freq data=Doc.data3;

table depart;

run;

/*75, 15%*/


/**************************************************************************/

/*** SGR SGR SGR **********************************************************/

proc logistic data=Doc.data3 outmodel=Doc.Model4 desc;

class sexe ansm ;

model sgr = sexe xp ansm / lackfit;

run;


proc logistic inmodel=Doc.Model4;

score data=Doc.data3 FITSTAT;

run;


/****************************************************************************/


proc logistic data=Doc.data3 outmodel=Doc.Model5 ;

class ansmedit alertlog sexe depart (param=ref ref='75') faculte (param=ref ref='Paris 6') exercice zone vidal prescrire net ansm doroz bcb;

model sgr = ansmedit faculte ansm age;

run;


proc logistic inmodel=Doc.Model5;

score data=Doc.data3 FITSTAT;

run;


/********************************************************************************/

/********************************************************************************/



/**************************************************************************/

/*** SGE SGE SGE **********************************************************/


proc logistic data=Doc.data3 outmodel=Doc.Model6;

class ansmedit alertlog sexe depart (param=ref ref='75') faculte (param=ref ref='Paris 6') exercice zone vidal prescrire net ansm doroz bcb;

model sge = ansmedit alertlog sexe age depart faculte exercice zone vidal prescrire net ansm doroz bcb / selection=forward;

run;


proc logistic inmodel=Doc.Model6;

score data=Doc.data3 FITSTAT;

run;


/***************************************************************************/

proc logistic data=Doc.data3 outmodel=Doc.Model7;

class ansmedit alertlog sexe depart (param=ref ref='75') faculte (param=ref ref='Paris 6') exercice zone vidal prescrire net ansm doroz bcb;

model sge = ansmedit age;

run;






/**********************************************************************************/


proc logistic data=Doc.data3 outmodel=Doc.Model7;

class ansmedit alertlog sexe depart (param=ref ref='75') faculte (param=ref ref='Paris 6') exercice zone vidal prescrire net ansm doroz bcb;

model sge = ansmedit faculte ansm age;

run;


proc logistic inmodel=Doc.Model7;

score data=Doc.data3 FITSTAT;

run;



/*http://www.ats.ucla.edu/stat/sas/dae/ologit.htm*/











ods pdf file="C:/doc2/m1.pdf" style=HTMLBLUE;

OPTIONS NODATE NONUMBER;

proc logistic data=Doc.data3 outmodel=Doc.Model4;

class sexe ansm ;

model sgr = sexe xp ansm / lackfit;

run;

ods pdf close;





/*********************************************************************************************************/

/*********************************************************************************************************/


/*77-78-91-92-93-94-95 departement exercice*/

/*fac parisienne (UVSQ*/


proc freq data=Doc.data3;

table faculte;

run;


proc means data=Doc.data3;

var xp;

run;


proc univariate data=Doc.data3;

var xp;

run;




/*

Paris 11

Paris 12

Paris 13

Paris 5

Paris 6

Paris 7

UVSQ*/









data doc.data4;

set doc.data3;


if depart = 75 then depart2 = "Paris";

else if depart = 77 then depart2 = "Paris";

else if depart = 78 then depart2 = "Paris";

else if depart = 91 then depart2 = "Paris";

else if depart = 92 then depart2 = "Paris";

else if depart = 93 then depart2 = "Paris";

else if depart = 94 then depart2 = "Paris";

else if depart = 95 then depart2 = "Paris";

else depart2 = "NonPa";


if faculte = "Paris 11" then faculte2 = "Paris";

else if faculte = "Paris 12" then faculte2 = "Paris";

else if faculte = "Paris 13" then faculte2 = "Paris";

else if faculte = "Paris 5" then faculte2 = "Paris";

else if faculte = "Paris 6" then faculte2 = "Paris";

else if faculte = "Paris 7" then faculte2 = "Paris";

else if faculte = "UVSQ" then faculte2 = "Paris";

else faculte2 = "NonPa";



if xp = 0 then xp2 = 0.5;

else xp2=xp;

run;


/**sgr**/

proc logistic data=Doc.data4 outmodel=Doc.Model8 desc;

class ansmedit alertlog sexe depart (param=ref ref='75') faculte (param=ref ref='Paris 6') exercice zone vidal prescrire net ansm doroz bcb depart2 faculte2;

model sgr = age ansmedit alertlog sexe depart2 faculte2 exercice zone vidal prescrire net ansm doroz bcb / selection=forward;

run;



proc logistic data=Doc.data4 outmodel=Doc.Model8 desc;

class ansmedit alertlog sexe depart (param=ref ref='75') faculte (param=ref ref='Paris 6') exercice zone vidal prescrire net ansm doroz bcb depart2 faculte2;

model sgr = age ansmedit faculte2 ansm;

run;



/**sge**/

proc logistic data=Doc.data4 outmodel=Doc.Model8 desc;

class ansmedit alertlog sexe depart (param=ref ref='75') faculte (param=ref ref='Paris 6') exercice zone vidal prescrire net ansm doroz bcb depart2 faculte2;

model sge = xp ansmedit alertlog sexe depart2 faculte2 exercice zone vidal prescrire net ansm doroz bcb / selection=forward;

run;



proc logistic data=Doc.data4 outmodel=Doc.Model8 desc;

class ansmedit alertlog sexe depart (param=ref ref='75') faculte (param=ref ref='Paris 6') exercice zone vidal prescrire net ansm doroz bcb depart2 faculte2;

model sge = xp ansmedit sexe;

run;


/*https://mchankins.wordpress.com/2013/04/21/still-not-significant-2/*/



/**************************************************************************************/

/***encore***********/

proc logistic data=Doc.data4 outmodel=Doc.Model8 desc;

class ansmedit (param=ref ref='Non') alertlog sexe (param=ref ref='Feminin') depart (param=ref ref='75') faculte (param=ref ref='Paris 6') exercice zone vidal prescrire net ansm (param=ref ref="0") doroz bcb depart2 faculte2;

model sgr = age ansmedit faculte2 ansm;

run;


proc logistic data=Doc.data4 outmodel=Doc.Model8 desc;

class ansmedit (param=ref ref='Non') alertlog sexe (param=ref ref='Feminin') depart (param=ref ref='75') faculte (param=ref ref='Paris 6') exercice zone vidal prescrire net ansm (param=ref ref="0") doroz bcb depart2 faculte2;

model sge = xp ansmedit sexe;

run;