DÚ 9

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.

LDA klasifikácia

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)

QDA klasifikácia

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)