Update

1 files changed, 9 insertions, 4 deletions

diff --git a/analysis/analysis.R b/analysis/analysis.R
index bad95ce..1a323a2 100644
--- a/analysis/analysis.R
+++ b/analysis/analysis.R
@@ -77,18 +77,23 @@ GG_VARIOUS=function(){
 
 
 GG_SIMULATIONS = function() {
+    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))+
+        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)")+
-        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")
+        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)
     ## Dimension Coverage
     ggplot(sim,aes(id,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")+
-        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")
+        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)
 
 }