setwd("/Users/kathrynnuss/Google Drive/Thesis/MapsTablesGraphs/PredictiveFiles")

install.packages("ggplot2")
install.packages("reshape")

library(ggplot2)
library(readxl)
library(reshape)

ThesisData <- read_excel("ThesisData.xlsx")
View(ThesisData)

attach(ThesisData)

# Creating a Relative Sea Level Curve for the study area 
plot(x = CalendarYearsBP, xlim = c(20000, 1000), y = rsl_m_below_modern, 
     ylim = c(180, 0),
     type = "o", xlab = "Calendar Years Before Present", 
     ylab = "Depth Below Modern Sea Level (m)", col = "blue", xaxt = 'n') 

axis(1, at = c(20000, 19000, 18000, 17000, 16000, 15000, 14000, 13000, 12000,
               11000, 10000,9000, 8000, 7000, 6000, 5000, 4000, 3000, 2000, 1000),
     labels = c(20000, 19000, 18000, 17000, 16000, 15000, 14000, 13000, 12000,
                11000, 10000,9000, 8000, 7000, 6000, 5000, 4000, 3000, 2000, 1000))

plot(x = CalendarYearsBP, xlim = c(20000, 1000), y = streams, ylim = c(20, 90),
     type = "o", xlab = "Calendar Years Before Present", 
     ylab = "Number of Streams", col = "blue", xaxt = 'n') 

axis(1, at = c(20000, 19000, 18000, 17000, 16000, 15000, 14000, 13000, 12000,
               11000, 10000,9000, 8000, 7000, 6000, 5000, 4000, 3000, 2000, 1000),
     labels = c(20000, 19000, 18000, 17000, 16000, 15000, 14000, 13000, 12000,
                11000, 10000,9000, 8000, 7000, 6000, 5000, 4000, 3000, 2000, 1000))
     

plot(x = CalendarYearsBP, xlim = c(20000, 1000), y = stream_length_km, 
     ylim = c(200, 700), type = "o", xlab = "Calendar Years Before Present", 
     ylab = "Length of Streams (km)", col = "blue", xaxt = 'n')

axis(1, at = c(20000, 19000, 18000, 17000, 16000, 15000, 14000, 13000, 12000,
               11000, 10000,9000, 8000, 7000, 6000, 5000, 4000, 3000, 2000, 1000),
     labels = c(20000, 19000, 18000, 17000, 16000, 15000, 14000, 13000, 12000,
                11000, 10000,9000, 8000, 7000, 6000, 5000, 4000, 3000, 2000, 1000))

plot(x = CalendarYearsBP, xlim = c(20000, 1000), y = total_area_kmsq, 
     ylim = c(500, 1650), type = "o", xlab = "Calendar Years Before Present", 
     ylab = "Total Area (km square)", col = "blue", xaxt = 'n')

axis(1, at = c(20000, 19000, 18000, 17000, 16000, 15000, 14000, 13000, 12000,
               11000, 10000,9000, 8000, 7000, 6000, 5000, 4000, 3000, 2000, 1000),
     labels = c(20000, 19000, 18000, 17000, 16000, 15000, 14000, 13000, 12000,
                11000, 10000,9000, 8000, 7000, 6000, 5000, 4000, 3000, 2000, 1000))

detach(ThesisData)