- Methodology
- Open Access

# Geographical, temporal and racial disparities in late-stage prostate cancer incidence across Florida: A multiscale joinpoint regression analysis

- Pierre Goovaerts
^{1}Email author and - Hong Xiao
^{2}

**10**:63

https://doi.org/10.1186/1476-072X-10-63

© Goovaerts and Xiao; licensee BioMed Central Ltd. 2011

**Received: **16 September 2011

**Accepted: **5 December 2011

**Published: **5 December 2011

## Abstract

### Background

Although prostate cancer-related incidence and mortality have declined recently, striking racial/ethnic differences persist in the United States. Visualizing and modelling temporal trends of prostate cancer late-stage incidence, and how they vary according to geographic locations and race, should help explaining such disparities. Joinpoint regression is increasingly used to identify the timing and extent of changes in time series of health outcomes. Yet, most analyses of temporal trends are aspatial and conducted at the national level or for a single cancer registry.

### Methods

Time series (1981-2007) of annual proportions of prostate cancer late-stage cases were analyzed for non-Hispanic Whites and non-Hispanic Blacks in each county of Florida. Noise in the data was first filtered by binomial kriging and results were modelled using joinpoint regression. A similar analysis was also conducted at the state level and for groups of metropolitan and non-metropolitan counties. Significant racial differences were detected using tests of parallelism and coincidence of time trends. A new disparity statistic was introduced to measure spatial and temporal changes in the frequency of racial disparities.

### Results

State-level percentage of late-stage diagnosis decreased 50% since 1981; a decline that accelerated in the 90's when Prostate Specific Antigen (PSA) screening was introduced. Analysis at the metropolitan and non-metropolitan levels revealed that the frequency of late-stage diagnosis increased recently in urban areas, and this trend was significant for white males. The annual rate of decrease in late-stage diagnosis and the onset years for significant declines varied greatly among counties and racial groups. Most counties with non-significant average annual percent change (AAPC) were located in the Florida Panhandle for white males, whereas they clustered in South-eastern Florida for black males. The new disparity statistic indicated that the spatial extent of racial disparities reached a peak in 1990 because of an early decline in frequency of late-stage diagnosis observed for black males.

### Conclusions

Analyzing temporal trends in cancer incidence and mortality rates outside a spatial framework is unsatisfactory, since it leads one to overlook significant geographical variation which can potentially generate new insights about the impact of various interventions. Differences observed among nested geographies in Florida show how the modifiable areal unit problem (MAUP) also impacts the analysis of temporal changes.

## Keywords

- Racial Disparity
- Black Male
- Prostate Cancer Mortality
- Modifiable Areal Unit Problem
- Disparity Statistic

## Background

Inequity in overall cancer survival by race is recognized and attributed to differences in the stage at which cancer is diagnosed, its treatment, and, to a lesser extent, in the aggressiveness of tumors. Given equal treatment, there is considerable evidence that African Americans and Whites could experience equal stage-specific survival [1–3]. Still, Whites are diagnosed at earlier stages than African Americans for thirty-one of the thirty-four tumor sites [4]. In particular, although the racial differences in participation in prostate cancer early detection programs are narrowing [5] elderly Blacks are substantially less likely to undergo PSA screening than elderly Whites, a difference that is not completely explained by differences in socioeconomic status and comorbid conditions [6]. Although cancer specialists remain deeply divided over the effectiveness of the PSA blood tests as a diagnostic tool for prostate cancer, some mathematical models projected that 45% to 70% of the observed decline in prostate cancer mortality could be plausibly attributed to the stage shift induced by PSA screening [7]. Other studies found however that racial disparity in PSA testing is probably not a major factor behind current racial differences in prostate cancer mortality rates and declines [8]. As stressed by Nancy Krieger [9] in her paper on social disparities in cancer, "*research is needed to improve monitoring of these disparities. Data on population trends not only reveal whether health inequalities are increasing or decreasing over time, but also stands as critical tests of our etiologic hypotheses. More specifically, if we cannot explain the observed patterns, our knowledge is likely incomplete and our interventions potentially misguided*".

