Archiwalne dane pomiarowe do roku 2015 znajdują się na stronie GIOŚ. W plikach .zip znajdują się arkusze kalkulacyjne z danymi dla poszczególnych związków chemicznych. Skrót 1g oznacza dane godzine, 24g - dobowe. Należy uważać z jednostkami: niektóre wyniki pomiarów są wyrażone w miligramach na metr sześcienny, a inne w mikrogramach na metr sześcienny. Wszystkie dane z 2016 roku są w tej drugiej jednostce.
Interesująca nas stacja nazywała się Wrocław Korzeniowskiego. Teraz nazywa się Wrocław Wybrzeże Conrada. Funkcje pobierające dane przyjmują dowolną z tych nazw i w razie potrzeby zmieniają ją na występującą w odpowiednim pliku.
Pakiety, które należy załadować, aby móc korzystać z pakietu to dplyr, tibble, readxl, lubridate i tidyr do pobierania danych oraz ggplot2 i EnvStats do pozostałych funkcji. W pakiecie EnvStats znajduje się funkcja rpareto potrzebna do symulowania symetrycznego rozkładu Pareto.
Funkcja pobierzXLSX służy do pobierania danych dla jednej stacji z surowego pliku ze strony GIOŚ. Informacje można uzyskać w pomocy do funkcji: ?pobierzXLSX. Jej ogólna składnia wygląda następująco.
pobierzXLSX(stacja, substancje, lata, liczbaGodzin = 1, sciezka = getwd(), pominiecie = 3)
Przykładowe wykorzystanie znajduje się poniżej.
pobierzXLSX("DsWrocKorzA", c("CO2", "NO2"), c("2014", "2015"))
Funkcja szuka plików o standardowych nazwach ze strony GIOŚ w folderze zadanym przez argument sciezka. Wynikiem takiego wywołania będzie ramka (dokładniej tibble) ze stężeniami CO2 i NO2 w latach 2014 i 2015. Ramka będzie miała kolumny stacja z nazwą stacji (dzięki czemu dane dla kilku stacji można łatwo pobrać, sklejając wyniki działania pobierzXLSX), dataPomiaru z datą pomiaru w formacie yyyy-mm-dd hh::mm:ss, substancja z nazwą związku i pomiar z wynikiem pomiaru (w miligramach na metr sześcienny). Funkcja pomija domyślnie pierwsze trzy wiersze w pliku .xls, bo znajdują się tam nazwy kolumn i inne informacje.
Funkcja pobierzCSV ma bardzo podobną składnię.
pobierzCSV(stacja, substancje, lata, liczbaGodzin = 1, sciezka = getwd())
Zwraca taką samą ramkę.
Odpowiednie statystyki dobowe (minimum, maksimum i średnią) można łatwo obliczyć, stosując funkcję obliczStatystykiDobowe.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.