使用 dplyr 和 rollapply 在数据框中滚动预测-6ren

使用 dplyr 和 rollapply 在数据框中滚动预测

转载作者：行者123 更新时间：2023-12-02 03:24:25

我的第一个问题:)

我的目标是:给定一个带有预测变量的数据框(每列一个预测变量/行观察值)使用 lm 拟合回归，然后使用滚动窗口使用最后一个观察值来预测值。

数据框如下所示:

> DfPredictor[1:40,]
           Y            X1           X2                 X3            X4               X5
1          3.2860       192.5115    2.1275              83381         11.4360           8.7440
2          3.2650       190.1462    2.0050              88720         11.4359           8.8971
3          3.2213       192.9773    2.0500              74130         11.4623           8.8380
4          3.1991       193.7058    2.1050              73930         11.3366           8.7536
5          3.2224       193.5407    2.0275              80875         11.3534           8.7555
6          3.2000       190.6049    2.0950              86606         11.3290           8.8555
7          3.1939       191.1390    2.0975              91402         11.2960           8.8433
8          3.1971       192.2921    2.2700              88181         11.2930           8.8681
9          3.1873       194.9700    2.3300             115959          1.9477           8.5245
10         3.2182       194.5396    2.4200             134754         11.3200           8.4990
11         3.2409       194.5396    2.2025             136685          1.9649           8.4192
12         3.2112       195.1362    2.1900             136316          1.9750           8.3752
13         3.2231       193.3560    2.2475             140295          1.9691           8.3546
14         3.2015       192.9649    2.2575             139474          1.9500           8.3116
15         3.1744       194.0154    2.1900             146202          1.8476           8.2225
16         3.1646       194.4423    2.2650             142983          1.8600           8.1948
17         3.1708       194.9473    2.2425             141377          1.8522           8.2589
18         3.1675       193.9788    2.2400             141377          1.8600           8.2600
19         3.1744       194.2563    2.3000             149875          1.8718           8.2899
20         3.1410       193.4316    2.2300             129561          1.8480           8.2395
21         3.1266       191.2633    2.2550             122636          1.8440           8.2396
22         3.1486       192.0354    2.3600             130996          1.8570           8.8640
23         3.1282       194.3351    2.4825              92430          1.7849           8.1291
24         3.1214       193.5196    2.4750              94814          1.7624           8.1991
25         3.1230       193.2017    2.3725              87590          1.7660           8.2310
26         3.1182       192.1642    2.4475              87715          1.6955           8.2414
27         3.1203       191.3744    2.3775              89857          1.6539           8.2480
28         3.1156       192.2646    2.3725              92159          1.5976           8.1676
29         3.1270       192.7555    2.3675              97425          1.5896           8.1162
30         3.1154       194.0375    2.3725              87598          1.5277           8.2640
31         3.1104       192.0596    2.3850              93236          1.5132           7.9999
32         3.0846       192.2792    2.2900              94608          1.4990           8.1600
33         3.0569       193.2573    2.3050              84663          1.4715           8.2200
34         3.0893       192.7632    2.2550              67149          1.4955           7.9590
35         3.0991       192.1229    2.3050              75519          1.4280           7.9183
36         3.0879       192.1229    2.3100              76756          1.3839           7.9133
37         3.0965       192.0502    2.2175              61748          1.3130           7.8750
38         3.0655       191.2274    2.2300              41490          1.2823           7.8656
39         3.0636       191.6342    2.1925              51049          1.1492           7.7447
40         3.1097       190.9312    2.2150              21934          1.1626           7.6895

例如，使用宽度 = 10 的滚动窗口，应估计回归，然后预测与 X1、X2、...、X5 相对应的“Y”。预测应包含在新列“Ypred”中。

有一些方法可以使用 rollapply + lm/predict + mudate 来做到这一点？

非常感谢!!

最佳答案

使用最后注释中的数据并假设在宽度为10的窗口中我们想要预测最后一个Y(即第10个)，那么:

library(zoo)

pred <- function(x) tail(fitted(lm(Y ~., as.data.frame(x))), 1)
transform(DF, pred = rollapplyr(DF, 10, pred, by.column = FALSE, fill = NA))