Joinpoint regression [10], also known as piecewise linear regression, is increasingly used to identify the timing and extent of changes in time series of health outcomes, thanks to a public-domain software developed at the US National Cancer Institute, NCI (http://srab.cancer.gov/joinpoint/). The basic idea is to model the time series using a few continuous linear segments. These segments are joined at points called joinpoints which represent the timing (i.e. year) for a statistically significant change in rate trend. The number of joinpoints, as well as the parameters of the piecewise linear regression, are estimated through an iterative procedure that tests whether models of increasing complexity (i.e. including more joinpoints) provide a significantly better goodness-of-fit than simpler models. The tests of significance use a Monte Carlo Permutation method. The approach yields estimates of average annual percent change (AAPC) which allows the summary and comparison of trends over a specified time interval [11]. Two segmented linear regressions (e.g. time trends for two ethnic groups) can also be compared and the parallelisms or identity of the two regression models can be tested [12].

There have been a few applications of joinpoint regression in cancer research, and it is now increasingly used to characterize long-terms trends in cancer mortality in the US [13] and foreign countries [14–16]. Interestingly, the approach was originally applied to prostate cancer incidence and mortality [10]. It helped determine when incidence started rising following the introduction of PSA test and estimate the amount of time diagnosis is advanced due to screening (lead time). On the other hand, a joinpoint model of cancer mortality was used to detect the possible benefits of PSA screening. Tests of comparison of joinpoint regression models were used by the same authors to compare female lung and breast cancer mortality rates between two registry areas and two states, respectively [12]. More recently, a similar methodology was applied to examine how socio-economic absolute and relative disparities in mammography use and associated changes in five breast cancer indicators varied over time [17]. A key difference was the use of rates that were smoothed using Bayesian hierarchical spatio-temporal methods, leading to more stable measures of disparities.

Most analyses of temporal trends have been aspatial and conducted at the National level or for a single cancer registry, with the implicit assumption that the trend parameters are constant across the study area. One exception is the study of temporal trends in breast cancer mortality by state and race from 1975 to 2004 [18]. Other authors also explored temporal trends in geographic disparities and classified 200 counties into priority groups based on changes in breast cancer incidence rates [19]. In a recent study on rates of low birth weight (LBW), joinpoint regression highlighted differences in temporal trends among the five geographical regions of Brazil [20]. Yet, no map was included in these studies, despite the benefit of mapping results of trend analysis.

Conducting a time series analysis through space raises several challenges, such as taking into account the instability of rates computed from smaller local populations, processing large amounts of results generated by multiple applications of joinpoint regression, or the inflated false discovery rate caused by multiple testing. This paper demonstrates how to apply the popular joinpoint regression approach at three different geographically nested levels: State, groups of metropolitan and non-metropolitan counties, and individual counties. Note that the analysis is conducted separately at each of these scales and does not involve a simultaneous evaluation of all three datasets.

The county-level analysis starts with a geostatistical pre-processing of the health data to filter their noise and ends with a post-processing that includes mapping of regression parameters and quantifying racial disparities in both space and time. In the first step, binomial kriging [21] allows one to capitalize on spatial autocorrelation and neighboring geographical units to filter the noise attached to health outcomes and to provide a measure of reliability (kriging variance) for the joinpoint regression. On the other hand, the combination of existing tests of comparison of time trends with multiple testing correction enables the application of joinpoint regression to the detection of geographical areas where time trends for the two races differ significantly. The disparity analysis is performed using a new statistic which measures the number of years where APC confidence intervals did not overlap. Summing up this statistic yearly over all geographical units provides an estimate of how the spatial extent of racial disparities changed with time. The approach is illustrated using times series of annual proportions of prostate cancer late-stage diagnosis available over the period 1981-2007 for each of the 67 counties in Florida.

## Methods

### Data

The data were downloaded from the Florida Cancer Data System website. They included the county-level incidence rates (used only to back-calculate population, see below) and number of cases of prostate cancer with associated stage at diagnosis that were recorded yearly from 1981 through 2007 for non-Hispanic white and black males. Proportions (rates) of late-stage diagnosis were computed over a 3-year moving window to reduce random fluctuations, yielding times series spanning 1982 through 2006. This computation was based only on cases 65 years and over to minimize the impact of disparities in age distribution across Florida and to attenuate the impact of variability in health coverage since all cases were covered by Medicare. The numbers of cases used in this study were: 206,993 white males and 19,442 black males.

Although a county-level analysis might seem rather crude and limits the interpretation of results because of potentially wide heterogeneity within a county, the present study represents a substantial improvement over most analyses of temporal trends which are usually aspatial and conducted at the National level or for a single cancer registry [22, 23]. In addition, county-level analysis allowed the use of a fine temporal resolution (i.e. year) which would not be possible for finer spatial resolutions because of rate instability caused by the small number problem.

*et al*[24] showed that non-metropolitan (non-metro) counties generally had higher death rates and incidence of late-stage disease and lower prevalence of prostate-specific antigen (PSA) screening test (53%) than metropolitan (metro) areas (58%). Following this study, results in this paper were interpreted on the basis of the US Department of Agriculture Rural-Urban Continuum Codes [25], which is most often referred to as the Beale codes, after its creator, Dr. Calvin Beale. This nine-part county codification distinguishes metro counties by the population size of their metro area, and non-metro counties by degree of urbanicity and adjacency to a metro area or areas (Table 1). This information was available for 1983, 1993 and 2003. For 1983 and 1993 codes 0 and 1 were combined to make these classifications comparable to the 2003's codification. These codes were linearly interpolated over the periods 1983-1993 and 1993-2003 to explore relationships between yearly health outcomes and urbanization. Figure 1A shows for each county the average urban code over the period 1981-2007. A metropolitan county is characterized by a Beale code of 3 or lower. The three metropolitan areas of Miami-Fort Lauderdale, Tampa and Orlando have the lowest Beale code. The more rural counties (highest Beale codes) are located in Northern Florida (Panhandle area) and in South Central Florida (Lake Okeechobee).

Definition of 2003 Rural-Urban Continuum Codes (From http://www.ers.usda.gov/data/RuralUrbanContinuumCodes/).

Code | Description |
---|---|

| |

| Counties in metro areas of 1 million population or more |

| Counties in metro areas of 250,000 to 1 million population |

| Counties in metro areas of fewer than 250,000 population |

| |

| Urban population of 20,000 or more, adjacent to a metro area |

| Urban population of 20,000 or more, not adjacent to a metro area |

| Urban population of 2,500 to 19,999, adjacent to a metro area |

| Urban population of 2,500 to 19,999, not adjacent to a metro area |

| Completely rural or less than 2,500 urban population, adjacent to a metro area |

| Completely rural or less than 2,500 urban population, not adjacent to a metro area |

Beale codes can be used as a proxy for provider's accessibility. Health provider information was obtained from the Florida Department of Health Division of Medical Quality Assurance to calculate provider to population ratios. The data included the name and address of each provider (family and internal medicine, urology), the county where they practise, the original date where their license was activated, their license status and expiration date. If the expiration date was missing (case for all urologists), we assumed that the provider was still in practice. The number of primary health providers and urologists was computed for each county and each year over the time period 1981-2007. For the same years, the county-level population size for white and black males was back calculated from the raw incidence rates of prostate cancer. Whenever the rate was zero (i.e. no case diagnosed that year), the population was computed by linear interpolation between years where cases were diagnosed. The provider ratio was then computed for each county and year as the number of primary health providers and urologists active within that county during that year divided by the corresponding yearly county-level population. To visualize the spatial distribution of provider ratios over the State of Florida, these annual county-level ratios were averaged over all 27 years and the results, multiplied by 1,000, are mapped in Figure 1B. Counties with the highest provider ratio are Miami-Dade and Alachua which hosts the main campus of the University of Florida (Gainesville). Both maps (Figures 1A and 1B) display similar patterns. Figure 1C incorporates the time dimension and indicates that the number of providers per 1,000 habitants has steadily increased with time in metropolitan counties, widening the gap between rural and non-rural areas.

### Joinpoint regression

*z*(

*r*;

*t*),

*t*= 1,...,

*T*} be the proportions or rates of late-stage diagnosis recorded for race

*r*at

*T*different time periods (e.g. years). Each observation

*z*(

*r*;

*t*) is computed as the ratio

*d*(

*r*;

*t*)/

*n*(

*r*;

*t*), where

*n*(

*r*;

*t*) is the total number of cases for race

*r*at time

*t*, and

*d*(

*r*;

*t*) is the number of late stages. Joinpoint regression [10] models each time series as a sequence of linear segments (Figure 2). In its log-linear version, the segmented regression model is written as:

where *ε*(*r*;*t*) is the residual for the *t*-th time, and the regression mean *μ*(*r*;*t*) over the entire time interval [*a, b*] is modelled as a succession of (*K*(*r*)+1*)* linear segments (e.g. 3 segments in Figure 2): [*a*,τ_{1}(*r*)] ... (τ_{k}(*r*),τ_{k+1}(*r*)] ... (τ_{K}(*r*), *b*]. The parameter *τ*_{k}(*r*) is the timing (joinpoint) for a statistically significant change in the slopes *β*_{k}(*r*) and *β*_{k+1}(*r*) of two successive segments.

#### Parameter estimation

*K*(

*r*) and values

*τ*

_{k}(

*r*) of the joinpoints, the intercept

*β*

_{0}(

*r*), as well as the slope

*β*

_{k}(

*r*) of each segment. Their estimation is performed in two steps: 1) a grid search method [26] is conducted over the set of possible joinpoints, and 2) at each step of the search the regression parameters and their standard errors are estimated by weighted least-square regression using the following criterion:

