script_descritiva.R

# PARTE I -----------------------------------------------------------------

# upload do conjunto de dados
# Metodo 1: ir ao botao 'Import Dataset'

# Metodo 2: via comandos abaixo
install.packages("readxl") # instalando pacote 'readxl'
library(readxl) # carregando pacote 'readxl'
dados_estud = read_excel("C:\caminho\para\arquivo\dados-exemplo-estudantes.xls")
View(dados_estud)

dados_estud
dados_estud$Sexo
dados_estud[, 3]

attach(dados_estud)

# slide 16
Sexo
table(Sexo)

n = length(Sexo)
n 

table(Sexo)/n
(table(Sexo)/n)*100 # em porcentagem

# slide 17
table(Toler)
table(Toler)/n
prop.table(table(Toler))
(table(Toler)/n)*100 

# slide 18
table(Idade)
table(Idade)/n
prop.table(table(Idade))
(table(Idade)/n)*100 # em porcentagem

cumsum(table(Idade)/n) 
(cumsum(table(Idade)/n))*100 # em porcentagem

# slide 19
tab_peso = hist(Peso, plot = FALSE, right = FALSE, 
                nclass = 6)
tab_peso
tab_peso$counts
tab_peso$counts/n  
cumsum(tab_peso$counts/n)

# slide 22
tab_tv = hist(TV, plot = FALSE, right = FALSE, nclass = 5)
tab_tv

tab_tv = hist(TV, breaks = c(0, 6, 12, 18, 24, 36), 
               plot = FALSE, right = FALSE)
tab_tv

# slide 23
barplot(table(Idade), xlab = "Idade", ylab = "Frequência",
        main = "Gráfico de barras para Idade",
        col = "grey") # frequencia
barplot(table(Idade)/n, xlab = "Idade", ylab = "Frequência Relativa",
        main = "Gráfico de barras para Idade",
        col = "grey") # frequencia # frequencia relativa

# slide 25
hist(Peso, right = FALSE, xlab = "Peso", ylab = "Frequência",
     main = "Histograma de Peso", col = "grey")
hist(Peso, freq = FALSE, right = FALSE)
hist(Peso, right = FALSE, plot = FALSE)


# PARTE II ----------------------------------------------------------------

# slide 31
median(Peso)

# slide 32
quantile(Peso, probs = 0.25) # Q1
quantile(Peso, probs = 0.50) # Q2
quantile(Peso, probs = 0.75, type = 9) # Q3

# slide 36
boxplot(Peso, ylab = "Peso (kg)")
boxplot(Peso, xlab = "Peso (kg)", horizontal = TRUE)

# slide 37
boxplot(Peso ~ Sexo, ylab = "Peso (kg)")
boxplot(Peso ~ Sexo, xlab = "Peso (kg)", horizontal = TRUE)

boxplot(Peso ~ Turma, ylab = "Peso (kg)")

boxplot(Peso ~ Fuma, ylab = "Peso (kg)")

# slide 42
mean(Idade)
summary(Idade)

# slide 43
mean(Peso)
sum(Peso)/n
min(Peso)
max(Peso)
summary(Peso)
max(Peso) - min(Peso)

# slide 55
quantile(Peso, probs = 0.75) - quantile(Peso, probs = 0.25)

# slide 56
var(Peso)
sd(Peso)

sd(Idade)

boxplot(Idade ~ Sexo)


# Parte III ---------------------------------------------------------------

# slide 63
table(Turma, Fuma)

# slide 65
table(Turma, Fuma)/n
prop.table(table(Turma, Fuma))
(table(Turma, Fuma)/n)*100

# slide 66
spineplot(as.factor(Fuma) ~ as.factor(Turma), xlab = "Turma",
          ylab = "Fuma")
spineplot(as.factor(Turma) ~ as.factor(Fuma), xlab = "Fuma",
          ylab = "Turma")
spineplot(as.factor(Fuma) ~ as.factor(Toler), xlab = "Tolerância",
          ylab = "Fuma")

# slide 68
plot(Alt, Peso, xlab = "Altura (m)", ylab = "Peso (kg)")
abline(v = mean(Alt), lty = 2, lwd = 2)
abline(h = mean(Peso), lty = 2, lwd = 2)

