# get.theta: Generalized version of 'get.theta' In IDSpatialStats: Estimate Global Clustering in Infectious Disease

## 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: for pairs that are (potentially) related for pairs that are unrelated 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

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.