Dane jakościowe

Katarzyna Smoter

kierunek: geoinformacja

wydział: Geodecji Górniczej i Inżynierii Środowiska

1. Wczytanie danych do R

Za pomocą poniższych poleceń wczytałam dane do R, a także potrzebne przy dalszej pracy pakiety.

library(readxl)
library(tidyverse)
library(lubridate)
library(ggplot2)
library(ggpmisc)

dane <- read_excel("C:/Users/Inni/Desktop/dane.xlsx")

Wczytane dane przedstawiłam w poniższej tabeli.

DT::datatable(dane, caption = "Tabela danych jakościowych", options = list(scrollX=TRUE))

2. Wykres stanu wszystkich zmiennych

Utworzyłam wykres wartości wszytkich zmiennych w zależności od lat i z podziałem na punkty pomiarowo kontrolne(P1,P2,P3,P4).

ggplot(dane, aes(Rok, Wartosc))+
  geom_point()+
  facet_wrap(vars(dane$Nr_pomiaru), scales="free")

3. Przetwarzanie danych

Swoje dane postanowiłam pogrupować ze względu na rodzaje zmiennych i zapisać je do nowych tabel.

BOD5 <- dane %>% filter(dane$Zmienna=="BOD5")
COD <- dane %>% filter(dane$Zmienna=="COD")
Ammonium_Nitrogen<- dane %>% filter(dane$Zmienna=="Ammonium Nitrogen")
Nitrite_Nitrogen<- dane %>% filter(dane$Zmienna=="Nitrite Nitrogen")
Nitrate_Nitrogen<- dane %>% filter(dane$Zmienna=="Nitrate Nitrogen")
Total_Nitrogen<- dane %>% filter(dane$Zmienna=="Total Nitrogen")
Total_Phosphorus<- dane %>% filter(dane$Zmienna=="Total Phosphorus")
Suspended_matter<- dane %>% filter(dane$Zmienna=="Suspended matter")

DT::datatable(BOD5, caption = "Tabela BOD5 ", options = list(scrollX=TRUE))

DT::datatable(COD, caption = "Tabela COD", options = list(scrollX=TRUE))

DT::datatable(Ammonium_Nitrogen, caption = "Tabela Ammonium Nitrogen", options = list(scrollX=TRUE))

DT::datatable(Nitrite_Nitrogen, caption = "Tabela Nitrite Nitrogen", options = list(scrollX=TRUE))

DT::datatable(Nitrate_Nitrogen, caption = "Tabela Nitrate Nitrogen", options = list(scrollX=TRUE))

DT::datatable(Total_Nitrogen, caption = "Tabela Total Nitrogen", options = list(scrollX=TRUE))

DT::datatable(Total_Phosphorus, caption = "Tabela Total Phosphorus", options = list(scrollX=TRUE))

DT::datatable(Suspended_matter, caption = "Tabela Suspended matter", options = list(scrollX=TRUE))

Następnie wartości dla każdej zmiennej z odpowiadających im tabel zapisałam w osobnej tabeli.

dane2 <- BOD5 %>%
  transmute(BOD5=BOD5$Wartosc,
            COD= COD$Wartosc,
            Ammonium_Nitrogen=Ammonium_Nitrogen$Wartosc,
            Nitrite_Nitrogen=Nitrite_Nitrogen$Wartosc,
            Nitrate_Nitrogen=Nitrate_Nitrogen$Wartosc,
            Total_Nitrogen=Total_Nitrogen$Wartosc,
            Total_Phosphorus=Total_Phosphorus$Wartosc,
            Suspended_matter=Suspended_matter$Wartosc)

DT::datatable(dane2, caption = "Tabela wartości zmiennych", options = list(scrollX=TRUE))

4. Podstawowe statystyki

Wykonałam także podstawowe statystki dla każdej zmiennej.

summary(dane2[1:8])

