Hi James,
There are multiple ways in which different scanners can affect the data -- signal and contrast to noise ratio, geometric distortion, field inhomogeneity, etc. Some of these differences will affect more or less certain steps of the processing, and often are difficult to predict. Some imaging modalities are more robust than others to scanner differences, and this also has an impact.
When it comes to the statistical analyses using the GLM, the options are essentially:
1) Add extra EVs to account for scanner(s).
2) Assume different variances for the different scanners.
Regarding (1), for two scanners, just one extra EV, coded for instance as 0 and 1, and mean-centered (make sure that the intercept is already included in the model), will take care of biases on global mean differences in signal intensity. This works for FEAT, randomise, and PALM.
Regarding (2), define one variance group for each scanner. In FEAT, this means that the design will need to be "separable" (see the documentation on how to accomplish this). In PALM, this means that permutations will happen among subjects acquired within scanner only, and a statistic that is robust to different variances will be computed. For randomise, also permutations will be within scanner, but in this case, unless the number of subjects in both scanners is the same, the results will tend to be conservative.
Hope this helps.