#############################################################################
#  EXEMPLOS AULA 12 - An. Descritiva de V. Quantitativas e Bidimensionais   #
#############################################################################

# LIMPA O R
rm(list = ls()); cat('\014')

# DIRETORIO
setwd(dirname(rstudioapi::getActiveDocumentContext()$path))

# PACOTES NECESSARIOS
library(texmex)


library(psych)
library(psychTools)


#############################################################################
# EXEMPLO 1 - slide 10 #
########################

# Conjunto de dados da aula 11 - slide 20
dados = read.csv("DadosS20.csv")
head(dados)

## Tabela por classe ##
fabs = table(cut(dados$Mcorporal, breaks = c(seq(49,98,7), 113), right = F))
fabs

Fa = cumsum(fabs); Fa

TAB = cbind(fabs, Fa); TAB
dimnames(TAB) = list(c("49|-57", "57|-65", "65|-73", "73|-81", 
                       "81|-89", "89|-97", "97|-105", "105|-113"),
                     c("f", "F"))
TAB

Media = mean(dados$Mcorporal); Media


#############################################################################
# EXEMPLO 2 - slide 38 #
########################

Cdados = c(2.1, 0.6, 0.5, 0.8, 1.1, 0.7, 1.3, 0.2, 
           2.5, 0.3, 1.6, 0.4, 2.8, 4.4, 1.4)

### Medidas de dispersao ###

# Amplitude
A = diff(range(Cdados)); A

# Variancia
s2 = var(Cdados); round(s2, 2)

# Desvio padrao
s = sd(Cdados); round(s, 2)


#############################################################################
# EXEMPLO 3 - slide 79 #
########################

data(winter); head(winter)

# O3   Daily maximum ozone in parts per billion.
# NO2  Daily maximum NO2 in parts per billion.
# NO   Daily maximum NO in parts per billion.
# SO2  Daily maximum SO2 in parts per billion.
# PM10 Daily maximum PM10 in micrograms/metre^3

describe(winter)
summary(winter)


# Grafico
pairs.panels(winter)