The weighting scheme takes into account the fact that the variance of the residuals *ε*(*r*;*t*) typically varies with time (heteroscedasticity) as the number of cases changes. It is also an important issue when assessing racial disparities (see below) since at any given time there are usually much fewer minority cases. These weights are the reciprocal of the variance that can be computed, inside the NCI's Joinpoint regression program, if the dependent variable counts follow a Poisson distribution. The binomial distribution is more appropriate in the present study since late-stage diagnosis is not a rare event, and the weights were thus computed as *n*(*r*;*t*)/[*z*(*r*;*t*)×(1-*z*(*r*;*t*))]. In addition to being heteroscedastic, the random errors in the regression model could be autocorrelated. For example, the average correlation among residuals computed over all 134 time series (i.e. 67 Florida counties and 2 races) is 0.37 for Δt = 1 year, yet only -0.01 for Δt = 2 years. Uncorrelated error models were considered here since these are the only models available in the NCI software for testing the hypothesis of coincidence or parallelisms of different trend models used in the racial disparity analysis.

The number *K*(*r*) of joinpoints is estimated through an iterative procedure that tests whether models of increasing complexity (i.e. including more joinpoints) provide a significantly better goodness-of-fit than simpler models [27]. The tests of significance use a Monte Carlo Permutation procedure described in [10]. A maximum number of joinpoints is typically specified (i.e. *K*_{max} = 3 here) to decrease the number of solutions and the computational time. To avoid that joinpoints get too close together or too close to either end of the time series, a minimum number of observations between joinpoints is also required and was set to 5 in the present application. This minimum number allows the computation of the standard error of the slope parameters, hence the calculation of their confidence intervals and the testing of whether these parameters are significantly different from zero.

#### Temporal trends

*r*over time interval [τ

_{k}(

*r*),

*τ*

_{k+1}(

*r*)] can be quantified by the annual percent change (APC) that is calculated from the slope of the regression model over that time interval as:

Like other regression parameters, confidence intervals can be computed for each APC and one can test whether an APC is significantly different from zero [10].

The trend over the entire time series [*a, b*] can be summarized by the average annual percent change (AAPC) that is the time-weighted average of the APC's from the joinpoint model: the weight of each APC is equal to the relative proportion of the time series [*a, b*] covered by the time interval [τ_{k}(*r*), *τ*_{
k+1
} (*r*)] . This measure is valid even if the joinpoint model indicates that there were changes in trends during those years [11]. Like for the APC, a (1-α) confidence interval can be computed and if it contains zero, then there is no evidence to reject the null hypothesis that the true AAPC is zero at the significance level of α.

#### Racial disparities

*r*and

*r'*can be detected by comparing the models fitted to their corresponding time series {

*z*(

*r*;

*t*),

*t*= 1, ...,

*T*}and {

*z*(

*r'*;

*t*),

*t*= 1, ...,

*T*}. Kim

*et al*. [12] proposed a permutation procedure to compare two segmented line regression functions and to test two types of hypothesis: 1) the two regression models are identical, or 2) the two mean functions are parallel allowing different intercepts. In addition to these global tests, we conducted a finer comparison by computing for every time period

*t*(i.e. year) the 95% confidence intervals of the APC for the two races and counting the number of times these two intervals did not overlap. This new racial disparity statistic can be expressed as:

where the indicator function I(.) = 1 if the following condition on the upper bounds (*U*) and lower bounds (*L*) of the two confidence intervals CI are met: *U*(*r*;*t*) < *L*(*r'*;*t*) or *L*(*r*;*t*) > *U*(*r'*;*t*). A large number indicates that rates of changes for the two races are consistently different over time. There is no statistical test associated with quantity *B*_{
rr'
} which is mainly descriptive.

#### Geographical disparities

*N*geographical units

*v*

_{α}; for example the 67 counties within the State of Florida. For any unit

*v*

_{α}, the regression model (1) is written as:

*K*(

*v*

_{α};

*r*), β

_{0}(

*v*

_{α};

*r*), and{(β

_{k}(

*v*

_{α};

*r*), τ

_{k}(

*v*

_{α};

*r*)),

*k*= 1,...,

*K*(

*v*

_{α};

*r*)}. Similarly, the estimation of temporal trend statistics (i.e. APC, AAPC) and the associated tests of hypothesis can be conducted for every geographical unit. Repeating the analysis across space allows exploring geographical disparities in temporal trends for each race separately but also investigating whether the magnitude of racial disparities changes across Florida. In particular, the racial disparity statistic (Equation 4) can now be computed for each geographical unit

*v*

_{α}as:

Because the tests of hypothesis are conducted for each geographical unit, there is a great likelihood that some tests will turn out significant by chance alone (i.e. false positives), even if the null hypothesis of absence of racial disparity is true in all cases. Multiple testing corrections reduce the significance level applied to each test so that the overall false positive rate is kept to less than or equal to the user-specified significance level α. We used the false discovery rate (FDR) approach which was proven to be less restrictive and more powerful than other approaches, such as the simple Bonferroni correction [28].

*z*(

*v*

_{ α };

*r*;

*t*) and corresponding weight

*ω*(

*v*

_{ α };

*r*;

*t*) in the weighted least-square criterion (Equation 2) by the binomial kriging estimate $\mathit{\u1e91}\left({v}_{\alpha};r;t\right)$ and the inverse of the kriging variance $1\u2215{\sigma}_{BK}^{2}$ (

*v*

_{α};

*r*;

*t*). Binomial kriging was preferred over Poisson kriging [29] to acknowledge the fact that late-stage cancer diagnosis is not a rare event. The kriging estimate, which is a noise-filtered version of the original rate, is computed as a linear combination of the kernel rate

*z*(

*v*

_{ α };

*r*;

*t*) and the rates observed in (

*n*-1) neighboring entities

*v*

_{i}at that time

*t*:

The weights *λ*_{
i
} (*r;t*) assigned to the *n* rates depend on the number of cases diagnosed within each county, the shape and size of the administrative units, as well as the pattern of spatial variability of the rates. These weights are solution of a system of linear equations, known as "binomial kriging" system; see [21, 30] for more details.

#### Temporal disparities

*t*is a measure of the geographical extent of racial disparities that existed at that time. This new geographical disparity statistic can be expressed as:

where the indicator function I(.) = 1 if the following condition on the upper bounds (*U*) and lower bounds (*L*) of the two confidence intervals CI are met: *U*(*v*_{α};*r*;*t*) < *L*(*v*_{α};*r'*;*t*) or *L*(*v*_{α};*r*;*t*) > *U*(*v*_{α};*r'*;*t*). A large number indicates that rates of changes for the two races were consistently different over a large part of Florida at that given point in time.

