get.theta {IDSpatialStats} | R Documentation |
Generalized version of get.theta
Description
Generalized version of the get.theta
function that takes in an arbitrary function and
returns the odds that a point within a particular range of a point of interest shares the relationship
specified by the passed in function with that point.
Usage
get.theta(posmat, fun, r = 1, r.low = rep(0, length(r)), data.frame = TRUE)
Arguments
posmat |
a matrix with columns x, y and any other named columns columns needed by fun |
fun |
a function that takes in two rows of posmat and returns:
Note that names from |
r |
the series of spatial distances (or there maximums) we are interested in |
r.low |
the low end of each range, 0 by default |
data.frame |
logical indicating whether to return results as a data frame (default = TRUE) |
Value
theta value for each distance range that we look at. Where:
\theta(d_1,d_2) = \frac{\sum \boldsymbol{1} d_{ij} \in [d_1,d_2)) \boldsymbol{1} (f(i,j)=1) }{\sum \sum \boldsymbol{1} d_{ij} \in [d_1,d_2)) \boldsymbol{1} (f(i,j)=2) }
Author(s)
Justin Lessler and Henrik Salje
See Also
Other get.theta:
get.theta.bootstrap()
,
get.theta.ci()
,
get.theta.permute()
,
get.theta.typed()
,
get.theta.typed.bootstrap()
,
get.theta.typed.permute()
Other spatialtau:
get.pi()
,
get.tau()
Examples
data(DengueSimR02)
r.max<-seq(20,1000,20)
r.min<-seq(0,980,20)
sero.type.func<-function(a,b,tlimit=20){
if(a[5]==b[5]&(abs(a[3]-b[3])<=tlimit)){rc=1}
else{rc=2}
return(rc)
}
sero.theta<-get.theta(DengueSimR02,sero.type.func,r=r.max,r.low=r.min)