1 files changed, 43 insertions, 4 deletions
diff --git a/analysis/analysis.R b/analysis/analysis.R
index 5b3c9fe..2e54763 100644
--- a/analysis/analysis.R
+++ b/analysis/analysis.R
@@ -7,13 +7,52 @@ library("tidyverse")
 ## seed           [1,200]
 ## node           on[0,12]
 ## isSender       {0,1}
+## dataSize       {1MB}
 
 ## Metrics:
 ## energy         [0,+inf)
 ## nDataRcv       [0,+inf)
 
-data=read_csv("../CCGRID2022.csv")
-data_seed=data%>%
-    group_by(simkey,wireless,wakeupfor,seed)%>%mutate(coverage=sum(nDataRcv))%>%
-    group_by(simkey,wireless,wakeupfor,seed,isSender)%>%summarize(energy=mean(energy))%>%
+nseed=200
+nwakeupfor=2
+nwireless=2
+nsimkey=4
+nsimulations=nseed*nwakeupfor*nwireless*nsimkey # Must be 3200
+
+## 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)
+data_seed_isSender=data%>%group_by(simkey,wireless,wakeupfor,seed,isSender)%>%summarize(energy_mean=mean(energy))%>%
+    left_join(tmp_data_coverage,by=c("simkey","wireless","wakeupfor","seed"))%>%
+    mutate(efficiency=energy_mean/coverage)%>%
+    ungroup()
+data_seed=data%>%group_by(simkey,wireless,wakeupfor,seed)%>%summarize(energy=sum(energy),coverage=sum(nDataRcv))%>%
+    mutate(efficiency=energy/coverage)%>%
     ungroup()
+
+
+
+## Perform various plots to understand the data
+GG_VARIOUS=function(){
+    ## Dimension Coverage
+    ggplot(data_seed%>%mutate(wakeupfor=as.character(wakeupfor)),
+           aes(simkey,coverage,color=wakeupfor))+geom_point()+facet_wrap(~wireless+wakeupfor)+ggtitle("Dimension Coverage")
+    ggsave("figures/dimension_coverage.pdf")
+    ## Dimension Energy
+    ggplot(data_seed%>%mutate(wakeupfor=as.character(wakeupfor)),
+           aes(simkey,energy,color=wakeupfor))+geom_point()+facet_wrap(~wireless+wakeupfor)+ggtitle("Dimension Energy")+ylab("Energy Consumption (J)")
+    ggsave("figures/dimension_energy.pdf")
+    ## Dimension Efficiency
+    ggplot(data_seed%>%mutate(wakeupfor=as.character(wakeupfor)),
+           aes(simkey,efficiency,color=wakeupfor))+
+        geom_point()+facet_wrap(~wireless+wakeupfor)+ggtitle("Dimension Efficiency")+ylab("Efficiency (J/Succ)")
+    ggsave("figures/dimension_efficiency.pdf")
+}
+
+## Plot to generate:
+GG_VARIOUS()
+
+
+
+system("rm -f figures/combined.pdf")
+system("pdfunite figures/*.pdf figures/combined.pdf")