Well, certainly communicating between packages means
it is hard to verify what you are doing from our end.
1. I don't know if the average squared residual (or
just the average, after sqrt()) is as appropriate as
the actual set of residual volumes, or if that is as
appropriate as the set of first level residuals (per
subject) that we tend to work with.
I would certainly want to compare those estimates with
those that came from the first level, before having
confidence in them.
2. Even from the group level, a "time" series would be
expected (one volume per subject). You are using the
mean of that, it seems.
3. If you have the ability to compute parameters from
the first level results, it would be a great idea to
compare them with those from the second level.
Perhaps Gang will chime in and confirm whether the
second level residuals are appropriate. He might
even prefer them.
- rick