DistObs: Calculates the set of close points and the mean distance in a...

In IndTestPP: Tests of Independence and Analysis of Dependence Between Point Processes in Time

Description

Given a set of two or three processes, this function calculates the set of close points and the mean distance for each point in the first process.

It calls the functions calcdist, not intended for the users, and uniongentri.

Usage

DistObs(posx, posy, posz=NULL, info = FALSE, PA = FALSE, procName=c('X','Y','Z'),...)

Arguments

posx	Numeric vector. Position of the occurrence points in the first process.
posy	Numeric vector. Position of the occurrence points in the second process.
posz	Optional. Numeric vector. Position of the occurrence points in the third process.
info	Optional logical flag. If it is TRUE, information about the generated points is showed on the screen and dotcharts and bivariate charts of the occurrence points of the processes are displayed.
PA	Optional logical flag. If it is TRUE, the close point relation is broadened by including the previous and the following points to the overlapping intervals.
procName	Vector of character strings. Labels for the first, second and third processes.
...	Further arguments to be passed to the function plot, if argument info="TRUE".

Details

Given a set of two or three point proccesses, for each point t_{x_i} in the first process of the set, this function calcultes its set of close points and the mean distance to its close points. The definition of set of close points can be found in Abaurrea et al. (2015), and the distances are defined as |t_{y_j}-t_{x_i}| if there are two processes, and as |t_{y_j}-t_{x_i}|+|t_{z_k}-t_{x_i}| if there are three.

Value

DistTri

The vector of the mean distances of points t_{x_i}.

References

Abaurrea, J. Asin, J. and Cebrian, A.C. (2015). A Bootstrap Test of Independence Between Three Temporal Nonhomogeneous Poisson Processes and its Application to Heat Wave Modeling. Environmental and Ecological Statistics, 22(1), 127-144.

See Also

TestIndNH, DistSim, uniongentri

Examples

data(TxBHZ)
dateT<-cbind(TxBHZ$year,TxBHZ$month,TxBHZ$day)  #year, month and day of the month
BivEv<-CPSPPOTevents(N1=TxBHZ$TxH,N2=TxBHZ$TxZ,thres1=36.4, thres2=37.8, date=dateT)
aux<-DistObs(BivEv$Px1, BivEv$Px2, BivEv$Px12,info = TRUE)

IndTestPP documentation built on Aug. 29, 2020, 1:06 a.m.