This archive contains ========================================================== data files from Kulldorff and Hjalmars (1999) on New Mexico lung cancer cases 1973-1991 Also available (06 Apr 2015) at http://www.satscan.org/datasets/nmlung ========================================================== $$ nm_l_readme.txt descriptions of file layouts $$ nm-km.geo.txt coordinates of New Mexico county seats $$ nm-m-pop.txt county-age-sex specific mid-year population estimates for 1973-1991 $$ lung.cas.txt lung cancer incidence by county,age,sex, and month ========================================================== R programs for the calculations in the paper ========================================================== $$ "1. big sim with New Mexico data.R" Simulates the null distribution of the standardized Knox statistic X* conditional on the number of events (but NOT on the marginal event counts), as in KH 1999. Compares the KH distribution to the true null under three different patterns of marginal space and time effects. $$ "2. Metropolis-Hastings with New Mexico data.R" Simulates the null distribution of X* conditional on marginal event counts, by Metropolis-Hastings sampling of event permutations $$ "3. Metropolis-Hastings COVARIATES with New Mexico data.R" Simulates the null distribution of X* conditional on marginal event counts, by Metropolis-Hastings sampling of event permutations INCLUDING SIMPLE COVARIATE EFFECTS $$ "4. plot KH vs Metropolis with and without covariates.R" Plots the results of the Markov chains from programs 2. and 3. $$ "5. complete enumeration 2x3 tables.R" Calculations for Appendix A ========================================================== output from R programs ========================================================== $$ "Fig1.eps" Plot of results from program 1. $$ "Fig2.eps" Plot of results from program 4. $$ "MH chains.eps" (Supplemental) plot showing the sequence of X* values generated by program 2. Black = standard Knox test (all permutations assumed equally likely) Red = new Knox test (prob. depends on exposure distrib/popn shifts) $$ "MH chains COVARIATES.eps" (Supplemental) plot showing the sequence of X* values generated by program 2. Black = standard Knox test (all permutations assumed equally likely) Red = new Knox test (prob. depends on exposure distrib/popn shifts AND on covariate effects) $$ "big sim with New Mexico data results Thu-24-Jul-2014-1211.dput" detailed output from program 1. $$ "Metropolis-Hastings chains.RData" detailed output from program 2. $$ "Metropolis-Hastings COVARIATES chains.RData" detailed output from program 3. $$ "NewMexicoLungCancer.png" An illustration showing total cases by county, for July 1985 +/-3 months All mapped cases are "close in time". Those that are "close in space" (<100 km between county seats) are connected with solid lines.