Title

library(ggplot2)
setwd("/home/rusty/Documents/INRIA/test")

dataset.td.lu.16.I<-read.table("dataframe.td.lu.16.I", header=TRUE)
dataset.td.lu.16.II<-read.table("dataframe.td.lu.16.I", header=TRUE)
dataset.td.lu.16.III<-read.table("dataframe.td.lu.16.I", header=TRUE)
dataset.td.lu.16.IV<-read.table("dataframe.td.lu.16.I", header=TRUE)
dataset.td.lu.16.V<-read.table("dataframe.td.lu.16.I", header=TRUE)
dataset.ti.lu.16.I<-read.table("dataframe.ti.lu.16.I", header=TRUE)
dataset.ti.lu.16.II<-read.table("dataframe.ti.lu.16.II", header=TRUE)
dataset.ti.lu.16.III<-read.table("dataframe.ti.lu.16.III", header=TRUE)
dataset.ti.lu.16.IV<-read.table("dataframe.ti.lu.16.IV", header=TRUE)
dataset.ti.lu.16.V<-read.table("dataframe.ti.lu.16.V", header=TRUE)

dI_I<-data.frame(xlabel=paste(dataset.td.lu.16.I$eventID, dataset.td.lu.16.I$subblock, sep="_"),
                 avgTD=dataset.td.lu.16.I$avg, minTD=dataset.td.lu.16.I$min, maxTD=dataset.td.lu.16.I$max,
                 avgTI=dataset.ti.lu.16.I$avg, minTI=dataset.ti.lu.16.I$min, maxTI=dataset.ti.lu.16.I$max,
                 ratioavg=dataset.td.lu.16.I$avg/dataset.ti.lu.16.I$avg,
                 ratiomin=dataset.td.lu.16.I$min/dataset.ti.lu.16.I$min,
                 ratiomax=dataset.td.lu.16.I$max/dataset.ti.lu.16.I$max)

dI_II<-data.frame(xlabel=paste(dataset.td.lu.16.I$eventID, dataset.td.lu.16.I$subblock, sep="_"),
                  avgTD=dataset.td.lu.16.I$avg, minTD=dataset.td.lu.16.I$min, maxTD=dataset.td.lu.16.I$max,
                  avgTI=dataset.ti.lu.16.II$avg, minTI=dataset.ti.lu.16.II$min, maxTI=dataset.ti.lu.16.II$max,
                  ratioavg=dataset.td.lu.16.I$avg/dataset.ti.lu.16.II$avg,
                  ratiomin=dataset.td.lu.16.I$min/dataset.ti.lu.16.II$min,
                  ratiomax=dataset.td.lu.16.I$max/dataset.ti.lu.16.II$max)

dI_III<-data.frame(xlabel=paste(dataset.td.lu.16.I$eventID, dataset.td.lu.16.I$subblock, sep="_"),
                   avgTD=dataset.td.lu.16.I$avg, minTD=dataset.td.lu.16.I$min, maxTD=dataset.td.lu.16.I$max,
                   avgTI=dataset.ti.lu.16.III$avg, minTI=dataset.ti.lu.16.III$min, maxTI=dataset.ti.lu.16.III$max,
                   ratioavg=dataset.td.lu.16.I$avg/dataset.ti.lu.16.III$avg,
                   ratiomin=dataset.td.lu.16.I$min/dataset.ti.lu.16.III$min,
                   ratiomax=dataset.td.lu.16.I$max/dataset.ti.lu.16.III$max)

dI_IV<-data.frame(xlabel=paste(dataset.td.lu.16.I$eventID, dataset.td.lu.16.I$subblock, sep="_"),
                  avgTD=dataset.td.lu.16.I$avg, minTD=dataset.td.lu.16.I$min, maxTD=dataset.td.lu.16.I$max,
                  avgTI=dataset.ti.lu.16.IV$avg, minTI=dataset.ti.lu.16.IV$min, maxTI=dataset.ti.lu.16.IV$max,
                  ratioavg=dataset.td.lu.16.I$avg/dataset.ti.lu.16.IV$avg,
                  ratiomin=dataset.td.lu.16.I$min/dataset.ti.lu.16.IV$min,
                  ratiomax=dataset.td.lu.16.I$max/dataset.ti.lu.16.IV$max)

