1 files changed, 14 insertions, 6 deletions
diff --git a/analysis/analysis.R b/analysis/analysis.R
index 110549f..bad95ce 100644
--- a/analysis/analysis.R
+++ b/analysis/analysis.R
@@ -32,10 +32,10 @@ data_seed=data%>%group_by(simkey,wireless,wakeupfor,seed)%>%summarize(energy=sum
 
 ## Load simulation results
 sim=read_csv("../simulations/results/results.csv")
-data_sim=sim%>%group_by(simkey,wireless,wakeupfor,seed)%>%summarize(energy=sum(energy),coverage=sum(nDataRcv))%>%
-    mutate(efficiency=energy/coverage)%>%
-    ungroup()
-
+sim=sim%>%group_by(simkey,wireless,wakeupfor,seed,energy_constraint,coverage_constraint,model)%>%summarize(energy=sum(energy),coverage=sum(nDataRcv))
+sim=sim%>%mutate(energy_delta=energy-energy_constraint,coverage_delta=coverage-coverage_constraint)
+sim=sim%>%group_by(simkey,wireless,wakeupfor,energy_constraint,coverage_constraint,model)%>%summarize(mean_energy_delta=mean(energy_delta),mean_coverage_delta=mean(coverage_delta),sd_energy_delta=sd(energy_delta),sd_coverage_delta=sd(coverage_delta))
+sim=sim%>%arrange(energy_constraint,coverage_constraint)%>%group_by(wireless,wakeupfor,model)%>%mutate(id=row_number())%>%ungroup() # Create id for each experiment
 
 ## Perform various plots to understand the data
 GG_VARIOUS=function(){
@@ -78,9 +78,17 @@ GG_VARIOUS=function(){
 
 GG_SIMULATIONS = function() {
     ## Dimension Energy
-    ggplot(sim%>%mutate(wakeupfor=as.character(wakeupfor)),
-           aes(simkey,energy,color=model))+geom_point()+facet_wrap(~wireless+wakeupfor)+ggtitle("Simulation Dimension Energy")+ylab("Energy Consumption (J)")
+    ggplot(sim,aes(id,mean_energy_delta,color=model))+
+        geom_ribbon(aes(ymin=mean_energy_delta-sd_energy_delta,ymax=mean_energy_delta+sd_energy_delta,fill=model),alpha=0.2,color=NA)+
+        geom_point(aes(shape=simkey),size=2)+geom_line()+facet_wrap(~wireless+wakeupfor,scale="free")+ggtitle("Simulation Dimension Energy")+ylab("Energy Relative Absolute Error (J)")+
+        scale_x_continuous("Experiments id",breaks=seq(1,max(sim$id)))+theme(panel.grid.minor = element_blank())+labs(shape="Predictions")
     ggsave("figures/sim_dimension_energy.pdf")
+    ## Dimension Coverage
+    ggplot(sim,aes(id,mean_coverage_delta,color=model))+
+        geom_ribbon(aes(ymin=mean_coverage_delta-sd_coverage_delta,ymax=mean_coverage_delta+sd_coverage_delta,fill=model),alpha=0.2,color=NA)+
+        geom_point(aes(shape=simkey),size=2)+geom_line()+facet_wrap(~wireless+wakeupfor,scale="free")+ggtitle("Simulation Dimension Coverage")+ylab("Coverage Relative Absolute Error")+
+        scale_x_continuous("Experiments id",breaks=seq(1,max(sim$id)))+theme(panel.grid.minor = element_blank())+labs(shape="Predictions")
+    ggsave("figures/sim_dimension_coverage.pdf")
 
 }