Hi AFNI team,
I was wondering if it would be possible to obtain model fit statistics (ideally log likelihood) at the voxel-level from 3dLME. The logic here being, as I add predictors to a multilevel model (e.g. gender, group status, etc.) I would like to identify regions (voxels) where not only does my predictor demonstrate a significant fixed effect, but its addition also creates a more parsimonious model (compared to a model without that predictor). To do this I would like to be able to compare the model fit statistics of multiple models (via chi-square or something of the sort) in a voxel-wise fashion. To my knowledge 3dLME doesn't contain an option to write out model fit statistics (please correct me if I'm wrong!), but I can image while it's running lme iteratively in R it would be possible to grab the logLik value along with the coefficients?
Thanks!
Scott
