Discrete 2 s periods of odorant (or sham) valve events were modeled in a within-subject general linear model, using as basis functions SPM's canonical hemodynamic response function (HRF). Initial testing showed that primary olfactory cortex (piriform) and olfactory association cortex (orbitofrontal) responses to odorants were maximized when the HRF onset was delayed by one second after the sniff instruction, with time and dispersion derivatives of the HRF accounting for slight variations in response onset and duration. Movement parameters from realignment were included as regressors to account for residual movement-induced effects. A high-pass filter with a cut-off of 1/128 Hz was applied to each voxel's time series to remove low frequency noise; auto-regression was not used with the long inter-stimulus interval (Della-Maggiore et al., 2002). This within-subject model yielded contrast images of activation within an odorant condition (AO, NApO, and ApCO), with each odorant set contrasted against sniffing of odorless control events (i.e., control valve opening without odorant delivery).
