use "E:\Projet OH\SNF\TREE1\TREE1_small_with_classification_20042022.dta" 


svyset _n [pweight=Pond], vce(linearized) singleunit(missing)


label define Classif_5 1 "Light" 2 "Heavy" 3 "Indécis" 4 "Abstinent" 5 "Light2heavy"
label val Classif_5 Classif_5

gen classheavy=.
replace classheavy=0 if (Classif_5==2 | Classif_5==5)
replace classheavy=1 if (Classif_5==4)

gen classlight=.
replace classlight=0 if (Classif_5==1 | Classif_5==3)
replace classlight=1 if (Classif_5==4)

gen t1drink=.
replace t1drink=0 if (t1drug1==1)
replace t1drink=1 if (t1drug1==2)
replace t1drink=2 if (t1drug1>2 & t1drug1<6)

gen t2drink=.
replace t2drink=0 if (t2drug1==1)
replace t2drink=1 if (t2drug1==2)
replace t2drink=2 if (t2drug1>2 & t2drug1<6)

gen age2=age/12

gen parents=.
replace parents=0 if (famstruc==2)
recode parents 0=0 else=1

gen parentsR=.
replace parentsR=0 if (parents==1)
replace parentsR=1 if (parents==0)

gen sibs=.
replace sibs=0 if (nsib==0)
replace sibs=1 if (nsib==1)
replace sibs=2 if (nsib>=2)

gen sibsord=.
replace sibsord=0 if (brthord<2)
replace sibsord=1 if (brthord>=2)

gen gymnasien=.
replace gymnasien=1 if (t1educ22==8)
recode gymnasien 1=1 else=0

gen ase1=.
replace ase1=1 if (t1grad==2 | t1grad==3)
recode ase1 1=1 else=0
label define ase1 0 "Autre" 1 "Oui"
label value ase1 ase1

gen ase2=.
replace ase2=1 if (t1crit06==2 | t1crit06==3)
recode ase2 1=1 else=0
label define ase2 0 "Autre" 1 "Oui"
label value ase2 ase2

gen ase3=.
replace ase3=1 if (t1crit07==2 | t1crit07==3)
recode ase3 1=1 else=0
label define ase3 0 "Autre" 1 "Oui"
label value ase3 ase3

gen ase4=.
replace ase4=1 if (t1crit11==2 | t1crit11==3)
recode ase4 1=1 else=0
label define ase4 0 "Autre" 1 "Oui"
label value ase4 ase4

gen ase5=.
replace ase5=1 if (t1crit12==2 | t1crit12==3)
recode ase5 1=1 else=0
label define ase5 0 "Autre" 1 "Oui"
label value ase5 ase5



gen asescale=ase1+ase2+ase3+ase4+ase5
alpha ase1 ase2 ase3 ase4 ase5 

alpha t1grad t1crit06 t1crit07 t1crit11 t1crit12

gen asegr=.
replace asegr=0 if (asescale==0)
replace asegr=1 if (asescale>0)
label define asegr 0 "Aucune" 1 "Au moins 1"
label value asegr asegr

gen family=.
replace family=1 if (famstruc==2)
recode family 1=1 else=0

gen track=.
replace track=1 if (t1educ22==2 |t1educ22==3 |t1educ22==4)
replace track=2 if (t1educ22==5 |t1educ22==6 |t1educ22==7)
replace track=3 if (t1educ22==8)
replace track=4 if (t1educ22==1 |(t1educ22>=9 & t1educ22<=22))
replace track=5 if (t1educ22>22)
label define track 1 "Apprentissage" 2 "Ecole supérieure" 3 "Gymnase" 4 "Autre" 5 "Rien"
label val track track

recode t1clev03 2=1 else=0
recode t1clev05 2=1 else=0
recode t1clev06 2=1 else=0
recode t1clev07 2=1 else=0
recode t1clev08 2=1 else=0
recode t1clev09 2=1 else=0
recode t1clev10 2=1 else=0

gen ASE= t1grad +t1crit06+ t1crit07 +t1crit11+ t1crit12
alpha t1grad t1crit06 t1crit07 t1crit11 t1crit12

gen stomac=.
replace stomac=0 if (t1heal1==1 | t1heal1==2)
replace stomac=1 if (t1heal1==3 | t1heal1==4)

gen appetit=.
replace appetit=0 if (t1heal2==1 | t1heal2==2)
replace appetit=1 if (t1heal2==3 | t1heal2==4)

gen concentre=.
replace concentre=0 if (t1heal3==1 | t1heal3==2)
replace concentre=1 if (t1heal3==3 | t1heal3==4)

gen dos=.
replace dos=0 if (t1heal4==1 | t1heal4==2)
replace dos=1 if (t1heal4==3 | t1heal4==4)

gen vertige=.
replace vertige=0 if (t1heal5==1 | t1heal5==2)
replace vertige=1 if (t1heal5==3 | t1heal5==4)

gen dormir=.
replace dormir=0 if (t1heal6==1 | t1heal6==2)
replace dormir=1 if (t1heal6==3 | t1heal6==4)

