library("gamlss2")
## load the data
data("storms", package = "gamlss2")
## yearly observations
plot(counts ~ year, data = storms)
## count distribution
barplot(table(storms$counts))
## NBI model
## model formula including spatial effect
s(year) + s(alt) + s(lon, lat)
## estimate model
b <- gamlss2(f, data = storms, family = NBI)
## estimated effects
plot(b)
## residual diagnostics
plot(b, which = "resid")Severe Storms in Germany
Description
According to the Beaufort scale, severe storms occur from a wind speed of 24.5-28.4 m/s. This dataset contains annual severe storm counts from weather stations in Germany from 1981 to 2021.
Usage
data("storms", package = "gamlss2")Format
A data frame containing 3494 observations on 8 variables.
- id
- Factor, the weather station id.
- county
- Character, the county in Germany where the weather station is located.
- state
- Character, the state in Germany where the weather station is located.
- year
- Integer, the year the observation was measured.
- counts
- Integer, the number of severe storms in this year.
- alt
- Numeric, the altitude in meters above sea level of the weather station.
- lon
- Numeric, the longitude coordinate of the weather station.
- lat
- Numeric, the latitude coordinate of the weather station.
Source
Severe Storms Data:
- Data Source:
- Deutscher Wetterdienst (DWD), Climate Data Center (CDC).
- Licence:
- CC BY 4.0
- URL:
- https://opendata.dwd.de/climate_environment/CDC/
- Coordinate Reference System:
- Longitude/latitude and the WGS84 datum.