给予:

        Y       X1     X2     X3      X4     X5     pred
1  3.2860 192.5115 2.1275  83381 11.4360 8.7440       NA
2  3.2650 190.1462 2.0050  88720 11.4359 8.8971       NA
3  3.2213 192.9773 2.0500  74130 11.4623 8.8380       NA
4  3.1991 193.7058 2.1050  73930 11.3366 8.7536       NA
5  3.2224 193.5407 2.0275  80875 11.3534 8.7555       NA
6  3.2000 190.6049 2.0950  86606 11.3290 8.8555       NA
7  3.1939 191.1390 2.0975  91402 11.2960 8.8433       NA
8  3.1971 192.2921 2.2700  88181 11.2930 8.8681       NA
9  3.1873 194.9700 2.3300 115959  1.9477 8.5245       NA
10 3.2182 194.5396 2.4200 134754 11.3200 8.4990 3.219764
11 3.2409 194.5396 2.2025 136685  1.9649 8.4192 3.241614
12 3.2112 195.1362 2.1900 136316  1.9750 8.3752 3.225423
13 3.2231 193.3560 2.2475 140295  1.9691 8.3546 3.217797
14 3.2015 192.9649 2.2575 139474  1.9500 8.3116 3.205856
15 3.1744 194.0154 2.1900 146202  1.8476 8.2225 3.177928
16 3.1646 194.4423 2.2650 142983  1.8600 8.1948 3.156405
17 3.1708 194.9473 2.2425 141377  1.8522 8.2589 3.176243
18 3.1675 193.9788 2.2400 141377  1.8600 8.2600 3.177165
19 3.1744 194.2563 2.3000 149875  1.8718 8.2899 3.177211
20 3.1410 193.4316 2.2300 129561  1.8480 8.2395 3.145533
21 3.1266 191.2633 2.2550 122636  1.8440 8.2396 3.127410
22 3.1486 192.0354 2.3600 130996  1.8570 8.8640 3.148792
23 3.1282 194.3351 2.4825  92430  1.7849 8.1291 3.124913
24 3.1214 193.5196 2.4750  94814  1.7624 8.1991 3.124992
25 3.1230 193.2017 2.3725  87590  1.7660 8.2310 3.117981
26 3.1182 192.1642 2.4475  87715  1.6955 8.2414 3.117679
27 3.1203 191.3744 2.3775  89857  1.6539 8.2480 3.119898
28 3.1156 192.2646 2.3725  92159  1.5976 8.1676 3.121039
29 3.1270 192.7555 2.3675  97425  1.5896 8.1162 3.123903
30 3.1154 194.0375 2.3725  87598  1.5277 8.2640 3.119438
31 3.1104 192.0596 2.3850  93236  1.5132 7.9999 3.113963
32 3.0846 192.2792 2.2900  94608  1.4990 8.1600 3.101229
33 3.0569 193.2573 2.3050  84663  1.4715 8.2200 3.076817
34 3.0893 192.7632 2.2550  67149  1.4955 7.9590 3.083266
35 3.0991 192.1229 2.3050  75519  1.4280 7.9183 3.089377
36 3.0879 192.1229 2.3100  76756  1.3839 7.9133 3.084225
37 3.0965 192.0502 2.2175  61748  1.3130 7.8750 3.075252
38 3.0655 191.2274 2.2300  41490  1.2823 7.8656 3.063025
39 3.0636 191.6342 2.1925  51049  1.1492 7.7447 3.068808
40 3.1097 190.9312 2.2150  21934  1.1626 7.6895 3.091819

注意:以可重现的形式输入 DF 为:

