如何访问R语言fable包中distribution对象的元素？

我正在使用fable包进行时间序列分析，通过以下代码生成了fable类的预测对象my_forecast，想要访问表格第三列中的元素（比如第一行的0.58和0.41），求可行的实现方法。

library(fable)
library(fabletools)
library(feasts)
library(tsibble)
library(dplyr)
library(lubridate)
library(fpp3)

fit <- us_change %>%
  model(trend_model = TSLM(Consumption ~ trend()))

my_forecast <- fit %>%
  forecast(h = "12 years")

my_forecast

解决方案 #

fable生成的预测对象中，目标变量列（这里是Consumption）是嵌套的列表列，每个元素包含预测值、分位数等信息。可以通过以下几种方式提取第三列的元素：

• 方法一：展开嵌套列后提取
  使用unnest()将嵌套的预测列展开为普通列，再直接提取对应值：

  # 展开嵌套的Consumption列
  forecast_df <- my_forecast %>% unnest(Consumption)
  
  # 查看展开后的列名，确认第三列名称
  names(forecast_df)
  
  # 提取第一行第三列的值，示例中第三列假设为`lower 80%`
  forecast_df$`lower 80%`[1]
  # 或用dplyr链式语法
  forecast_df %>% slice(1) %>% pull(`lower 80%`)
  

• 方法二：直接访问嵌套列表元素
  直接通过索引访问嵌套列表中的对应位置：

  # 获取第一行预测结果里的第三列值
  my_forecast$Consumption[[1]][[3]]
  

• 方法三：转换为普通tibble后处理
  将fable对象转为普通tibble，再操作嵌套列：

  forecast_tbl <- as_tibble(my_forecast)
  # 提取第一行第三列的分位数
  forecast_tbl$Consumption[[1]][[3]]
  

注意 #

可以先执行names(my_forecast$Consumption[[1]])查看嵌套列的具体名称，根据实际列名调整提取代码，确保精准获取目标值。

内容的提问来源于stack exchange，提问作者dkent