use "\\file3.intranet.chuv\data3\UNISANTE_DESS\S_SCS_GRSA\ADO\Documents de JC Suris\SNF OH\GenFRee\GenFRee_original_data_and_classification.dta" 
svyset _n [pweight=v1], vce(linearized) singleunit(missing)

recode Classif_5 1=3 2=2 3=1 4=4 5=5

label define Classif_5 1 "Light" 2 "Heavy" 3 "Abstinent" 4 "Indécis"  5 "Lateonset light"
label val Classif_5 Classif_5

gen classif3=.
replace classif3=1 if (Classif_5==1 | Classif_5==4 | Classif_5==5)
replace classif3=2 if (Classif_5==2)
replace classif3=3 if (Classif_5==3)
label define classif3 1 "Light" 2 "Heavy" 3 "Abstinent"
label val classif3 classif3

gen classheavy=.
replace classheavy=0 if (classif3==2)
replace classheavy=1 if (classif3==3)

gen classlight=.
replace classlight=0 if (classif3==1)
replace classlight=1 if (classif3==3)

gen lightgroup=.
replace lightgroup=1 if (Classif_5==1)
replace lightgroup=2 if (Classif_5==4)
replace lightgroup=3 if (Classif_5==5)
label define lightgroup 1 "Light" 2  "Indécis"  3"Lateonset light"
label val lightgroup lightgroup

gen T1_filiere3=.
replace T1_filiere3=0 if (T1_Q5==1)
replace T1_filiere3=1 if (T1_Q5==2 | T1_Q5==4)
replace T1_filiere3=2 if (T1_Q5==3 | T1_Q5==5 | T1_Q5==6 | T1_Q5==7)
label define T1_filiere3 0 "Gymnase" 1 "École" 2 "FP"
label val T1_filiere3 T1_filiere3


gen T1_filiere4=.
replace T1_filiere4=1 if (T1_Q5==1)
replace T1_filiere4=2 if (T1_Q5==2 | T1_Q5==4)
recode T1_filiere4 1=1 2=2 else=3
label define T1_filiere4 1 "Gymnase" 2 "École" 3 "FP"
label val T1_filiere4 T1_filiere4

gen T1_filiere5=.
replace T1_filiere5=3 if (T1_Q5==1)
replace T1_filiere5=2 if (T1_Q5==2 | T1_Q5==4)
replace T1_filiere5=1 if (T1_Q5==3)
replace T1_filiere5=4 if (T1_Q5==5 | T1_Q5==6 | T1_Q5==7)
label define T1_filiere5 3 "Gymnase" 2 "École" 1 "FP" 4"Autre"
label val T1_filiere5 T1_filiere5

gen parents2=.
replace parents2=0 if (T1_Q13==1)
replace parents2=1 if (T1_Q13>1)

gen parents2R=.
replace parents2R=0 if (parents2==1)
replace parents2R=1 if (parents2==0)


gen tabac2=.
replace tabac2=0 if (T1_Q46==1 | T1_Q46==2)
replace tabac2=1 if (T1_Q46>2)

gen maria302=.
replace maria302=0 if (T1_Q49_2==1)
replace maria302=1 if (T1_Q49_2>1)

**Variable T1health
gen T1health=.
replace T1health =0 if (T1_Q31<4)
replace T1health =1 if (T1_Q31>3)

**Variable AVENIR
gen T1avenir=T1_Q45_1+T1_Q45_2+T1_Q45_3+T1_Q45_4+T1_Q45_5+T1_Q45_6
alpha T1_Q45_1 T1_Q45_2 T1_Q45_3 T1_Q45_4 T1_Q45_5 T1_Q45_6

gen socialsup=T1_Q61_3+ T1_Q61_4+ T1_Q61_5+ T1_Q61_10
alpha T1_Q61_3 T1_Q61_4 T1_Q61_5 T1_Q61_10 

gen parentsup=T1_Q61_1 + T1_Q61_2
alpha T1_Q61_1  T1_Q61_2

gen sociallife=T1_Q62_1+ T1_Q62_2+ T1_Q62_3+ T1_Q62_4
alpha T1_Q62_1 T1_Q62_2 T1_Q62_3 T1_Q62_4


gen abstinent=.
replace abstinent=1 if (Classif_5==1)
replace abstinent=0 if (Classif_5>1)

gen ses=.
replace ses=0 if (T1_Q15<3)
replace ses=1 if (T1_Q15==3)
label define ses 0 "Autre" 1 "Dessous"
label val ses ses

gen eleve=.
replace eleve=0 if (T1_Q10!=3)
replace eleve=1 if (T1_Q10==3)
label define eleve 0 "Autre" 1 "Dessous"
label val eleve eleve