Lines <- "           Y            X1           X2                 X3            X4               X5
1          3.2860       192.5115    2.1275              83381         11.4360           8.7440
2          3.2650       190.1462    2.0050              88720         11.4359           8.8971
3          3.2213       192.9773    2.0500              74130         11.4623           8.8380
4          3.1991       193.7058    2.1050              73930         11.3366           8.7536
5          3.2224       193.5407    2.0275              80875         11.3534           8.7555
6          3.2000       190.6049    2.0950              86606         11.3290           8.8555
7          3.1939       191.1390    2.0975              91402         11.2960           8.8433
8          3.1971       192.2921    2.2700              88181         11.2930           8.8681
9          3.1873       194.9700    2.3300             115959          1.9477           8.5245
10         3.2182       194.5396    2.4200             134754         11.3200           8.4990
11         3.2409       194.5396    2.2025             136685          1.9649           8.4192
12         3.2112       195.1362    2.1900             136316          1.9750           8.3752
13         3.2231       193.3560    2.2475             140295          1.9691           8.3546
14         3.2015       192.9649    2.2575             139474          1.9500           8.3116
15         3.1744       194.0154    2.1900             146202          1.8476           8.2225
16         3.1646       194.4423    2.2650             142983          1.8600           8.1948
17         3.1708       194.9473    2.2425             141377          1.8522           8.2589
18         3.1675       193.9788    2.2400             141377          1.8600           8.2600
19         3.1744       194.2563    2.3000             149875          1.8718           8.2899
20         3.1410       193.4316    2.2300             129561          1.8480           8.2395
21         3.1266       191.2633    2.2550             122636          1.8440           8.2396
22         3.1486       192.0354    2.3600             130996          1.8570           8.8640
23         3.1282       194.3351    2.4825              92430          1.7849           8.1291
24         3.1214       193.5196    2.4750              94814          1.7624           8.1991
25         3.1230       193.2017    2.3725              87590          1.7660           8.2310
26         3.1182       192.1642    2.4475              87715          1.6955           8.2414
27         3.1203       191.3744    2.3775              89857          1.6539           8.2480
28         3.1156       192.2646    2.3725              92159          1.5976           8.1676
29         3.1270       192.7555    2.3675              97425          1.5896           8.1162
30         3.1154       194.0375    2.3725              87598          1.5277           8.2640
31         3.1104       192.0596    2.3850              93236          1.5132           7.9999
32         3.0846       192.2792    2.2900              94608          1.4990           8.1600
33         3.0569       193.2573    2.3050              84663          1.4715           8.2200
34         3.0893       192.7632    2.2550              67149          1.4955           7.9590
35         3.0991       192.1229    2.3050              75519          1.4280           7.9183
36         3.0879       192.1229    2.3100              76756          1.3839           7.9133
37         3.0965       192.0502    2.2175              61748          1.3130           7.8750
38         3.0655       191.2274    2.2300              41490          1.2823           7.8656
39         3.0636       191.6342    2.1925              51049          1.1492           7.7447
40         3.1097       190.9312    2.2150              21934          1.1626           7.6895"

DF <- read.table(text = Lines, header = TRUE)

关于使用 dplyr 和 rollapply 在数据框中滚动预测，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/46694555/

文章推荐： java - 在字节前面添加零

文章推荐： sql-server - 如何在 SQL Server 中使用 NOT IN 创建过滤索引

文章推荐： maven - 无法在启用 LDAP 的情况下使用加密密码部署 Artifact

