# HG changeset patch
# User deepakjadmin
# Date 1483635813 18000
# Node ID d56ce97c4985f08b642536b00c99348233acdb96
# Parent  b84589b7c014ba1e9dadd02bd7d625da353c4df5
Uploaded

diff -r b84589b7c014 -r d56ce97c4985 feature_selection.R
--- a/feature_selection.R	Tue Jan 03 02:30:32 2017 -0500
+++ b/feature_selection.R	Thu Jan 05 12:03:33 2017 -0500
@@ -9,7 +9,9 @@
 arg7 <- args[7]
 arg8 <- args[8]
 arg9 <- args[9]
+arg10 <- args[10]
 library(caret)
+library(doMC)
 load(arg1)
 
 #RAWDATA <- dataX
@@ -133,7 +135,7 @@
 
 
 
- 
+if (as.numeric(arg10) == 1){ 
 Profile <- rfe(dx, dy,sizes = subsets,rfeControl = ctrl1)
 
 pred11 <- predictors(Profile)
@@ -144,4 +146,19 @@
 save(dataX,dataY,file=arg3)
 rm(dataX)
 rm(dataY)
+} else if (as.numeric(arg10) > 1){
+registerDoMC(cores = as.numeric(arg10))
 
+Profile <- rfe(dx, dy,sizes = subsets,rfeControl = ctrl1)
+
+pred11 <- predictors(Profile)
+save(Profile,file=arg2)
+dataX <- rawData[,pred11]
+dataY <- rawData$outcome
+
+save(dataX,dataY,file=arg3)
+rm(dataX)
+rm(dataY)
+} else { stop("something went wrong. please see the parameters")}  
+
+