alpha T0_Q41_1 T0_Q41_2 T0_Q41_3 T0_Q41_4 T0_Q41_5




****ANALYSES T1/2

svy linearized: tab Classif_5 T1_Q1, count col row
svy linearized: mean T1_Q2, over (Classif_5)
svy linearized: regress T1_Q2 Classif_5
svy linearized: tab Classif_5 parents2, count col row
svy linearized: tab Classif_5 T1_filiere5, count col row
svy linearized: tab Classif_5 eleve, count col row
svy linearized: tab Classif_5 T1_Q11, count col row
svy linearized: mean T1_Q14M, over (Classif_5)
svy linearized: regress T1_Q14M Classif_5
svy linearized: mean T1_Q14P, over (Classif_5)
svy linearized: regress T1_Q14P Classif_5
svy linearized: tab Classif_5 ses, count col row
svy linearized: tab Classif_5 T1_Q16, count col row
svy linearized: tab Classif_5 WHO5, count col row
svy linearized: mean T1_Q26, over (Classif_5)
svy linearized: regress T1_Q26 Classif_5
svy linearized: tab T1health Classif_5,count col row
svy linearized: mean T1_Q39, over (Classif_5)
svy linearized: regress T1_Q39 Classif_5
svy linearized: mean T1avenir, over (Classif_5)
svy linearized: regress T1avenir Classif_5
svy linearized: tab tabac2 Classif_5,count col row
svy linearized: tab maria302 Classif_5,count col row
svy linearized: mean socialsup, over (Classif_5)
svy linearized: regress socialsup Classif_5
svy linearized: mean parentsup, over (Classif_5)
svy linearized: regress parentsup Classif_5
svy linearized: mean sociallife, over (Classif_5)
svy linearized: regress sociallife Classif_5


***ANALYSE MULTIVARIEE BACKWARD

svy linearized: mlogit Classif_5 T1_Q1 T1_Q2 parents2R i.T1_filiere5 T1_Q11 T1_Q16 tabac2 maria302, rrr base(3)

**BETWEEN LIGHTS
svy linearized: tab lightgroup T1_Q1, count col row
svy linearized: mean T1_Q2, over (lightgroup)
svy linearized: regress T1_Q2 lightgroup
svy linearized: tab lightgroup parents2R, count col row
svy linearized: tab lightgroup T1_filiere5, count col row
svy linearized: tab lightgroup eleve, count col row
svy linearized: tab lightgroup T1_Q11, count col row
svy linearized: mean T1_Q14M, over (lightgroup)
svy linearized: regress T1_Q14M lightgroup
svy linearized: mean T1_Q14P, over (lightgroup)
svy linearized: regress T1_Q14P lightgroup
svy linearized: tab lightgroup ses, count col row
svy linearized: tab lightgroup T1_Q16, count col row
svy linearized: tab lightgroup WHO5, count col row
svy linearized: mean T1_Q26, over (lightgroup)
svy linearized: regress T1_Q26 lightgroup
svy linearized: tab T1health lightgroup,count col row
svy linearized: mean T1_Q39, over (lightgroup)
svy linearized: regress T1_Q39 lightgroup
svy linearized: mean T1avenir, over (lightgroup)
svy linearized: regress T1avenir lightgroup
svy linearized: tab tabac2 lightgroup,count col row
svy linearized: tab maria302 lightgroup,count col row
svy linearized: mean socialsup, over (lightgroup)
svy linearized: regress socialsup lightgroup
svy linearized: mean parentsup, over (lightgroup)
svy linearized: regress parentsup lightgroup
svy linearized: mean sociallife, over (lightgroup)
svy linearized: regress sociallife lightgroup


***ABSTINENT vs. HEAVY
svy linearized: tab classheavy T1_Q1, count col row
svy linearized: mean T1_Q2, over (classheavy)
svy linearized: regress T1_Q2 classheavy
svy linearized: tab classheavy parents2R, count col row
svy linearized: tab classheavy T1_filiere5, count col row
svy linearized: tab classheavy eleve, count col row
svy linearized: tab classheavy T1_Q11, count col row
svy linearized: mean T1_Q14M, over (classheavy)
svy linearized: regress T1_Q14M classheavy
svy linearized: mean T1_Q14P, over (classheavy)
svy linearized: regress T1_Q14P classheavy
svy linearized: tab classheavy ses, count col row
svy linearized: tab classheavy T1_Q16, count col row
svy linearized: tab classheavy WHO5, count col row
svy linearized: mean T1_Q26, over (classheavy)
svy linearized: regress T1_Q26 classheavy
svy linearized: tab T1health classheavy,count col row
svy linearized: mean T1_Q39, over (classheavy)
svy linearized: regress T1_Q39 classheavy
svy linearized: mean T1avenir, over (classheavy)
svy linearized: regress T1avenir classheavy
svy linearized: tab tabac2 classheavy,count col row
svy linearized: tab maria302 classheavy,count col row
svy linearized: mean socialsup, over (classheavy)
svy linearized: regress socialsup classheavy
svy linearized: mean parentsup, over (classheavy)
svy linearized: regress parentsup classheavy
svy linearized: mean sociallife, over (classheavy)
svy linearized: regress sociallife classheavy