dI_V<-data.frame(xlabel=paste(dataset.td.lu.16.I$eventID, dataset.td.lu.16.I$subblock, sep="_"),
                 avgTD=dataset.td.lu.16.I$avg, minTD=dataset.td.lu.16.I$min, maxTD=dataset.td.lu.16.I$max,
                 avgTI=dataset.ti.lu.16.V$avg, minTI=dataset.ti.lu.16.V$min, maxTI=dataset.ti.lu.16.V$max,
                 ratioavg=dataset.td.lu.16.I$avg/dataset.ti.lu.16.V$avg,
                 ratiomin=dataset.td.lu.16.I$min/dataset.ti.lu.16.V$min,
                 ratiomax=dataset.td.lu.16.I$max/dataset.ti.lu.16.V$max)

dII_I<-data.frame(xlabel=paste(dataset.td.lu.16.II$eventID, dataset.td.lu.16.II$subblock, sep="_"),
                  avgTD=dataset.td.lu.16.II$avg, minTD=dataset.td.lu.16.II$min, maxTD=dataset.td.lu.16.II$max,
                  avgTI=dataset.ti.lu.16.I$avg, minTI=dataset.ti.lu.16.I$min, maxTI=dataset.ti.lu.16.I$max,
                  ratioavg=dataset.td.lu.16.II$avg/dataset.ti.lu.16.I$avg,
                  ratiomin=dataset.td.lu.16.II$min/dataset.ti.lu.16.I$min,
                  ratiomax=dataset.td.lu.16.II$max/dataset.ti.lu.16.I$max)

dII_II<-data.frame(xlabel=paste(dataset.td.lu.16.II$eventID, dataset.td.lu.16.II$subblock, sep="_"),
                   avgTD=dataset.td.lu.16.II$avg, minTD=dataset.td.lu.16.II$min, maxTD=dataset.td.lu.16.II$max,
                   avgTI=dataset.ti.lu.16.II$avg, minTI=dataset.ti.lu.16.II$min, maxTI=dataset.ti.lu.16.II$max,
                   ratioavg=dataset.td.lu.16.II$avg/dataset.ti.lu.16.II$avg,
                   ratiomin=dataset.td.lu.16.II$min/dataset.ti.lu.16.II$min,
                   ratiomax=dataset.td.lu.16.II$max/dataset.ti.lu.16.II$max)

dII_III<-data.frame(xlabel=paste(dataset.td.lu.16.II$eventID, dataset.td.lu.16.II$subblock, sep="_"),
                    avgTD=dataset.td.lu.16.II$avg, minTD=dataset.td.lu.16.II$min, maxTD=dataset.td.lu.16.II$max,
                    avgTI=dataset.ti.lu.16.III$avg, minTI=dataset.ti.lu.16.III$min, maxTI=dataset.ti.lu.16.III$max,
                    ratioavg=dataset.td.lu.16.II$avg/dataset.ti.lu.16.III$avg,
                    ratiomin=dataset.td.lu.16.II$min/dataset.ti.lu.16.III$min,
                    ratiomax=dataset.td.lu.16.II$max/dataset.ti.lu.16.III$max)

dII_IV<-data.frame(xlabel=paste(dataset.td.lu.16.II$eventID, dataset.td.lu.16.II$subblock, sep="_"),
                   avgTD=dataset.td.lu.16.II$avg, minTD=dataset.td.lu.16.II$min, maxTD=dataset.td.lu.16.II$max,
                   avgTI=dataset.ti.lu.16.IV$avg, minTI=dataset.ti.lu.16.IV$min, maxTI=dataset.ti.lu.16.IV$max,
                   ratioavg=dataset.td.lu.16.II$avg/dataset.ti.lu.16.IV$avg,
                   ratiomin=dataset.td.lu.16.II$min/dataset.ti.lu.16.IV$min,
                   ratiomax=dataset.td.lu.16.II$max/dataset.ti.lu.16.IV$max)

dII_V<-data.frame(xlabel=paste(dataset.td.lu.16.II$eventID, dataset.td.lu.16.II$subblock, sep="_"),
                  avgTD=dataset.td.lu.16.II$avg, minTD=dataset.td.lu.16.II$min, maxTD=dataset.td.lu.16.II$max,
                  avgTI=dataset.ti.lu.16.V$avg, minTI=dataset.ti.lu.16.V$min, maxTI=dataset.ti.lu.16.V$max,
                  ratioavg=dataset.td.lu.16.II$avg/dataset.ti.lu.16.V$avg,
                  ratiomin=dataset.td.lu.16.II$min/dataset.ti.lu.16.V$min,
                  ratiomax=dataset.td.lu.16.II$max/dataset.ti.lu.16.V$max)