### Software

Joinpoint regression was conducted using the public-domain Joinpoint Regression Program 3.5.1 July 2011 [10] developed at the US National Cancer Institute, NCI (http://surveillance.cancer.gov/joinpoint/). Binomial kriging and multiple testing correction of *p*-values computed by Joinpoint Regression were performed using the commercial software SpaceStat 2.2 [31]. The three-dimensional display of county-level time series was created using SGeMS (Stanford Geostatistical Modelling Software [32]) 3D visualization panel and FORTRAN programs developed to format the data. All other computations, including the calculation of disparity statistics, were accomplished using FORTRAN programs developed by the first author.

## Results and discussion

### State-level analysis

*T*= 27) for non-Hispanic white and black males in Florida. Joinpoint regression models were fitted to each time series and the parameters are listed in Table 2. For all parameters, the confidence intervals are wider for black males because of the fewer cases: 19,442 black males versus 206,993 white males.

Joinpoint regression analysis of state-level time series of proportion of prostate cancer late-stage diagnosis.

Parameter | White males | Black males | ||||
---|---|---|---|---|---|---|

Estimate | CI | CI | Estimate | CI | CI | |

| ||||||

τ | 1989 | 1988 | 1990 | 1988 | 1986 | 1993 |

τ | 2000 | 1999 | 2002 | 1999 | 1996 | 2002 |

| ||||||

