multinom.stat: Multinomial Scan Statistic

Description Usage Arguments Details Value Author(s) References Examples

Description

The multinomial scan statistic evaluates the statistic which compares the node attribute within the subgraph with that outside the subgraph while the node attribute follows the multinomial distribution.

Usage

1
multinom.stat(obs, pop = 1, zloc)

Arguments

obs

Numeric vector of observation values.

pop

Numeric vector of population values (default is 1).

zloc

Numeric vector of selected nodes.

Details

A network with interested attributes is denoted as G=(V,E,X), where X=(x_1,…,x_{|V|}) follows a defined distribution. Suppose a subgraph, Z, is selected. Suppose there are k categories in an interested data, the multinomial test statistic is expressed as

λ_A(Z)=∑_k \{ n_{zk} \ln\big(\frac{n_{zk}}{n_{z}}\big) + (n_k-n_{zk}) \ln\big(\frac{n_k-n_{zk}}{n-n_{z}}\big) -n_k \ln \big(\frac{n_k}{n}\big) \},

where n is the total number of observations (nodes), n_{z} is total number of observations in z, n_{zk} is total number of k category in z, and n_{k} is total number of k category in all data.

Value

Three values will be returned. The first value is test statistic. The second is the estimated means which estimated outside the selected nodes. The third is the estimated means estimated within the selected nodes.

Author(s)

Taichi Wang <taichi43@stat.sinica.edu.tw>

References

Jung, I., Kulldorff, M., & Richard, O. J. (2010). A spatial scan statistic for multinomial data. Statistics in medicine, 29(18), 1910–1918.

Examples

1
multinom.stat(obs=rep(1:5,each=10),zloc=1:5)

SNscan documentation built on May 2, 2019, 8:26 a.m.