V tejto domácej úlohe opäť pracujeme s dátami carc
obsahujúcimi 66 pozorovaní a 13 premenných. Budeme rozlišovať sídla automobilky podľa toho, či ležia v Európe alebo nie.
Obrázok nižšie zobrazuje, ako bola klasifikované odozva. X-ová os zobrazuje hodnoty koeficientov 1. smeru z LDA modelu. Všimnime si, že dáta sa rozlišujú približne v okolí 2.
data_new = data
data_new$C = as.numeric(data_new$C)
data_new <- transform(data_new,
C = factor(C, levels = 1:3,
labels = c("nonEurope", "nonEurope", "Europe")))
lda.model <- lda(C ~ P + M + R78 + R77 + H + R + Tr + W + L + T + D + G, data = data_new)
pred <- predict(lda.model)
par(mar=c(4,4,4,4))
ldahist(data = pred$x[,1], g=data_new$C)
Na obrázku nižšie môžeme vidieť pozorovania všetkých dvojích uvažovaných premenných.
par(mar = c(1, 1, 1, 1))
partimat(C ~ P + M + H + R + Tr,
method = "lda", data = data_new)
Nasledujúci obrázok ponúka iný spôsob delenia ako pri LDA klasifikácii.
partimat(C ~ P + M + H + R + Tr,
method = "qda", data = data_new)