rm(list=ls())
setwd('/Users/mpetric/Desktop/Grinn_WSN')
require(xlsx)## Loading required package: xlsx## Loading required package: rJava## Loading required package: xlsxjarsrequire(ggplot2)## Loading required package: ggplot2out.del<-read.table("outside delta.csv", sep="")
room.del<-read.table("room delta.csv",sep="")
sens.names<-paste('sens',c(1:25))summary for outdoor delT
summary(out.del)## V1 V2 V3 V4
## Min. :-0.40500 Min. :-0.8020 Min. :-0.7740 Min. :-0.1480
## 1st Qu.: 0.08975 1st Qu.:-0.3850 1st Qu.:-0.3123 1st Qu.: 0.2160
## Median : 0.28150 Median :-0.2380 Median :-0.0995 Median : 0.4000
## Mean : 0.40670 Mean :-0.1413 Mean :-0.0363 Mean : 0.4869
## 3rd Qu.: 0.81925 3rd Qu.: 0.1405 3rd Qu.: 0.2660 3rd Qu.: 0.7240
## Max. : 1.24000 Max. : 0.6800 Max. : 0.8300 Max. : 1.2800
## V5 V6 V7 V8
## Min. :-0.5770 Min. :-0.90000 Min. :-0.2870 Min. :-0.6080
## 1st Qu.:-0.1150 1st Qu.:-0.61250 1st Qu.: 0.2110 1st Qu.:-0.2165
## Median : 0.1330 Median :-0.32800 Median : 0.5155 Median : 0.0350
## Mean : 0.2469 Mean :-0.24770 Mean : 0.5735 Mean : 0.1229
## 3rd Qu.: 0.5473 3rd Qu.:-0.00825 3rd Qu.: 0.8692 3rd Qu.: 0.4365
## Max. : 1.2380 Max. : 0.62700 Max. : 1.4700 Max. : 1.0300
## V9 V10 V11 V12
## Min. :-0.3540 Min. :-1.0340 Min. :-0.81000 Min. :-0.6350
## 1st Qu.: 0.1103 1st Qu.:-0.5497 1st Qu.:-0.52250 1st Qu.:-0.2795
## Median : 0.3405 Median :-0.2810 Median :-0.24300 Median : 0.0040
## Mean : 0.4097 Mean :-0.2163 Mean :-0.17670 Mean : 0.0658
## 3rd Qu.: 0.7385 3rd Qu.: 0.1542 3rd Qu.: 0.03675 3rd Qu.: 0.3465
## Max. : 1.3100 Max. : 0.7400 Max. : 0.65700 Max. : 0.9900
## V13 V14 V15 V16
## Min. :-0.53500 Min. :-0.7670 Min. :-0.5100 Min. :-0.2840
## 1st Qu.:-0.05075 1st Qu.:-0.3175 1st Qu.:-0.1850 1st Qu.: 0.1978
## Median : 0.16900 Median :-0.0820 Median : 0.1070 Median : 0.4690
## Mean : 0.32470 Mean :-0.0031 Mean : 0.2063 Mean : 0.5327
## 3rd Qu.: 0.83300 3rd Qu.: 0.2447 3rd Qu.: 0.4793 3rd Qu.: 0.8968
## Max. : 1.17000 Max. : 0.8580 Max. : 1.1870 Max. : 1.4900
## V17 V18 V19 V20
## Min. :-0.6820 Min. :-0.6720 Min. :-0.3240 Min. :-0.7170
## 1st Qu.:-0.2525 1st Qu.:-0.2200 1st Qu.: 0.1477 1st Qu.:-0.2625
## Median :-0.0930 Median :-0.0180 Median : 0.3590 Median : 0.0030
## Mean : 0.0197 Mean : 0.1017 Mean : 0.4147 Mean : 0.1049
## 3rd Qu.: 0.3405 3rd Qu.: 0.4680 3rd Qu.: 0.7192 3rd Qu.: 0.4098
## Max. : 0.9000 Max. : 1.0800 Max. : 1.2600 Max. : 1.0780
## V21 V22 V23 V24
## Min. :-0.7400 Min. :-0.4260 Min. :-0.7670 Min. :-0.3250
## 1st Qu.:-0.4250 1st Qu.:-0.1960 1st Qu.:-0.3275 1st Qu.: 0.0450
## Median :-0.1130 Median :-0.0050 Median :-0.0820 Median : 0.2765
## Mean :-0.0077 Mean : 0.0596 Mean : 0.0319 Mean : 0.3952
## 3rd Qu.: 0.2943 3rd Qu.: 0.1442 3rd Qu.: 0.3272 3rd Qu.: 0.7400
## Max. : 1.0070 Max. : 0.8170 Max. : 1.0080 Max. : 1.2140
## V25
## Min. :-0.3200
## 1st Qu.:-0.0125
## Median : 0.2570
## Mean : 0.3443
## 3rd Qu.: 0.5667
## Max. : 1.2470summary for indoor delT
summary(room.del)## V1 V2 V3 V4
## Min. :-0.1810 Min. :-0.5310 Min. :-0.4310 Min. :-0.0910
## 1st Qu.:-0.0490 1st Qu.:-0.4205 1st Qu.:-0.3230 1st Qu.: 0.0105
## Median : 0.1070 Median :-0.3170 Median :-0.2570 Median : 0.0760
## Mean : 0.1659 Mean :-0.2195 Mean :-0.1773 Mean : 0.1608
## 3rd Qu.: 0.3215 3rd Qu.:-0.0010 3rd Qu.:-0.0385 3rd Qu.: 0.2315
## Max. : 1.0120 Max. : 0.5520 Max. : 0.8820 Max. : 1.2520
## V5 V6 V7 V8
## Min. :-0.8200 Min. :-0.2510 Min. :-0.29100 Min. :-0.7410
## 1st Qu.:-0.3760 1st Qu.:-0.1265 1st Qu.:-0.20250 1st Qu.:-0.6595
## Median :-0.3270 Median :-0.0140 Median :-0.10100 Median :-0.6010
## Mean :-0.2972 Mean : 0.1208 Mean :-0.01706 Mean :-0.4125
## 3rd Qu.:-0.2430 3rd Qu.: 0.4425 3rd Qu.: 0.16950 3rd Qu.:-0.1605
## Max. : 0.5720 Max. : 1.0920 Max. : 0.47200 Max. : 0.7620
## V9 V10 V11 V12
## Min. :-0.8930 Min. :-0.4910 Min. :-0.37100 Min. :-0.40100
## 1st Qu.:-0.4245 1st Qu.:-0.3840 1st Qu.:-0.26600 1st Qu.:-0.27850
## Median :-0.3820 Median :-0.2730 Median :-0.18400 Median :-0.14500
## Mean :-0.3855 Mean :-0.1518 Mean :-0.07568 Mean :-0.03194
## 3rd Qu.:-0.3510 3rd Qu.: 0.1475 3rd Qu.: 0.18350 3rd Qu.: 0.27250
## Max. : 0.5520 Max. : 0.7620 Max. : 0.75200 Max. : 0.93200
## V13 V14 V15 V16
## Min. :-0.42100 Min. :-0.37100 Min. :-0.9980 Min. :-0.0710
## 1st Qu.:-0.30550 1st Qu.:-0.27750 1st Qu.:-0.4845 1st Qu.: 0.0355
## Median :-0.17700 Median :-0.16100 Median :-0.4480 Median : 0.1630
## Mean :-0.07471 Mean :-0.06061 Mean :-0.4080 Mean : 0.2534
## 3rd Qu.: 0.10400 3rd Qu.: 0.13950 3rd Qu.:-0.3495 3rd Qu.: 0.4155
## Max. : 1.07200 Max. : 0.97200 Max. : 0.7420 Max. : 1.3920
## V17 V18 V19 V20
## Min. :-0.6210 Min. :-0.6750 Min. :-0.1710 Min. :-0.3920
## 1st Qu.:-0.5625 1st Qu.:-0.4015 1st Qu.:-0.0715 1st Qu.:-0.2890
## Median :-0.5250 Median :-0.3700 Median : 0.0370 Median :-0.2550
## Mean :-0.4664 Mean :-0.3205 Mean : 0.1086 Mean :-0.2035
## 3rd Qu.:-0.4100 3rd Qu.:-0.2900 3rd Qu.: 0.2645 3rd Qu.:-0.1750
## Max. :-0.0080 Max. : 0.6920 Max. : 1.2720 Max. : 0.4820
## V21 V22 V23 V24
## Min. :-0.925 Min. :-3.007 Min. :-0.4610 Min. :-0.3310
## 1st Qu.:-0.470 1st Qu.:-1.349 1st Qu.:-0.3745 1st Qu.:-0.2645
## Median :-0.441 Median :-0.669 Median :-0.3070 Median :-0.1930
## Mean : 5.742 Mean :-1.014 Mean :-0.1718 Mean :-0.1281
## 3rd Qu.:11.739 3rd Qu.:-0.574 3rd Qu.: 0.0975 3rd Qu.:-0.0420
## Max. :23.632 Max. : 0.642 Max. : 0.7720 Max. : 0.8120
## V25
## Min. :-0.2010
## 1st Qu.:-0.0640
## Median : 0.0760
## Mean : 0.1521
## 3rd Qu.: 0.2865
## Max. : 1.2920out.rmse<-apply(out.del,2,function(x) sqrt(mean(x^2)))
room.rmse<-apply(room.del,2,function(x) sqrt(mean(x^2)))
q.out<-data.frame(cbind(c(1:25),out.rmse))
colnames(q.out)<-c("sensor","rmse")
q.room<-data.frame(cbind(c(1:25),room.rmse))
colnames(q.room)<-c("sensor",'rmse')
qplot(sensor,rmse,data=q.out,geom=c('line',"point")) + geom_vline(xintercept=c(1:25),col="seashell3")
summary(out.del[,c(7,16)])## V7 V16
## Min. :-0.2870 Min. :-0.2840
## 1st Qu.: 0.2110 1st Qu.: 0.1978
## Median : 0.5155 Median : 0.4690
## Mean : 0.5735 Mean : 0.5327
## 3rd Qu.: 0.8692 3rd Qu.: 0.8968
## Max. : 1.4700 Max. : 1.4900qplot(sensor,rmse,data=q.room,geom=c('line',"point")) + geom_vline(xintercept=c(1:25),col="seashell3")
comp.room<-read.xlsx("ED_Sensor_comparison.xlsx",2,fill=TRUE)
comp.out<-read.xlsx("ED_Sensor_comparison.xlsx",1,fill=TRUE)summary(room.del[,21])## Min. 1st Qu. Median Mean 3rd Qu. Max.
## -0.925 -0.470 -0.441 5.742 11.740 23.630