[1981, τ | -0.5 | -1.6 | 0.6 | 0.5 | -3.0 | 4.1 |

[τ | -9.7* | -10.4 | -9.0 | -10.2* | -12.2 | -8.1 |

[τ | 1.8 | -0.2 | 3.8 | -0.2 | -4.3 | 3.9 |

| ||||||

1981-2007 | -3.9* | -4.5 | -3.3 | -4.4* | -6.0 | -2.7 |

Two joinpoints were fitted to each curve and the timing of significant changes for the two races was fairly similar. Both ethnic groups displayed a significant decline in the 90's when PSA screening test was introduced: Annual Percentage Change=-9.69% and -10.20%. The slopes of the two curves are non-significantly different since their confidence intervals overlap. During the two other time periods (late eighties and early 2000), time trends for white and black males were of opposite signs, albeit not significant. Over the entire time period, the average annual rate of change was slighter greater for black (AAPC = -4.4%) than for white (AAPC = -3.9%) although once again the difference was not statistically significant. The racial disparity statistic *B*_{
rr'
} (Equation 4) equalled 2 yrs since the APC confidence intervals did not overlap in 1999 and 2000.

### Urban versus rural areas

The impact of urbanization on temporal trends was explored by grouping, every year, the sixty seven counties based on whether their interpolated Beale index exceeded 3 (non-metropolitan or rural group) or was below 4 (metropolitan or urban group). For both races, the percentage of cases in rural counties decreased with time but these cases tend to be diagnosed at later stages than in more urbanized counties, which confirms results by Jemal *et al*. [24].

*p*-value = 0.461) and the racial disparity statistic equalled 1 since the APC confidence intervals did not overlap in 1989.

Joinpoint regression analysis of proportions of prostate cancer cases diagnosed late in metropolitan and non-metropolitan counties.

Parameter | White males | Black males | ||||
---|---|---|---|---|---|---|

Estimate | CI | CI | Estimate | CI | CI | |

| ||||||

| ||||||

τ | 1989 | 1988 | 1990 | 1988 | 1985 | 1992 |

τ | 2000 | 1999 | 2002 | 2000 | 1999 | 2002 |

| ||||||

[1981, τ | -0.5 | -1.6 | 0.6 | 0.6 | -2.9 | 4.3 |

[τ | -9.9* | -10.7 | -9.2 | -10.2* | -12.0 | -8.3 |

[τ | 2.1* | 0.1 | 4.1 | 1.6 | -3.6 | 6.9 |

| ||||||

1981-2007 | -4.0* | -4.6 | -3.3 | -4.3* | -6.0 | -2.5 |

| ||||||

| ||||||

τ | 1992 | 1989 | 1994 | - | - | - |

τ | 1996 | 1994 | 2003 | - | - | - |

| ||||||

[1981, τ | -1.6* | -2.9 | -0.2 | - | - | - |

[τ | -12.9* | -23.4 | -0.8 | - | - | - |

[τ | -2.1 | -4.3 | 0.1 | - | - | - |

| ||||||

1981-2007 | -3.6* | -5.7 | -1.5 | -5.1* | -6.2 | -3.9 |

Since fewer cases were diagnosed in non-metro counties, their time series were more irregular and confidence intervals for a few parameters widened greatly. No joinpoint was estimated for black males whose time series was fitted using a single exponential curve. For white males, all the linear segments had negative slopes which were significantly different from zero over the first two time periods. This continuous decline in non-metro counties was a significant departure from the recent increase observed in metropolitan counties. This was confirmed by the rejection at α = 0.05 of the hypothesis of parallelism of joinpoint regression models fitted to rates for white males in metro and non-metro counties. This hypothesis was not rejected for any comparison involving black males because of the uncertainty caused by smaller population sizes. The racial disparity statistic (Equation 4) was much larger in rural counties (*B*_{
rr'
} = 12 yrs) than in urban counties (*B*_{
rr'
} = 1 yr).

### County-level analysis

The above analysis of temporal trends in metro and non-metro areas revealed the existence of geographical and racial disparities that could not be detected at the State level. An even finer spatial analysis could be performed by examining the county-level time series of proportions of late-stage diagnosis. Such an analysis needed however to be conducted after smoothing using binomial kriging because of the larger rate instability observed when moving to smaller geographical units. Another reason for the application of smoothing techniques was the existence of missing values (i.e. years where no case was diagnosed) that needed to be replaced by rate estimates in order to run NCI joinpoint regression program. Missing values were only observed for black males and represented 6.27% of all rate-years. Binomial kriging was conducted using as neighbours the counties that share a common border or vertex with the county being smoothed (1^{st} order Queen's adjacency). A population-weighted variogram of percentage of late-stage diagnosis was computed every year for both races. The fitted variogram model had on average a longer range of autocorrelation for white males (110 km) than for black males (85.5 km).

#### Visualization of geographical and temporal disparities

#### Analysis of heavily populated counties

Racial disparities within these same geographical units were investigated by conducting a joinpoint regression analysis and testing the hypothesis of parallelism of regression models for both races. Figure 7 showcased some of the geographical and racial disparities observed over Florida. Although all geographical units and races experienced a significant decline in percentage of late-stage diagnosis during the nineties, joinpoint regression models differed greatly. For some counties, like Alachua or Duval County, this decline followed a significant increase in the eighties while for other counties the rates remained stable or slightly decreased during that time period. Similarly, a wide spectrum of temporal trends was observed for the most recent years; for example for white males a significant decline in Escambia County contrasted with a significant increase in Miami-Dade County. The hypothesis of parallelism was rejected at α = 0.05 for all six counties and two groups of rural counties. The use of the disparity statistic *B*_{
rr'
} (Equation 6) allowed discriminating the different geographical units based on the frequency of racial disparities in annual percent changes (APC). Figure 7 indicates that disparities were the smallest for Miami-Dade County (*B*_{
rr'
} = 2 yrs) where both curves had similar joinpoints and slopes of similar sign. Racial disparities in temporal trends are much larger for Escambia and Orange Counties (*B*_{
rr'
} = 16 yrs), and reached a maximum for the group of very rural counties (*B*_{
rr'
} = 25 yrs). In the later case no joinpoint was estimated for black males whose time series was fitted using a single exponential curve with significant APC. For white males, the decline was only significant between 1991 and 1997. Note that since the racial disparities in APC mostly (19 years out of 25 years) took the form of larger declines for black males compared to white males, the racial disparities in percentage of late-stage diagnosis at the end of the time period are actually one of the smallest among all geographical units.

#### Analysis of all Florida counties

For white males, most counties with non-significant AAPC were located in the Florida Panhandle, which explained the largest proportion of late-stage diagnosis observed on average over the entire time period. For black males, counties with non-significant AAPC clustered in southeast Florida. Both races experienced the largest overall decline in percentage of late-stage diagnosis in Miami-Dade County and Central-west Florida (Tampa Bay); this decline started earlier for Miami-Dade County, in particular for white males (Figure 8B-D). Among the three metropolitan areas with the lowest Beale codes (Figure 1A), Orlando area (Orange county) showed the least progress in lowering the frequency of late-stage diagnosis, yet the initial proportion of late-stage cases was smaller in that county. Surprisingly, Alachua County (University of Florida) that has the second largest provider ratio in the State (Figure 1B) has a very low AAPC, which is not significantly different from zero for black males. The Florida Panhandle encompassed large differences in onset years for white males: significant declines in frequency of late-stage diagnosis started much later in the Big Bend region than in the ten counties west of it.

#### Racial disparities

_{rr'}statistic equalled 12.8 yrs. This result confirms the conclusions drawn from Figure 4 where time series for non metro counties had much different shapes (B

_{rr'}= 12 yrs) than metro counties (B

_{rr'}= 1 yr).

### Sensitivity analysis

The second comparison criteria were the two summary statistics mapped in Figure 8: the average annual percent change (AAPC) and the joinpoint corresponding to the first significant decline in proportion of late-stage diagnosis (i.e. onset years). For both parameters, the use of raw rates yields a wider range of values because of the lack of reliability of time series recorded for sparsely populated counties. When looking at counties with at least 5 cases per year on average, the AAPC computed from joinpoint regression of raw rates and noise-filtered rates are relatively similar: the linear correlation coefficient is 0.62 for white males (54 counties) and 0.69 for black males (23 counties). Smaller correlation coefficients were observed for onset years, in particular for black males: 0.45 versus 0.68 for white males. Although the underlying temporal trends are unknown, simulation studies [21, 29] have demonstrated the benefit of noise-filtering by kriging to estimate the "true" cancer risks relatively to the use of raw estimates. We can thus hypothesize that the greater prediction accuracy of kriging translates into a more accurate modelling of temporal trends.

The impact of noise filtering on the detection of racial disparities by the new disparity statistic *D*_{
rr'
} (Equation 8) was also explored. This statistic measures how the geographical extent of significant racial disparities changed yearly over the period 1981-2007. Figure 10C showed that the number of counties with significant racial disparities peaked around 1990, because the decline started earlier for black males. Both approaches lead to similar conclusions although the 1990's peak was more apparent when analyzing raw rates. The time series in Figure 10D revealed that the percentage of counties with more favorable changes for black males (i.e. significantly smaller APC) sharply dropped in the early nineties when PSA screening was introduced. In this case, the use of noise filtered rates enhanced the 1990's peak. For both statistics, using raw rates led on average to the detection of fewer significant disparities because of the uncertainty attached to these rates.

### The MAUP effect

_{rr'}(Equation 6), greatly increased when the analysis was conducted on cases diagnosed in non-metropolitan Florida (Figure 11): the number of years with non-overlapping confidence intervals for APC jumped from 2 to 12. Dividing this geographical unit based on whether the county-level Beale index exceeds 6 or not enhanced these disparities even more, with the disparity statistic reaching a maximum of 25 years for the most rural counties (Table 4). On the other hand, this stratification of the State based on the Beale index revealed a greater similarity between races in urban areas (B

_{rr'}= 1); an interesting result was the recent increase in the frequency of late-stage diagnosis shared by both races. The county-level analysis however showed a great heterogeneity among cities. For example, racial disparities in Miami-Dade (B

_{rr'}= 2) were much smaller than in Orange county (B

_{rr'}= 16).

Impact of geographical scale on parameters of the joinpoint regression models and detection of racial disparities.

Geographical | AAPC | Onset year | Racial disparities | |||
---|---|---|---|---|---|---|

Units | WM | BM | WM | BM | H | B |

| ||||||

Florida | -3.9 | -4.4 | 1989 | 1988 | NR | 2 |

| ||||||

Urban counties | -4.0 | -4.3 | 1989 | 1988 | NR | 1 |

Rural counties | -3.6 | -5.1 | 1982 | 1982 | NR | 12 |

| ||||||

Miami-Dade | -5.2 | -6.4 | 1982 | 1986 | R | 2 |

Duval | -3.0 | -3.3 | 1989 | 1989 | R | 9 |

Hillsborough | -4.8 | -7.0 | 1990 | 1988 | R | 6 |

Orange | -3.1 | -3.0 | 1986 | 1992 | R | 16 |

Escambia | -4.2 | -3.1 | 1989 | 1982 | R | 16 |

Alachua | -3.0 | -3.3 | 1982 | 1988 | R | 6 |

Beale 3-6 | -3.8 | -3.7 | 1987 | 1987 | R | 7 |

Beale > 6 | -2.6 | -4.9 | 1991 | 1982 | R | 25 |

In their study on breast and prostate cancer survival in Michigan over the period 1985-2002, Meliker *et al*. [35] used changes in the magnitude of absolute and relative disparity statistics across geographic scales to evaluate the relative importance of innate and societal-level factors in explaining racial disparities. Unlike in the present study, racial disparities diminished and virtually disappeared in smaller geographic units (state House districts and urban neighborhoods) where the population at risk is more uniform with regard to modifiable risk factors, such as socio-economic status and proximity to screening facilities. The authors concluded that modifiable factors were responsible for apparent racial disparities observed at larger scales. The opposite trend was observed in another study [36] on the detection of disparities in breast cancer mortality among three ethnic groups in Texas (period 1995-2005). The frequency of racial disparities increased when moving down from the county level to the ZIP code and census tract levels. According to the authors, this may indicate that contextual and environmental risk factors exert different roles on health at different aggregation levels. Another culprit was the attenuation of health difference within larger geographic regions where the impact of population concentration and racial residential segregation of minorities into small and specific areas is diluted. A similar reasoning could apply here where the mixing of patients diagnosed in rural and urban areas tend to blur differences between races at the State level. However the analysis of time series distinguishes the present study from other analyses of the impact of spatial scale on racial disparities that were static and conducted for a single time period.

One solution to the MAUP effect is the use of disaggregated data [37]. Census tract level rates are the most spatially detailed data available for the present application. Thus, the creation of a continuous map of proportion of late-stage diagnosis would require a disaggregation of these data, for example using area-to-point binomial kriging [38]. For Oppenshaw and Rao [39], the answer to the MAUP effect resides in the creation of zones of approximately equal population size, or tailored to standardize the results of specific analyses. This has led to the development of automated zone matching (AZM) methodology [40] for automated zone design. For example, in a recent study on low birth weight and infant mortality in Michigan [41] AZM was used to create aggregates of ZIP codes that meet a series of constraints, such as a minimum number of cases per unit, spatial compactness and maximum intra-area correlation to ensure homogeneity in terms of race and educational level. A similar approach was used here since counties were aggregated based on their Rural-Urban Continuum Codes which proved to be an important factor in explaining the magnitude of racial disparities.

The county-level joinpoint regression analysis revealed large differences among counties within the same zone, for example urban counties (Figure 7 and Table 4), which reflects the spatial non-stationarity of the parameters of the time trend models. Geographically-weighted regression [37] was developed to account for the spatial variability of regression parameters. In this approach the regression analysis is conducted within local windows where each observation is weighted according to its proximity to the centre of the window. Such locally varying models may not be influenced by MAUP issues to the same extent as the global regression models that ignore the spatial location of the data [37]. Future research should explore the generalization of this concept to joinpoint regression models.

## Conclusions

Too often racial disparities in health outcomes are evaluated simply by computing the difference or ratio between crude rates, ignoring the lack of reliability of rates recorded for small minority populations, as well as the spatial and temporal dimensions of the data.

This paper went beyond earlier methodological work on quantifying racial disparities [35, 36, 42] by incorporating the temporal dimension through the comparison of time series instead of rates aggregated over a given time period. This comparison relied on the innovative application of joinpoint regression to rates that were stabilized using binomial kriging. A second innovation was the introduction of a new disparity statistic to supplement the tests of parallelism and coincidence available in joinpoint regression for the comparison of time series. A major advantage of this disparity statistic, which measures the number of years where APC confidence intervals did not overlap, is that it can be summed up yearly over all geographical units, providing an estimate of how the spatial extent of racial disparities changed with time. It is noteworthy that unlike previous disparity measures this statistic quantifies disparities in the pace of changes instead of the values of the health outcomes themselves. Within the context of cancer control and surveillance, this statistic facilitates the quantification of how health outcomes for different racial groups changed following strategies to improve cancer prevention and early detection, which should help better understand the causes underlying observed racial disparities in cancer incidence, mortality and morbidity.

A major challenge when working in both the spatial and temporal domains is the issue of scale or resolution. An accurate determination of when the slopes of the linear models fitted to time series undergo statistically significant changes (i.e. number and location of joinpoints) requires data that are reliable and with a good temporal resolution. Joinpoint regression has been traditionally applied to yearly time series at the State or Federal level. Keeping the same temporal resolution while zooming into the county level for example enhanced the small number problem; in particular as the focus is here on minority populations that represent on average only 10% of cases. The issue of rate instability was tackled by using both the rate estimates and standard errors provided by binomial kriging as input to joinpoint regression. A sensitivity analysis showed that kriging-based noise-filtering improved the fit by the joinpoint regression models (i.e. lower residual variability) compared to the modelling of raw rates. Another benefit of using noise-filtered rates was the clearer detection of the drop in the percentage of counties with more favorable changes for black males (i.e. significantly smaller APC) that occurred in the early nineties when PSA screening was introduced.

Another issue associated with the spatialization of joinpoint regression results is the repetition of tests of hypothesis that need to be conducted for each geographical unit, increasing the risk of false positives. Multiple testing correction was here applied using the traditional implementation of the false discovery rate (FDR) approach which is based on the underlying assumption of independence of tests. This assumption might not be appropriate for adjacent geographical units whose kriged estimates are based on common neighbors. Several techniques were proposed recently to account for highly correlated test statistics in the FDR approach [43, 44]. These approaches might however be too conservative since tests for geographical units that are further apart are independent and only adjacent counties were used in binomial kriging. In addition, some authors [28] hypothesized that the spatial dependence could be controlled by the traditional FDR correction given its statistical properties. More research is needed on this issue of multiple testing correction [45].

The case-study illustrated very well how the proportion of late-stage diagnosis for a common disease, such a prostate cancer, can change dramatically over time (i.e. 50% decline over 20 years) and display striking geographical and racial disparities within a single State. Thus, a comprehensive picture of the burden of cancer and the impact of various interventions can only be achieved through the simultaneous incorporation of the spatial and temporal dimensions in the visualization and analysis of health outcomes and putative covariates. State-level percentage of late-stage diagnosis decreased 50% since 1981; a decline that started slightly earlier for black males which had historically high rates of late-stage diagnosis. This decline accelerated in the 90's when PSA screening was introduced and stopped abruptly in the early 2000. Analysis at the metropolitan and non-metropolitan levels revealed a recent urban increase in the frequency of late-stage diagnosis for both races, and this trend was significant for white males. This result has important public health implications since it might cause an increase in prostate cancer mortality in the future. Non-metropolitan counties, despite displaying a steady decline in percentage of late-stage diagnosis, are still behind urban areas that generally have higher providers to cases ratio. Larger differences in temporal trends for Whites and Blacks were also observed in these rural counties, as measured by the new disparity statistic. In particular, the group of counties with a Beale index above 6 (i.e. most rural counties) received the maximum value for the racial disparity statistic, indicating that the yearly rate of change was significantly different between races for each of the 25 years.

The annual rate of decrease in prostate cancer late-stage diagnosis and the onset years for significant declines varied greatly among counties and racial groups. This spatial heterogeneity reflects the non-stationarity of the parameters of the time trend model even within the metropolitan and non-metropolitan strata and might warrant the development of a geographically-weighted version of joinpoint regression. Most counties with non-significant average annual percent change (AAPC) were located in the Florida Panhandle for white males, whereas they clustered in south-eastern Florida for black males. The new disparity statistic indicated that the spatial extent of racial disparities reached a peak in 1990 because of an early decline in frequency of late-stage diagnosis observed for black males. This result suggests the existence of racial disparities in the application or availability of PSA testing, in particular as the new screening procedure was introduced.

The present study was mainly methodological and the interpretation of the results suffers from limitations typically associated with ecological studies. As discussed before, the modelling of temporal trends requires some level of spatial aggregation in order to capture enough cases for a reliable estimation of percentages of late-stage diagnosis on an annual basis, which is the main culprit for the MAUP effect. In addition, the use of cases 65 year old and older allowed controlling for one source of individual-level heterogeneity because the entire study population had Medicare coverage. Individual-level data available for the same period are being analyzed to explore the impact of individual characteristics, area-level census measures of education, income, and environmental exposure on prostate cancer mortality, incidence and stage at diagnosis. These data will help conduct joinpoint regression at the sub-county level for heavily populated areas where enough cases are available for geographically detailed analysis and modelling.

## Declarations

### Acknowledgements

This research was funded by grants R43CA150496-01 and R44CA132347-02 from the National Cancer Institute, as well as grant #RSGT-10-082-01-CPHPS from the American Cancer Society. The views stated in this publication are those of the authors and do not necessarily represent the official views of the NCI and ACS. The authors thank four anonymous reviewers for their comments and suggestions to improve the presentation of the methodology and discussion of results.

## Authors’ Affiliations

## References

- Bradley CJ, Given CW, Roberts C: Race, socioeconomic status, and breast cancer treatment and survival. Journal of the National Cancer Institute. 2002, 94: 490-496.View ArticlePubMedGoogle Scholar
- Rogers SO, Ray WA, Smalley WE: A population-based study of survival among elderly persons diagnosed with colorectal cancer: Does race matter if all are insured?. Cancer Causes and Control. 2004, 15: 193-199.View ArticlePubMedGoogle Scholar
- Du W, Simon MS: Racial disparities in treatment and survival of women with stage I-III breast cancer at a large academic medical center in Metropolitan Detroit. Breast Cancer Research and Treatment. 2005, 91: 243-248. 10.1007/s10549-005-0324-9.View ArticlePubMedGoogle Scholar
- Virnig BA, Baxter NN, Habermann EB, Feldman RD, Bradley CJ: A matter of race: early- versus late-stage cancer diagnosis. Health Affairs. 2009, 28 (1): 160-168. 10.1377/hlthaff.28.1.160.PubMed CentralView ArticlePubMedGoogle Scholar
- Mohler JL: Racial Differences in Prostate Cancer Mortality. Prostate Cancer: Biology, Genetics, and New Therapeutics. Edited by: Simons JW, Chung LWK, Isaacs WB. 2006, Humana Press Inc, 355-376. SecondGoogle Scholar
- Gilligan T, Wang PS, Levin R, Kantoff PW, Avorn J: Racial differences in screening for prostate cancer in the elderly. Archives of Internal Medicine. 2004, 164: 1858-1864. 10.1001/archinte.164.17.1858.View ArticlePubMedGoogle Scholar
- Etzioni R, Tsodikov A, Mariotto A, Szabo A, Falcon S, Wegelin J, diTommaso D, Karnofski K, Gulati R, Penson DFG, Feuer EJ: Quantifying the Role of PSA Screening in the US Prostate Cancer Mortality Decline. Cancer Causes and Control. 2008, 1: 175-181.View ArticleGoogle Scholar
- Mariotto AB, Etzioni R, Krapcho M, Feuer EJ: Reconstructing PSA testing patterns between black and white men in the US from Medicare claims and the National Health Interview Survey. Cancer. 2007, 109: 1877-1886. 10.1002/cncr.22607.View ArticlePubMedGoogle Scholar
- Krieger N: Defining and investigating social disparities in cancer: critical issues. Cancer Causes and Control. 2005, 16: 5-14. 10.1007/s10552-004-1251-5.View ArticlePubMedGoogle Scholar
- Kim HJ, Fay MP, Feuer EJ, Midthune DN: Permutation tests for joinpoint regression with applications to cancer rates. Statistics in Medicine. 2000, 19: 335-351. 10.1002/(SICI)1097-0258(20000215)19:3<335::AID-SIM336>3.0.CO;2-Z. (correction: 2001;20:655)View ArticlePubMedGoogle Scholar
- Clegg LX, Hankey BF, Tiwari R, Feuer EJ, Edwards BK: Estimating average annual percent change in trend analysis. Statistics in Medicine. 2009, 28: 3670-3682. 10.1002/sim.3733.PubMed CentralView ArticlePubMedGoogle Scholar
- Kim HJ, Fay MP, Yu B, Barrett MJ, Feuer EJ: Comparability of segmented line regression models. Biometrics. 2004, 1005-1014.Google Scholar
- Jemal A, Thun MJ, Ries LA, Howe HL, Weir HK, Center MM, Ward E, Wu X-C, Eheman C, Anderson R, Ajani UA, Kohler B, Edwards BK: Annual report to the nation on the status of cancer, 1975-2005, featuring trends in lung cancer, tobacco use, and tobacco control. Journal of the National Cancer Institute. 2008, 100: 1672-1694. 10.1093/jnci/djn389.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang L, Parkin DM, Li L, Chen Y: Time trends in cancer mortality in China: 1987-1999. International Journal of Cancer. 2003, 106: 771-783. 10.1002/ijc.11300.View ArticleGoogle Scholar
- Qiu D, Katanoda K, Marugame T, Sobue T: A Joinpoint regression analysis of long-term trends in cancer mortality in Japan (1958-2004). International Journal of Cancer. 2009, 124: 443-448. 10.1002/ijc.23911.View ArticleGoogle Scholar
- La Vecchia C, Bosetti C, Lucchini F, Bertuccio P, Negri E, Boyle P, Levi F: Cancer mortality in Europe, 2000-2004, and an overview of trends since 1975. Annals of Oncology. 2010, 21 (6): 1323-1360. 10.1093/annonc/mdp530.View ArticlePubMedGoogle Scholar
- Schootman M, Lian M, Deshpande AD, Baker EA, Pruitt SL, Aft R, Jeffe DB: Temporal trends in area socioeconomic disparities in breast-cancer incidence and mortality, 1988-2005. Breast Cancer Research and Treatment. 2010, 122 (2): 533-543. 10.1007/s10549-009-0729-y.PubMed CentralView ArticlePubMedGoogle Scholar
- DeSantis C, Jemal A, Ward E, Thun MJ: Temporal trends in breast cancer mortality by state and race. Cancer Causes Control. 2008, 19: 537-545. 10.1007/s10552-008-9113-1.View ArticlePubMedGoogle Scholar
- Schootman M, Lian M, Deshpande AD, Baker EA, Pruitt SL, Aft R, Jeffe DB: Temporal trends in geographic disparities in small-area breast cancer incidence and mortality, 1988 to 2005. Cancer Epidemiology Biomarkers and Prevention. 2010, 19 (4): 1122-1131. 10.1158/1055-9965.EPI-09-0966.View ArticleGoogle Scholar
- Silva AAM, Silva LM, Barbieri MA, Bettiol H, Carvalho , Ribeiro VS, Goldani MZ: The epidemiologic paradox of low birth weight in Brazil. Rev Saúde Pública. 2010, 44: 767-775.View ArticlePubMedGoogle Scholar
- Goovaerts P: Combining area-based and individual-level data in the geostatistical mapping of late-stage cancer incidence. Spatial and Spatio-temporal Epidemiology. 2009, 1: 61-71. 10.1016/j.sste.2009.07.001.PubMed CentralView ArticlePubMedGoogle Scholar
- Coldman AG, Phillips N, Pickles TA: Trends in prostate cancer incidence and mortality: an analysis of mortality change by screening intensity. Canadian Medical Association Journal. 2003, 168 (1): 31-5.PubMed CentralPubMedGoogle Scholar
- Shaw PA, Etzioni R, Zeliadt SB, Mariotto A, Karnofski K, Penson DF, Weiss NS, Feuer EJ: An ecologic study of prostate-specific antigen screening and prostate cancer mortality in nine geographic areas of the United States. American Journal of Epidemiology. 2004, 160: 1059-1069. 10.1093/aje/kwh336.View ArticlePubMedGoogle Scholar
- Jemal E, Ward E, Wu X, Martin HJ, McLaughlin CC, Thun MJ: Geographic patterns of prostate cancer mortality and variations in access to medical care in the United States. Cancer Epidemiology, Biomarkers & Prevention. 2005, 14: 590-595. 10.1158/1055-9965.EPI-04-0522.View ArticleGoogle Scholar
- USDA: Measuring rurality: rural-urban continuum codes: Economic Research Service: US Department of Agriculture. 2004, Accessed September 14, 2011, http://www.ers.usda.gov/briefing/Rurality/RuralUrbCon/http://www.ers.usda.gov/briefing/Rurality/RuralUrbCon/Google Scholar
- Lerman PM: Fitting segmented regression models by grid search. Applied Statistics. 1980, 29: 77-84. 10.2307/2346413.View ArticleGoogle Scholar
- Kim HJ, Yu B, Feuer EJ: Selecting the number of change-points in segmented line regression. Statistica Sinica. 2009, 19 (2): 597-609.PubMed CentralPubMedGoogle Scholar
- Castro MC, Singer BH: Controlling the false discovery rate: a new application to account for multiple and dependent tests in local statistics of spatial association. Geographical Analysis. 2006, 38: 180-208. 10.1111/j.0016-7363.2006.00682.x.View ArticleGoogle Scholar
- Goovaerts P: Geostatistical analysis of disease data: estimation of cancer mortality risk from empirical frequencies using Poisson kriging. International Journal of Health Geographics. 2005, 4: 31-10.1186/1476-072X-4-31.PubMed CentralView ArticlePubMedGoogle Scholar
- Walker E, Monestiez P, Renard D, Bez N: Kriging of the latent probability of a binomial variable: application to fish statistics. Geostatistics 2008. Edited by: Ortiz J, Emery X. 2008, Santiago, Chile, GECAMIN Ltd, 981-990.Google Scholar
- BioMedware, Inc: SpaceStat User Manual version 2.2. 2011Google Scholar
- Remy N, Boucher A, Wu J: Applied Geostatistics with SGeMS: A User's Guide. 2008, New-York: Cambridge University PressGoogle Scholar
- Krieger N, Chen JT, Waterman PD, Soobader M-J, Subramanian SV, Carson R: Geocoding and monitoring of US socio-economic inequalities in mortality and cancer incidence: Does the choice of area-based measure and geographic level matter? -- The public health disparities geocoding project. American Journal of Epidemiology. 2002, 156: 471-482. 10.1093/aje/kwf068.View ArticlePubMedGoogle Scholar
- Soobader M, LeClere FB, Hadden W, Maury B: Using aggregate geographic data to proxy individual socioeconomic status: does size matter?. American Journal of Public Health. 2001, 91: 632-636.PubMed CentralView ArticlePubMedGoogle Scholar
- Meliker JR, Goovaerts P, Jacquez GM, AvRuskin GA Copeland G: Breast and prostate cancer survival in Michigan: Can geographic analyses assist in understanding racial disparities?. Cancer. 2009, 115: 2212-2221. 10.1002/cncr.24251.PubMed CentralView ArticlePubMedGoogle Scholar
- Tian N, Goovaerts P, Zhan FB, Wilson Wilson: Identification of racial disparities in breast cancer mortality: does scale matter?. International Journal of Health Geographics. 2010, 9: 35-10.1186/1476-072X-9-35.PubMed CentralView ArticlePubMedGoogle Scholar
- Fotheringham AS, Brunsdon C, Charlton ME: Geographically weighted regression: the analysis of spatially varying relationships. 2002, United Kingdom: WileyGoogle Scholar
- Goovaerts P: Geostatistical analysis of health data with different levels of spatial aggregation. Spatial and Spatio-temporal Epidemiology. 2011, 2012Google Scholar
- Openshaw S, Rao L: Algorithms for reengineering 1991 Census geography. Environment and Planning A. 1995, 27: 425-446. 10.1068/a270425.View ArticlePubMedGoogle Scholar
- Martin D, Nolan A, Trammer M: The application of zone-design methodology in the 2001 UK Census. Environment and Planning A. 2001, 33: 1949-1962. 10.1068/a3497.View ArticleGoogle Scholar
- Grady SC, Enander H: Geographic analysis of low birthweight and infant mortality in Michigan using automated zoning methodology. International Journal of Health Geographics. 2009, 8 (10):Google Scholar
- Goovaerts P, Meliker J, Jacquez GM: A comparative analysis of aspatial statistics for detecting racial disparities in cancer mortality rates. International Journal of Health Geographics. 2007, 6: 32-10.1186/1476-072X-6-32.PubMed CentralView ArticlePubMedGoogle Scholar
- Benjamini Y, Yekutieli D: The control of the false discovery rate in multiple testing under dependency. The Annals of Statistics. 2001, 29: 1165-1188. 10.1214/aos/1013699998.View ArticleGoogle Scholar
- Romano JP, Shaikh AM, Wolf M: Control of the false discovery rate under dependence using the bootstrap and subsampling. TEST. 2008, 17: 417-442. 10.1007/s11749-008-0126-6.View ArticleGoogle Scholar
- Goovaerts P: Visualizing and testing the impact of place on late-stage breast cancer incidence: A non-parametric geostatistical approach. Health and Place. 2010, 16: 321-330. 10.1016/j.healthplace.2009.10.017.PubMed CentralView ArticlePubMedGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.