*** BANCO 3 - LOW ***

* Variáveis quantitativas continuas
* AGE - continua
tab age low, row chi
logit low age
* AGE - categorica
generate age1 = age
recode age1 1/19=0 20/45=1
tab age1 low, row chi
logit low i.age1

* LWT - continua
tab lwt low, row chi
logit low lwt
* LWT - categorica
generate lwt1 = lwt
recode lwt1 80/131=0 132/150=1 151/250=2
tab lwt1 low, row chi
logit low i.lwt1
*LWT - categorica
generate lwt2 = lwt
recode lwt2 80/130=0 131/150=1 151/250=2
tab lwt2 low, row chi
logit low i.lwt2

* Variaveis qualitativas categorica (3 categorias)
tab race low, row chi
logit low i.race
tab low race, chi
prop low race
tabulate low race, cchi2 cell chi2 column expected row rowsort colsort

* Variaveis qualitativas dicotomicas (0 e 1)
* SMOKE
tab smoke low, row chi
logit low i.smoke
cc smoke low
cs low smoke

* PTL
tab ptl low, row chi
generate ptl1 = ptl
recode ptl1 1/3=1
tab ptl1
logit low ptl1
cc ptl1 low
cs low ptl1

* HT
tab ht low, row chi
logit low i.ht
cc ht low
cs low ht

* UI
tab ui low, row chi
logit low i.ui
cc ui low
cs low ui

* Variavel quantitativa discreta
* FTV - discreta
tab ftv low, row chi
logit low i.ftv
* FTV - categorica
generate float ftv1 = ftv
recode ftv1 0=0 1/3=1 4/6=2
tab ftv1 low, row chi
logit low i.ftv1

* Modelo final
logit low i.ptl1 i.ui i.smoke i.ht i.race
logistic low i.ptl1 i.ui i.smoke i.ht i.race

* Teste de Hosmer-Lemeshow
* Digitar o comando abaixo após o modelo final, o valor de p deve rejeitar Ho (> ou = 0,05)
estat gof, group (10) table


*** BANCO 4 - STA ***

* AGE continuo
tab age sta, row chi
logit sta age
logistic sta age

* AGE categorizado
gen age1 = age
recode age1 16/60=0 61/92=1
tab age1 sta, row chi
logit sta i.age1
logistic sta i.age1

* SEX
tab sex sta, row chi
logit sta i.sex
logistic sta i.sex

* RACE
tab race sta, row chi
logit sta i.race
logistic sta i.race

* SER 
tab ser sta, row chi
logit sta i.ser
logistic sta i.ser

* CAN
tab can sta, row chi
logit sta i.can
logistic sta i.can

* CRN
tab crn sta, row chi
logit sta i.crn
logistic sta i.crn

* INF
tab inf sta, row chi
logit sta i.inf
logistic sta i.inf

* CPR
tab cpr sta, row chi
logit sta i.cpr
logistic sta i.cpr

* SYS - continua
tab sys sta, row chi
logit sta sys
logistic sta sys
* SYS1 - categorica
gen sys1 = sys
recode sys1 36/120=0 121/256=1
tab sys1 sta, row chi
logit sta i.sys1
logistic sta i.sys1

* HRA - continua
tab hra sta, row chi
logit sta hra
logistic sta hra
* HRA - categórica
gen hra1=hra
recode hra1 39/80= 0 81/192=1
tab hra1 sta, row chi
logit sta i.hra1
logistic sta i.hra1

* PRE
tab pre sta, row chi
logit sta i.pre
logistic sta i.pre

* TYP
tab typ sta, row chi
logit sta i.typ
logistic sta i.typ

* FRA
tab fra sta, row chi
logit sta i.fra
logistic sta i.fra

* PO2
tab po2 sta, row chi
logit sta i.po2
logistic sta i.po2

* PH
tab ph sta, row chi
logit sta i.ph
logistic sta i.ph

* PCO
tab pco sta, row chi
logit sta i.pco
logistic sta i.pco

* BIC
tab bic sta, row chi
logit sta i.bic
logistic sta i.bic

* CRE
tab cre sta, row chi
logit sta i.cre
logistic sta i.cre

* LOC
tab loc sta, row chi
logit sta i.loc
logistic sta i.loc
* LOC - recodificado (tinha casela vazia)
gen loc1 = loc
recode loc1 2=1
tab loc1 sta, row chi
logit sta i.loc1
logistic sta i.loc1

* MODELO - OR
logit sta i.loc1
logit sta i.loc1 i.typ
logit sta i.loc1 i.typ i.cpr
logit sta i.loc1 i.typ i.cpr i.ser
logit sta i.loc1 i.typ i.cpr i.inf
logit sta i.loc1 i.typ i.cpr i.crn
logit sta i.loc1 i.typ i.cpr i.cre
logit sta i.loc1 i.typ i.cpr i.age1

* MODELO - OR
logistic sta i.loc1
logistic sta i.loc1 i.typ
logistic sta i.loc1 i.typ i.cpr
logistic sta i.loc1 i.typ i.cpr i.ser
logistic sta i.loc1 i.typ i.cpr i.inf
logistic sta i.loc1 i.typ i.cpr i.crn
logistic sta i.loc1 i.typ i.cpr i.cre
logistic sta i.loc1 i.typ i.cpr i.age1

* Teste de Hosmer-Lemeshow – imediatamente após o modelo
estat gof, group (10) table

* DICAS

* 1. Para odds ratio, pode ser usado um dos dois comando:
logit v.dependente v.independente
* ou
logistic v.dependente v.independente

* 2. Para valores do X2:
tab v.independente v.dependente, row all exact

*3. Para criar variáveis utilizando o commando replace
gen ptl1=.
replace ptl1=0 if ptl<=0