# slide 72
cor(Alt, Peso)

# slide 75
aggregate(Peso ~ Sexo, dados_estud, mean)
aggregate(Peso ~ Sexo, dados_estud, sd)
aggregate(Peso ~ Sexo, dados_estud, min)
aggregate(Peso ~ Sexo, dados_estud, max)
aggregate(Peso ~ Sexo, dados_estud, summary)

# slide 76
boxplot(Peso ~ Sexo, ylab = "Peso (kg)", xlab = "Sexo")

# slide 77
cdplot(as.factor(Sexo) ~ Peso, xlab = "Peso (kg)", ylab = "Sexo")
cdplot(as.factor(Fuma) ~ Peso, xlab = "Peso (kg)", ylab = "Fuma")
cdplot(as.factor(Toler) ~ Peso, xlab = "Peso (kg)", ylab = "Tolerância")


# Parte IV ----------------------------------------------------------------


# slide 84
table(Sexo, Fuma)
mosaicplot(table(Sexo, Fuma), color = TRUE, xlab = "Sexo", 
           ylab = "Fuma", main = "")

# slide 85
table(Fuma, Turma, Sexo)
mosaicplot(table(Sexo, Fuma, Turma), color = TRUE, xlab = "Sexo", 
           ylab = "Fuma", main = "")

# slide 86
coplot(Peso ~ Alt | Idade, pch = 20, xlab = "Altura")
coplot(Peso ~ Alt | Idade, panel = panel.smooth, pch = 20,
       lty = 2, xlab = "Altura")

# slide 87
coplot(Peso ~ Alt | Idade * Turma, panel = panel.smooth, pch = 20,
       lty = 2, xlab = "Altura")

# slide 88
coplot(Peso ~ Alt | Idade * Turma, panel = panel.smooth, 
       col = c("blue", "green")[as.factor(Sexo)], lty = 2, 
       pch = 20, xlab = "Altura")
# Sexo: blue para 'F', green para 'M'

# slide 89
coplot(Peso ~ Alt | Sexo, panel = panel.smooth, lty = 2,
       pch = 20, xlab = "Altura")

# slide 90
coplot(Peso ~ Alt | Sexo * Turma, panel = panel.smooth,
       lty = 2, pch = 20, xlab = "Altura")

#- 
coplot(Peso ~ Alt | Sexo * Turma, panel = panel.smooth, lty = 2, 
       pch = 20, col = c("blue", "green")[as.factor(Fuma)])
# Fuma: blue para 'NAO', green para 'SIM'

# slide 91
pairs(~ Alt + Peso + Idade, pch = 20)
pairs(~ Alt + Peso + Idade, pch = 20, panel = panel.smooth, lty = 2)

# slide 92
panel.cor <- function(x, y, digits = 2, prefix = "", cex.cor, ...) {
        usr <- par("usr"); on.exit(par(usr))
        par(usr = c(0, 1, 0, 1))
        r <- abs(cor(x, y))
        txt <- format(c(r, 0.123456789), digits = digits)[1]
        txt <- paste0(prefix, txt)
        if(missing(cex.cor)) cex.cor <- 0.8/strwidth(txt)
        text(0.5, 0.5, txt, cex = cex.cor * r)
}

pairs(~ Alt + Peso + Idade, pch = 20, lower.panel = panel.cor,
      upper.panel = panel.smooth, lty = 2)

# slide 93
panel.hist <- function(x, ...) {
        usr <- par("usr"); on.exit(par(usr))
        par(usr = c(usr[1:2], 0, 1.5) )
        h <- hist(x, plot = FALSE)
        breaks <- h$breaks; nB <- length(breaks)
        y <- h$counts; y <- y/max(y)
        rect(breaks[-nB], 0, breaks[-1], y, col = "grey", ...)
}
pairs(~ Alt + Peso + Idade, pch = 20, lower.panel = panel.cor,
      upper.panel = panel.smooth, diag.panel = panel.hist,
      lty = 2)

#-
pairs(~ Alt + Peso + Idade, pch = 20, panel = panel.smooth,
      lty = 2, col = c("blue", "magenta")[as.factor(Sexo)])
# Sexo: blue para 'F', magenta para 'M'