get.theta: Generalized version of 'get.theta'

Description Usage Arguments Value Author(s) See Also Examples

View source: R/spatialfuncs.r

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

1
2
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:

θ(d_1,d_2) = \frac{∑ \boldsymbol{1} d_{ij} \in [d_1,d_2)) \boldsymbol{1} (f(i,j)=1) }{∑ ∑ \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.bootstrap, get.theta.typed.permute, get.theta.typed

Other spatialtau: get.pi, get.tau

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
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)

IDSpatialStats documentation built on Nov. 16, 2019, 1:08 a.m.