dIII_I<-data.frame(xlabel=paste(dataset.td.lu.16.III$eventID, dataset.td.lu.16.III$subblock, sep="_"),
                   avgTD=dataset.td.lu.16.III$avg, minTD=dataset.td.lu.16.III$min, maxTD=dataset.td.lu.16.III$max,
                   avgTI=dataset.ti.lu.16.I$avg, minTI=dataset.ti.lu.16.I$min, maxTI=dataset.ti.lu.16.I$max,
                   ratioavg=dataset.td.lu.16.III$avg/dataset.ti.lu.16.I$avg,
                   ratiomin=dataset.td.lu.16.III$min/dataset.ti.lu.16.I$min,
                   ratiomax=dataset.td.lu.16.III$max/dataset.ti.lu.16.I$max)

dIII_II<-data.frame(xlabel=paste(dataset.td.lu.16.III$eventID, dataset.td.lu.16.III$subblock, sep="_"),
                    avgTD=dataset.td.lu.16.III$avg, minTD=dataset.td.lu.16.III$min, maxTD=dataset.td.lu.16.III$max,
                    avgTI=dataset.ti.lu.16.II$avg, minTI=dataset.ti.lu.16.II$min, maxTI=dataset.ti.lu.16.II$max,
                    ratioavg=dataset.td.lu.16.III$avg/dataset.ti.lu.16.II$avg,
                    ratiomin=dataset.td.lu.16.III$min/dataset.ti.lu.16.II$min,
                    ratiomax=dataset.td.lu.16.III$max/dataset.ti.lu.16.II$max)

dIII_III<-data.frame(xlabel=paste(dataset.td.lu.16.III$eventID, dataset.td.lu.16.III$subblock, sep="_"),
                     avgTD=dataset.td.lu.16.III$avg, minTD=dataset.td.lu.16.III$min, maxTD=dataset.td.lu.16.III$max,
                     avgTI=dataset.ti.lu.16.III$avg, minTI=dataset.ti.lu.16.III$min, maxTI=dataset.ti.lu.16.III$max,
                     ratioavg=dataset.td.lu.16.III$avg/dataset.ti.lu.16.III$avg,
                     ratiomin=dataset.td.lu.16.III$min/dataset.ti.lu.16.III$min,
                     ratiomax=dataset.td.lu.16.III$max/dataset.ti.lu.16.III$max)

dIII_IV<-data.frame(xlabel=paste(dataset.td.lu.16.III$eventID, dataset.td.lu.16.III$subblock, sep="_"),
                    avgTD=dataset.td.lu.16.III$avg, minTD=dataset.td.lu.16.III$min, maxTD=dataset.td.lu.16.III$max,
                    avgTI=dataset.ti.lu.16.IV$avg, minTI=dataset.ti.lu.16.IV$min, maxTI=dataset.ti.lu.16.IV$max,
                    ratioavg=dataset.td.lu.16.III$avg/dataset.ti.lu.16.IV$avg,
                    ratiomin=dataset.td.lu.16.III$min/dataset.ti.lu.16.IV$min,
                    ratiomax=dataset.td.lu.16.III$max/dataset.ti.lu.16.IV$max)

dIII_V<-data.frame(xlabel=paste(dataset.td.lu.16.III$eventID, dataset.td.lu.16.III$subblock, sep="_"),
                   avgTD=dataset.td.lu.16.III$avg, minTD=dataset.td.lu.16.III$min, maxTD=dataset.td.lu.16.III$max,
                   avgTI=dataset.ti.lu.16.V$avg, minTI=dataset.ti.lu.16.V$min, maxTI=dataset.ti.lu.16.V$max,
                   ratioavg=dataset.td.lu.16.III$avg/dataset.ti.lu.16.V$avg,
                   ratiomin=dataset.td.lu.16.III$min/dataset.ti.lu.16.V$min,
                   ratiomax=dataset.td.lu.16.III$max/dataset.ti.lu.16.V$max)