gen nerveux=.
replace nerveux=0 if (t1heal7==1 | t1heal7==2)
replace nerveux=1 if (t1heal7==3 | t1heal7==4)

gen fatigue=.
replace fatigue=0 if (t1heal8==1 | t1heal8==2)
replace fatigue=1 if (t1heal8==3 | t1heal8==4)

gen céphale=.
replace céphale=0 if (t1heal9==1 | t1heal9==2)
replace céphale=1 if (t1heal9==3 | t1heal9==4)

gen t1tabac=.
replace t1tabac=1 if (t1drug2>1 & t1drug2<6)
recode t1tabac 1=1 else=0

gen t1maria=.
replace t1maria=1 if (t1drug4>1 & t1drug4<6)
recode t1maria 1=1 else=0

gen sex2=.
replace sex2=1 if (sex==2)
replace sex2=2 if (sex==1)
label define sex2 1 "Male" 2 "Female"
label val sex2 sex2

gen swiss=.
replace swiss=0 if (st16n01==1)
replace swiss=1 if (st16n01>1)

gen swiss2=.
replace swiss2=1 if (st16n01==1)
replace swiss2=0 if (st16n01>1)

gen notesinsuf=.
replace notesinsuf=0 if (t1grad==1)
replace notesinsuf=1 if (t1grad>1)

gen notesinsuf2=.
replace notesinsuf2=1 if (t1grad==2 | t1grad==3)
recode notesinsuf2 1=1 else=0

gen estsoi=t1sele1+t1sele2+t1sele3+t1sele4+t1sele5

gen depre=t1seld1+t1seld2+t1seld3+t1seld4+t1seld5

gen positif=t1posl1+t1posl2+t1posl3+t1posl4+t1posl5

gen soutiensoc= t5pasi1+t5pasi2+t5pasi3+t5pasi4



svy linearized: tab Classif_5 sex2, count col row
svy linearized: mean age2, over (Classif_5)
svy linearized: regress age2 Classif_5
svy linearized: tab reg_ling Classif_5, count col row
svy linearized: tab Classif_5 agglo, count col row
svy linearized: tab Classif_5 parents, count col row
svy linearized: tab Classif_5 sibs, count col row
svy linearized: tab Classif_5 sibsord, count col row
svy linearized: tab Classif_5 swiss, count col row
svy linearized: tab Classif_5 notesinsuf2, count col row
svy linearized: tab asegr Classif_5, count col row
svy linearized: tab track Classif_5, count col row
svy linearized: tab t1clev03 Classif_5, count col row
svy linearized: tab t1clev05 Classif_5, count col row
svy linearized: tab t1clev06 Classif_5, count col row
svy linearized: tab t1clev07 Classif_5, count col row
svy linearized: tab t1clev08 Classif_5, count col row
svy linearized: tab t1clev09 Classif_5, count col row
svy linearized: tab t1clev10 Classif_5, count col row
svy linearized: tab asegr Classif_5, count col row
svy linearized: tab stomac Classif_5, count col row
svy linearized: tab appetit Classif_5, count col row
svy linearized: tab concentre Classif_5, count col row
svy linearized: tab dos Classif_5, count col row
svy linearized: tab vertige Classif_5, count col row
svy linearized: tab dormir Classif_5, count col row
svy linearized: tab nerveux Classif_5, count col row
svy linearized: tab fatigue Classif_5, count col row
svy linearized: tab céphale Classif_5, count col row
svy linearized: tab t1tabac Classif_5, count col row
svy linearized: tab t1maria Classif_5, count col row
svy linearized: mean estsoi, over (Classif_5)
svy linearized: regress estsoi Classif_5
svy linearized: mean depre, over (Classif_5)
svy linearized: regress depre Classif_5
svy linearized: mean positif, over (Classif_5)
svy linearized: regress positif Classif_5
svy linearized: mean soutiensoc, over (Classif_5)
svy linearized: regress soutiensoc Classif_5



****MULTIVARIE/3 THE GOOD ONE
svy linearized: mlogit Classif_5 sex2 age2 agglo swiss2 i.track t1clev03 t1clev07 t1tabac t1maria, rrr base(4)
svy linearized: mlogit Classif_5 sex2 age2 agglo swiss2 i.track t1clev03 t1tabac t1maria, rrr base(4)
svy linearized: mlogit Classif_5 sex2 age2 agglo swiss2 i.track t1clev03 t1tabac, rrr base(4)
svy linearized: mlogit Classif_5 sex2 agglo swiss2 i.track t1clev03 t1tabac, rrr base(4)


