---
title: "basicr on data.luftkval.oslo10"
output: rmarkdown::html_vignette
vignette: >
  %\VignetteIndexEntry{luft}
  %\VignetteEngine{knitr::rmarkdown}
  %\VignetteEncoding{UTF-8}
---

```{r setup, include = FALSE}
library(basicr)
knitr::opts_chunk$set(collapse = T, comment = "#>")
```

## Basic analysis on 'data.luftkval.oslo10' df.

### install
```R
library("devtools")
devtools::install_git("https://git.nilu.no/rdiagnosis/basicr.git")
devtools::install_git("https://git.nilu.no/rfriendlyr/friendlyr.git")
devtools::install_git("https://git.nilu.no/rfriendlyr/d.luft.oslyr.git")

```

### Get the 'data.luftkval.oslo10' dataset with friendlyr.

The dataset represents real observations, meaning that it contains errors.
It is possible to know the quality of the value by requesting the qa flags with 'friendlyr::qa()', cf this [vignette](https://git.nilu.no/rfriendlyr/friendlyr/blob/master/vignettes/data.luftkval.oslo10.Rmd)

```R

fromDate <- "2004/01/01 00:00:00"
toDate <- "2018/12/31 23:00:00"
timeResolution <- "hourly"
precision <- "hourly"

element_id <- c("NO2")

# station ids
station_id <- c("7","848","464","827","665","9","11","163","504","809")


# luft
luft <- friendlyr::ts(service="d.luft.oslyr",
                 fromDate=fromDate,toDate=toDate,precision=timeResolution,timeResolution=timeResolution,
                 element_id=element_id,station_id=station_id)

# filter negatives values, used as NA
varnames <- setdiff(names(luft), friendlyts::get_name_date_col(precision))

library(dplyr)
library(tidyr)
tmp <- luft %>% select(names(luft)) %>%
                gather(variable, value, -c(friendlyts::get_name_date_col(precision))) %>%
                mutate(value = replace(value, value < 0, NA)) %>%
                spread(variable, value) %>% friendlyts::as_tbl_friendlyts(precision=precision,date_type="friendlyts")

```

### basics characteristics

```R

basicr::basic_info(df=tmp,target="NO2_7",group="none")
basicr::basic_info(df=tmp,target="NO2_7",group="YEAR")
basicr::basic_info(df=tmp,target="NO2_7",group="MONTH")
basicr::basic_info(df=tmp,target="NO2_7",group="WDAY")
basicr::basic_info(df=tmp,target="NO2_7",group="HOUR")
basicr::basic_info(df=tmp,target="NO2_7",group=c("YEAR","MONTH"))
res <- basicr::basic_info(df=tmp,target="NO2_7",group=c("MONTH","WDAY","HOUR"))

basicr::plot_tbl_basic(tmp,target="NO2_7",group="HOUR",y_name="[ug/m3]")
basicr::plot_tbl_basic(tmp,target="NO2_7",group="WDAY",y_name="[ug/m3]")
basicr::plot_tbl_basic(tmp,target="NO2_7",group="MONTH",y_name="[ug/m3]")

```

### Timeseries of basic characteristics and box-whisker timeseries plot with observation timeseries

```R

res <- basicr::basic_ts(df=tmp,target="NO2_7",group=c("WDAY","HOUR"),precision="hourly")
basicr::plot_boxwhisk_ts(x=res,fromDate="2004/01/05 00:00:00",toDate="2004/01/06 00:00:00",timeResolution="hourly",DESKTOP=TRUE)

```

### L-moment

```R

lmom_ <- basicr::samlmu2(df=tmp,target=varnames)

basicr::plot_lmrd(x=lmom_,DESKTOP=TRUE)


```