R/etc/dosvm.r

   1 #!/usr/bin/env littler
   2
   3 #install.packages('e1071')
   4 require('e1071')
   5 #mdat <- matrix(c(21,.1,.11,.12,.14, 23,.15,.16,.17,.15, 27,.2,.19,.18,.3, 32,.3,.33,.32,.34, 36,.4,.41,.42,.4, 42,.5,.51,.52,.54), 5, 5, TRUE, dimnames = list(c(),c("Age", "Pos1", "Pos2","Pos3","Pos4")))
   6 mdat <- read.delim('dosvm.txt',comment.char='#')
   7 x <- subset(mdat, select = -Age)
   8 y <- mdat[,1]
   9 model <- svm(Age ~ ., data = mdat, type='eps-regression', gamma=0.1, cost=2, epsilon = 0.1)
  10 model2 <- svm(Age ~ a11044875 + a11044880, data = mdat, type='eps-regression', gamma=0.1, cost=2, epsilon = 0.1)
  11 model3 <- svm(Age ~ a11044875 + a11044880 + b11044877 + b11044888 + b11044894, data = mdat, type='eps-regression', gamma=0.1, cost=2, epsilon = 0.1)
  12 pred_result <- predict(model, x)
  13 pred_result2 <- predict(model2, x)
  14 pred_result3 <- predict(model3, x)
  15
  16 #print(table(pred_result,y))
  17
  18 print(model)
  19 print(y)
  20 print(pred_result)
  21 print(pred_result2)
  22 print(pred_result3)
  23
  24 #plot(model, mdat, Pos1 ~ Pos2,color.palette = terrain.colors)