Hello Susanna,

The perl script that I use to calculate the motion parameter summary statistics can be downloaded here: http://www.nitrc.org/frs/downloadlink.ph... or from the Athena wiki page.

I did use the max _absolute_ value for each of the columns, maybe that is the difference between your script and mine?  

To calculate max motion, I first compute the geometric mean of the translation at each TR and then take the maximum of these statistics.  I did the same for rotation.

This might not actually be the best way to summarize motion for quality control.  Another method measures the amount of motion between each pair of TRs and takes the maximum of this value.  The justification is that even a very large amount of motion can be corrected if it is very slow motion, i.e. no quick large transitions.

Cheers,
Cameron