Summary of “A Spatial Scan Statistic” by M. Kulldorff

1. Summary of “A Spatial Scan Statistic” by M. Kulldorff Presented by Gauri S. Datta gauri@stat.uga.edu SAMSI September 29, 2005

2. Background • Scan Statistic • A tool to detect cluster in a Point Process • Naus (1965 JASA) studied in one dimension • tests if a 1-dim point process is purely random • Point Process • Consider a time interval [a,b] and a window A=[t,t+w] of fixed width w • (A)= # of e-mails arrived in the time window A • n(A) ´ nA = # of junk e-mails = number of “points” • Arrival times of junk e-mails define a “Point Process”

3. Main Idea in Scan Statistic • Move a window [t,t+w] of size w < b-a over a time interval [a,b] • Over all possible values of t, record the maximum number of points in the window • Compare this number with cut off points under the the hypothesis of a purely Poisson Process

4. Building block of Scan Test • Repeated use of tests to test equality of two Binomial or Poisson populations • Two populations are defined by the scanning window A and its complement Ac • As in multiple comparison, these tests are dependent as one moves the scanning window

5. Spatial Scan Statistic (SSS) • Kulldorff (1997) used SSS to detect clusters in spatial process • SSS can be used • In multi-dim point process • With variable window size • With baseline process an inhomogeneous Poisson process or Bernoulli Process

6. SSS (continued) • Scanning window can be any predefined shape • SSS is on a geographical space G with a measure  • In traditional point process, G is a line,  is a uniform measure • In 2-dim, G is a plane,  a Lebesgue measure

7. More Examples • Forestry: • Spatial clustering of trees. • Want to see for clusters of a specific kind of trees after adjusting fro uneven spatial distribution of all trees • (A)=Total # of trees in region A • nA=# of trees in A of specific kind

8. More Examples (continued) • Epidemiology • Interest in detecting geographical clusters of disease • Need to adjust for uneven population density • Rural vs. urban population • For data aggregated into census districts, measure is concentrated at the central coordinates of districts

9. More Examples (continued) • If interest is in space-time clusters of a disease, the measure will still be concentrated in the geographical region as in the prior example • Adjusting for uneven population distribution is not always enough. Should take confounding factors into account. E.g., in epidemiology measure can reflect standardized expected incidence rate

10. SS = LR statistic • For a fixed size window, scan statistic is the maximum # of points in the window at any given time/geographical region • Test Stat is equivalent to LR test statistic for testing H0:1=2 vs. Ha:1>2 • Generalization to LR test is important for variable window

11. Generalized SS: Notation • G= Geographical area / study space • A= Window ½ G • N(A)= Random # of points in A • A spatial point process • Goal to find a zone Z, the prominent cluster ={Z: Z ½ G} = collection of zones over G

12. Standard Models for SS • Two useful models for point process • (a) Bernoulli model • (b) Poisson model For Bernoulli model, measure  is such that (A) is an integer for all subsets A of G • Two states (disease “point” or no disease) for each unit • Location of the points define a point process

13. More on Bernoulli model • There is exactly one zone Z ½ G such that p=probability of “point” in Z q=probability of “point” in Zc Test H0: p=q vs. Ha: p>q, Z 2 Under H0, n(A)» bin((A), p) for all A ½ G, Under Ha, n(A)» bin((A), p) if A ½ Z; n(A) » bin ((A), q) if A ½ ZC

14. Poisson Model • In Poisson model, points are generated by inhomogeneous Poisson process • n(A)= # of cases in zone A • Exactly one zone Z½ G such that n(A)» Po(p(A Z) + q(A  Zc)) for all A ½ G • H0: p=q vs. Ha: p>q, Z 2 • Under H0, n(A)» Po((A)) for all A, free from Z

15. Choice of Zones • How is  selected? Possibilities: • All circular subsets • All circles centered at anyof several foci on a fixed grid, with a possible upper limit on size • Same as (2) but with a fixed size • All rectangles of fixed size and shape • If looking for space-time clusters, use “cylinders” scanning circular geographical areas over variable time intervals

16. Bernoulli vs. Posson Model • Choice between a Bernoulli or Poisson model does not matter much if n(G) << (G) In other cases, use the model most appropriate for application

17. Likelihood Ratio Test • Bernoulli model Likelihood L(Z,p,q) = pn(z)(1-p)(z)-n(z)£ qn(G)-n(Z)(1-q)(G)-(Z)-n(G)+n(Z) L(Z)= {sup L(Z,p,q): p>q} \hat Z is such that: L(Z^\prime) \le L(\hat Z) for all Z^\prime in \xi

18. LR Test (continued) • L0={sup L(Z,p,q): p=q} = sup pn(G)(1-p)(G)-n(G) free from Z • Likelihood ratio: • = L(\hat Z)/L_0

19. A Useful Result • For LR test for Poisson model, see the paper • An important result on most likely cluster based on these models is given in the paper. It states that as long as the points within the zone constituting the most likely cluster are located where they are, H_0 will be rejected irrespective of the other points in G. If a cluster is located in Seattle, locations of the points in the east coast of U.S. do not matter (Theorem 1)

20. Application of SSS to SIDS • Bernoulli and Poisson models are illustrated using the SIDS data from NC • For 100 counties in NC, total # of live births and # of SIDS cases for 1974-84. • Live births range from 567 to 52345 • Location of county seats are the coordinates. Measure is the # of live births in a county

21. Application to SIDS (continued) • Zones for scanning window are circles centered at a county coordinate point including at most half of the total population • Zones are circular only wrt the aggregated data. As circles around a county seat are drawn, other counties are will either be completely part of a zone or else not at all, depending on whether its county seat is within the circle or not

22. Bernoulli model for SIDS • Bernoulli model is very natural. Each birth can correspond to at most one SID. Table 1 summarizes the results of the analysis. • From Figure 1, the most likely cluster A, consists of Bladen, Columbus, Hoke, Robeson, and Scotland. • Using a conservative test, a secondary cluster is B, consists of Halifax, Hartford and Northampton counties.

23. Poisson model for SIDS • For a rare disease SIDS, Poisson model gives a close approximation to Bernoulli. Results are reported in Table 1 • Both models detect the same cluster • P-values for the primary cluster are same for both the models; p-values for the secondary cluster are very close

24. Application to SIDS (continued)

25. Two significant clusters based on SSS

26. SSS adjusted for Race • For SIDS one useful covariate is race • Race is related to SIDS through unobserved covariates such as quality of housing, access to health care • Overall incidence of SIDS for white children is 1.512 per 1000 and for black children is 2.970 per 1000.

27. SSS: race-adjusted (continued) • Racial distribution differs widely among the counties in NC • This analysis leads to the same primary cluster (see Figure 2) • Previous secondary cluster disappeared but a third secondary cluster C emerges. Cluster C consists of a bunch of counties in the western part of the state

28. Application to SIDS (continued)

29. SSS to SIDS adjusted for race