deletions | additions
diff --git a/Appendix.tex b/Appendix.tex
index 5a7e9a6..4bd42a6 100644
--- a/Appendix.tex
+++ b/Appendix.tex
...
library(xlsx)
file <- '20140410-Set-Progress.xlsx'
sheets <- c('Broad', 'JHU', 'PNNL', 'VUMC-OVelos', 'VUMC-Elite')
data <- lapply(sheets, function(x)
{read.xlsx2(file, {
read.xlsx2(file,
sheetName=x,
colClasses=c('character', 'character', rep('numeric', 44)),
stringsAsFactors=FALSE)}) stringsAsFactors=FALSE)
})
data <- mapply(function(x, sheet) {
x$lab <- sheet;
x[[2]] <- as.POSIXct(x[[2]], format='%Y-%m-%dT%H:%M:%SZ');
...
trainingSet <- which(data[[4]]$StartTimeStamp > as.POSIXct('2013/02/25') &
data[[4]]$StartTimeStamp < as.POSIXct('2013/04/15'))
testSet <- unlist(standardInd)[!unlist(standardInd)
\%in\% %in% trainingSet]
### PCA
pcaChoice <- pca(dataMatrix[[4]][trainingSet,], 'nipals', 5, scale='uv', center=T, cv='q2')
...
### One-class SVM
SVM <- ksvm(dataMatrix[[4]][trainingSet,], scaled=T, type='one-svc', cross=10, nu=0.1, kernel='vanilladot')
out <- predict(SVM, dataMatrix[[4]][testSet,])
testRes$outlierSVM <- testRes$time
\%in\% %in% data[[4]]$StartTimeStamp[testSet][out]
pdf('svmPCA.pdf', width=7, height=9)
qplot(time, value, data=testRes, geom='point', colour=outlierSVM) +
facet_grid(variable~., scales = 'free_y')+theme_bw() +