PHOS graphs scratch

Exploratory graphs, testing first upload of TMDL data 2015.

library(RODBC)
library(reshape2)
library(plyr)
library(lubridate)
library(dplyr)
library(pander)
library(ggplot2)
library(dplyr)
library(knitr)
library(tidyr)

setwd("\\\\psf/Home/Dropbox/MysticDB")
source("./Rcode/Sandbox/Jeff/wq-precip-scripts/load_wq.R")
wq <- load_wq()

## Loading wq database: \\psf\Home\Dropbox\MysticDB\MysticDB_20151224.accdb
## Fetching tables...done
## Merging tables...done
## Excluding field blanks, duplicates...done
## Excluding flagged results...done

source('./Rcode/Sandbox/Jeff/wq-precip-scripts/load_precip.R')
precip <- load_precip()

## Loading precip file: \\psf\Home\Dropbox\MysticDB\Processed\Precip\LoganPrecip.xlsx

wq2 <- append_weather(wq,precip)

## Computing antecedent precip...done
## Computing DateHour column in wq dataframe...done
## Merging wq and precip...done

wq3 <- wq2 %>%
  filter(ProjectID == "TMDL")%>%
  dcast(VisitID + LocationID ~ CharacteristicID, value.var = "ResultValue") 

TURB vs TP

p<- wq3%>%
  ggplot(aes(TURB, TP)) + geom_point() + xlim(c(0,20)) + ylim(c(0,75)) 
p 

By location

p<- wq3%>%
  ggplot(aes(TURB, TP)) + geom_point() + xlim(c(0,20)) + ylim(c(0,75)) + facet_wrap( ~ LocationID)
p 

TP vs CHLA

p<- wq3%>%
  ggplot(aes(TP, CHLA)) + geom_point(aes(color=LocationID), size =3) + xlim(c(0,60)) + ylim(c(0,75)) 
p