# EXPLORAČNÍ FA

installed.packages("readxl")        # Instalace knihovny načítající Excel
library(readxl)                     # Aktivace knihovny

dat <- as.data.frame(read_excel("FA.15.7.20.xlsx",range = "A1:Q101"))

names(dat)

pca <- princomp(dat[,5:17], cor=T)         # Metoda hlavních komponent, velmi základní, pro všechny řádky a sloupce 5-17 datového souboru "dat"
plot(pca)                           # Graf rozptylů hlavních komponent
summary(pca)
loadings(pca)

install.packages("GPArotation")   # Instalace knihovny pro rotace
library(GPArotation)              # Aktivace knihovny

pca2 <- principal(dat[,5:17], nfactors=3, rotate="varimax")

??principal    # Funkce principal havarovala - resp. program ji nenašel, proto ji hledáme, v Helpu je nápověda na knihovnu "psych"

install.packages("psych")
library(psych)

pca2 <- principal(dat[,5:17], nfactors=3, rotate="varimax")
pca2

?psych

# KONFIRMAČNÍ FA

data()      # uKÁŽE DATA VE VŠECH DOSTUPNÝCH (AKTIVNÍCH) KNIHOVNÁCH
data(package = "psych")   # Ukáže data pouze pro knihovnu "psych"

bfi
bfi.dat <- as.data.frame(bfi)   # Načtení dat

install.packages("lavaan")      # Knihovna pro konfirmační FA (LAtent VAriable ANalysis)
library(lavaan)

bfi.model <- '
agree =~ A1 + A2 + A3 + A4 + A5
consc =~ C1 + C2 + C3 + C4 + C5
extra =~ E1 + E2 + E3 + E4 + E5
neuro =~ N1 + N2 + N3 + N4 + N5
opene =~ O1 + O2 + O3 + O4 + O5
'

bfi.model     # Kontrola modelu

describe(bfi.dat)     # Popisné statistiky dat

bfi.cfa <- cfa(bfi.model, data=bfi.dat[,1:25])    # Spuštění konfirmační analýzy

bfi.cfa     # Velmi základní pohled na objekt - výsledek analýzy funkcí "cfa"
summary(bfi.cfa, fit.measures = T)      # Jeden ze způsobů prohlídky výsledků      

summary(bfi.cfa, fit.measures = T, standardized = T)