Comprendre l'analyse Coincidence de l'événement du package CoinCalc R

Question

J'ai deux événements binarisées (EVENTA et eventB), je veux savoir s'il y a une coïncidence dans ces deux événements. Donc, je vais utiliser le nouveau paquet CoinCalc pour étudier la relation potentielle entre ces deux.

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 

deux événements binaires

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) 

analyse initiée ECA

ca.out < - CC.eca.ts (Eventa, eventB, DELT = 2, tau = 2)

ce rendement:

$ NH precursor

1 TRUE

NH trigger $

1 FAUX

$ p-value precursor

1 0,2544052

$ p-value trigger

1 0,003287963

$ precursor coincidence rate

1 0,8243243

$ trigger coincidence rate

1 0,9285714

Je veux vous assurer que je comprends cela correctement. D'après les résultats, l'hypothèse nulle ne peut être rejetée pour le déclencheur, ce qui est statistiquement significatif au niveau 0,003, et le taux de coïncidence est 0,92 (très haut, est-ce équivalent à R2?). Cela peut-il être interprété que eventB a une forte influence sur eventA, mais pas le contraire?

Ensuite, je peux tracer ces deux événements en utilisant la fonction CC.plot:

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

Ce qui donne:

Est-il possible de modifier les paramètres graphiques dans CC.plot? Les années fictives ne sont pas visibles dans cette intrigue. Je voudrais changer les polices, taille, couleurs, etc. Est-il possible de tirer la même figure en appelant la sortie du modèle (ca.out)?

Merci d'avance!

Answer 1

Je vais essayer de répondre à vos questions:

Question n ° 1: Le problème le plus important que je vois dans votre exemple est que vos événements ne sont pas « rares ». Par conséquent, la pré-condition la plus importante du test de signification analytique que vous avez utilisé par défaut (sigtest = "poisson") n'est pas remplie. Un autre "problème" est que les événements des deux séries semblent être groupés (cela peut aussi être un effet du nombre élevé d'événements). Je recommande d'utiliser sigtest = "shuffle.surrogate" qui est plus approprié pour ce cas. Plus d'informations sur le test de signification peuvent être trouvées chez Siegmund et al. 2017 (http://www.sciencedirect.com/science/article/pii/S0098300416305489) L'exécution de cette opération révèle que les deux taux de coïncidence ne sont pas significatifs. Soit dit en passant: avec un tel nombre d'événements, il est extrêmement improbable que vous obteniez un «taux de coïncidence significatif», parce que la probabilité que les simultanéités se produisent au hasard est très très élevée. Néanmoins, si le taux de coïncidence de déclenchement est significatif et le précurseur non, votre interprétation est possible. Question n ° 2: Le problème avec le tracé est encore une fois, qu'il y a trop d'événements (par rapport à ce à quoi la méthode a été conçue à l'origine). C'est pourquoi tout semble si désordonné. La fonction était plutôt d'aider à expliquer comment la méthode fonctionne et ce que vous avez fait. Si vous effectuez par ex. seule parcelle, par ex. 20 ans de vos données

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

vous obtiendrez une bien meilleure image que pourtant, en raison de la haute densité d'événements de près de 50%, ne sont pas très agréable. CoinCalc plot

Pour l'instant, il n'y a pas d'options pour modifier les paramètres de tracé. Cela pourrait venir pour une future version du paquet.

J'espère que cela vous aide un peu!