dIV_I<-data.frame(xlabel=paste(dataset.td.lu.16.IV$eventID, dataset.td.lu.16.IV$subblock, sep="_"),
                  avgTD=dataset.td.lu.16.IV$avg, minTD=dataset.td.lu.16.IV$min, maxTD=dataset.td.lu.16.IV$max,
                  avgTI=dataset.ti.lu.16.I$avg, minTI=dataset.ti.lu.16.I$min, maxTI=dataset.ti.lu.16.I$max,
                  ratioavg=dataset.td.lu.16.IV$avg/dataset.ti.lu.16.I$avg,
                  ratiomin=dataset.td.lu.16.IV$min/dataset.ti.lu.16.I$min,
                  ratiomax=dataset.td.lu.16.IV$max/dataset.ti.lu.16.I$max)

dIV_II<-data.frame(xlabel=paste(dataset.td.lu.16.IV$eventID, dataset.td.lu.16.IV$subblock, sep="_"),
                   avgTD=dataset.td.lu.16.IV$avg, minTD=dataset.td.lu.16.IV$min, maxTD=dataset.td.lu.16.IV$max,
                   avgTI=dataset.ti.lu.16.II$avg, minTI=dataset.ti.lu.16.II$min, maxTI=dataset.ti.lu.16.II$max,
                   ratioavg=dataset.td.lu.16.IV$avg/dataset.ti.lu.16.II$avg,
                   ratiomin=dataset.td.lu.16.IV$min/dataset.ti.lu.16.II$min,
                   ratiomax=dataset.td.lu.16.IV$max/dataset.ti.lu.16.II$max)

dIV_III<-data.frame(xlabel=paste(dataset.td.lu.16.IV$eventID, dataset.td.lu.16.IV$subblock, sep="_"),
                    avgTD=dataset.td.lu.16.IV$avg, minTD=dataset.td.lu.16.IV$min, maxTD=dataset.td.lu.16.IV$max,
                    avgTI=dataset.ti.lu.16.III$avg, minTI=dataset.ti.lu.16.III$min, maxTI=dataset.ti.lu.16.III$max,
                    ratioavg=dataset.td.lu.16.IV$avg/dataset.ti.lu.16.III$avg,
                    ratiomin=dataset.td.lu.16.IV$min/dataset.ti.lu.16.III$min,
                    ratiomax=dataset.td.lu.16.IV$max/dataset.ti.lu.16.III$max)

dIV_IV<-data.frame(xlabel=paste(dataset.td.lu.16.IV$eventID, dataset.td.lu.16.IV$subblock, sep="_"),
                   avgTD=dataset.td.lu.16.IV$avg, minTD=dataset.td.lu.16.IV$min, maxTD=dataset.td.lu.16.IV$max,
                   avgTI=dataset.ti.lu.16.IV$avg, minTI=dataset.ti.lu.16.IV$min, maxTI=dataset.ti.lu.16.IV$max,
                   ratioavg=dataset.td.lu.16.IV$avg/dataset.ti.lu.16.IV$avg,
                   ratiomin=dataset.td.lu.16.IV$min/dataset.ti.lu.16.IV$min,
                   ratiomax=dataset.td.lu.16.IV$max/dataset.ti.lu.16.IV$max)

dIV_V<-data.frame(xlabel=paste(dataset.td.lu.16.IV$eventID, dataset.td.lu.16.IV$subblock, sep="_"),
                  avgTD=dataset.td.lu.16.IV$avg, minTD=dataset.td.lu.16.IV$min, maxTD=dataset.td.lu.16.IV$max,
                  avgTI=dataset.ti.lu.16.V$avg, minTI=dataset.ti.lu.16.V$min, maxTI=dataset.ti.lu.16.V$max,
                  ratioavg=dataset.td.lu.16.IV$avg/dataset.ti.lu.16.V$avg,
                  ratiomin=dataset.td.lu.16.IV$min/dataset.ti.lu.16.V$min,
                  ratiomax=dataset.td.lu.16.IV$max/dataset.ti.lu.16.V$max)

