### ---- Representação de variáveis qualitativas

### ---- Exemplo

x<-c("SE","SE","NE","NE","S","CO","S","NE","SE","SE","SE","NE","SE","SE","NE","SE","SE","SE","NE","NE","S","SE","SE","NE","SE","SE","S","S","SE","SE","SE","NE","SE","SE","SE","N","S","S","NE","NE","S","CO","S","N","CO","NE","S","SE","S","CO","NE","NE","SE","NE","SE","NE","SE","NE","SE","SE","N","S","SE","S","CO","S","SE","S","SE","SE","SE","NE","NE","SE","NE","S","S","S","NE","SE")

n = length(x)


### ---- Tabela de frequências

table(x) / n

table(x) / n * 100

sort(table(x) / n) * 100

freqa =  table(x)

freqr =  table(x) / n * 100

sort(freqr, decreasing = TRUE)


### ---- Restringindo valores

freqr[which(freqr > 10)]

freqr[freqr > 10]

sum(freqr[c("CO", "SE")])

which.max(freqr)   

sort(table(x), decreasing = TRUE) / n * 100

names(freqr)

freqr[2]

freqr[setdiff(names(freqr), c("N", "S"))]

freqr[which.max(freqr)]



#### Gráficos de barras

freqa =  table(x)

freqr =  sort(table(x)) / n * 100

barplot(freqa)



#### Gráficos de barras - Acrescentando valores das frequências

posb = barplot(freqr, xlab = "Região de origem",  ylab = "Percentagem")

text(posb, freqr / 2, freqr)

box()


#### Gráficos de barras - formato horizontal

posb = barplot(freqr, ylab = "Região de origem",  xlab = "Percentagem", horiz = TRUE)

text(freqr / 2, posb, freqr)

box()

posb = barplot(freqr, xlab = "Região de origem",  ylab = "Percentagem", angle = 45, density = 15, col = "black")

text(posb, freqr / 2, freqr, col = “red”)

box()


#### Gráficos de barras - cores

posb = barplot(freqr, ylab = "Região de origem",  xlab = "Percentagem", horiz = TRUE, col = "cyan")

text(freqr / 2, posb, freqr)

box()

mcores = rainbow(length(freqr))

posb = barplot(freqr, xlab = "Região de origem",  ylab = "Percentagem", col = mcores)

text(posb, freqr / 2, freqr)

box()



### Gráfico de Pareto

library(qcc)
pareto.chart(freqa, main = "", xlab = "Região de origem", ylab = "Frequencia")



### Gráfico de Pareto - cores

pareto.chart(freqa, main = "", xlab = "Região de origem", ylab = "Frequencia",  col = “lightgray”)




### Gráfico de setores

pie(freqr)



### Gráfico de setores - cores

pie(freqr, col=mcores)



### Gráficos de setores - com labels

pie(freqr, col = mcores, labels = paste(names(freqr), freqr))


### Gráficos de setores - em tons de cinza






### Dados FUVEST 2009

dados <- c(
rep("Pais",1476), 
rep("trabalho + família",974), 
rep("trabalho para família",108), 
rep("por conta própria",193), 
rep("bolsa/crédito",422), 
rep("bolsa/crédito + tr + fam",653), 
rep("outros",74))

library(qcc)

pareto.chart(table(dados))