在传递给 R 中 Arima() 的 xreg 参数之前,我们是否需要对外生变量进行差分?

2024-02-28

我正在尝试在 R 中使用 ARIMAX 构建预测模型,并需要一些关于如何在 xreg 参数中处理协变量的指导。

据我了解,auto.arima 函数在拟合模型(来自训练期数据)时负责协变量的差异,并且我也不需要差异协变量来生成测试期(未来值)的预测。 然而,当使用 R 中的 Arima() 与自定义 (p, d, q) 和 (P, D, Q)[m] 值(d 或 D 大于 0)拟合模型时,我们是否需要手动对协变量? 如果我进行差分,我会得到这样的问题:差分协变量矩阵的长度小于因变量的数据点的数量。

应该如何处理这个问题呢?

  • 我应该按原样发送协变量矩阵吗?即没有差异?
  • 我是否应该进行差分但忽略没有差分协变量数据的前几个观察结果?
  • 我是否应该保留前几行的实际值,其中差异协变量值不可用,而其余行具有差异值?
  • 如果我必须将标志变量 (1/0) 传递给 xreg 矩阵,我是否应该对这些变量进行差分,或者将标志变量的实际值与剩余变量的差分值进行绑定?

另外,在生成未来时期的预测时,如何传递协变量值(按原样或差分后)?

我正在使用以下代码:

ndiff <- ifelse(((pdq_order == "auto") || (PDQ_order == "auto")), ndiffs(ts_train_PowerTransformed), pdq_order[2])
nsdiff <- ifelse(((pdq_order == "auto") || (PDQ_order == "auto")), nsdiffs(ts_train_PowerTransformed), PDQ_order$order[2])

# Creating the appropriate covariates matrix after doing differencing

ifelse(nsdiff >= 1
      , ifelse(ndiff >= 1
                , xreg_differenced <- diff(diff(ts_CovariatesData_TrainingPeriod, lag =  PDQ_order$period, differences = nsdiff),  lag = 1, differences = ndiff)
                , xreg_differenced <- diff(ts_CovariatesData_TrainingPeriod , lag =  PDQ_order$period, differences = nsdiff)
                )
      , ifelse(ndiff >= 1
               , xreg_differenced <- diff( ts_CovariatesData,  lag = 1, differences = ndiff)
               , xreg_differenced <- ts_CovariatesData
 )

# Fitting the model
model_arimax <- Arima(ts_train_PowerTransformed, order = pdq_order, seasonal = PDQ_order, xreg = xreg_differenced)) 

