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


R AirPassengers 1949 年至 1960 年每月航空公司乘客数量


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

说明

经典的 Box & Jenkins 航空公司数据。 1949 年至 1960 年每月国际航空旅客总数。

用法

AirPassengers

格式

每月时间序列,以千为单位。

例子

## Not run: 
## These are quite slow and so not run by example(AirPassengers)

## The classic 'airline model', by full ML
(fit <- arima(log10(AirPassengers), c(0, 1, 1),
              seasonal = list(order = c(0, 1, 1), period = 12)))
update(fit, method = "CSS")
update(fit, x = window(log10(AirPassengers), start = 1954))
pred <- predict(fit, n.ahead = 24)
tl <- pred$pred - 1.96 * pred$se
tu <- pred$pred + 1.96 * pred$se
ts.plot(AirPassengers, 10^tl, 10^tu, log = "y", lty = c(1, 2, 2))

## full ML fit is the same if the series is reversed, CSS fit is not
ap0 <- rev(log10(AirPassengers))
attributes(ap0) <- attributes(AirPassengers)
arima(ap0, c(0, 1, 1), seasonal = list(order = c(0, 1, 1), period = 12))
arima(ap0, c(0, 1, 1), seasonal = list(order = c(0, 1, 1), period = 12),
      method = "CSS")

## Structural Time Series
ap <- log10(AirPassengers) - 2
(fit <- StructTS(ap, type = "BSM"))
par(mfrow = c(1, 2))
plot(cbind(ap, fitted(fit)), plot.type = "single")
plot(cbind(ap, tsSmooth(fit)), plot.type = "single")

## End(Not run)

来源

Box, G. E. P.、Jenkins, G. M. 和 Reinsel, G. C. (1976) 时间序列分析、预测和控制。第三版。 Holden-Day。 G 系列。

相关用法


注:本文由纯净天空筛选整理自R-devel大神的英文原创作品 Monthly Airline Passenger Numbers 1949-1960。非经特殊声明,原始代码版权归原作者所有,本译文未经允许或授权,请勿转载或复制。