Logistic discrimination


In [1]:
load("PLS.rda")
load("DP.rda")
suppressMessages(library(caret))
set.seed(201703)

In [2]:
options(warn=-1)
# Logistic Discrimination
pls_ld_s = train(response~., data = pls_train, method = "glm", family = "binomial", trControl = trainControl(method = "LOOCV"))
pls_ld_te = predict(pls_ld_s, data.frame(pls_test_s))
pls_ld_ac = mean(pls_ld_te == golub_test_r)
pls_ld_re = c(LOOCV = pls_ld_s$results$Accuracy, Test = pls_ld_ac)

In [3]:
pls_ld_re


LOOCV
1
Test
0.970588235294118