DistObs {IndTestPP} | R Documentation |
Calculates the set of close points and the mean distance in a vector of processes, for each point in the first process
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 |
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 |
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)