Title

library(ggplot2)

dataset.td.lu.16.I<-read.table("dataframe.comp.td.I", header=TRUE)
dataset.td.lu.16.II<-read.table("dataframe.comp.td.II", header=TRUE)
dataset.td.lu.16.III<-read.table("dataframe.comp.td.III", header=TRUE)
dataset.td.lu.16.IV<-read.table("dataframe.comp.td.IV", header=TRUE)
dataset.td.lu.16.V<-read.table("dataframe.comp.td.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.td.lu.16.I$avg, minTI=dataset.td.lu.16.I$min, maxTI=dataset.td.lu.16.I$max,
                 ratioavg=dataset.td.lu.16.I$avg/dataset.td.lu.16.I$avg,
                 ratiomin=dataset.td.lu.16.I$min/dataset.td.lu.16.I$min,
                 ratiomax=dataset.td.lu.16.I$max/dataset.td.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.td.lu.16.II$avg, minTI=dataset.td.lu.16.II$min, maxTI=dataset.td.lu.16.II$max,
                  ratioavg=dataset.td.lu.16.I$avg/dataset.td.lu.16.II$avg,
                  ratiomin=dataset.td.lu.16.I$min/dataset.td.lu.16.II$min,
                  ratiomax=dataset.td.lu.16.I$max/dataset.td.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.td.lu.16.III$avg, minTI=dataset.td.lu.16.III$min, maxTI=dataset.td.lu.16.III$max,
                   ratioavg=dataset.td.lu.16.I$avg/dataset.td.lu.16.III$avg,
                   ratiomin=dataset.td.lu.16.I$min/dataset.td.lu.16.III$min,
                   ratiomax=dataset.td.lu.16.I$max/dataset.td.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.td.lu.16.IV$avg, minTI=dataset.td.lu.16.IV$min, maxTI=dataset.td.lu.16.IV$max,
                  ratioavg=dataset.td.lu.16.I$avg/dataset.td.lu.16.IV$avg,
                  ratiomin=dataset.td.lu.16.I$min/dataset.td.lu.16.IV$min,
                  ratiomax=dataset.td.lu.16.I$max/dataset.td.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.td.lu.16.V$avg, minTI=dataset.td.lu.16.V$min, maxTI=dataset.td.lu.16.V$max,
                 ratioavg=dataset.td.lu.16.I$avg/dataset.td.lu.16.V$avg,
                 ratiomin=dataset.td.lu.16.I$min/dataset.td.lu.16.V$min,
                 ratiomax=dataset.td.lu.16.I$max/dataset.td.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.td.lu.16.I$avg, minTI=dataset.td.lu.16.I$min, maxTI=dataset.td.lu.16.I$max,
                  ratioavg=dataset.td.lu.16.II$avg/dataset.td.lu.16.I$avg,
                  ratiomin=dataset.td.lu.16.II$min/dataset.td.lu.16.I$min,
                  ratiomax=dataset.td.lu.16.II$max/dataset.td.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.td.lu.16.II$avg, minTI=dataset.td.lu.16.II$min, maxTI=dataset.td.lu.16.II$max,
                   ratioavg=dataset.td.lu.16.II$avg/dataset.td.lu.16.II$avg,
                   ratiomin=dataset.td.lu.16.II$min/dataset.td.lu.16.II$min,
                   ratiomax=dataset.td.lu.16.II$max/dataset.td.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.td.lu.16.III$avg, minTI=dataset.td.lu.16.III$min, maxTI=dataset.td.lu.16.III$max,
                    ratioavg=dataset.td.lu.16.II$avg/dataset.td.lu.16.III$avg,
                    ratiomin=dataset.td.lu.16.II$min/dataset.td.lu.16.III$min,
                    ratiomax=dataset.td.lu.16.II$max/dataset.td.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.td.lu.16.IV$avg, minTI=dataset.td.lu.16.IV$min, maxTI=dataset.td.lu.16.IV$max,
                   ratioavg=dataset.td.lu.16.II$avg/dataset.td.lu.16.IV$avg,
                   ratiomin=dataset.td.lu.16.II$min/dataset.td.lu.16.IV$min,
                   ratiomax=dataset.td.lu.16.II$max/dataset.td.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.td.lu.16.V$avg, minTI=dataset.td.lu.16.V$min, maxTI=dataset.td.lu.16.V$max,
                  ratioavg=dataset.td.lu.16.II$avg/dataset.td.lu.16.V$avg,
                  ratiomin=dataset.td.lu.16.II$min/dataset.td.lu.16.V$min,
                  ratiomax=dataset.td.lu.16.II$max/dataset.td.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.td.lu.16.I$avg, minTI=dataset.td.lu.16.I$min, maxTI=dataset.td.lu.16.I$max,
                   ratioavg=dataset.td.lu.16.III$avg/dataset.td.lu.16.I$avg,
                   ratiomin=dataset.td.lu.16.III$min/dataset.td.lu.16.I$min,
                   ratiomax=dataset.td.lu.16.III$max/dataset.td.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.td.lu.16.II$avg, minTI=dataset.td.lu.16.II$min, maxTI=dataset.td.lu.16.II$max,
                    ratioavg=dataset.td.lu.16.III$avg/dataset.td.lu.16.II$avg,
                    ratiomin=dataset.td.lu.16.III$min/dataset.td.lu.16.II$min,
                    ratiomax=dataset.td.lu.16.III$max/dataset.td.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.td.lu.16.III$avg, minTI=dataset.td.lu.16.III$min, maxTI=dataset.td.lu.16.III$max,
                     ratioavg=dataset.td.lu.16.III$avg/dataset.td.lu.16.III$avg,
                     ratiomin=dataset.td.lu.16.III$min/dataset.td.lu.16.III$min,
                     ratiomax=dataset.td.lu.16.III$max/dataset.td.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.td.lu.16.IV$avg, minTI=dataset.td.lu.16.IV$min, maxTI=dataset.td.lu.16.IV$max,
                    ratioavg=dataset.td.lu.16.III$avg/dataset.td.lu.16.IV$avg,
                    ratiomin=dataset.td.lu.16.III$min/dataset.td.lu.16.IV$min,
                    ratiomax=dataset.td.lu.16.III$max/dataset.td.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.td.lu.16.V$avg, minTI=dataset.td.lu.16.V$min, maxTI=dataset.td.lu.16.V$max,
                   ratioavg=dataset.td.lu.16.III$avg/dataset.td.lu.16.V$avg,
                   ratiomin=dataset.td.lu.16.III$min/dataset.td.lu.16.V$min,
                   ratiomax=dataset.td.lu.16.III$max/dataset.td.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.td.lu.16.I$avg, minTI=dataset.td.lu.16.I$min, maxTI=dataset.td.lu.16.I$max,
                  ratioavg=dataset.td.lu.16.IV$avg/dataset.td.lu.16.I$avg,
                  ratiomin=dataset.td.lu.16.IV$min/dataset.td.lu.16.I$min,
                  ratiomax=dataset.td.lu.16.IV$max/dataset.td.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.td.lu.16.II$avg, minTI=dataset.td.lu.16.II$min, maxTI=dataset.td.lu.16.II$max,
                   ratioavg=dataset.td.lu.16.IV$avg/dataset.td.lu.16.II$avg,
                   ratiomin=dataset.td.lu.16.IV$min/dataset.td.lu.16.II$min,
                   ratiomax=dataset.td.lu.16.IV$max/dataset.td.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.td.lu.16.III$avg, minTI=dataset.td.lu.16.III$min, maxTI=dataset.td.lu.16.III$max,
                    ratioavg=dataset.td.lu.16.IV$avg/dataset.td.lu.16.III$avg,
                    ratiomin=dataset.td.lu.16.IV$min/dataset.td.lu.16.III$min,
                    ratiomax=dataset.td.lu.16.IV$max/dataset.td.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.td.lu.16.IV$avg, minTI=dataset.td.lu.16.IV$min, maxTI=dataset.td.lu.16.IV$max,
                   ratioavg=dataset.td.lu.16.IV$avg/dataset.td.lu.16.IV$avg,
                   ratiomin=dataset.td.lu.16.IV$min/dataset.td.lu.16.IV$min,
                   ratiomax=dataset.td.lu.16.IV$max/dataset.td.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.td.lu.16.V$avg, minTI=dataset.td.lu.16.V$min, maxTI=dataset.td.lu.16.V$max,
                  ratioavg=dataset.td.lu.16.IV$avg/dataset.td.lu.16.V$avg,
                  ratiomin=dataset.td.lu.16.IV$min/dataset.td.lu.16.V$min,
                  ratiomax=dataset.td.lu.16.IV$max/dataset.td.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.td.lu.16.I$avg, minTI=dataset.td.lu.16.I$min, maxTI=dataset.td.lu.16.I$max,
                 ratioavg=dataset.td.lu.16.V$avg/dataset.td.lu.16.I$avg,
                 ratiomin=dataset.td.lu.16.V$min/dataset.td.lu.16.I$min,
                 ratiomax=dataset.td.lu.16.V$max/dataset.td.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.td.lu.16.II$avg, minTI=dataset.td.lu.16.II$min, maxTI=dataset.td.lu.16.II$max,
                  ratioavg=dataset.td.lu.16.V$avg/dataset.td.lu.16.II$avg,
                  ratiomin=dataset.td.lu.16.V$min/dataset.td.lu.16.II$min,
                  ratiomax=dataset.td.lu.16.V$max/dataset.td.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.td.lu.16.III$avg, minTI=dataset.td.lu.16.III$min, maxTI=dataset.td.lu.16.III$max,
                   ratioavg=dataset.td.lu.16.V$avg/dataset.td.lu.16.III$avg,
                   ratiomin=dataset.td.lu.16.V$min/dataset.td.lu.16.III$min,
                   ratiomax=dataset.td.lu.16.V$max/dataset.td.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.td.lu.16.IV$avg, minTI=dataset.td.lu.16.IV$min, maxTI=dataset.td.lu.16.IV$max,
                  ratioavg=dataset.td.lu.16.V$avg/dataset.td.lu.16.IV$avg,
                  ratiomin=dataset.td.lu.16.V$min/dataset.td.lu.16.IV$min,
                  ratiomax=dataset.td.lu.16.V$max/dataset.td.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.td.lu.16.V$avg, minTI=dataset.td.lu.16.V$min, maxTI=dataset.td.lu.16.V$max,
                 ratioavg=dataset.td.lu.16.V$avg/dataset.td.lu.16.V$avg,
                 ratiomin=dataset.td.lu.16.V$min/dataset.td.lu.16.V$min,
                 ratiomax=dataset.td.lu.16.V$max/dataset.td.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.0000 1.009 1.0017 1.0035 1.0077
## [2,] 0.9953 1.000 0.9944 0.9951 0.9989
## [3,] 1.0108 1.018 1.0000 1.0105 1.0161
## [4,] 1.0025 1.009 1.0008 1.0000 1.0067
## [5,] 0.9982 1.003 0.9973 0.9979 1.0000

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.000 1.004 0.9956 1.0002 1.0010
## [2,] 1.001 1.000 0.9938 0.9991 0.9992
## [3,] 1.016 1.017 1.0000 1.0139 1.0146
## [4,] 1.004 1.006 0.9975 1.0000 1.0034
## [5,] 1.004 1.005 0.9971 1.0023 1.0000

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,] 1.0000 1.078 1.070 1.0398 1.072
## [2,] 0.9809 1.000 1.026 0.9871 1.018
## [3,] 1.0002 1.047 1.000 1.0076 1.041
## [4,] 1.0236 1.060 1.071 1.0000 1.061
## [5,] 1.0054 1.037 1.043 1.0113 1.000