In [1]:

    

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


    

In [2]:

    

options(warn=-1)
# QDA
pca_qda_s = train(response~., data = pca_train, method = "qda", trControl = trainControl(method = "LOOCV"))
pca_qda_te = predict(pca_qda_s, data.frame(pca_test_s))
pca_qda_ac = mean(pca_qda_te == golub_test_r)
pca_qda_re = c(LOOCV = pca_qda_s$results$Accuracy, Test = pca_qda_ac)


    

In [3]:

    

pca_qda_re


    

    






	
LOOCV
		
1
	
Test
		
0.941176470588235