1. If your f-statistic is "labeled properly" (according to the help) it should
figure it out itself. If you generated it with 3dAnova, I'd guess it should
be properly labeled. If it doesn't work, let me know, and I'll look through
the code.

2. For multiple hotspots, I'd just run a different cluster analysis, with
a higher z-value threshold. I don't think there's a way to do this in
3dclust or 3dmerge automatically. You'd have to set the z-threshold
high enough to separate the hotspots into different clusters. Then there
would probably be some really cool way to recombine the outputs of
large cluster and small cluster/hotspot data. I'd have to think about that,
depending on what you wanted for a final result.

rich