Questão 3 trabalho - multivariada

library("ggraph")

## Loading required package: ggplot2

library("factoextra")

## Welcome! Want to learn more? See two factoextra-related books at https://goo.gl/ve3WBa

alimentos = read.table("alimentos.txt",header=T)
alimentos_padronizados = scale(alimentos[,-1])
dados <- cbind(alimentos$pais, alimentos_padronizados)
dados = as.data.frame(dados)
names(dados)[names(dados) == "V1"] <- "pais"
row.names(dados) = dados$pais
dados$aveg = as.numeric(dados$aveg)
dados$aani = as.numeric(dados$aani)
dados$fleg = as.numeric(dados$fleg)
dados$carn = as.numeric(dados$carn)
dados$pesc = as.numeric(dados$pesc)
dados$lech = as.numeric(dados$lech)
dados$card = as.numeric(dados$card)
dados$cancer = as.numeric(dados$cancer)
dados$lmas = as.numeric(dados$lmas)
dados$lfem = as.numeric(dados$lfem)

Letra A

matriz_dissimi  <- dist(dados[, c("lech", "fleg")], method = "euclidean")
complete_linkage <- hclust(matriz_dissimi, method = "complete")
#complete_linkage$height <- complete_linkage$height - min(complete_linkage$height)
plot(complete_linkage, cex = 0.8, main = "Dendrograma de Consumos")

Letra B

matriz_dissimi2 <- dist(dados[, -1], method = "euclidean")
complete_linkage2 <- hclust(matriz_dissimi2, method = "complete")
#complete_linkage2$height <- complete_linkage2$height - min(complete_linkage2$height)
plot(complete_linkage2, cex = 0.8, main = "Dendrograma de Consumo de Alimentos")

Letra c

average_linkage <- hclust(matriz_dissimi2, method = "average")
#average_linkage$height <- average_linkage$height - min(average_linkage$height)
plot(average_linkage, cex = 0.8, main = "Dendrograma - Average Linkage")

Letra d

clustersFL_2 <- cutree(complete_linkage, k = 2)
dados$clustersFL_2 = clustersFL_2
clustersFL_3 <- cutree(complete_linkage, k = 3)
dados$clustersFL_3 = clustersFL_3

grupos_b_2 <- cutree(complete_linkage2, k = 2)
grupos_b_3 <- cutree(complete_linkage2, k = 3)

Letra e

alim_2 = dados[, c("lech","fleg")]; alim_2

fviz_cluster(list(data = alim_2[, c("lech", "fleg")], cluster = clustersFL_2), 
             ellipse.type = "convex", 
             main = "Clusters de acordo com consumo de lech e fleg para K=2")

fviz_cluster(list(data = alim_2[, c("lech", "fleg")], cluster = clustersFL_3), 
             ellipse.type = "convex", 
             main = "Clusters de acordo com consumo de lech e fleg para K=3")

fviz_cluster(list(data = dados[,-1], cluster = grupos_b_2), 
             ellipse.type = "convex", 
             main = "Clusters de acordo com consumo de alimentos para K=2")

fviz_cluster(list(data = dados[,-1], cluster = grupos_b_3), 
             ellipse.type = "convex", 
             main = "Clusters de acordo com consumo de alimentos para K=3")