javascript - 阻止 body 滚动，但阻止触摸设备上的任何父 div 滚动
如果附加了 'not-scroll' 类，我希望我的 body 不滚动，否则它应该正常工作。我已经搜索这个问题两天了，但找不到任何适合我的解决方案。我想要的是向 body 添加一个 class，并
iOS Safari - 滚动 iFrame 导致下面的 div 滚动
我发现似乎是 iOS Safari 中的错误(我正在 iOS 8 上进行测试)。当绝对定位的 iFrame 漂浮在一段可滚动内容上方时，滚动 iFrame 也会滚动下面的内容。以下 HTML (ava
javascript - 如何启用 Safari 滚动 SVG(滚动 div 中的对象标签)？
我有以下代码来显示一系列投资组合图片，这些图片以 SVG 格式存储在滚动 div 中: 在 Safari 中滚动使用两根手指或鼠标滚轮当光标位于 SVG 之一上时不起作用。该页
Javascript 滚动
我想用 javascript 做的是: 一旦你向下滚动页面，将#sidebar-box-fixed 的位置从 position: relative; 更改为定位:固定；。改回position:rela
Elasticsearch 滚动
我对 Elasticsearch 的滚动功能有点困惑。在 elasticsearch 中，每当用户在结果集上滚动时，是否可以每次调用搜索 API？来自文档 "search_type" => "scan
Javascript 滚动
我试图做到这一点，以便当我向上或向下滚动页面时，它会运行不同的相应功能。我发现了一个类似的问题here但我已经尝试了他们的答案并且没有运气。注意:此页面没有正常显示的滚动条。没有地方可以滚动。 bo
gtk - 滚动 GtkDrawingArea
(C语言，GTK库) 在我的表单上，我有一个 GtkDrawingArea 小部件，我在上面使用 Cairo 绘制 GdkPixbufs(从文件加载)。我想要完成的是能够在窗口大小保持固定的情况下使用
SwiftUI 滚动/列表滚动事件
最近我一直在尝试创建一个拉到(刷新，加载更多)swiftUI ScrollView !!，灵感来自 https://cocoapods.org/pods/SwiftPullToRefresh 我正在努
jquery - 拖放 + 滚动
我正在开发一个应用程序，其中有两个带有可放置区域的列表和一个带有可拖动项目的侧面菜单。当我滚动屏幕时，项目的位置困惑。我试图在谷歌上寻找一些东西，最后得到了这个问题:jQuery draggabl
iphone - UIWebView 滚动
我在 UIWebView 中加载了一个 HTML 表单，而我的 UIWebView 恰好从 View 的中间开始并扩展。我必须锁定此 webView 不滚动并将其放在 ScrollView 之上以允许
javascript - 每个元素上的淡入淡出 - 滚动
如何在每个元素而不是整个元素上应用淡入淡出(与其高度相比)？ HTML: CSS: * { padding: 0; margin: 0; box-sizing: border
flutter - 滚动 PageView
我想使用带有垂直轴的 PageView 并使用鼠标滚动在页面之间移动，但是当我使用鼠标滚动时页面不滚动...仅页面单击并向上/向下滑动时滚动。有什么办法吗？我想保留属性 pageSnapping:
java - JTextField 滚动
我制作这个程序是为了好玩，但我被卡住了，因为程序在屏幕外运行。如何在不完全更改代码的情况下实现滚动条。 public static void main(String args[]) throws IO
flutter - 滚动 PageView
我想使用带有垂直轴的 PageView 并使用鼠标滚动在页面之间移动，但是当我使用鼠标滚动时页面不滚动...仅页面单击并向上/向下滑动时滚动。有什么办法吗？我想保留属性 pageSnapping:
java - 滚动()和添加()
很难说出这里问的是什么。这个问题是含糊的、模糊的、不完整的、过于宽泛的或修辞性的，无法以目前的形式得到合理的回答。如需帮助澄清此问题以便重新打开它，visit the help center 。已关
jquery - 滚动/跟随侧边栏帮助
使用 jquery 技术从 css-tricks.com 获得滚动/跟随侧边栏，如果您不知道我在说什么，这里是代码: $(function() { var $sidebar = $
Jquery-手机滑动菜单面板+滚动
我是 jQuery Mobile 新手。我需要向我的应用程序添加 Facebook 滑动面板功能。我经历了 sliding menu panel ，它工作正常，但我在菜单面板中的内容超出了窗口大小，
javascript - 如何在js或jQuery中检测ctrl +滚动？
有没有办法在 js 或 jQuery 或任何其他工具中检测 ctrl + 滚动。我正在尝试执行一些动态布局代码，我需要检测不同分辨率下的屏幕宽度，我通过使用 setTimeout() 的计时器实现了这
jquery 滚动 + 可拖动项目
我有一部分html代码:
c# - 滚动 RichTextBox
我想控制 RichTextBox 滚动，但在控件中找不到任何方法来执行此操作。这样做的原因是我希望当鼠标光标位于 RichTextBox 控件上时鼠标滚轮滚动有效(它没有事件焦点:鼠标滚轮事件由表单

行者123

个人简介

我是一名优秀的程序员,十分优秀！

作者热门文章

滴滴打车优惠券免费领取

全站热门文章

首页

博学

6Ren·AI

商城

使用 dplyr 和 rollapply 在数据框中滚动预测