---

---
title: "PallaviSaitu_ANYL512_Dashboard_Lab2"
author: "Pallavi Saitu"
date: "04/04/2019"
output: 
  flexdashboard::flex_dashboard:
    social: menu
    source_code: embed
---


```{r setup, message=FALSE, warning=FALSE}
library(flexdashboard)
library(dygraphs)
library(plotly)

library(ggplot2)
```




-------------------------------------------------------------------
#Carbon Dioxide Through the Years
```{r message=FALSE, warning=FALSE}

library(readxl)
data1 <- read_excel("/Users/pallavisaitu/Downloads/Climate_Data_2018.xls")
#data1
library(ggplot2)
library(plotly)

#Get information
Carbon_Dioxide = read.table("ftp://aftp.cmdl.noaa.gov/products/trends/co2/co2_mm_mlo.txt", header = FALSE, col.names = c("Year","Month","Decimal_Date","Average","Interpolated","Trend","Number_Days"), skip = 70)

# reformat date
Carbon_Dioxide$Average = replace(Carbon_Dioxide$Average, Carbon_Dioxide$Average == -99.99, NA)
Carbon_Dioxide$Number_Days = replace(Carbon_Dioxide$Days, Carbon_Dioxide $Number_Days == -1, NA)
scaled_Carbon_Dioxide  = Carbon_Dioxide 

# plot
scaled_Carbon_Dioxide_time_series2 = ts(Carbon_Dioxide$Average, frequency = 12, start = c(1958,3))
dygraph(scaled_Carbon_Dioxide_time_series2, main = "Carbon Dioxide Direct Measurements", xlab = "Year", ylab="CO2 (parts per million)") %>%
  dyRangeSelector() %>%
  dyLegend(width = 400, show = "onmouseover") %>%
  dyOptions(drawGrid = FALSE) %>%
  dyOptions(fillGraph = TRUE, fillAlpha = 0.5)

```


-------------------------------------------------------------------
#Harrisburg Precipitation in 2018
```{r message=FALSE, warning=FALSE}

library(devtools)
#install.packages("lawn")
library(lawn)
#install.packages("rnoaa")
library('rnoaa')

options(noaakey = "ouWMpKQfVPZcIlqJpnVxZBfpiTBPEOOP")


out <- ncdc(datasetid='GHCND', stationid='GHCND:USC00363698', datatypeid='PRCP', startdate = '2018-01-01', enddate = '2018-12-31', limit=500)
ncdc_plot(out, breaks="1 month", dateformat="%d/%m")

```

-------------------------------------------------------------------
#Global Temperature Change 1880~2017

```{r message=FALSE, warning=FALSE}
### Global Temperature Change 1880~2017

globalTemp = read.table("https://data.giss.nasa.gov/gistemp/graphs/graph_data/Global_Mean_Estimates_based_on_Land_and_Ocean_Data/graph.txt", header = FALSE, col.names = c("Year","No_Smoothing","Lowess(5)"),skip = 5)

Lowess_smoothing = ts(globalTemp$Lowess.5., frequency = 1, start=c(1880))
Annual_Mean = ts(globalTemp$No_Smoothing, frequency = 1, start=c(1880))
Temperatures <- cbind(Lowess_smoothing, Annual_Mean)

dygraph(Temperatures, main = "Global Land-Ocean Temperature Index", xlab = "Year", ylab="Temperature Anomaly (C)") %>%
dyRangeSelector() %>%
dyLegend(width = 500, show = "onmouseover") %>%
dyOptions(drawGrid = FALSE) %>%
dyOptions(colors = RColorBrewer::brewer.pal(4, "Set1"))

```