# Generating Forecast for the test period
fit.test <- model_arimax %>% forecast(h=length(ts_test), 
                                              xreg = as.data.frame(diff(diff(ts_CovariatesData_TestPeriod, lag =  PDQ_order$period, differences = nsdiff),  lag = 1, differences = ndiff))

请建议。


Arima将按照顺序和季节性参数中指定的方式对响应变量和 xreg 变量进行差异。您永远不需要自己进行差异化。

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

在传递给 R 中 Arima() 的 xreg 参数之前,我们是否需要对外生变量进行差分? 的相关文章

  • LaTex 中与 knit 和 xtable 交叉引用的问题

    我目前正在与 R Studio 合作 使用 LaTex 中的 R knitr 生成 PDF 文档 在这些文档中 我想在文本中引用的表格中展示我的部分结果 我使用 R 中的 xtable 包生成这些表 它运行良好并为我提供了正确的表 到目前为
  • R::bigmemory - 如何创建角色big.matrix?

    我尝试使用bigmemory封装在R我一开始就陷入困境 我愿意 temp lt matrix paste a 1 10 5 2 并得到一个字符矩阵 没关系 但后来我尝试 x lt as big matrix temp type char 我
  • ggplot2可以在一个图例中分别控制点大小和线大小(线宽)吗?

    一个使用的例子ggplot2绘制数据点组和连接每组均值的线 并使用相同的映射aes for shape并为linetype p lt ggplot mtcars aes gear mpg shape factor cyl linetype
  • 将 Instagram/youtube 嵌入 Shiny R 应用程序

    我想通过点击图表来播放 Instagram 或 Youtube 视频 例如显示异常值等 到目前为止 明确告诉 Shiny 视频内容是有效的 require shiny require ggplot2 data df lt data fram
  • 基于另一个数据集获取数据集的子集

    假设我有一个数据集 即 dat1 ID block plot SPID TotHeight 1 1 1 4 44 5 2 1 1 4 51 3 1 1 4 28 7 4 1 1 4 24 5 5 1 1 4 27 3 6 1 1 4 20
  • R 中具有稳健回归的异常值

    我正在使用lmrobR 中的函数使用robustbase用于稳健回归的库 我会把它用作 rob reg lt lmrob y 0 dat method MM control a1 当我想返回我使用的摘要时summary rob reg 稳健
  • 使用 broom 和 tidyverse 总结 r 平方游戏

    我发布了一个问题here https stackoverflow com questions 48627287 getting adjusted r squared value for each line in a geom smooth
  • 如何用外部图像填充地图边界?

    我正在创建一张带有州边界的巴西地图 这可以直接使用ggplot2 and geom sf 然而 这一次 我不想用数据填充每个状态的颜色 而是想用外部图像 png 填充每个状态的边界 类似于this https online olivet e
  • 计算每个唯一值出现的次数

    假设我有 v rep c 1 2 2 2 25 现在 我想计算每个唯一值出现的次数 unique v 返回唯一值是什么 但不返回它们的数量 gt unique v 1 1 2 我想要一些能给我的东西 length v v 1 1 25 le
  • 在shiny中过滤传单地图数据

    我在用传单地图设置这个闪亮的东西时遇到了麻烦 我的原帖 https stackoverflow com questions 50111566 applying leaflet map bounds to filter data within
  • R foreach问题(某些进程返回NULL)

    我遇到了问题foreach我正在 R 中使用的程序的一部分 该程序用于运行不同参数的模拟 然后将结果返回到单个列表 然后用于生成报告 当并非所有分配的模拟运行都在报告上实际可见时 就会出现问题 从各方面来看 似乎只有分配的运行的一个子集实际
  • 如何删除箱线图上的刻度线

    我试图从箱线图中删除 x 轴刻度线 但保留与刻度线关联的标签 这在基础 R 中可能吗 colors lt c lightskyblue3 gray78 gold1 wheat1 boxplot avgscore module data mi
  • 如何计算R中移动窗口内的平均斜率

    我的数据集包含2个变量y 和 t 05s y 每 05 秒测量一次 我正在尝试计算移动中的平均坡度20秒窗口 即计算第一个 20 秒斜率值后 窗口向前移动一个时间单位 05 秒 并计算下一个 20 秒窗口 在以下位置生成连续 20 秒斜率值
  • 选择 R 中的数据表中隐藏时(在绿色加号下方)列的显示顺序

    Context 使用 DataTables 库制作交互式表格时 当屏幕宽度对于列的数量和宽度来说太窄时 列将隐藏在绿色 号下 我有一个非常宽的表格 有 20 多列 其中一些内容非常冗长 因此某些列在所有屏幕宽度下总是隐藏的 每次隐藏新列时
  • 绘制点之间的所有线

    我有以下 R 代码 x lt c 0 01848598 0 08052353 0 06741172 0 11652034 y lt c 0 4177541 0 4042247 0 3964025 0 4074685 d lt data fr
  • 使用 R 选择第一个非 NA 值

    df lt data frame ID c 1 1 1 2 3 3 3 test c NA 5 5 6 4 NA 7 3 NA 10 9 我想创建一个名为 value 的变量 它是每个单独 ID 测试的第一个非 NA 值 对于只有NA的个体
  • 在 R 中使用 lapply 绘制多个数据帧

    我正在尝试使用 lapply 函数绘制多个数据帧 每个数据帧一个图 但是尽管有关此主题的所有帖子我都找不到答案 因为我不断收到错误 图的输出列表为空 我的数据结构如下 df1 lt mtcars gt group by cyl gt tal
  • 朴素贝叶斯分类器仅基于先验概率做出决策

    我试图根据推文的情绪将推文分为三类 买入 持有 卖出 我正在使用 R 和包 e1071 我有两个数据框 一个训练集和一组需要预测情绪的新推文 训练集数据框 text sentiment this stock is a good buy Bu
  • 将数据框中重叠的范围合并到唯一的组中

    我有一个 n 行 3 的数据框 df lt data frame start c 178 400 983 1932 33653 end c 5025 5025 5535 6918 38197 group c 1 1 2 2 3 df sta
  • 将阴影区域添加到五分位数之间的直方图中

    All 我有一个包含 2 个直方图的图表 其中我还绘制了代表第 20 40 60 和 80 个百分位数的线条 下面的代码使用虚拟数据重现了类似的图表 data lt rbind data frame x rnorm 1000 0 1 g o

随机推荐

  • Yeoman 无法识别全球安装的发电机

    I ran sudo npm install g generator flask它被保存了 我能够通过运行来确认npm list g 然后我跑yo flask 我以为这会起作用 但收到一条错误消息 Error flask You don t
  • Oracle在Java中的RETURNING INTO用法(JDBC、Prepared Statement)

    我正在使用 JDBC 执行 Oracle 语句 如下所示 INSERT INTO MYTABLE MYDATA VALUES RETURNING MY CALCULATED DATA INTO MYTABLE s def makes MY
  • Intellij 完成贡献者

    我正在为 intellij 开发一个插件 我想向基于 xsd 的 xml 编辑器添加自定义建议 到目前为止 我可以从 xsd 文件中获取所需的建议 我已经实现了 xml 的完成贡献者 如下所示 import com intellij cod
  • 将两个元素对齐在同一行

    我正在尝试对齐两个不同的 p 同一行上的元素 我让它工作了 但它完全弄乱了我的标志 他们俩 p p 元素是您在徽标两侧看到的 TEST 对齐的代码是这样的 HTML代码 p div class container div align cen
  • php 解析错误不会显示

    我正在运行 ubuntu 10 04 nginx php fpm 5 4 如果我设置display errors On在 php ini 中 所有错误都会被打印 如果相反 我将其设置为关闭 然后使用 ini set display erro
  • ggplotly 上的主标题和 legend.position 问题

    我在如何在 ggplotly 中定位主标题和图例时遇到问题 我希望我的主标题位于图表顶部并左对齐 我还希望我的图例位于图表的底部中心 这是我的代码 library ggplot2 library dplyr library tidyr li
  • 如何更改bookdown中的图形标题格式

    使用 bookdown 单个文档 时 图会自动编号为 图 1 图标题文本 在化学中 惯例是将主要数字标记为 图1 图标题的文本 对于支持信息文件 图S1 图标题的文本 另外 在文本中的图形参考中 我们需要 如图 1 所示 所以参考文本不应该
  • Windows8 中的首选字体大小需要计算器吗?

    在查看一些官方 Windows8 Metro 材料时 我看到了这行 xaml
  • Ken Burns 效果与 nivo 滑块

    有没有人设置一个尼沃滑块 http nivo dev7studios com 平移每个图像 又名肯伯恩斯效应 http en wikipedia org wiki Ken Burns effect 我正在尝试实现它 但有点棘手 事实上 我的
  • 需要对 Docpad 持久性进行解释

    我对 Docpad 中数据持久化背后的架构感到非常困惑 从博客和论坛中 我了解到内存中 和 或输出目录 用于生成的内容 但Docpad的卖点之一是 完全基于文件 从表面上看 将其托管在 Heroku 或任何临时文件系统上似乎不合逻辑 谁能给
  • JavaScript indexOf 忽略大小写

    我正在尝试查找图像的源名称中是否包含noPic可以是大写或小写 var noPic largeSrc indexOf nopic 我应该写 var noPic largeSrc toLowerCase indexOf nopic 但这个解决
  • Django-[Errno 111]使用 smtp 时连接被拒绝

    我正在开发一个 Django 应用程序 该应用程序具有向用户发送电子邮件的功能 在测试用例中 我有一个联系表单 它使用 smtp 通过 gmail 将表单数据作为电子邮件的一部分提交 我按照建议禁用了验证码该视频播放列表 https www
  • IIS 7.5 应用程序初始化(预热)和 HTTPS

    我们已经使用 IIS 7 5 的应用程序初始化模块有一段时间了 它总是工作得很好 然而 我们刚刚开始实施 SSL 这似乎与热身产生了冲突 我已经做了很多研究 但到目前为止还没有解决方案 基本上问题是初始化模块不遵循重定向 我们必须为网站的某
  • mod_wsgi - 致命 Python 错误:initfsencoding:无法加载文件系统编解码器

    使用 Red Hat apache 2 4 6 worker mpm mod wsgi 4 6 5 和 Python 3 7 当我启动 httpd 时 出现上述错误 ModuleNotFoundError No module named e
  • Webstorm 中的 Sails.js 智能感知

    Webstorm for Sails js 应用程序有可用的智能感知吗 在我的所有控制器中 我收到模型未定义的消息 即使它工作得很好 呼叫服务也是如此 我知道这有点旧 但我最终所做的是首先使用与文件相同的名称 如类名 命名所有 module
  • 返回标识值时的 ExecuteScalar 与 ExecuteNonQuery

    试图弄清楚是否最好使用ExecuteScalar or ExecuteNonQuery如果我想返回新插入行的标识列 我读过了这个问题 https stackoverflow com questions 2974154 what is the
  • 将照片保存到 iOS 8 中的自定义相册

    我在这里需要一点帮助 我有一个方法可以将 UIImage 保存到相机胶卷中 在 iOS 8 中没有问题 该方法如下 PHPhotoLibrary sharedPhotoLibrary performChanges PHAssetChange
  • 如何在 IntelliJ 中永久启用行号?

    如何在 IntelliJ IDEA 中永久启用行号 IntelliJ 14 X 及以上版本 从14 0版本开始 设置对话框的路径略有不同 General子菜单已添加到Editor and 外貌如下所示 IntelliJ 8 1 2 13 X
  • 背景附件:已修复在 iPad 上不起作用

    是否有 CSS Modernizr 方法来了解浏览器是否支持 background attachment fixed 我同时使用背景大小和背景附件 background size cover background attachment fi
  • 在传递给 R 中 Arima() 的 xreg 参数之前,我们是否需要对外生变量进行差分?

    我正在尝试在 R 中使用 ARIMAX 构建预测模型 并需要一些关于如何在 xreg 参数中处理协变量的指导 据我了解 auto arima 函数在拟合模型 来自训练期数据 时负责协变量的差异 并且我也不需要差异协变量来生成测试期 未来值