Advertisement
Guest User

Failing code

a guest
Dec 21st, 2014
143
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
R 1.45 KB | None | 0 0
  1. library(e1071)
  2. library(kernlab)
  3.  
  4. errors <- c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0)
  5.  
  6. for(i in seq(1,50)) {
  7.  
  8. enquesta<-read.csv("APA-TEMA7-P1.csv", header=TRUE, sep=",",dec=",")
  9. enquesta["Salary.3"]<-rep(nrow(enquesta), 0)
  10. enquesta$Salary.3[enquesta$Salary == "35k"] <- 2
  11. enquesta$Salary.3[enquesta$Salary %in% c("45k", ">45k")] <- 3
  12. enquesta$Salary.3[is.na(enquesta$Salary.3)] <- 1
  13. enquesta$Salary.3 <- as.factor(enquesta$Salary.3)
  14. L <- 5
  15. K <- 5
  16. tr <- sample(rep(1:nrow(enquesta)),round(.8*nrow(enquesta)))
  17. dades<-enquesta[tr,]
  18. test<-enquesta[-tr,]
  19. vectorC<-c(.25,.5,.75,1)
  20. errorsC <- rep(0, length(vectorC))
  21. for(cross in seq(1,L)) {
  22.   folds <- sample(rep(1:K, length=nrow(dades)), nrow(dades), replace=FALSE)
  23.   for(k in seq(1:K)){
  24.     training <- dades[folds != k,]
  25.     validation <- dades[folds == k,]
  26.     for(i in seq(1,length(vectorC))){
  27.       training.svm<-ksvm(training$Salary.3~., training, C=vectorC[i])
  28.       p <- predict(training.svm, newdata=validation)
  29.       error <- 1 - (sum(diag(table(p, validation$Salary.3)))/length(p))
  30.       errorsC[i] <- errorsC[i] + error
  31.     }
  32.   }
  33. }
  34. errorsC <- errorsC/(L*K)
  35. bestC.index <- which(errorsC == min(errorsC))
  36. ksvm.result<-ksvm(dades$Salary.3~., dades, C=vectorC[bestC.index])
  37. p.result <- predict(ksvm.result, newdata=test)
  38. (error.result <- 1 - (sum(diag(table(p, validation$Salary.3)))/length(p)))
  39.  
  40. errors[i] <- error.result
  41.  
  42. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement