In [1]:

    

library("dplyr")
library("readr")


    

    




Attaching package: 'dplyr'

The following objects are masked from 'package:stats':

    filter, lag

The following objects are masked from 'package:base':

    intersect, setdiff, setequal, union

In [2]:

    

medians <- readr::read_csv("NSQD_Res_TSS.csv") %>% 
    dplyr::group_by(parameter, units, season) %>% 
    summarize(median(res))

medians


    

    
Out[2]:





parameter
units
season
median(res)

	
1
Total Suspended Solids
mg/L
autumn
55
	
2
Total Suspended Solids
mg/L
spring
100
	
3
Total Suspended Solids
mg/L
summer
93
	
4
Total Suspended Solids
mg/L
winter
92