R 从分布到置信区间的寓言

我设法用寓言来做预测，然后得到结果

对于如何将此分布更改为80%95%置信区间，我可以提供一些指导吗？谢谢大家!

您可以使用这里的示例代码来获取分发

result <–USAccDeaths %>% as_tsibble %>% 
  model(arima = ARIMA(log(value) ~ pdq(0,1,1) + PDQ(0,1,1)))%>%
  forecast(h=12)

结果%as\u tsibble%>%
模型（arima=arima（对数（值）~pdq（0,1,1）+pdq（0,1,1）））%>%
预测（h=12）

使用

hilo（）

函数可以从预测分布中提取置信区间。它既可以用于分布向量，也可以用于寓言本身

库（tidyverse）
图书馆（寓言）
结果%
模型（arima=arima（对数（值）~pdq（0,1,1）+pdq（0,1,1）））%>%
预测（h=12）
结果%>%
突变（`80%`=hilo（值，80））
#>#寓言：12 x 5[1M]
#>#键：。型号[1]
#>.模型指数值.平均值'80%`
#>                                          
#>1 arima 1979年1月t（N（9，0.0014））8290。[ 7899.082,  8689.169]80
#>2 arima 1979年2月t（N（8.9,0.0018））7453。[ 7055.860,  7859.100]80
#>3阿里玛1979年3月t（N（9，0.0022））8276。[ 7789.719,  8774.054]80
#>4 arima 1979年4月t（N（9.1,0.0025））8584。[ 8036.304,  9144.752]80
#>5阿里玛1979年5月t（N（9.2,0.0029））9499。[ 8849.860, 10166.302]80
#>6阿里玛1979年6月t（N（9.2,0.0033））9900。[ 9180.375, 10639.833]80
#>7阿里玛1979年7月t（N（9.3，0.0037））10988。[10145.473, 11857.038]80
#>8 arima 1979年8月t（N（9.2,0.0041））10132。[ 9315.840, 10974.140]80
#>9 arima 1979年9月t（N（9.1,0.0045））9138。[ 8368.585,  9933.124]80
#>10 arima 1979年10月t（N（9.1,0.0049））9391。[ 8567.874, 10243.615]80
#>11ARIMA1979年11月t（N（9.1,0.0052））8863。[ 8056.754,  9699.824]80
#>12 arima 1979年12月t日（N（9.1,0.0056））9356。[ 8474.732, 10271.739]80
结果%>%
hilo（等级=c（80,95））
#>#一个气泡：12 x 6[1M]
#>#键：。型号[1]
#>.模型指数值.平均值'80%`
#>                                          
#>1 arima 1979年1月t（N（9，0.0014））8290。[ 7899.082,  8689.169]80
#>2 arima 1979年2月t（N（8.9,0.0018））7453。[ 7055.860,  7859.100]80
#>3阿里玛1979年3月t（N（9，0.0022））8276。[ 7789.719,  8774.054]80
#>4 arima 1979年4月t（N（9.1,0.0025））8584。[ 8036.304,  9144.752]80
#>5阿里玛1979年5月t（N（9.2,0.0029））9499。[ 8849.860, 10166.302]80
#>6阿里玛1979年6月t（N（9.2,0.0033））9900。[ 9180.375, 10639.833]80
#>7阿里玛1979年7月t（N（9.3，0.0037））10988。[10145.473, 11857.038]80
#>8 arima 1979年8月t（N（9.2,0.0041））10132。[ 9315.840, 10974.140]80
#>9 arima 1979年9月t（N（9.1,0.0045））9138。[ 8368.585,  9933.124]80
#>10 arima 1979年10月t（N（9.1,0.0049））9391。[ 8567.874, 10243.615]80
#>11ARIMA1979年11月t（N（9.1,0.0052））8863。[ 8056.754,  9699.824]80
#>12 arima 1979年12月t日（N（9.1,0.0056））9356。[ 8474.732, 10271.739]80
#>#…还有一个变量：`95%`

要从

对象中提取数值，可以使用

unpack\u hilo（）

函数

结果%>%
hilo（水平=c（80,95））%>%
打开包装（“80%”）
#>#一个气泡：12 x 7[1M]
#>#键：。型号[1]
#>.模型指数值。平均值'80%'下限'80%'上限`
#>                                     
#>1 arima 1979年1月t（N（9，0.0014））8290。78998689
#>2 arima 1979年2月t（N（8.9,0.0018））7453。70567859
#>3阿里玛1979年3月t（N（9，0.0022））8276。77908774
#>4 arima 1979年4月t（N（9.1,0.0025））8584。80369145
#>5阿里玛1979年5月t（N（9.2,0.0029））9499。885010166
#>6阿里玛1979年6月t（N（9.2,0.0033））9900。918010640
#>7阿里玛1979年7月t（N（9.3，0.0037））10988。1014511857
#>8 arima 1979年8月t（N（9.2,0.0041））10132。931610974
#>9 arima 1979年9月t（N（9.1,0.0045））9138。83699933
#>10 arima 1979年10月t（N（9.1,0.0049））9391。856810244
#>11ARIMA1979年11月t（N（9.1,0.0052））8863。80579700
#>12 arima 1979年12月t日（N（9.1,0.0056））9356。847510272
#>#…还有一个变量：`95%`

---

由（v0.3.0）于2020-04-08创建。

您好，米切尔，谢谢您的输入。我在尝试第二种方法。但仍然如此。我怎样才能把数字中的两个数字7899.082和8689.169从帧中取出？我需要使用这两个数字来画一条线，或者将结果分为两列推送到数据库中。非常感谢。能够使用summary（级别=80）获取对象。但无法以数据帧格式获取数字。创建一个函数。（解包hilo%Unest（）%%>%t（）%%>%matrix（nrow=3）%%>%t（）%%>%as.data.frame（）%%>%rename（lwr=V1，upr=V2，level=V3）}；使用t\p%summary（level=85）%%>%as_tible（）%%>%select（'85%'）%%>%能够获取信息，但它们是“列表”格式，无法将其更改为数据帧，无法将其保存在CSV中。此操作非常有效，谢谢！我以前使用了一个mutate（lower_95=map_dbl（

95%

，~.[1]]）、upper_95=map_dbl（

95%

，~.[2]]）来获得下限和上限，而随着最近的软件包更新，这段代码不再起作用了（也许它从来没有打算这样做）现在获取上下值更容易了，您可以使用

$

提取

下

、

上

和

级别

（带有代码完成提示）。例如，

fbl$hilo$upper

将从

fbl

的

hilo

列中获取

上

绑定。