当前位置: 首页>>代码示例 >>用法及示例精选 >>正文


R UKDriverDeaths 1969-84 年英国道路伤亡情况


R语言 UKDriverDeaths 位于 datasets 包(package)。

说明

UKDriverDeaths 是一个时间序列,显示了 1969 年 1 月至 1984 年 12 月期间英国每月死亡或重伤的汽车司机总数。1983 年 1 月 31 日开始强制佩戴安全带。

Seatbelts 是有关同一问题的更多信息。

用法

UKDriverDeaths
Seatbelts

格式

Seatbelts 是一个多时间序列,包含列

DriversKilled

汽车司机被杀。

drivers

UKDriverDeaths 相同。

front

front-seat 乘客死亡或重伤。

rear

rear-seat 乘客死亡或重伤。

kms

行驶距离。

PetrolPrice

汽油价格。

VanKilled

货车(“轻型货车”)司机的数量。

law

0/1:该法律当月生效吗?

例子

require(stats); require(graphics)
## work with pre-seatbelt period to identify a model, use logs
work <- window(log10(UKDriverDeaths), end = 1982+11/12)
par(mfrow = c(3, 1))
plot(work); acf(work); pacf(work)
par(mfrow = c(1, 1))
(fit <- arima(work, c(1, 0, 0), seasonal = list(order = c(1, 0, 0))))
z <- predict(fit, n.ahead = 24)
ts.plot(log10(UKDriverDeaths), z$pred, z$pred+2*z$se, z$pred-2*z$se,
        lty = c(1, 3, 2, 2), col = c("black", "red", "blue", "blue"))

## now see the effect of the explanatory variables
X <- Seatbelts[, c("kms", "PetrolPrice", "law")]
X[, 1] <- log10(X[, 1]) - 4
arima(log10(Seatbelts[, "drivers"]), c(1, 0, 0),
      seasonal = list(order = c(1, 0, 0)), xreg = X)

来源

哈维,A.C. (1989)。预测、结构时间序列模型和卡尔曼滤波器。剑桥大学出版社,第 519-523 页。

Durbin, J. 和 Koopman, S. J. (2001)。通过状态空间方法进行时间序列分析。牛津大学出版社。

参考

Harvey, A. C. and Durbin, J. (1986). The effects of seat belt legislation on British road casualties: A case study in structural time series modelling. Journal of the Royal Statistical Society series A, 149, 187-227. doi:10.2307/2981553.

相关用法


注:本文由纯净天空筛选整理自R-devel大神的英文原创作品 Road Casualties in Great Britain 1969–84。非经特殊声明,原始代码版权归原作者所有,本译文未经允许或授权,请勿转载或复制。