dV_I<-data.frame(xlabel=paste(dataset.td.lu.16.V$eventID, dataset.td.lu.16.V$subblock, sep="_"),
                 avgTD=dataset.td.lu.16.V$avg, minTD=dataset.td.lu.16.V$min, maxTD=dataset.td.lu.16.V$max,
                 avgTI=dataset.ti.lu.16.I$avg, minTI=dataset.ti.lu.16.I$min, maxTI=dataset.ti.lu.16.I$max,
                 ratioavg=dataset.td.lu.16.V$avg/dataset.ti.lu.16.I$avg,
                 ratiomin=dataset.td.lu.16.V$min/dataset.ti.lu.16.I$min,
                 ratiomax=dataset.td.lu.16.V$max/dataset.ti.lu.16.I$max)

dV_II<-data.frame(xlabel=paste(dataset.td.lu.16.V$eventID, dataset.td.lu.16.V$subblock, sep="_"),
                  avgTD=dataset.td.lu.16.V$avg, minTD=dataset.td.lu.16.V$min, maxTD=dataset.td.lu.16.V$max,
                  avgTI=dataset.ti.lu.16.II$avg, minTI=dataset.ti.lu.16.II$min, maxTI=dataset.ti.lu.16.II$max,
                  ratioavg=dataset.td.lu.16.V$avg/dataset.ti.lu.16.II$avg,
                  ratiomin=dataset.td.lu.16.V$min/dataset.ti.lu.16.II$min,
                  ratiomax=dataset.td.lu.16.V$max/dataset.ti.lu.16.II$max)

dV_III<-data.frame(xlabel=paste(dataset.td.lu.16.V$eventID, dataset.td.lu.16.V$subblock, sep="_"),
                   avgTD=dataset.td.lu.16.V$avg, minTD=dataset.td.lu.16.V$min, maxTD=dataset.td.lu.16.V$max,
                   avgTI=dataset.ti.lu.16.III$avg, minTI=dataset.ti.lu.16.III$min, maxTI=dataset.ti.lu.16.III$max,
                   ratioavg=dataset.td.lu.16.V$avg/dataset.ti.lu.16.III$avg,
                   ratiomin=dataset.td.lu.16.V$min/dataset.ti.lu.16.III$min,
                   ratiomax=dataset.td.lu.16.V$max/dataset.ti.lu.16.III$max)

dV_IV<-data.frame(xlabel=paste(dataset.td.lu.16.V$eventID, dataset.td.lu.16.V$subblock, sep="_"),
                  avgTD=dataset.td.lu.16.V$avg, minTD=dataset.td.lu.16.V$min, maxTD=dataset.td.lu.16.V$max,
                  avgTI=dataset.ti.lu.16.IV$avg, minTI=dataset.ti.lu.16.IV$min, maxTI=dataset.ti.lu.16.IV$max,
                  ratioavg=dataset.td.lu.16.V$avg/dataset.ti.lu.16.IV$avg,
                  ratiomin=dataset.td.lu.16.V$min/dataset.ti.lu.16.IV$min,
                  ratiomax=dataset.td.lu.16.V$max/dataset.ti.lu.16.IV$max)

dV_V<-data.frame(xlabel=paste(dataset.td.lu.16.V$eventID, dataset.td.lu.16.V$subblock, sep="_"),
                 avgTD=dataset.td.lu.16.V$avg, minTD=dataset.td.lu.16.V$min, maxTD=dataset.td.lu.16.V$max,
                 avgTI=dataset.ti.lu.16.V$avg, minTI=dataset.ti.lu.16.V$min, maxTI=dataset.ti.lu.16.V$max,
                 ratioavg=dataset.td.lu.16.V$avg/dataset.ti.lu.16.V$avg,
                 ratiomin=dataset.td.lu.16.V$min/dataset.ti.lu.16.V$min,
                 ratiomax=dataset.td.lu.16.V$max/dataset.ti.lu.16.V$max)

mean.avg<-matrix(0, 5, 5)
mean.avg[1,]<-c(mean(dI_I$ratioavg), mean(dI_II$ratioavg), mean(dI_III$ratioavg), mean(dI_IV$ratioavg), mean(dI_V$ratioavg))
mean.avg[2,]<-c(mean(dII_I$ratioavg), mean(dII_II$ratioavg), mean(dII_III$ratioavg), mean(dII_IV$ratioavg), mean(dII_V$ratioavg))
mean.avg[3,]<-c(mean(dIII_I$ratioavg), mean(dIII_II$ratioavg), mean(dIII_III$ratioavg), mean(dIII_IV$ratioavg), mean(dIII_V$ratioavg))
mean.avg[4,]<-c(mean(dIV_I$ratioavg), mean(dIV_II$ratioavg), mean(dIV_III$ratioavg), mean(dIV_IV$ratioavg), mean(dIV_V$ratioavg))
mean.avg[5,]<-c(mean(dV_I$ratioavg), mean(dV_II$ratioavg), mean(dV_III$ratioavg), mean(dV_IV$ratioavg), mean(dV_V$ratioavg))
print(mean.avg)

