Hi Gang,
Thanks for the helpful reply! To be frank, I hadn't thought about 3dREMLfit. My near-term goal was to get a volume of beta-weight estimates and corresponding t-scores for each stimulus event in an event-related design (for later use with RSA). Before 3dLSS existed the only way I knew to do that was to run 3dDeconvolve once for each stimulus, modeling all other stimuli separately. Just using 3dDeconvovle to generate the design matrix and then feeding that to 3dLSS to make the maps for each stimulus is orders of magnitude faster, so that's what I was referring to.
But now that you mention it, if I read the program help correct, it looks like 3dREMLfit will also generate a map per stimulus without having to run it separately for each stimulus, AND it will output nice t-maps as well. If that's right, then it definitely seems worth a shot!
-Will