BrainVoyager v23.0

Multi-Subject RFX GLM

Introduction

For second-level random effects (RFX) analyses, statistical maps containing estimated effects (beta values) separately for each subject are required as input. These beta values can be computed with the RFX-GLM function, which is an optimized version of a standard "separate-subject predictors" GLM. While a standard multi-subject GLM can also be used to prepare RFX analyses, the RFX-GLM routine operates substantially faster and requires much less working memory by running internally a series of individual single-subject GLMs instead of one monolithic GLM (for details, see below). In order to use the RFX-GLM, the RFX GLM option must be checked in the General Linear Model: Multi-Study, Multi-Subject dialog. The multi-subject design matrix (MDM) is specified as usual by adding funtional runs (VTCs or MTCs) and associated single-run (RTC) design matrix files.

Detailed Description

Random effects analysis can be performed in BrainVoyager by running a multi-subject GLM with predictors separated for each included subject. For the resulting GLM the random effects option can be turned on to test effects across subjects using the summary statistics approach. To perform the RFX analysis, the same contrast is specified for each subject and the resulting mean value across the subjects is tested against zero using a t-test. In a similar way, two groups can be compared by computing the mean of the summary statistic (contrast) for each group followed by a t-test to compare the group means. In case of multi-factorial designs with within and between factors and potentially covariates, the ANCOVA dialog can be used.

Why RFX-GLM?

A limiting factor for all random effects analyses is that the standard separate-subject predictor GLM does not scale well to cases with a large number (e.g. more than 50) subjects. This property of the standard multi-subject GLM rests in the used approach to construct an overall design matrix for all included studies of all subjects. The resulting huge matrix is then further processed to estimate simultaneously all beta values for all subjects. While this approach constitutes the most powerful approach (allowing to formulate e.g. interactions between subjects and conditions), the possibilities of a monolithic design matrix are not needed for hierarchical random effects analyses. For RFX analyses the monolithic design matrix approach is inefficient since the created overall design matrix contains mostly zero values: Only diagonal sections where time courses and subject predictors are crossed are filled in the design matrix. In order to speed up the computation of random effects analyses, the new RFX-GLM has been introduced in version 1.3 of BrainVoyager QX. The major difference to the standard multi-subject GLM is that an overall design matrix is not created. Instead, the beta values are estimated individually for each included study and the results from the studies belonging to the same subject are pooled. This approach results in the same beta values as the standard multi-subject GLM approach but the computations are obtained much faster with a substantially reduced memory load. This approach internally runs a series of single-study GLMs resulting in linear scaling behavior with respect to both memory requirements and computation time. While running a RFX-GLM, the Log pane presents information about each conducted single-study GLM. The resulting GLM file can not be used to test fixed effects contrasts since the standard errors of the estimated betas are not computed (in the future the first-level standard errors may be integraated to calculate mixed effects analyses). It can be used, however, for second-level (summary statistics) random-effects analysies using the Overlay GLM Contrasts and ANCOVA dialog. Because of its advantages, it is recommended to use the RFX-GLM for random effects analyses. Since the RFX-GLM estimates beta values as input for a second level random effects analysis, it is also not strictly necessary to turn on the "correct for serial correlations" option since the estimate of betas is unbiased even in the case of serial correlations and the fixed-effects first-level standard error values are not used (at present) in the second level analysis; turning off the correct for serial correlations option further increases the speed of GLM calculations The RFX-GLM uses as default the percent signal change time course normalization option, which is recommended for random effects analyses because it appears to reflect differences in effect size between subjects better than z normalization. When running a RFX-GLM, estimated beta values of all subjects are shown in the Voxel Beta Plot window while moving the mouse cursor over individual voxels.


Copyright © 2023 Rainer Goebel. All rights reserved.