##       BOD5             COD        Ammonium_Nitrogen Nitrite_Nitrogen
##  Min.   : 1.000   Min.   :17.20   Min.   :0.1500    Min.   :0.0300  
##  1st Qu.: 4.000   1st Qu.:21.50   1st Qu.:0.3675    1st Qu.:0.1775  
##  Median : 5.000   Median :22.50   Median :0.4800    Median :0.2450  
##  Mean   : 5.371   Mean   :22.57   Mean   :0.5507    Mean   :0.2547  
##  3rd Qu.: 7.000   3rd Qu.:23.73   3rd Qu.:0.6925    3rd Qu.:0.2925  
##  Max.   :12.000   Max.   :27.10   Max.   :1.1000    Max.   :0.5700  
##  Nitrate_Nitrogen Total_Nitrogen  Total_Phosphorus Suspended_matter
##  Min.   :0.0600   Min.   :1.500   Min.   :0.1400   Min.   :10.00   
##  1st Qu.:0.1200   1st Qu.:1.837   1st Qu.:0.2000   1st Qu.:12.00   
##  Median :0.1750   Median :2.335   Median :0.2750   Median :13.00   
##  Mean   :0.1983   Mean   :2.546   Mean   :0.2614   Mean   :15.32   
##  3rd Qu.:0.2700   3rd Qu.:3.502   3rd Qu.:0.3200   3rd Qu.:16.00   
##  Max.   :0.5400   Max.   :3.730   Max.   :0.3500   Max.   :31.00

5. Wykresy rozrzutu

Utworzyłam wykresy rozrzutu dla każdej zmiennej z każdą.