***ABSTINENT vs. HEAVY
svy linearized: tab classheavy sex2, count col row
svy linearized: mean age2, over (classheavy)
svy linearized: regress age2 classheavy
svy linearized: tab reg_ling classheavy, count col row
svy linearized: tab classheavy agglo, count col row
svy linearized: tab classheavy parentsR, count col row
svy linearized: tab classheavy sibs, count col row
svy linearized: tab classheavy sibsord, count col row
svy linearized: tab classheavy swiss2, count col row
svy linearized: tab classheavy notesinsuf2, count col row
svy linearized: tab track classheavy, count col row
svy linearized: tab t1clev03 classheavy, count col row
svy linearized: tab t1clev05 classheavy, count col row
svy linearized: tab t1clev06 classheavy, count col row
svy linearized: tab t1clev07 classheavy, count col row
svy linearized: tab t1clev08 classheavy, count col row
svy linearized: tab t1clev09 classheavy, count col row
svy linearized: tab t1clev10 classheavy, count col row
svy linearized: tab asegr classheavy, count col row
svy linearized: tab stomac classheavy, count col row
svy linearized: tab appetit classheavy, count col row
svy linearized: tab concentre classheavy, count col row
svy linearized: tab dos classheavy, count col row
svy linearized: tab vertige classheavy, count col row
svy linearized: tab dormir classheavy, count col row
svy linearized: tab nerveux classheavy, count col row
svy linearized: tab fatigue classheavy, count col row
svy linearized: tab céphale classheavy, count col row
svy linearized: tab t1tabac classheavy, count col row
svy linearized: tab t1maria classheavy, count col row
svy linearized: mean estsoi, over (classheavy)
svy linearized: regress estsoi classheavy
svy linearized: mean depre, over (classheavy)
svy linearized: regress depre classheavy
svy linearized: mean positif, over (classheavy)
svy linearized: regress positif classheavy
svy linearized: mean soutiensoc, over (classheavy)
svy linearized: regress soutiensoc classheavy


svy linearized: logit classheavy sex2 age2 agglo parentsR swiss2 i.track t1clev05 t1clev09 fatigue t1tabac t1maria, or
svy linearized: logit classheavy sex2 age2 agglo parentsR swiss2 i.track t1clev09 fatigue t1tabac t1maria, or
svy linearized: logit classheavy sex2 agglo parentsR swiss2 i.track t1clev09 fatigue t1tabac t1maria, or
svy linearized: logit classheavy sex2 agglo parentsR swiss2 i.track fatigue t1tabac t1maria, or
svy linearized: logit classheavy sex2 agglo parentsR swiss2 i.track fatigue t1tabac, or
svy linearized: logit classheavy sex2 parentsR swiss2 i.track fatigue t1tabac, or
svy linearized: logit classheavy sex2 swiss2 i.track fatigue t1tabac, or



***ABSTINENT vs. LIGHT 
svy linearized: tab classlight sex2, count col row
svy linearized: mean age2, over (classlight)
svy linearized: regress age2 classlight
svy linearized: tab reg_ling classlight, count col row
svy linearized: tab classlight agglo, count col row
svy linearized: tab classlight parentsR, count col row
svy linearized: tab classlight sibs, count col row
svy linearized: tab classlight sibsord, count col row
svy linearized: tab classlight swiss2, count col row
svy linearized: tab classlight notesinsuf2, count col row
svy linearized: tab track classlight, count col row
svy linearized: tab t1clev03 classlight, count col row
svy linearized: tab t1clev05 classlight, count col row
svy linearized: tab t1clev06 classlight, count col row
svy linearized: tab t1clev07 classlight, count col row
svy linearized: tab t1clev08 classlight, count col row
svy linearized: tab t1clev09 classlight, count col row
svy linearized: tab t1clev10 classlight, count col row
svy linearized: tab asegr classlight, count col row
svy linearized: tab stomac classlight, count col row
svy linearized: tab appetit classlight, count col row
svy linearized: tab concentre classlight, count col row
svy linearized: tab dos classlight, count col row
svy linearized: tab vertige classlight, count col row
svy linearized: tab dormir classlight, count col row
svy linearized: tab nerveux classlight, count col row
svy linearized: tab fatigue classlight, count col row
svy linearized: tab céphale classlight, count col row
svy linearized: tab t1tabac classlight, count col row
svy linearized: tab t1maria classlight, count col row
svy linearized: mean estsoi, over (classlight)
svy linearized: regress estsoi classlight
svy linearized: mean depre, over (classlight)
svy linearized: regress depre classlight
svy linearized: mean positif, over (classlight)
svy linearized: regress positif classlight
svy linearized: mean soutiensoc, over (classlight)
svy linearized: regress soutiensoc classlight

svy linearized: logit classlight sex2 age2 parentsR swiss2 i.track t1clev03 t1clev05 t1clev09 fatigue, or
svy linearized: logit classlight age2 parentsR swiss2 i.track t1clev03 t1clev05 t1clev09 fatigue, or
svy linearized: logit classlight age2 parentsR swiss2 i.track t1clev05 t1clev09 fatigue, or
svy linearized: logit classlight parentsR swiss2 i.track t1clev05 t1clev09 fatigue, or
svy linearized: logit classlight parentsR swiss2 i.track t1clev09 fatigue, or
svy linearized: logit classlight parentsR swiss2 i.track fatigue, or




***T1 & T2

svy linearized: tab t1drink, count col
svy linearized: tab t2drink, count col
svy linearized: tab t1drink t2drink, count col row
