Impact of inconsistent resolution on VBM studies João M.S. Pereira, Peter J. Nestor, and Guy B. Williams Department of Clinical Neurosciences, University of Cambridge School of Clinical Medicine, Cambridge, UK Received 13 September 2007; revised 16 January 2008; accepted 22 January 2008 Available online 1 February 2008 This paper considers the effects of using magnetic resonance scans with different voxel dimensions in voxel-based morphometry studies. This is of potential relevance to many longitudinal studies or any ad-hoc study that relies on pre-existing databases of subjects. In order to study this effect, a group of controls were contrasted with a group of semantic dementia as well as with a group of Alzheimer's disease patients using a mixture of different voxel dimensions scans on each side of the sta- tistical test. Scans were interpolated using a sinc function in order to obtain a different voxel depth. The effects were measured by com- paring the output of each analysis to the benchmark in which all scans had the original depth (and highest resolution), both visually and through the computation of the rootmeansquare error difference between the resulting t-maps. It was shown that the impact is highly dependent on the scan itself, with some images showing more ro- bustness to the interpolation process, and hence yielding fewer differences. A measure of robustness is proposed, which may be used in order to understand the impact of mixing different dimensions or adjusting them for each scan. Indiscriminate use of voxel dimensions on both groups was found to produce more errors (false positives/false negatives) than does an approach involving the use of balanced groups and a voxel dimension nuisance covariate. © 2008 Elsevier Inc. All rights reserved. Keywords: VBM; Voxel dimensions; Interpolation; Degradation Introduction Voxel based morphometry (VBM) is a popular tool for com- parison of structural changes, such as regional grey matter density, in patient cohorts. Ideally, raw data for VBM processing should be as coherent as possible, i.e., scans should be acquired in the same scanner with the same imaging parameters (Good et al., 2001). In general, studies only use consistently acquired data (Mechelli et al., 2005; Gitelman et al., 2001) but this may be an issue for longitudinal or multi-centre studies. For instance, it is conceivable that in rare diseases, acquisition of a cohort may take several years, spanning scanner upgrades or changes to acquisition protocols. In such circumstances, it would be potentially advantageous to include such mixed scans in a single analysis. Since the usual VBM method involves reslicing the data to a lower resolution and application of a smoothing kernel, it is feasible that such an analysis might be unbiased. The following experiments explore the impact of one source of scan heterogeneityvariation in slice thickness. This paper addressed the impact of interpolation on results. The modelling of voxel differences through a nuisance covariate was also explored. Methods Subjects Subjects were chosen from two different diseased cohorts with different patterns of atrophy: Alzheimer s disease (ADmild diffuse atrophy) and semantic dementia (SDsevere focal atrophy). Each cohort had 14 subjects and was contrasted with a sample of 14 normal controls. Demographic details are provided in the Supplementary material. These dementia cohorts present very distinct atrophic patterns: SD is characterised by severe, typically asymmetric, atrophy of the anterior temporal lobes (Williams et al., 2005), whereas AD shows a more diffuse pattern of cortical atrophy more so posteriorly as well as an atrophic hippocampal region (Nestor et al., 2006; Lerch et al., 2005). Imaging Subjects were scanned using a 1.5-T GE MRI scanner. The images were acquired using a coronal T1-weighted 3D spoiled gradient echo sequence with echo time of 4.2 ms, inversion time of 650 ms and flip angle of 20°. All scans had original voxel dimensions of 0.86 × 0.86 × 1.5 mm 3 and number of slices 256 × 256 × FoV in which FoV is the field of view on the coronal direction, in this case ranging from 110 to 124. Two additional T1 3D MPRAGE scans were acquired from a single healthy volunteer using a Siemens Trio 3-T MRI scanner, with voxel dimensions 1×1×1 mm 3 (matrix size 240 × 256 × 176) and 1 × 1 × 1.5 mm 3 (matrix size 240 × 256 × 256) with echo time 5.9 ms, inversion time 900 ms and flip angle 9°. www.elsevier.com/locate/ynimg NeuroImage 40 (2008) 1711 1717 Corresponding author. E-mail address: jmsp2@wbic.cam.ac.uk (J.M.S. Pereira). Available online on ScienceDirect (www.sciencedirect.com). 1053-8119/$ - see front matter © 2008 Elsevier Inc. All rights reserved. doi:10.1016/j.neuroimage.2008.01.031