ggplot(data=dane2, aes(x=dane2$BOD5, y=dane2$COD))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("BOD5") +
  ylab("COD") +
  ggtitle("Wykres rozrzutu BOD5 od COD") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$BOD5, y=dane2$Ammonium_Nitrogen))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("BOD5") +
  ylab("Ammonium Nitrogen") +
  ggtitle("Wykres rozrzutu BOD5 od Ammonium Nitrogen") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$BOD5, y=dane2$Nitrite_Nitrogen))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("BOD5") +
  ylab("Nitrite Nitrogen") +
  ggtitle("Wykres rozrzutu BOD5 od Nitrite Nitrogen") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$BOD5, y=dane2$Nitrate_Nitrogen))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("BOD5") +
  ylab("Nitrate Nitrogen") +
  ggtitle("Wykres rozrzutu BOD5 od Nitrate Nitrogen") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$BOD5, y=dane2$Total_Nitrogen))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("BOD5") +
  ylab("Total Nitrogen") +
  ggtitle("Wykres rozrzutu BOD5 od Total Nitrogen") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$BOD5, y=dane2$Total_Phosphorus))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("BOD5") +
  ylab("Total Phosphorus") +
  ggtitle("Wykres rozrzutu BOD5 od Total Phosphorus") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$BOD5, y=dane2$Suspended_matter))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("BOD5") +
  ylab("Suspended matte") +
  ggtitle("Wykres rozrzutu BOD5 od Suspended matte") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$COD, y=dane2$Ammonium_Nitrogen))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("COD") +
  ylab("Ammonium Nitrogen") +
  ggtitle("Wykres rozrzutu COD od Ammonium Nitrogen") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$COD, y=dane2$Nitrite_Nitrogen))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("COD") +
  ylab("Nitrite Nitrogen") +
  ggtitle("Wykres rozrzutu COD od Nitrite Nitrogen") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$COD, y=dane2$Nitrate_Nitrogen))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("COD") +
  ylab("Nitrate Nitrogen") +
  ggtitle("Wykres rozrzutu COD od Nitrate Nitrogen") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$COD, y=dane2$Total_Nitrogen))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("COD") +
  ylab("Total Nitrogen") +
  ggtitle("Wykres rozrzutu COD od Total Nitrogen") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$COD, y=dane2$Total_Phosphorus))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("COD") +
  ylab("Total Phosphorus") +
  ggtitle("Wykres rozrzutu COD od Total Phosphorus") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$COD, y=dane2$Suspended_matter))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("COD") +
  ylab("Suspended matter") +
  ggtitle("Wykres rozrzutu COD od Suspended matter") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$Ammonium_Nitrogen, y=dane2$Nitrite_Nitrogen))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("Ammonium Nitrogen") +
  ylab("Nitrite Nitrogen") +
  ggtitle("Wykres rozrzutu Ammonium Nitrogen od Nitrite Nitrogen") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$Ammonium_Nitrogen, y=dane2$Nitrate_Nitrogen))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("Ammonium Nitrogen") +
  ylab("Nitrate Nitrogen") +
  ggtitle("Wykres rozrzutu Ammonium Nitrogen od Nitrate Nitrogen") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$Ammonium_Nitrogen, y=dane2$Total_Nitrogen))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("Ammonium Nitrogen") +
  ylab("Total Nitrogen") +
  ggtitle("Wykres rozrzutu Ammonium Nitrogen od Total Nitrogen") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$Ammonium_Nitrogen, y=dane2$Total_Phosphorus))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("Ammonium Nitrogen") +
  ylab("Total Phosphorus") +
  ggtitle("Wykres rozrzutu Ammonium Nitrogen od Total Phosphorus") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$Ammonium_Nitrogen, y=dane2$Suspended_matter))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("Ammonium Nitrogen") +
  ylab("Suspended matter") +
  ggtitle("Wykres rozrzutu Ammonium Nitrogen od Suspended matter") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$Nitrite_Nitrogen, y=dane2$Nitrate_Nitrogen))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("Nitrite Nitrogen") +
  ylab("Nitrate Nitrogen") +
  ggtitle("Wykres rozrzutu Nitrite Nitrogen od Nitrate Nitrogen") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$Nitrite_Nitrogen, y=dane2$Total_Nitrogen))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("Nitrite Nitrogen") +
  ylab("Total Nitrogen") +
  ggtitle("Wykres rozrzutu Nitrite Nitrogen od Total Nitrogen") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$Nitrite_Nitrogen, y=dane2$Total_Phosphorus))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("Nitrite Nitrogen") +
  ylab("Total Phosphorus") +
  ggtitle("Wykres rozrzutu Nitrite Nitrogen od Total Phosphorus") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$Nitrite_Nitrogen, y=dane2$Suspended_matter))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("Nitrite Nitrogen") +
  ylab("Suspended matter") +
  ggtitle("Wykres rozrzutu Nitrite Nitrogen od Suspended matter") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$Nitrate_Nitrogen, y=dane2$Total_Nitrogen))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("Nitrate Nitrogen") +
  ylab("Total Nitrogen") +
  ggtitle("Wykres rozrzutu Nitrate Nitrogen od Total Nitrogen") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$Nitrate_Nitrogen, y=dane2$Total_Phosphorus))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("Nitrate Nitrogen") +
  ylab("Total Phosphorus") +
  ggtitle("Wykres rozrzutu Nitrate Nitrogen od Total Phosphorus") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$Nitrate_Nitrogen, y=dane2$Suspended_matter))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("Nitrate Nitrogen") +
  ylab("Suspended matter") +
  ggtitle("Wykres rozrzutu Nitrate Nitrogen od Suspended matter") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$Total_Nitrogen, y=dane2$Total_Phosphorus))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("Total Nitrogen") +
  ylab("Total Phosphorus") +
  ggtitle("Wykres rozrzutu Total Nitrogen od Total Phosphorus") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$Total_Nitrogen, y=dane2$Suspended_matter))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("Total Nitrogen") +
  ylab("Suspended matter") +
  ggtitle("Wykres rozrzutu Total Nitrogen od Suspended matter") +
  theme_minimal()

ggplot(data=dane2, aes(x=dane2$Total_Phosphorus, y=dane2$Suspended_matter))+
  geom_point(alpha=0.6, color="darkred") +
  geom_smooth(method="lm", formula = y~x-1) +
  geom_abline(intercept=0, col="grey", linetype="dashed", size=1) +
  stat_poly_eq(aes(label = paste(..eq.label.., ..rr.label.., sep = "~~~")), 
               label.x.npc = "right", label.y.npc = 0,
               formula = y~x-1, parse = TRUE, size = 5) +
  xlab("Total Phosphorus") +
  ylab("Suspended matter") +
  ggtitle("Wykres rozrzutu Total Phosphorus od Suspended matter") +
  theme_minimal()