clear
set more off
set mem 400m

*World Values Survey
use "C:\DCP854B\wvs1981_2008.dta", clear

*Mantendo observações em que índice sobre
*"valores racionais: tradicional/secular" (tradrat5)
*são diferentes de missing
keep if tradrat5!=.

*Salvando banco
saveold "C:\DCP854B\tradrat5.dta", replace

*Mas vamos abrir o banco original
use "C:\DCP854B\wvs1981_2008.dta", clear

*Explorando índice de "valores racionais: tradicional/secular"

*Tamanho do banco
count

*Histograma
hist tradrat5, frequency
graph export "C:\DCP854B\hist_tradrat5.wmf", replace

*Diagrama de dispersão entre índice (tradrat5) e idade (x003)
scatter tradrat5 x003

*Estatísticas descritivas
sum tradrat5
sum tradrat5, d

*Entendendo o peso
*Peso "s017a" é o homogeneizado
browse s002 s003 s003a s017 s017a

*Estatísticas descritivas com peso
sum tradrat5 [aweight=s017a], d

*Para calcular moda, devemos instalar módulo "mmodes"
do "C:\DCP854B\modulos\mmodes.ado"
mmodes tradrat5 if tradrat5!=.

*Ponto médio = (MAX + MIN) / 2
di (3.844721 + (- 0.9999222)) / 2

*Amplitude = MAX - MIN
di (3.844721 - (- 0.9999222))

*Regra empírica da amplitude
*Desvio padrão = amplitude / 4
di (3.844721 - (- 0.9999222)) / 4 

*Valor real do desvio padrão
sum tradrat5, d

*Regra empírica da amplitude
*Valor amostral mínimo = média - 2*desvio padrão
di 0.2250021 - (2*1.2111608)

*Valor amostral máximo = média + 2*desvio padrão
di 0.2250021 + (2*1.2111608)

*Coeficiente de variação
*Idade
sum x003
di 16.00398 / 39.98154 * 100

*Índice de "valores racionais: tradicional/secular"
sum tradrat5
di 0.8854222 / 0.2250021 * 100

*Para calcular escores "z", devemos instalar módulo "zscore"
do "C:\DCP854B\modulos\zscore.ado"
zscore (tradrat5)
sum z_tradrat5, d

*Valores comuns
count if (z_tradrat5>=-2) & (z_tradrat5<=2)

*Valores não-usuais
count if (z_tradrat5>2) | (z_tradrat5<-2)

*Estatísticas definidas por quartis e percentis
sum tradrat5 [aweight=s017a], d

*Intervalo interquartil (IIQ) = Q3–Q1
di 0.7238476 - (-0.4261761)

*Intervalo semi-interquartil = (Q3–Q1) / 2
di (0.7238476 - (-0.4261761)) /2

*Ponto médio dos quartis= (Q3+Q1) / 2
di (0.7238476 + (-0.4261761)) /2

*Intervalo percentílico 10–90 = P90–P10
di 1.49589 - (-0.8049945)

*Diagrama de caixa ("boxplot"), excluindo "outliers"
graph hbox tradrat5, noout

*Diagrama de caixa, indicando "outliers" com pontos
graph hbox tradrat5

*Diagrama de caixa de "tradrat5" por sexo (x001)
graph hbox tradrat5, over(x001)

*Diagrama de caixa de "tradrat5" por estado civil (x007)
graph hbox tradrat5, over(x007)