---
title: "Lab2 - Data Exploration and Analysis Laboratory"
output: 
  flexdashboard::flex_dashboard:
      
    orientation: columns
    vertical_layout: fill
    source_code: embed

---

```{r setup, include=FALSE}
library(flexdashboard)
```


Tornado Data
===========================================

Column {data-width=650}
-----------------------------------------------------------------------


### Tornado Trend from 1980 to 2018

```{r}
tornado <- read.csv(url("https://www.spc.noaa.gov/wcm/data/1950-2017_actual_tornadoes.csv"))
##install.packages("usmap")
library(ggplot2)
trend <- as.data.frame(table(tornado$yr))
trend$Var1 <- as.Date(trend$Var1,format="%Y")

trend$ma <- zoo::rollmean(trend$Freq,5,fill=NA)
plotchart <- ggplot(trend,aes(Var1,Freq))+geom_line()+labs(x="Date",y="Frequency Count")+ggtitle("Tornado Trend: 1980-2018")
plotly::ggplotly(plotchart)
```

Column {data-width=350}
-----------------------------------------------------------------------

### Tornado Frequency by State

```{r}

library(usmap)
library(ggplot2)
freq <- as.data.frame(table(tornado$st))
colnames(freq) <- c("state","count")
 plot_usmap(data=freq,values="count",lines="blue")+scale_fill_continuous(low="white",high="blue",name="Count",label=scales::comma)+ggtitle("Tornado Frequency by State: 1950-2018")+theme(legend.position="right",plot.title=element_text(hjust=0.5,size=15,face="bold"))

```

### Tornado Tensity

```{r}

sta <- as.data.frame(table(tornado$yr,tornado$mag),stringsAsFactors=F)
sta$Var1 <- as.integer(sta$Var1)
sta$Decade <- dplyr::case_when(sta$Var1>=1950&sta$Var1<=1959~"1950-1959",sta$Var1>=1960&sta$Var1<=1969~"1960-1969",sta$Var1>=1970&sta$Var1<=1979~"1970-1979",sta$Var1>=1980&sta$Var1<=1989~"1980-1989",sta$Var1>=1990&sta$Var1<=1999~"1990-1999",sta$Var1>=2000&sta$Var1<=2009~"2000-2009",TRUE~"2010-")

p2 <- ggplot(sta,aes(Var2,Freq,fill=Decade))+geom_bar(stat="identity")+scale_fill_brewer(palette="GREEN")+labs(x="(E)F-scale",y="Frequency Count")+ggtitle("Tornado Intensity")
plotly::ggplotly(p2)

```


Precitipation Data
===========================================
Column {data-width=650}
-----------------------------------------------------------------------
### Precitipation comparison between 2009, 2010 and 2011 summer 
```{r}
library("rnoaa")
options(noaakey = "DnQuBjdgsReVpOaRpAbERINIpPODRShS")
out1 <- ncdc(datasetid='GHCND', stationid='GHCND:USC00363698', datatypeid='PRCP',
   startdate = '2009-06-01', enddate = '2009-08-31', limit=500)
out2 <- ncdc(datasetid='GHCND', stationid='GHCND:USC00363698', datatypeid='PRCP',
   startdate = '2010-06-01', enddate = '2010-08-31', limit=500)
out3 <- ncdc(datasetid='GHCND', stationid='GHCND:USC00363698', datatypeid='PRCP',
   startdate = '2011-06-01', enddate = '2011-08-31', limit=500)
ncdc_plot(out1, out2, out3,breaks = "1 month")
```
Column {data-width=350}
-----------------------------------------------------------------------

### Highest Temperature comparison between 2009, 2010 and 2011 summer
```{r}
library("rnoaa")
options(noaakey = "DnQuBjdgsReVpOaRpAbERINIpPODRShS")
out4 <- ncdc(datasetid='GHCND', stationid='GHCND:USC00363698', datatypeid='TMAX',
   startdate = '2009-06-01', enddate = '2009-08-31', limit=500)
out5 <- ncdc(datasetid='GHCND', stationid='GHCND:USC00363698', datatypeid='TMAX',
   startdate = '2010-06-01', enddate = '2010-08-31', limit=500)
out6 <- ncdc(datasetid='GHCND', stationid='GHCND:USC00363698', datatypeid='TMAX',
   startdate = '2011-06-01', enddate = '2011-08-31', limit=500)
ncdc_plot(out4, out5, out6,breaks = "1 month")
```