##       [,1]   [,2]   [,3]   [,4]   [,5]
## [1,] 1.113 0.7104 0.9504 0.6915 0.9507
## [2,] 1.113 0.7104 0.9504 0.6915 0.9507
## [3,] 1.113 0.7104 0.9504 0.6915 0.9507
## [4,] 1.113 0.7104 0.9504 0.6915 0.9507
## [5,] 1.113 0.7104 0.9504 0.6915 0.9507

mean.min<-matrix(0, 5, 5)
mean.min[1,]<-c(mean(dI_I$ratiomin), mean(dI_II$ratiomin), mean(dI_III$ratiomin), mean(dI_IV$ratiomin), mean(dI_V$ratiomin))
mean.min[2,]<-c(mean(dII_I$ratiomin), mean(dII_II$ratiomin), mean(dII_III$ratiomin), mean(dII_IV$ratiomin), mean(dII_V$ratiomin))
mean.min[3,]<-c(mean(dIII_I$ratiomin), mean(dIII_II$ratiomin), mean(dIII_III$ratiomin), mean(dIII_IV$ratiomin), mean(dIII_V$ratiomin))
mean.min[4,]<-c(mean(dIV_I$ratiomin), mean(dIV_II$ratiomin), mean(dIV_III$ratiomin), mean(dIV_IV$ratiomin), mean(dIV_V$ratiomin))
mean.min[5,]<-c(mean(dV_I$ratiomin), mean(dV_II$ratiomin), mean(dV_III$ratiomin), mean(dV_IV$ratiomin), mean(dV_V$ratiomin))
print(mean.min)

##       [,1]  [,2]  [,3]  [,4]  [,5]
## [1,] 1.074 1.216 0.899 1.106 1.038
## [2,] 1.074 1.216 0.899 1.106 1.038
## [3,] 1.074 1.216 0.899 1.106 1.038
## [4,] 1.074 1.216 0.899 1.106 1.038
## [5,] 1.074 1.216 0.899 1.106 1.038

mean.max<-matrix(0, 5, 5)
mean.max[1,]<-c(mean(dI_I$ratiomax), mean(dI_II$ratiomax), mean(dI_III$ratiomax), mean(dI_IV$ratiomax), mean(dI_V$ratiomax))
mean.max[2,]<-c(mean(dII_I$ratiomax), mean(dII_II$ratiomax), mean(dII_III$ratiomax), mean(dII_IV$ratiomax), mean(dII_V$ratiomax))
mean.max[3,]<-c(mean(dIII_I$ratiomax), mean(dIII_II$ratiomax), mean(dIII_III$ratiomax), mean(dIII_IV$ratiomax), mean(dIII_V$ratiomax))
mean.max[4,]<-c(mean(dIV_I$ratiomax), mean(dIV_II$ratiomax), mean(dIV_III$ratiomax), mean(dIV_IV$ratiomax), mean(dIV_V$ratiomax))
mean.max[5,]<-c(mean(dV_I$ratiomax), mean(dV_II$ratiomax), mean(dV_III$ratiomax), mean(dV_IV$ratiomax), mean(dV_V$ratiomax))
print(mean.max)

##       [,1]   [,2]  [,3]   [,4]  [,5]
## [1,] 2.261 0.7099 1.579 0.7091 1.414
## [2,] 2.261 0.7099 1.579 0.7091 1.414
## [3,] 2.261 0.7099 1.579 0.7091 1.414
## [4,] 2.261 0.7099 1.579 0.7091 1.414
## [5,] 2.261 0.7099 1.579 0.7091 1.414

#p <- ggplot(d, aes(d$xlabel, d$ratio))
#ggplot(data=mtcars, aes(x=wt, y=qsec, color=factor(cyl), shape=factor(gear))) + geom_point() + theme_bw() + geom_smooth (method="lm") + facet_wrap (~ cyl, scale="free_y")