results/analysis.R


1
2
3
4
5
6
7
8
9
10

library("tidyverse")


data=read_csv("results.csv")

final=data%>%group_by(simkey)%>%summarize(energy=mean(energy),sd=sd(energy))

ggplot(final,aes(x=simkey,y=energy))+geom_bar(stat="identity")+
    xlab("Scenario")+ylab("Energy (J)")+
     geom_errorbar(aes(ymin=energy-sd, ymax=energy+sd),width=.2,size=1.5)