Brain metastases: Insights from statistical modeling of size distribution

Abstract

BACKGROUND AND PURPOSE: Brain metastases are a common finding on brain MRI. However, the factors that dictate their size and distribution are incompletely understood. Our aim was to discover a statistical model that can account for the size distribution of parenchymal metastases in the brain as measured on contrast-enhanced MR imaging. MATERIALS AND METHODS: Tumor volumes were calculated on the basis of measured tumor diameters from contrast-enhanced T1-weighted spoiled gradient-echo images in 68 patients with untreated parenchymal metastatic disease. Tumor volumes were then placed in rank-order distributions and compared with 11 different statistical curve types. The resultant R2 values to assess goodness of fit were calculated. The top 2 distributions were then compared using the likelihood ratio test, with resultant R values demonstrating the relative likelihood of these distributions accounting for the observed data. RESULTS: Thirty-nine of 68 cases best fit a power distribution (mean R2 ¼ 0.938 6 0.050), 20 cases best fit an exponential distribution (mean R2 ¼ 0.957 6 0.050), and the remaining cases were scattered among the remaining distributions. Likelihood ratio analysis revealed that 66 of 68 cases had a positive mean R value (1.596 6 1.316), skewing toward a power law distribution. CONCLUSIONS: The size distributions of untreated brain metastases favor a power law distribution. This finding suggests that metastases do not exist in isolation, but rather as part of a complex system. Furthermore, these results suggest that there may be a relatively small number of underlying variables that substantially influence the behavior of these systems. The identification of these variables could have a profound effect on our understanding of these lesions and our ability to treat them.