- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我有 5 个地方 15 天的时间序列数据(以日格式)存储为 matrix
。数据结构为
meter_daywise<-structure(c(24.4745528484842, 21.5936510486629, 58.9120896540103,
49.4188338105575, 568.791971631185, 27.1682608244523, 23.3482757939878,
74.710966227615, 82.6947717673258, 704.212340152625, 23.7581651139442,
21.154634543401, 64.9680107059625, 420.903181621575, 672.629513512841,
128.22871420984, 601.521395359887, 74.6606087800009, 335.87599588534,
576.451039365565, 641.329910104503, 1010.78497435794, 72.6159099850862,
225.153924410613, 582.652388366075, 529.082673064516, 1151.87208010484,
76.9939865858514, 198.567927906582, 641.511944831027, 280.685806121688,
998.647413766557, 73.2033388656998, 337.966543898629, 847.24874747014,
76.7357959402453, 1065.75153722813, 220.286408574643, 301.120955096701,
552.703945876515, 206.496034127105, 1053.49582469841, 206.187963352323,
219.791668265415, 655.496754449233, 172.87981151456, 1018.01514547636,
544.551001017031, 227.116788647859, 656.566145328213, 373.484460701849,
1503.65562864399, 117.732932835236, 251.383369528816, 802.871808716031,
150.471195301885, 1414.88799728991, 14.6490905509617, 203.429955747521,
622.731792495107, 548.093577186778, 1076.5618643676, 15.5135269483705,
256.581499048612, 644.572474965446, 63.2304035656636, 1538.07906461011,
15.0980567507389, 261.513768642083, 622.17970609429, 210.786387991582,
996.998005580537, 15.8138368515615, 157.390773346978, 573.477606081416
), .Dim = c(5L, 15L), .Dimnames = list(c("apFac_401", "apFac_403",
"apFac_501", "apFac_503", "apFac_601"), c("D1", "D2", "D3", "D4",
"D5", "D6", "D7", "D8", "D9", "D10", "D11", "D12", "D13", "D14",
"D15")))
早些时候,我在计算不同系列之间的相关性,使用
library(corrplot)# for plotting correlation matrix
corrplot(cor(t(meter_daywise)),method = "number",type="lower")# have taken transpose of above structure
因此,有了这个,我得到了一个很好的相关矩阵,显示了不同系列之间的相关性。
但是,在观察相关值时我发现有些地方不对劲,在搜索时我发现了这个 link ,其中提到我们需要计算互相关。因此,现在我需要像上面那样计算互相关矩阵。因此,我发现了一些功能,例如
1. ccf() #in base packages
2. diss(meter_daywise,METHOD = "CORT",deltamethod = "DTW")#in TSclust package
我在使用上述功能时遇到了两个问题:
ccf
不将全矩阵作为输入diss()
采用输入矩阵并生成一些矩阵,但在观察值时我发现它不是互相关矩阵,因为值不在 -1
和 1
之间。 所以问题是我们如何计算 R 中不同时间序列值的互相关矩阵?
最佳答案
您可以将 matrix
转换为 ts
对象的 list
,然后使用 do.call
。
meter_daywise <- structure(c(24.4745528484842, 21.5936510486629, 58.9120896540103,
49.4188338105575, 568.791971631185, 27.1682608244523, 23.3482757939878,
74.710966227615, 82.6947717673258, 704.212340152625, 23.7581651139442,
21.154634543401, 64.9680107059625, 420.903181621575, 672.629513512841,
128.22871420984, 601.521395359887, 74.6606087800009, 335.87599588534,
576.451039365565, 641.329910104503, 1010.78497435794, 72.6159099850862,
225.153924410613, 582.652388366075, 529.082673064516, 1151.87208010484,
76.9939865858514, 198.567927906582, 641.511944831027, 280.685806121688,
998.647413766557, 73.2033388656998, 337.966543898629, 847.24874747014,
76.7357959402453, 1065.75153722813, 220.286408574643, 301.120955096701,
552.703945876515, 206.496034127105, 1053.49582469841, 206.187963352323,
219.791668265415, 655.496754449233, 172.87981151456, 1018.01514547636,
544.551001017031, 227.116788647859, 656.566145328213, 373.484460701849,
1503.65562864399, 117.732932835236, 251.383369528816, 802.871808716031,
150.471195301885, 1414.88799728991, 14.6490905509617, 203.429955747521,
622.731792495107, 548.093577186778, 1076.5618643676, 15.5135269483705,
256.581499048612, 644.572474965446, 63.2304035656636, 1538.07906461011,
15.0980567507389, 261.513768642083, 622.17970609429, 210.786387991582,
996.998005580537, 15.8138368515615, 157.390773346978, 573.477606081416
), .Dim = c(5L, 15L), .Dimnames = list(c("apFac_401", "apFac_403",
"apFac_501", "apFac_503", "apFac_601"), c("D1", "D2", "D3", "D4",
"D5", "D6", "D7", "D8", "D9", "D10", "D11", "D12", "D13", "D14",
"D15")))
tss <- unlist(apply(meter_daywise, 1 , function(x) list(ts(x))), recursive = FALSE)
tssu <- do.call(ts.union, tss)
plot(acf(tssu))
关于r - R中不同时间序列数据值的互相关,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/33537687/
您好,我是使用 xampp 的 PHPmyadmin 新手,没有 MYSQL 背景。当我喜欢研究它是如何工作的时,我的脑海中浮现出一个想法,它让我一周都无法休眠,因为我似乎无法弄清楚如何使用 MIN(
Go docs say (强调): Programs using times should typically store and pass them as values, not pointers.
我有一组用户在 8 月 1 日有一个条目。我想找到在 8 月 1 日有条目但在 8 月 2 日没有做任何事情的用户。 现在是 10 月,所以事件已经过去很久了。 我有限的知识说: SELECT * F
我有以下代码,主要编码和取消编码时间结构。这是代码 package main import ( "fmt" "time" "encoding/json" ) type chec
您能详细解释一下“用户 CPU 时间”和“系统 CPU 时间”吗?我读了很多,但我不太理解。 最佳答案 区别在于时间花在用户空间还是内核空间。用户 CPU 时间是处理器运行程序代码(或库中的代码)所花
应用程序不计算东西,但做输入/输出、读取文件、使用网络。我希望探查器显示它。 我希望像 callgrind 中的东西一样,在每个问题中调用 clock_gettime。 或者像 oprofile 那样
目前我的 web 应用程序接收 websocket 数据来触发操作。 这会在页面重新加载时中断,因此我需要一个能够触发特定事件的客户端解决方案。 这个想法可行吗? 假设你有 TimeX = curre
很难说出这里问的是什么。这个问题是含糊的、模糊的、不完整的、过于宽泛的或修辞性的,无法以目前的形式得到合理的回答。如需帮助澄清此问题以便重新打开它,visit the help center 。 已关
我有一个 Instant (org.joda.time.Instant) 的实例,我在一些 api 响应中得到它。我有另一个来自 (java.time.Instant) 的实例,这是我从其他调用中获得
如何集成功能 f(y) w.r.t 时间;即 'y'是一个包含 3000 个值和值 time(t) 的数组从 1 到 3000 不等。所以,在整合 f(y) 后我需要 3000 个值. 积分将是不确定
可以通过 CLI 创建命名空间,但是如何使用 Java SDK 来创建命名空间? 最佳答案 它以编程方式通过 gRPC API 完成由服务公开。 在 Java 中,生成的 gRPC 客户端可以通过 W
我有一个函数,它接受 2 组日期(开始日期和结束日期),这些日期将用于我的匹配引擎 我必须知道start_date1和end_date1是否在start_date2和end_date2内 快进:当我在
我想从 Python 脚本运行“time”unix 命令,以计算非 Python 应用程序的执行时间。我会使用 os.system 方法。有什么方法可以在Python中保存这个输出吗?我的目标是多次运
我正在寻找一种“漂亮的数字”算法来确定日期/时间值轴上的标签。我熟悉 Paul Heckbert's Nice Numbers algorithm . 我有一个在 X 轴上显示时间/日期的图,用户可以
在 PowerShell 中,您可以格式化日期以返回当前小时,如下所示: Get-Date -UFormat %H 您可以像这样在 UTC 中获取日期字符串: $dateNow = Get-Date
我正在尝试使用 Javascript 向父子窗口添加一些页面加载检查功能。 我的目标是“从父窗口”检测,每次子窗口完全加载然后执行一些代码。 我在父窗口中使用以下代码示例: childPage=wi
我正在尝试设置此 FFmpeg 命令的 drawtext 何时开始,我尝试使用 start_number 但看起来它不会成功。 ffmpeg -i 1.mp4 -acodec aac -keyint_
我收到了一个 Excel (2010) 电子表格,它基本上是一个文本转储。 单元格 - J8 具有以下信息 2014 年 2 月 4 日星期二 00:08:06 EST 单元格 - L8 具有以下信息
我收到的原始数据包含一列具有以下日期和时间戳格式的数据: 2014 年 3 月 31 日凌晨 3:38 单元格的格式并不一致,因为有些单元格有单个空格,而另一些单元格中有两个或三个字符之间的空格。所以
我想知道是否有办法在我的 Grails 应用程序顶部显示版本和构建日期。 编辑:我应该说我正在寻找构建应用程序的日期/时间。 最佳答案 在您的主模板中,或任何地方。 Server version:
我是一名优秀的程序员,十分优秀!