Commit ce35cd60 authored by Jana Ulrich's avatar Jana Ulrich
Browse files

changed filnames of images created by README.Rmd

parent 05fc2f51
......@@ -48,14 +48,14 @@ devtools::install_git("https://gitlab.met.fu-berlin.de/Rpackages/IDF")
Here are a few examples to illustrate the order in which the functions are intended to be used.
* Step 0: sample 20 years of example hourly 'precipitation' data
```{r sample data}
```{r sampledata}
dates <- seq(as.POSIXct("2000-01-01 00:00:00"),as.POSIXct("2019-12-31 23:00:00"),by = 'hour')
sample.precip <- rgamma(n = length(dates), shape = 0.05, rate = 0.4)
precip.df <- data.frame(date=dates,RR=sample.precip)
```
* Step 1: get annual maxima
```{r annual max}
```{r annualmax}
library(IDF)
durations <- 2^(0:6) # accumulation durations [h]
......
......@@ -31,7 +31,7 @@ install.packages("IDF")
or from gitlab using:
``` r
devtools::install_git("https://gitlab.met.fu-berlin.de/Rpackages/idf_package")
devtools::install_git("https://gitlab.met.fu-berlin.de/Rpackages/IDF")
```
## Example
......@@ -62,7 +62,7 @@ ann.max <- IDF.agg(list(precip.df),ds=durations,na.accept = 0.1)
plot(ann.max$ds,ann.max$xdat,log='xy',xlab = 'Duration [h]',ylab='Intensity [mm/h]')
```
<img src="man/figures/README-annual max-1.png" width="100%" />
<img src="man/figures/README-annualmax-1.png" width="100%" />
- Step 2: fit d-GEV to annual maxima
......@@ -74,13 +74,13 @@ fit <- gev.d.fit(xdat = ann.max$xdat,ds = ann.max$ds,sigma0link = make.link('log
#> [1] 0
#>
#> $nllh
#> [1] 73.81918
#> [1] 62.01441
#>
#> $mle
#> [1] 5.384207e+00 6.400833e-01 -9.082818e-03 5.752210e-09 8.153861e-01
#> [1] 5.983501e+00 4.844650e-01 -1.860657e-02 2.126704e-08 7.908172e-01
#>
#> $se
#> [1] 3.531081e-01 7.670363e-02 5.422432e-02 2.000074e-06 1.279425e-02
#> [1] 3.974995e-01 8.085273e-02 8.330635e-02 2.000063e-06 1.230878e-02
# checking the fit
gev.d.diag(fit,pch=1,)
```
......@@ -92,7 +92,7 @@ gev.d.diag(fit,pch=1,)
params <- gev.d.params(fit)
print(params)
#> mut sigma0 xi theta eta
#> 1 5.384207 1.896639 -0.009082818 5.75221e-09 0.8153861
#> 1 5.983501 1.623306 -0.01860657 2.126704e-08 0.7908172
# plotting the probability density for a single duration
q.min <- floor(min(ann.max$xdat[ann.max$ds%in%1:2]))
......
man/figures/README-fit-1.png

43.9 KB | W: | H:

man/figures/README-fit-1.png

45.8 KB | W: | H:

man/figures/README-fit-1.png
man/figures/README-fit-1.png
man/figures/README-fit-1.png
man/figures/README-fit-1.png
  • 2-up
  • Swipe
  • Onion skin
man/figures/README-fit-2.png

22 KB | W: | H:

man/figures/README-fit-2.png

23.5 KB | W: | H:

man/figures/README-fit-2.png
man/figures/README-fit-2.png
man/figures/README-fit-2.png
man/figures/README-fit-2.png
  • 2-up
  • Swipe
  • Onion skin
man/figures/README-idf-1.png

38 KB | W: | H:

man/figures/README-idf-1.png

35.4 KB | W: | H:

man/figures/README-idf-1.png
man/figures/README-idf-1.png
man/figures/README-idf-1.png
man/figures/README-idf-1.png
  • 2-up
  • Swipe
  • Onion skin
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment