Hi. I would like a voxelwise effect size measure to report alongside (i.e., in a separate map) my eff Z of my within-subject quantitative variable from my 3dLME results. To give a little background, I am determining auditory arousal thresholds about eight times throughout the night during all-night functional MRI sleep studies (https://doi.org/10.1016/j.jneumeth.2018.09.019). Therefore, I have one within-subject factor of condition with up to eight levels (Arousal 1, Arousal 2, etc.) and one within-subject quantitative variable (the auditory arousal threshold for each arousal).

Hi, how are you?! I ran 3dLME two times. In reality, I am waiting on the fetal fat_lme_prep.py program to be birthed by Paul for the first one, but let me write in past tense for both.

For the first, I used the Fisher-transformed correlations from 3dNetCorr as input. For the second, I used the unstandardized regression coefficients from afni_proc.py (via 3dDeconvolve) from a single seed region correlation/connectivity analysis. That is, I used stats.â€ś$label_fmriâ€ťâ€ś$regionâ€ťâ€ś$regsizâ€ťâ€ś$preproâ€ť.nii"[1]" from the afni_proc.py results directory where the stats file was created with a single predictor entered via -regress_stim_files. See below for the gory details.

For your second analysis, I assume that youâ€™re interested in the cross-brain correlation with the seed, not the scaling factors relative to the response at the second. So, would it better to use the Fisher-transformed correlation values from each subject as input for population-level analysis? That way, your effect magnitude at the population level would be directly related to correlation.

Hi again. â€śI assume that youâ€™re interested in the cross-brain correlation[s] with the seed [at the subject level of the analysis].â€ť Yes, I am interested in that.

â€śNot the scaling factors relative to the response at the [group level of the analysis].â€ť I am interested in that, also. My ultimate question is: does the response during each arousal correlate with the cross-brain correlations with the seed?

â€śSo, would it better to use the Fisher-transformed correlation values from each subject as input for population-level analysis? That way, your effect magnitude at the population level would be directly related to correlation.â€ť I am working on updating my analyses to use Fisher-transformed correlation values.

What are your thoughts on how I might calculate and show an effect size for the group level analysis? I would very much like to hear them.

Hi. This approach would work if I only have Fisher's z transformed r values in the model. Then the population-level coefficients would be in units of Fisher's z. However, I am interested in the correlation between Fisher's z transformed r values and another variable: auditory arousal threshold (AAT). This other variable is entered in the linear mixed model as a within-subject quantitative variable.

Therefore, the units are no longer Fisher's z transformed r values. Going back to our discussion in February, perhaps I should have said that I am most interested in "the scaling factors relative to the response at the second," if I understand what you said then correctly.

How might I use AFNI to calculate and show a voxel-wise effect size for the population-level analysis of my within-subject quantitative variable? What do you think of the articles in my initial post?

Hi again. I have an update. The approach in those references is not what I want. Specifically, the method in the Jaeger et al. is for calculating R^2 of the fixed effects in linear mixed effects.

My desire is to obtain the effect size for my one within-subject quantitative fixed effects variable. The 3dLME model I use is effectively Fisher's z transformed r ~ 1 + AAT + cond + AAT by cond + (1 + AAT | subjects). The three fixed effect variables are: AAT, cond, and AAT by cond. The method in Jaeger et al. can only measure the sum of all three fixed effects variables (AAT, cond, AAT by cond) together. The effect size measure that I want is how ordinary R^2 changes when adding my one within-subject quantitative fixed effects variable (AAT) to the model (Selya et al., 2012).

Can this be implemented in 3dLME? Alternatively, can I get R^2 of the whole equation (fixed and random effects) so that I can calculate Cohen's f^2? Please say yes.

Sincerely,

Dante

References
Selya, A. S., Rose, J. S., Dierker, L. C., Hedeker, D., & Mermelstein, R. J. (2012). A practical guide to calculating Cohen's f(2), a measure of local effect size, from PROC MIXED. Frontiers in Psychology, 3, 111. https://doi.org/10.3389/fpsyg.2012.00111

I want is how ordinary R^2 changes when adding my one within-subject quantitative fixed effects variable (AAT) to the model

Adding a variable to a model requires careful consideration. In your case, the relationship between 'cond' and 'AAT' introduces complications. Does 'AAT' differ across conditions? If so, extreme caution must be exercised.

The beta values combined with their precision in the output represent a form of effect size. However, I am not aware of any implementation for calculating R^2 or Cohen's f^2 based on the R packages nlme and lme4.

Can the marginal R-square (fixed only) and conditional R-square (fixed + random) be obtained by using r.squaredGLMM(fit1.lme) from the MuMIn R package? See

for my source of this potentially useful information. Would adding the MuMIn R package to 3dLME's capabilities help all AFNI users and not just me in the Intramural Research Program? Yes, beta values are a form of effect size, but they are not comparable across studies.

I have cautiously considered this linear mixed effects analysis and its assumptions. I have read the linear mixed effects chapter in what I consider to be the regression/correlation Bible (Cohen et al., 2003) and performed all the corresponding exercises. I have performed your hands-on exercise

for 3dLME and adapted it to my current analysis. I read the 3dLME -help output in its entirety and cautiously adapted Example 2 in my currently analysis. That is, reaction time (RT) in that example is exactly analogous to auditory arousal threshold (AAT) in my adaption of it.

Sincerely,

Dante

References
Cohen, J., Cohen, P., West, S. G., & Aiken, L. S. (2003). Random coefficient regression and multilevel models. In Applied multiple regression/correlation analysis for the behavioral sciences (3rd ed., pp. 536-567). Mahwah, NJ: Lawrence Erlbaum Associates.

Hi. How are you? How are you coping with the chaos of the impending holidays? Based on your advice from other posts, I have switched from 3dLME to 3dLMEr, so the current thread is now regarding 3dLMEr.

Do you now have time to add the MuMIn R package to 3dLMEr's capabilities? It is my hope that this or something similar will benefit all AFNI users, and I very much need the above-described effect size measure for my manuscript. I am now finished with the figures with the exception of the effect-size figures, so there is little else that I can do without them.

I've incorporated the -R2 option into 3dLMEr to offer both conditional and marginal coefficients of determination (R^2), utilizing the MuMIn R package. This enhancement will be part of the upcoming official AFNI release. In the interim, you can test it by replacing your current 3dLMEr.R file in the AFNI binary directory with this temporary file. Please let me know if it functions as intended.

The
National Institute of Mental Health (NIMH) is part of the National Institutes of
Health (NIH), a component of the U.S. Department of Health and Human
Services.