autohrf-an R package for generating data-informed event models for general linear modeling of task-based fMRI data

The analysis of task-related fMRI data at the level of individual participants is commonly based on general linear modeling (GLM), which allows us to estimate the extent to which the BOLD signal can be explained by the task response predictors specified in the event model. The predictors are constru...

Full description

Bibliographic Details
Main Authors: Nina Purg, Jure Demšar, Alan Anticevic, Grega Repovš
Format: Article
Language:English
Published: Frontiers Media S.A. 2022-12-01
Series:Frontiers in Neuroimaging
Subjects:
Online Access:https://www.frontiersin.org/articles/10.3389/fnimg.2022.983324/full
_version_ 1797902577290969088
author Nina Purg
Jure Demšar
Jure Demšar
Alan Anticevic
Alan Anticevic
Grega Repovš
author_facet Nina Purg
Jure Demšar
Jure Demšar
Alan Anticevic
Alan Anticevic
Grega Repovš
author_sort Nina Purg
collection DOAJ
description The analysis of task-related fMRI data at the level of individual participants is commonly based on general linear modeling (GLM), which allows us to estimate the extent to which the BOLD signal can be explained by the task response predictors specified in the event model. The predictors are constructed by convolving the hypothesized time course of neural activity with an assumed hemodynamic response function (HRF). However, our assumptions about the components of brain activity, including their onset and duration, may be incorrect. Their timing may also differ across brain regions or from person to person, leading to inappropriate or suboptimal models, poor fit of the model to actual data, and invalid estimates of brain activity. Here, we present an approach that uses theoretically driven models of task response to define constraints on which the final model is computationally derived using actual fMRI data. Specifically, we developed autohrf–an R package that enables the evaluation and data-driven estimation of event models for GLM analysis. The highlight of the package is the automated parameter search that uses genetic algorithms to find the onset and duration of task predictors that result in the highest fitness of GLM based on the fMRI signal under predefined constraints. We evaluated the usefulness of the autohrf package on two original datasets of task-related fMRI activity, a slow event-related spatial working memory study and a mixed state-item study using the flanker task, and on a simulated slow event-related working memory data. Our results suggest that autohrf can be used to efficiently construct and evaluate better task-related brain activity models to gain a deeper understanding of BOLD task response and improve the validity of model estimates. Our study also highlights the sensitivity of fMRI analysis with GLM to precise event model specification and the need for model evaluation, especially in complex and overlapping event designs.
first_indexed 2024-04-10T09:19:47Z
format Article
id doaj.art-64da6d9f921f4646b94f7bb85561b8a5
institution Directory Open Access Journal
issn 2813-1193
language English
last_indexed 2024-04-10T09:19:47Z
publishDate 2022-12-01
publisher Frontiers Media S.A.
record_format Article
series Frontiers in Neuroimaging
spelling doaj.art-64da6d9f921f4646b94f7bb85561b8a52023-02-20T15:55:16ZengFrontiers Media S.A.Frontiers in Neuroimaging2813-11932022-12-01110.3389/fnimg.2022.983324983324autohrf-an R package for generating data-informed event models for general linear modeling of task-based fMRI dataNina Purg0Jure Demšar1Jure Demšar2Alan Anticevic3Alan Anticevic4Grega Repovš5Department of Psychology, Faculty of Arts, University of Ljubljana, Ljubljana, SloveniaDepartment of Psychology, Faculty of Arts, University of Ljubljana, Ljubljana, SloveniaFaculty of Computer and Information Science, University of Ljubljana, Ljubljana, SloveniaDepartment of Psychiatry, Yale University School of Medicine, New Haven, CT, United StatesDepartment of Psychology, Yale University School of Medicine, New Haven, CT, United StatesDepartment of Psychology, Faculty of Arts, University of Ljubljana, Ljubljana, SloveniaThe analysis of task-related fMRI data at the level of individual participants is commonly based on general linear modeling (GLM), which allows us to estimate the extent to which the BOLD signal can be explained by the task response predictors specified in the event model. The predictors are constructed by convolving the hypothesized time course of neural activity with an assumed hemodynamic response function (HRF). However, our assumptions about the components of brain activity, including their onset and duration, may be incorrect. Their timing may also differ across brain regions or from person to person, leading to inappropriate or suboptimal models, poor fit of the model to actual data, and invalid estimates of brain activity. Here, we present an approach that uses theoretically driven models of task response to define constraints on which the final model is computationally derived using actual fMRI data. Specifically, we developed autohrf–an R package that enables the evaluation and data-driven estimation of event models for GLM analysis. The highlight of the package is the automated parameter search that uses genetic algorithms to find the onset and duration of task predictors that result in the highest fitness of GLM based on the fMRI signal under predefined constraints. We evaluated the usefulness of the autohrf package on two original datasets of task-related fMRI activity, a slow event-related spatial working memory study and a mixed state-item study using the flanker task, and on a simulated slow event-related working memory data. Our results suggest that autohrf can be used to efficiently construct and evaluate better task-related brain activity models to gain a deeper understanding of BOLD task response and improve the validity of model estimates. Our study also highlights the sensitivity of fMRI analysis with GLM to precise event model specification and the need for model evaluation, especially in complex and overlapping event designs.https://www.frontiersin.org/articles/10.3389/fnimg.2022.983324/fullfMRIGLMassumed modelingtask-related activityautohrfR
spellingShingle Nina Purg
Jure Demšar
Jure Demšar
Alan Anticevic
Alan Anticevic
Grega Repovš
autohrf-an R package for generating data-informed event models for general linear modeling of task-based fMRI data
Frontiers in Neuroimaging
fMRI
GLM
assumed modeling
task-related activity
autohrf
R
title autohrf-an R package for generating data-informed event models for general linear modeling of task-based fMRI data
title_full autohrf-an R package for generating data-informed event models for general linear modeling of task-based fMRI data
title_fullStr autohrf-an R package for generating data-informed event models for general linear modeling of task-based fMRI data
title_full_unstemmed autohrf-an R package for generating data-informed event models for general linear modeling of task-based fMRI data
title_short autohrf-an R package for generating data-informed event models for general linear modeling of task-based fMRI data
title_sort autohrf an r package for generating data informed event models for general linear modeling of task based fmri data
topic fMRI
GLM
assumed modeling
task-related activity
autohrf
R
url https://www.frontiersin.org/articles/10.3389/fnimg.2022.983324/full
work_keys_str_mv AT ninapurg autohrfanrpackageforgeneratingdatainformedeventmodelsforgenerallinearmodelingoftaskbasedfmridata
AT juredemsar autohrfanrpackageforgeneratingdatainformedeventmodelsforgenerallinearmodelingoftaskbasedfmridata
AT juredemsar autohrfanrpackageforgeneratingdatainformedeventmodelsforgenerallinearmodelingoftaskbasedfmridata
AT alananticevic autohrfanrpackageforgeneratingdatainformedeventmodelsforgenerallinearmodelingoftaskbasedfmridata
AT alananticevic autohrfanrpackageforgeneratingdatainformedeventmodelsforgenerallinearmodelingoftaskbasedfmridata
AT gregarepovs autohrfanrpackageforgeneratingdatainformedeventmodelsforgenerallinearmodelingoftaskbasedfmridata