2 つの 2 値化されたイベント (eventA と eventB) があります。これら 2 つのイベントに偶然の一致があるかどうかを知りたいです。そこで、新しいパッケージ CoinCalc を使用して、これら 2 つの潜在的な関係を調査します。
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
2 つのバイナリ イベント
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 関数を使用してこれら 2 つのイベントをプロットできます。
CC.plot(eventA,eventB,dates=c(1900:2040),delT=2, tau=2, seriesAname = 'EventA', seriesBname = 'EventB')
どちらが得られますか:
CC.plot のグラフィカル パラメータを変更する方法はありますか? このプロットでは、ダミーの年は表示されません。フォント、サイズ、色などを変更したいのですが、モデル出力(ca.out)を呼び出して同じ図を描く方法はありますか?
前もって感謝します！