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.