diff --git a/Analysis.r b/Analysis.r
new file mode 100644
index 0000000..4a18fa9
--- /dev/null
+++ b/Analysis.r
@@ -0,0 +1,20 @@
+library(tidyverse)
+library(RcppRoll)
+RES <- readRDS("Results/Storage_Values_by_Facility_and_Variable_Discounts.Rds")
+RES[[1]]
+names(RES)
+read_csv("Data/Curie_Spent_Fuel_Site_Totals.csv")
+CDATA <- RES$Per_5  %>% left_join(read_csv("Data/Curie_Spent_Fuel_Site_Totals.csv")) %>%  mutate(Revenue=parse_number(Revenue),Year=parse_number(Year),Marginal_Value=Revenue/Total_Tons) %>% select(Facility,Year,Revenue,Marginal_Value,Total_Tons)
+CDATA <- CDATA %>% group_by(Year) %>% arrange(Year,desc(Marginal_Value)) %>% mutate(CIFS_total=cumsum(Total_Tons)) %>% ungroup 
+CDATA <- CDATA %>% group_by(Year) %>% mutate(CIFS_20=CIFS_total<=20000,CIFS_40=CIFS_total<=40000,CIFS_80=CIFS_total<=80000,CIFS_150=CIFS_total<=150000) %>% ungroup
+SUMMARY <- rbind(CDATA %>% filter(CIFS_20)  %>% group_by(Year) %>% summarize(Revenue=sum(Revenue),Marginal_Value=min(Marginal_Value),Total_tons=sum(Total_Tons),CIFS_Size=20000),
+		 CDATA %>% filter(CIFS_40)  %>% group_by(Year) %>% summarize(Revenue=sum(Revenue),Marginal_Value=min(Marginal_Value),Total_tons=sum(Total_Tons),CIFS_Size=40000),
+CDATA %>% filter(CIFS_80)  %>% group_by(Year) %>% summarize(Revenue=sum(Revenue),Marginal_Value=min(Marginal_Value),Total_tons=sum(Total_Tons),CIFS_Size=80000),
+CDATA %>% filter(CIFS_150)  %>% group_by(Year) %>% summarize(Revenue=sum(Revenue),Marginal_Value=min(Marginal_Value),Total_tons=sum(Total_Tons),CIFS_Size=150000))
+
+###Graph the waiting value (value added by higher need, vs cost of waiting)
+COST <- 500
+SUMMARY <- SUMMARY %>% group_by(CIFS_Size) %>% mutate(Waiting_Rev=lead(Revenue,1)-Revenue,Waiting_Rev=ifelse(is.na(Waiting_Rev),0,Waiting_Rev),Waiting_Cost=0.05*(Revenue-COST),Waiting_Profit=Waiting_Rev-Waiting_Cost)
+ggsave(ggplot(SUMMARY ,aes(x=Year,y=Waiting_Profit,color=as.factor(CIFS_Size)))+geom_point()+ scale_x_continuous(breaks = seq(1960,2083, by=5))+theme_bw(),filename="test.png",width=20)
+
+ggsave(ggplot(SUMMARY ,aes(x=Year,y=Revenue,color=as.factor(CIFS_Size)))+geom_point()+geom_smooth()+ scale_x_continuous(breaks = seq(1960,2083, by=5)),filename="test.png",width=20)