In [1]:

    

load("PLS.rda")
load("DP.rda")
suppressMessages(library(caret))
set.seed(201703)


    

In [2]:

    

options(warn=-1)
# QDA
pls_qda_s = train(response~., data = pls_train, method = "qda", trControl = trainControl(method = "LOOCV"))
pls_qda_te = predict(pls_qda_s, data.frame(pls_test_s))
pls_qda_ac = mean(pls_qda_te == golub_test_r)
pls_qda_re = c(LOOCV = pls_qda_s$results$Accuracy, Test = pls_qda_ac)


    

In [3]:

    

pls_qda_re


    

    






	
LOOCV
		
1
	
Test
		
0.970588235294118

In [ ]: