PlotMCPSP: Plotting the occurrence points of the marginal processes in a CPSP

Description

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

Usage

PlotMCPSP(pos1,pos2, T, date=NULL, axispoints=NULL, ...)

Arguments

pos1

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

pos2

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

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

The points in the two marginal processes \(N_{1}\), \(N_{2}\) of a bivariate CPSP are plotted versus the time index. The simultaneous points (points of the indicator process \(N_{(12)}\)) are drawn in red.

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.

Examples

Run this code

# 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:T)[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)
PlotMCPSP(pos1=union(BivEv$Px1, BivEv$Px12),pos2=union(BivEv$Px2,BivEv$Px12), T=T)

marca<- c(1:T)[c(1,diff(dateT[,1]))==1]
PlotMCPSP(pos1=union(BivEv$Px1, BivEv$Px12),pos2=union(BivEv$Px2,BivEv$Px12), T=T, 
	date=dateT[,1], axispoints=marca)

# }