0

我有两个二值化事件（eventA 和 eventB），我想知道这两个事件是否有任何巧合。所以我将使用新的包 CoinCalc 来研究这两者之间的潜在关系。

library(CoinCalc) #note that the package is not visible (at least for) me in CRAN. I got it from GitHub https://github.com/JonatanSiegmund/CoinCalc

两个二元事件

eventA= c(0,1,0,0,1,1,0,0,1,1,0,0,1,0,1,0,1,0,1,1,1,1,0,0,0,1,1,1,0,0,1,1,0,1,1,0,1,0,0,0,1,1,0,0,0,1,1,0,1,1,1,1,1,1,0,1,0,0,0,1,1,0,0,0,0,0,1,1,0,0,1,1,1,0,0,1,0,1,1,1,0,0,0,1,1,0,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,1,0,0,1,1,0,0,1,1,0,1,1,0,0,0,1,0,0,0,0,0,1,1,0,1,0,1,1,0,1,0,0,0,1,0,0,1,0,1)

eventB = c(0,1,0,0,0,0,1,0,0,0,1,1,1,1,1,0,0,0,0,1,0,1,0,0,0,1,0,1,0,1,0,1,0,1,1,0,1,1,1,0,0,1,1,1,0,0,0,1,1,0,1,1,1,1,1,1,0,1,1,1,0,0,1,0,1,1,1,1,1,1,0,0,1,1,1,0,1,1,1,1,0,1,1,0,0,0,0,0,0,1,1,1,1,0,0,0,1,1,1,0,1,0,0,0,0,0,1,0,1,1,1,1,0,0,0,0,1,0,0,0,1,0,0,0,0,1,1,0,1,0,0,1,0,1,0,1,1,0,0,0)

运行 ECA 分析

ca.out <- CC.eca.ts(eventA, eventB,delT=2,tau=2)

这产生：

$NH precursor

1对

$NH trigger

1错误

$p-value precursor

1 0.2544052

$p-value trigger

1 0.003287963

$precursor coincidence rate

1 0.8243243

$trigger coincidence rate

1 0.9285714

我想确保我正确理解了这一点。根据结果​​，只能拒绝触发器的零假设，在0.003水平上具有统计学意义，符合率是0.92（非常高，这相当于R2吗？）。这是否可以解释为 eventB 对 eventA 的影响很大，而不是相反？

然后我可以使用 CC.plot 函数绘制这两个事件：

CC.plot(eventA,eventB,dates=c(1900:2040),delT=2, tau=2, seriesAname = 'EventA', seriesBname = 'EventB')

产生： CC.plot 输出使用二进制数据

有什么方法可以修改 CC.plot 中的图形参数吗？虚拟年份在该图中不可见。我想更改字体、大小、颜色等。有没有办法通过调用模型输出（ca.out）来绘制相同的图形？

提前致谢！

4

1 回答 1

1

我将尝试回答您的问题：

问题 #1：我在您的示例中看到的最重要的问题是您的事件并非“罕见”。因此，您默认使用的分析显着性检验 ( sigtest="poisson" ) 中最重要的先决条件未得到满足。另一个“问题”是，两个系列中的事件似乎是聚集在一起的（也可能是大量事件的影响）。我建议使用更适合这种情况的sigtest="shuffle.surrogate" 。有关显着性检验的更多信息，请参见 Siegmund 等人。2017 ( http://www.sciencedirect.com/science/article/pii/S0098300416305489) 执行此操作表明，两个重合率都不显着。顺便说一句：在如此多的事件中，您极不可能获得“显着的符合率”，因为随机同时发生的可能性非常高。

然而，如果触发重合率很重要而前兆不重要，那么您的解释是可能的。

问题 #2：情节的问题再次是，事件太多（与最初设计的方法相比）。这就是为什么一切看起来如此凌乱的原因。该功能更像是帮助解释该方法的工作原理以及您所做的工作。如果您仅绘制例如 20 年的数据

CC.plot(eventA[120:140],eventB[120:140],dates=c(2020:2040),delT=2, tau=2, seriesAname = 'EventA', seriesBname = 'EventB')

你会得到一个更好的图像，然而，由于事件密度接近 50%，这不是很好。 CoinCalc 图

目前，没有选项可以更改绘图参数。这可能会出现在该软件包的未来版本中。

我希望这对您有所帮助！

于 2016-11-01T16:21:44.440 回答