*use "C:\Documents and Settings\Verlato\Desktop\didattic\OtherCourses\Dottorato\Corso Stata 2011\FreqDrug.dta" 
*use "FreqDrug2.dta" 
*use "C:\Users\Giuseppe\Desktop\didattica\Stata6\Stata 2021\1_dati_FreqDrug.dta" 

use 1_dati_FreqDrug.dta 

sum _all
list

* STATISTICA DESCRITTIVA DELLA VARIABILI QUALITATIVA *.
label define farmaci 1 "diuretico" 2 "beta-bloccante"
label values drug farmaci
tab drug

* STATISTICA DESCRITTIVA DELLA VARIABILE QUANTITATIVA *.
sum FreqCard,detail
graph box FreqCard
hist FreqCard,width(4) norm
*sktest FreqCard

* CONFRONTO FRA I DUE GRUPPI DI TRATTAMENTO: STATISTICA DESCRITTIVA *
tabstat FreqCard,by(drug) stat(n mean sd)
tabstat FreqCard,by(drug) stat(min p25 p50 p75 max)
graph box FreqCard,by(drug)
hist FreqCard,width(4) norm by(drug)
tabstat FreqCard,by(drug) stat(n mean sd skewness kurtosis)


* CONFRONTO FRA I DUE GRUPPI DI TRATTAMENTO: STATISTICA INFERENZIALE *

* CONTROLLO DELLE ASSUNZIONI DEI TEST PARAMETRICI *
* Test per l'omogeneità della varianza *
sdtest FreqCard,by(drug)

* Test per la normalità delle distribuzioni *
*LOW POWER: swilk can be used with 4<=n<=2000 observations. sfrancia can be used with 10<=n<=5000 observations *
*swilk FreqCard
sort drug
by drug:swilk FreqCard

* Test t di Student (varianza omogenea) *
ttest FreqCard, by(drug)

* Test t di Student (varianza disomogenea) *
*ttest FreqCard, by(drug) unequal

* Test non-parametrico: test di Wilcoxon-Mann-Whitney *
ranksum FreqCard, by(drug)