svy linearized: logit classheavy T1_Q1 T1_Q2 parents2R i.T1_filiere5 T1_Q11 T1_Q16 T1_Q26 T1_Q39 tabac2 maria302 sociallife, or
svy linearized: logit classheavy T1_Q1 T1_Q2 parents2R i.T1_filiere5 T1_Q11 T1_Q16 T1_Q26 tabac2 maria302 sociallife, or
svy linearized: logit classheavy T1_Q1 T1_Q2 parents2R T1_Q11 T1_Q16 T1_Q26 tabac2 maria302 sociallife, or



***ABSTINENT vs. LIGHT
svy linearized: tab classlight T1_Q1, count col row
svy linearized: mean T1_Q2, over (classlight)
svy linearized: regress T1_Q2 classlight
svy linearized: tab classlight parents2R, count col row
svy linearized: tab classlight T1_filiere5, count col row
svy linearized: tab classlight eleve, count col row
svy linearized: tab classlight T1_Q11, count col row
svy linearized: mean T1_Q14M, over (classlight)
svy linearized: regress T1_Q14M classlight
svy linearized: mean T1_Q14P, over (classlight)
svy linearized: regress T1_Q14P classlight
svy linearized: tab classlight ses, count col row
svy linearized: tab classlight T1_Q16, count col row
svy linearized: tab classlight WHO5, count col row
svy linearized: mean T1_Q26, over (classlight)
svy linearized: regress T1_Q26 classlight
svy linearized: tab T1health classlight,count col row
svy linearized: mean T1_Q39, over (classlight)
svy linearized: regress T1_Q39 classlight
svy linearized: mean T1avenir, over (classlight)
svy linearized: regress T1avenir classlight
svy linearized: tab tabac2 classlight,count col row
svy linearized: tab maria302 classlight,count col row
svy linearized: mean socialsup, over (classlight)
svy linearized: regress socialsup classlight
svy linearized: mean parentsup, over (classlight)
svy linearized: regress parentsup classlight
svy linearized: mean sociallife, over (classlight)
svy linearized: regress sociallife classlight

svy linearized: logit classlight T1_Q1 T1_Q2 T1_Q11 T1_Q26 maria302, or
svy linearized: logit classlight T1_Q2 T1_Q11 T1_Q26 maria302, or
svy linearized: logit classlight T1_Q11 T1_Q26 maria302, or



***ABSTINENT vs. LIGHT vs. HEAVY
svy linearized: tab classif3 T1_Q1, count col row
svy linearized: mean T1_Q2, over (classif3)
svy linearized: regress T1_Q2 i.classif3
svy linearized: tab classif3 parents2R, count col row
svy linearized: tab classif3 T1_filiere, count col row
svy linearized: tab classif3 eleve, count col row
svy linearized: tab classif3 T1_Q11, count col row
svy linearized: mean T1_Q14M, over (classif3)
svy linearized: regress T1_Q14M classif3
svy linearized: mean T1_Q14P, over (classif3)
svy linearized: regress T1_Q14P classif3
svy linearized: tab classif3 ses, count col row
svy linearized: tab classif3 T1_Q16, count col row
svy linearized: tab classif3 WHO5, count col row
svy linearized: mean T1_Q26, over (classif3)
svy linearized: regress T1_Q26 classif3
svy linearized: tab T1health classif3,count col row
svy linearized: mean T1_Q39, over (classif3)
svy linearized: regress T1_Q39 classif3
svy linearized: mean T1avenir, over (classif3)
svy linearized: regress T1avenir classif3
svy linearized: tab tabac2 classif3,count col row
svy linearized: tab maria302 classif3,count col row
svy linearized: mean socialsup, over (classif3)
svy linearized: regress socialsup classif3
svy linearized: mean parentsup, over (classif3)
svy linearized: regress parentsup classif3
svy linearized: mean sociallife, over (classif3)
svy linearized: regress sociallife classif3

***MULTIVARIE

svy linearized: mlogit classif3 T1_Q1 T1_Q2 parents2R i.T1_filiere5 T1_Q11 T1_Q16 tabac2 maria302, rrr base (1)