get.theta {IDSpatialStats}

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:

1. for pairs that are (potentially) related

2. for pairs that are unrelated

3. for pairs that should be ignored all together

Note that names from posmat are not preserved in calls to fun, so the columns of the matrix should be referenced numerically so this is not available to the fun

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)