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;