PlotMCPSP {IndTestPP}

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}

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.

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.

Value

A plot.

See Also

CPSPpoints, PlotICPSP

Examples

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)