1 files changed, 25 insertions, 11 deletions
diff --git a/analysis/analysis.R b/analysis/analysis.R
index 1a323a2..f2343a7 100644
--- a/analysis/analysis.R
+++ b/analysis/analysis.R
@@ -36,6 +36,8 @@ sim=sim%>%group_by(simkey,wireless,wakeupfor,seed,energy_constraint,coverage_con
 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
+rand_id=sample(seq(1,NROW(sim)/8))
+sim=sim%>%mutate(id_rand=rand_id[id])
 
 ## Perform various plots to understand the data
 GG_VARIOUS=function(){
@@ -73,31 +75,43 @@ GG_VARIOUS=function(){
 }
 
 ## Plot to generate:
-#GG_VARIOUS()
+GG_VARIOUS()
 
 
-GG_SIMULATIONS = function() {
+GG_SIMULATIONS = function(file,key="NOKEY") {
+    ## Load simulation results
+    sim=read_csv(file)
+    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
+    rand_id=sample(seq(1,NROW(sim)/8))
+    sim=sim%>%mutate(id_rand=rand_id[id])
+
     sim_means=sim%>%group_by(wireless,wakeupfor,model)%>%summarize(mean_energy=mean(mean_energy_delta),mean_coverage=mean(mean_coverage_delta))
     ## Dimension Energy
-    ggplot(sim,aes(id,mean_energy_delta,color=model))+
+    ggplot(sim,aes(id_rand,mean_energy_delta,color=model))+
         geom_hline(data=sim_means,aes(yintercept=mean_energy,color=model),linetype="dashed")+
         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)")+
+        geom_point(aes(shape=simkey),size=2)+geom_line()+facet_wrap(~wireless+wakeupfor,scale="free")+ggtitle(paste0("Simulation Dimension Energy (",key,")"))+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")+
-        guides(fill=guide_legend(title="Models"),color=guide_legend(title="Models"))
-    ggsave("figures/sim_dimension_energy.pdf",width=11)
+        guides(fill=guide_legend(title="Models"),color=guide_legend(title="Models"))+
+        theme(axis.text.x=element_blank(),axis.ticks.x=element_blank())
+    ggsave(paste0("figures/sim_dimension_energy_",key,".pdf"),width=20)
     ## Dimension Coverage
-    ggplot(sim,aes(id,mean_coverage_delta,color=model))+
+    ggplot(sim,aes(id_rand,mean_coverage_delta,color=model))+
         geom_hline(data=sim_means,aes(yintercept=mean_coverage,color=model),linetype="dashed")+
         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")+
+        geom_point(aes(shape=simkey),size=2)+geom_line()+facet_wrap(~wireless+wakeupfor,scale="free")+ggtitle(paste0("Simulation Dimension Coverage (",key,")"))+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")+
-        guides(fill=guide_legend(title="Models"),color=guide_legend(title="Models"))
-    ggsave("figures/sim_dimension_coverage.pdf",width=11)
+        guides(fill=guide_legend(title="Models"),color=guide_legend(title="Models"))+
+        theme(axis.text.x=element_blank(),axis.ticks.x=element_blank())
+    ggsave(paste0("figures/sim_dimension_coverage_",key,".pdf"),width=20)
 
 }
 
-GG_SIMULATIONS()
+GG_SIMULATIONS("../simulations/results/results_NoHintIsFALSE.csv",key="WITH_HINT")
+GG_SIMULATIONS("../simulations/results/results_NoHintIsTRUE.csv",key="NO_HINT")
 
 system("rm -f figures/combined.pdf")
 system("pdfunite figures/*.pdf figures/combined.pdf")