## @knitr leggi_librerie
# carico i package --------------------------------------------------------
library(tidyverse)
library(janitor)
library(corrr)
library(FactoMineR)
library(factoextra)
library(ggthemes)

## @knitr prepara_dati
# preparo la tabella dati -------------------------------------------------
iris_data <- clean_names(iris)

## @knitr matrice_correlazione
# matrice di correlazione -------------------------------------------------
correlate(iris_data, diagonal = 0)

## @knitr esegui_pca
# PCA ---------------------------------------------------------------------
iris_pca <- PCA(select(iris_data, -species), graph = FALSE)

## @knitr grafico_autovalori
# andamento autovalori ----------------------------------------------------
fviz_eig(iris_pca) +
  theme_light() +
  labs(x = "Dimensioni",
       y = "Percentuale di varianza spiegata",
       title = "Screeplot",
       subtitle = "Andamento della variabilità spiegata dalle dimensioni di sintesi",
       caption = "[Grafico ottenuto usando FactoMiner e factoextra]")

## @knitr cerchio_correlazioni
# cerchio delle correlazioni ----------------------------------------------
fviz_pca_var(iris_pca)

## @knitr grafico_unita
# grafico delle unità -----------------------------------------------------
fviz_pca_ind(iris_pca, habillage = iris_data$species, label = "none")
