For analyses of effects of illness duration, medication, episode, and symptom severity, we calculated a linear mixed-effects models (LMM) by using 3dLMMr in AFNI: y ∼1 + HAMD + Episode + Medication + Duration + Education + Age + Sex + Motion + (1| Site). This LMM describes the relationship between a response variable (e.g., VMHC) and independent variables (here, HAMD, Episode, Medication, Duration, and covariates of age, sex, education, and head motion), with coefficients that can vary with respect to the site.
By setting a threshold of p<0.001, we found significant HAMD effect and Episode effects(i.e., 2 clusters: Clust_HAMDEffect_p0.001_10_mask+tlrc and Clust_EpisodeEffect_p0.001_10_mask+tlrc). The current reviewer wants us to show HAMD effect with scatter charts and Episode effect with bar charts.
How to deal with this problem effectively? Can the Cluster Explorer feature in AFNI handle this? Maybe there’s a better way? Thanks!
PS: I have: (more specifics below):
An output dataset from 3dLMEr.
All of my individual subject datasets.
i.e., 2000 VMHC maps that went into the group level analysis (.nii).
A table(.txt) with my subject ID, HAMD, Episode, Medication, Duration, Education, Age, Sex, Motion, Site, InputFile(xxx.nii).
The master template to which the input datasets were aligned (i.e., ch2bet.nii.gz). This will serve as an underlay for viewing results.