library(MASS) # library for the 'stepAIC()' function

TOMCAT <- read.csv("cohesionMetricsFinal_Tomcat.csv") # reading the data (could be also the tabular file)
TOMCAT2 <- na.omit(TOMCAT) # omiting N/A measures
linear.TOMCAT <- lm(CC ~ LCOM2+LCOM3+LCOM5+TCC+LCbC+LCOM4, data=TOMCAT2) # running the multivariate linear regression
step.linear.TOMCAT <- stepAIC(linear.TOMCAT, direction="both") # running the stepwise function for the linear regression
step.linear.TOMCAR$anova # display results 

RHINO <- read.csv("cohesionMetricsFinal_Rhino.csv")
RHINO2 <- na.omit(RHINO)
linear.RHINO <- lm(CC ~ LCOM2+LCOM3+LCOM5+TCC+LCbC+LCOM4, data=RHINO2)
step.linear.RHINO <- stepAIC(linear.RHINO, direction="both")
step.linear.RHINO$anova

JFREECHART <- read.csv("cohesionMetricsFinal_JFreeChart.csv")
JFREECHART2 <- na.omit(JFREECHART)
linear.JFREECHART <- lm(CC ~ LCOM2+LCOM3+LCOM5+TCC+LCbC+LCOM4, data=JFREECHART2)
step.linear.JFREECHART <- stepAIC(linear.JFREECHART, direction="both")
step.linear.JFREECHART$anova

JEDIT <- read.csv("cohesionMetricsFinal_jEdit.csv")
JEDIT2 <- na.omit(JEDIT)
linear.JEDIT <- lm(CC ~ LCOM2+LCOM3+LCOM5+TCC+LCbC+LCOM4, data=JEDIT2)
step.linear.JEDIT <- stepAIC(linear.JEDIT, direction="both")
step.linear.JEDIT$anova

FINDBUGS <- read.csv("cohesionMetricsFinal_Findbugs.csv")
FINDBUGS2 <- na.omit(FINDBUGS)
linear.FINDBUGS <- lm(CC ~ LCOM2+LCOM3+LCOM5+TCC+LCbC+LCOM4, data=FINDBUGS2)
step.linear.FINDBUGS <- stepAIC(linear.FINDBUGS, direction="both")
step.linear.FINDBUGS$anova

FREECOL <- read.csv("cohesionMetricsFinal_Freecol.csv")
FREECOL2 <- na.omit(FREECOL)
linear.FREECOL <- lm(CC ~ LCOM2+LCOM3+LCOM5+TCC+LCbC+LCOM4, data=FREECOL2)
step.linear.FREECOL <- stepAIC(linear.FREECOL, direction="both")
step.linear.FREECOL$anova