1 files changed, 25 insertions, 1 deletions
diff --git a/results/analysis.R b/results/analysis.R
index d2baea7..2a946fb 100644
--- a/results/analysis.R
+++ b/results/analysis.R
@@ -1,7 +1,9 @@
 library("tidyverse")
 library("knitr")
 
-data=read_csv("results.csv")
+data=read_csv("results.csv")    
+ccgrid=read_csv("ccgrid.csv")
+
 
 r_=function(x){round(x,digits=2)}
 
@@ -30,11 +32,32 @@ receiverUptime=data%>%filter(isSender==0)%>%group_by(simkey,wireless,wakeupfor)%
 finalStats=finalStats%>%left_join(senderUptime,by=c("simkey","wireless","wakeupfor"))
 finalStats=finalStats%>%left_join(receiverUptime,by=c("simkey","wireless","wakeupfor"))
 
+# Comparison
+finalStats=finalStats%>%left_join(ccgrid,by=c("simkey","wireless","wakeupfor"),suffix = c("","_ccgrid"))
+finalStats=finalStats%>%mutate(senderEnergy_err=abs(r_((senderEnergy-senderEnergy_ccgrid)/senderEnergy)*100))
+finalStats=finalStats%>%mutate(receiverEnergy_err=abs(r_((receiverEnergy-receiverEnergy_ccgrid)/receiverEnergy)*100))
+finalStats=finalStats%>%mutate(senderUptime_err=abs(r_((senderUptime-senderUptime_ccgrid)/senderUptime)*100))
+finalStats=finalStats%>%mutate(receiverUptime_err=abs(r_((receiverUptime-receiverUptime_ccgrid)/receiverUptime)*100))
 
 # Add nSendFail to stats
 sendFail=data%>%filter(node=="on0")%>%group_by(simkey,wireless,wakeupfor)%>%summarise(nSendFail=mean(nSendFail))
 finalStats=finalStats%>%left_join(sendFail,by=c("simkey","wireless","wakeupfor"))
 
+
+ggplot(finalStats,aes(x=simkey,senderEnergy_err,color=wireless,group=wireless))+geom_point()+geom_line()+facet_wrap(~wakeupfor)+xlab("Scenarios") + ylab("Sender Energy Relative Error (%)") + expand_limits(x = 0, y = 0)
+ggsave("err_sender_energy.png",width=7,height = 5)
+
+ggplot(finalStats,aes(x=simkey,receiverEnergy_err,color=wireless,group=wireless))+geom_point()+geom_line()+facet_wrap(~wakeupfor)+xlab("Scenarios") + ylab("Receiver Energy Relative Error (%)") + expand_limits(x = 0, y = 0)
+ggsave("err_receiver_energy.png",width=7,height = 5)
+
+ggplot(finalStats,aes(x=simkey,senderUptime_err,color=wireless,group=wireless))+geom_point()+geom_line()+facet_wrap(~wakeupfor)+xlab("Scenarios") + ylab("Sender AccUptime Energy Relative Error (%)") + expand_limits(x = 0, y = 0)
+ggsave("err_sender_accuptime.png",width=7,height = 5)
+
+ggplot(finalStats,aes(x=simkey,receiverUptime_err,color=wireless,group=wireless))+geom_point()+geom_line()+facet_wrap(~wakeupfor)+xlab("Scenarios") + ylab("Receiver AccUptime Energy Relative Error (%)") + expand_limits(x = 0, y = 0)
+ggsave("err_rcv_accuptime.png",width=7,height = 5)
+
+
+stopifnot(1)
 # Generate plots
 ggplot(finalEnergy,aes(x=simkey,y=energy,fill=wireless))+
     geom_bar(stat="identity",position=position_dodge())+
@@ -57,6 +80,7 @@ ggplot(finalStats,aes(x=simkey,y=success,fill=wireless))+
     facet_wrap(~ wakeupfor) +
 ggsave("success.png")
 
+##### Save tables #####
 tf="tables-energy.org"
 write("",file=tf)
 for(wake in unique(data$wakeupfor)){