Cleaning analysis code

1 files changed, 4 insertions, 5 deletions

diff --git a/analysis/analysis.R b/analysis/analysis.R
index 90874ed..46d8788 100644
--- a/analysis/analysis.R
+++ b/analysis/analysis.R
@@ -20,7 +20,7 @@ nwireless=2
 nsimkey=4
 nsimulations=nseed*nwakeupfor*nwireless*nsimkey # Must be 3200
 
-if(F){ # Toggle
+if(!exists("data_seed")){ # Auto-load data only one time
 ## Load data
 data=read_csv("../CCGRID2022.csv")%>%distinct() # Note that in the data experiment wireless=="lora",seed==1,wakeupfor==60,simkey=="baseline" is present 2 times in the CSV file
 tmp_data_coverage=data%>%group_by(simkey,wireless,wakeupfor,seed)%>%mutate(coverage=sum(nDataRcv))%>%ungroup()%>%filter(isSender==1)%>%select(simkey,wireless,wakeupfor,seed,coverage)
@@ -52,10 +52,9 @@ GG_VARIOUS=function(){
         guides(colour = guide_legend(override.aes = list(alpha = 1)))+
         facet_wrap(~wakeupfor+wireless,scale="free")+scale_color_viridis(discrete=TRUE,option="B")+
         theme(legend.position = c(0.15,0.87),
-                         legend.background = element_rect(fill = "white", color = "black",size=0.8))+
+                         legend.background = element_rect(fill = "white", color = "black",linewidth=0.8))+
         scale_x_continuous(breaks=seq(0,12))
     ggsave("figures/dimension_energy-coverage.pdf")
-    stopifnot(1)
     ## Dimension Energy/Coverage Wakeupfor
     ggplot(data_seed%>%mutate(wakeupfor=as.character(wakeupfor)),
            aes(coverage,energy,color=wakeupfor))+geom_point(alpha=0.5)+labs(color="Policy")+theme(legend.position="top")+
@@ -67,7 +66,7 @@ GG_VARIOUS=function(){
         ggtitle("Dimension Energy/Coverage Policy")+xlab("Coverage")+ylab("Nodes energy consumption (J)")
     ggsave("figures/dimension_energy-coverage-policy.pdf")
 
-    
+
 }
 
 ## Plot to generate:
@@ -86,7 +85,7 @@ GG_SIMULATIONS = function(file,key="NOKEY") {
 
     sim_means=sim%>%group_by(wireless,wakeupfor,model)%>%summarize(mean_energy=mean(mean_energy_delta),mean_coverage=mean(mean_coverage_delta),sd_energy=sd(sd_energy_delta),sd_coverage=sd(sd_coverage_delta))
     write.csv(sim_means%>%mutate_if(is.numeric, ~round(.,digits=2)),paste0("figures/sim_results_",key,".csv"),row.names=FALSE, quote=FALSE)
-    
+
     ## Dimension Energy
     ggplot(sim,aes(id_rand,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)+