IndTestPP (version 3.0)

PlotICPSP: Plotting the occurrence points of the indicator processes in a CPSP

Description

This function plots the points in the three indicator processes \(N_{(1)}\), \(N_{(2)}\) and \(N_{(12)}\) of a bivariate Common Poisson shock process (CPSP).

Usage

PlotICPSP(posi1,posi2, posi12, T, date=NULL,axispoints=NULL, ...)

Arguments

posi1

Numeric vector of the points in \(N_{(1)}\)

posi2

Numeric vector of the points in \(N_{(2)}\)

posi12

Numeric vector of the points in \(N_{(12)}\)

T

Numeric value. The length of the observed period of the CPSP.

date

Optional. A vector indicating the date of each observation to be used in the axis of the plot.

axispoints

Optional. Numeric vector with the points in the time index in which axis ticks and labels (from date) will be drawn.

...

Further arguments to be passed to the function plot.

Value

A plot.

Details

A CPSP \(N\) can be decomposed into three independent indicator processes: \(N_{(1)}\), \(N_{(2)}\) and \(N_{(12)}\), the processes of the points occurring only in the first marginal process, only in the second and in both of them (simultaneous points).

The points in the three indicator processes are plotted versus the time index. If one of the arguments date and axispoints is NULL, default axis are used. Otherwise, the values in axispoints are used as the points in the time index in which axis ticks and labels from date are drawn.

See Also

CPSPpoints, PlotMCPSP, PlotMargP

Examples

# NOT RUN {
data(TxBHZ)
T<-length(TxBHZ$TxH)
dateT<-cbind(TxBHZ$year,TxBHZ$month,TxBHZ$day)  #year, month and day of the month
marca<- c(1:length(TxBHZ$TxH))[c(1,diff(dateT[,1]))==1]  #points at first day of the year
BivEv<-CPSPPOTevents(N1=TxBHZ$TxH,N2=TxBHZ$TxZ,thres1=37.8, thres2=36.4, date=dateT, 
                     axispoints=marca)
PlotICPSP(posi1=BivEv$Px1,posi2=BivEv$Px2, posi12=BivEv$Px12, T=T)
PlotICPSP(posi1=BivEv$Px1,posi2=BivEv$Px2, posi12=BivEv$Px12, T=T, date=dateT[,1], 
          axispoints=marca)
# }