diff options
| author | Loic Guegan <manzerbredes@mailbox.org> | 2021-05-11 08:53:11 +0200 |
|---|---|---|
| committer | Loic Guegan <manzerbredes@mailbox.org> | 2021-05-11 08:53:11 +0200 |
| commit | 1477290330630bdf2d2ccf019b6bc416c8866dcd (patch) | |
| tree | a5b608bcc809b2171f6387e6070ec335e2a152cd /results/analysis.R | |
| parent | 9c9be6752d32106d18226e4eed7b7d68da79d178 (diff) | |
Improve data analysis and simulation
Diffstat (limited to 'results/analysis.R')
| -rw-r--r-- | results/analysis.R | 14 |
1 files changed, 10 insertions, 4 deletions
diff --git a/results/analysis.R b/results/analysis.R index 837daf6..c3eb093 100644 --- a/results/analysis.R +++ b/results/analysis.R @@ -7,7 +7,13 @@ finalEnergy=data%>%group_by(simkey,wireless,wakeupfor)%>%summarize(wireless=firs finalStats=data%>%group_by(simkey,wireless,wakeupfor,seed)%>% summarise(success=sum(nSend),totalUptime=mean(totalUptime))%>% ungroup()%>%group_by(simkey,wireless,wakeupfor)%>% - summarise(sd=sd(success),success=mean(success),totalUptime=first(totalUptime)) + summarise(sd_success=sd(success),success=mean(success),totalUptime=first(totalUptime)) + +# Sender energy +senderEnergy=data%>%filter(isSender!=0)%>%group_by(simkey,wireless,wakeupfor)%>%summarise(senderEnergy=mean(energy)) +receiverEnergy=data%>%filter(isSender==0)%>%group_by(simkey,wireless,wakeupfor)%>%summarise(receiverEnergy=mean(energy)) +finalStats=finalStats%>%left_join(senderEnergy,by=c("simkey","wireless","wakeupfor")) +finalStats=finalStats%>%left_join(receiverEnergy,by=c("simkey","wireless","wakeupfor")) # Add nSendFail to stats sendFail=data%>%filter(node=="on0")%>%group_by(simkey,wireless,wakeupfor)%>%summarise(nSendFail=mean(nSendFail)) @@ -28,9 +34,9 @@ ggsave("energy.png") ggplot(finalStats,aes(x=simkey,y=success,fill=wireless))+ geom_bar(stat="identity",position=position_dodge())+ xlab("Scenarios")+ylab("Energy (J)")+ - geom_errorbar(position=position_dodge(0.9),aes(ymin=success-sd, ymax=success+sd),width=0.5,size=1.5,) + - geom_text(aes(label = round(success+sd,digits=1), y= success+sd),vjust=-0.4,position=position_dodge(0.9))+ - geom_text(aes(label = round(success-sd,digits=1), y= success-sd),vjust=+1.5,position=position_dodge(0.9))+ + geom_errorbar(position=position_dodge(0.9),aes(ymin=success-sd_success, ymax=success+sd_success),width=0.5,size=1.5,) + + geom_text(aes(label = round(success+sd_success,digits=1), y= success+sd_success),vjust=-0.4,position=position_dodge(0.9))+ + geom_text(aes(label = round(success-sd_success,digits=1), y= success-sd_success),vjust=+1.5,position=position_dodge(0.9))+ geom_text(aes(label = round(success,digits=1), y= success/2),fontface="bold",colour="white",position=position_dodge(0.9))+ facet_wrap(~ wakeupfor) + ggsave("success.png") |