Odpływ jednostkowy

Wczytanie i uporządkowanie danych

dane <- read.xlsx("odpływ-jednostkowy.xlsx")

dane$data <- convertToDate(dane$data, format="%Y-%m-%d")

colnames(dane)[2] <- "Odpływ jednostkowy"

Utworzenie nowych kolumn

dane_rh <- mutate(dane, "Rok hydrologiczny" = data)

dane_rh <- filter(dane_rh, dane_rh$`Rok hydrologiczny` > "2017-10-31")

dane_nz <- mutate(dane_rh, "Niżówki zwykłe" = dane_rh$`Odpływ jednostkowy`)

dane_ng <- mutate(dane_rh, "Niżówki głębokie" = dane_rh$`Odpływ jednostkowy`)

dane_all <- mutate(dane, "Rok hydrologiczny" = data,"Niżówki zwykłe" = dane$`Odpływ jednostkowy`,"Niżówki głębokie" = dane$`Odpływ jednostkowy`)

DT::datatable(dane_all, filter ='top', options = list(scrollX=TRUE))

Sortowanie danych

Niżówki zwykłe w czasie roku hydrologicznego

dane_nz_sort <- arrange(dane_nz, dane_nz$`Niżówki zwykłe`)

dane_nz_sort <- filter(dane_nz_sort, dane_nz_sort$`Niżówki zwykłe`<10)

#dane_nz_sort <- arrange(dane_nz, dane_nz$`Rok hydrologiczny`)

nz <- select(dane_nz_sort, "Rok hydrologiczny","Niżówki zwykłe")

DT::datatable(nz, filter ='top', options = list(scrollX=TRUE))

Niżówki głębokie w czasie roku hydrologicznego

dane_ng_sort <- arrange(dane_ng,dane_ng$`Niżówki głębokie`)

dane_ng_sort <- filter(dane_ng_sort, dane_ng_sort$`Niżówki głębokie`<5)

ng <- select(dane_ng_sort, "Rok hydrologiczny","Niżówki głębokie")

DT::datatable(ng, filter ='top', options = list(scrollX=TRUE))

Tworzenie wykresów uzyskanych danych

# Wykres niżówek zwykłych

nz%>%
  ggplot(., mapping = aes(x = `Rok hydrologiczny` ,y = `Niżówki zwykłe`))+
  geom_point(shape = 21,fill = "white", size = 0.5)+
  geom_rug(sides = "b",aes(colour = `Niżówki zwykłe` ))+
  geom_line( colour = "red",alpha = 1/2, size = 0.5)+
  labs(title = "Wykres niżówek zwykłych",
       subtitle = "2017-2018")

# Wykres niżówek głębokich

ng%>%
  ggplot(., mapping = aes(x= `Rok hydrologiczny` ,y = `Niżówki głębokie`))+
  geom_point(shape = 21, size = 0.75)+
  geom_line(colour = "darkgreen",alpha=1/2,size = 0.5 )+
  geom_rug(sides = "b", aes(colour = `Niżówki głębokie` ))+
  labs(title = "Wykres niżówek głębokich",
       subtitle = "2017-2018")

dane_all_rh <- mutate(dane_nz, "Nizówki głębokie" = dane_ng$`Niżówki głębokie`)

dane_all_rh <- select(dane_all_rh, "Rok hydrologiczny", "Odpływ jednostkowy")

dane_all_rh%>%
  filter(.,`Rok hydrologiczny`> "2017-10-31" & `Rok hydrologiczny`< "2017-12-01") -> listopad 
dane_all_rh%>%
  filter(.,`Rok hydrologiczny`> "2017-11-30" & `Rok hydrologiczny`< "2018-01-01") -> grudzien 
dane_all_rh%>%
  filter(.,`Rok hydrologiczny`> "2017-12-31" & `Rok hydrologiczny`< "2018-02-01") -> styczen
dane_all_rh%>%
  filter(.,`Rok hydrologiczny`> "2018-01-31" & `Rok hydrologiczny`< "2018-03-01") -> luty 
dane_all_rh%>%
  filter(.,`Rok hydrologiczny`> "2018-02-28" & `Rok hydrologiczny`< "2018-04-01") -> marzec
dane_all_rh%>%
  filter(.,`Rok hydrologiczny`> "2018-03-31" & `Rok hydrologiczny`< "2018-05-01") -> kwiecien
dane_all_rh%>%
  filter(.,`Rok hydrologiczny`> "2018-04-30" & `Rok hydrologiczny`< "2018-06-01") -> maj
dane_all_rh%>%
  filter(.,`Rok hydrologiczny`> "2018-05-31" & `Rok hydrologiczny`< "2018-07-01") -> czerwiec
dane_all_rh%>%
  filter(.,`Rok hydrologiczny`> "2018-06-30" & `Rok hydrologiczny`< "2018-08-01") -> lipiec
dane_all_rh%>%
  filter(.,`Rok hydrologiczny`> "2018-07-31" & `Rok hydrologiczny`< "2018-09-01") -> sierpien
dane_all_rh%>%
  filter(.,`Rok hydrologiczny`> "2018-08-31" & `Rok hydrologiczny`< "2018-10-01") -> wrzesien
dane_all_rh%>%
  filter(.,`Rok hydrologiczny`> "2018-09-30" & `Rok hydrologiczny`< "2018-11-01") -> pazdziernik

sr_listopad <- mean(listopad$`Odpływ jednostkowy`)
sr_grudzien <- mean(grudzien$`Odpływ jednostkowy`)
sr_styczen <- mean(styczen$`Odpływ jednostkowy`)
sr_luty <- mean(luty$`Odpływ jednostkowy`)
sr_marzec <- mean(marzec$`Odpływ jednostkowy`)
sr_kwiecien <- mean(kwiecien$`Odpływ jednostkowy`)
sr_maj <- mean(maj$`Odpływ jednostkowy`)
sr_czerwiec <- mean(czerwiec$`Odpływ jednostkowy`)
sr_lipiec <- mean(lipiec$`Odpływ jednostkowy`)
sr_sierpien <- mean(sierpien$`Odpływ jednostkowy`)
sr_wrzesien <- mean(wrzesien$`Odpływ jednostkowy`)
sr_pazdziernik <- mean(pazdziernik$`Odpływ jednostkowy`)

rok_hydrologiczny <- cbind(sr_listopad,sr_grudzien,sr_styczen,sr_luty,sr_marzec,sr_kwiecien,sr_maj,sr_czerwiec,sr_lipiec,sr_sierpien,sr_wrzesien,sr_pazdziernik)
as.data.frame(rok_hydrologiczny)

##   sr_listopad sr_grudzien sr_styczen  sr_luty sr_marzec sr_kwiecien   sr_maj
## 1    11.24922    10.60318   8.236757 5.427664   5.65653    9.418655 4.778383
##   sr_czerwiec sr_lipiec sr_sierpien sr_wrzesien sr_pazdziernik
## 1    3.929554  6.415801      7.3652    3.926617       1.933962

barplot(rok_hydrologiczny,
        main = "Wykres średniomiesięcznych wartości odpływu jednostkowego",
        horiz = T,
        las = 1,
        col = "darkred")