fit_to_constraints

A helper function for fitting a model to constraints.

Analysis of task-related functional magnetic resonance imaging (fMRI) activity at the level of individual participants is commonly based on general linear modelling (GLM) that allows us to estimate to what extent the blood oxygenation level dependent (BOLD) signal can be explained by task response predictors specified in the GLM model. The predictors are constructed by convolving the hypothesised timecourse of neural activity with an assumed hemodynamic response function (HRF). To get valid and precise estimates of task response, it is important to construct a model of neural activity that best matches actual neuronal activity. The construction of models is most often driven by predefined assumptions on the components of brain activity and their duration based on the task design and specific aims of the study. However, our assumptions about the onset and duration of component processes might be wrong and can also differ across brain regions. This can result in inappropriate or suboptimal models, bad fitting of the model to the actual data and invalid estimations of brain activity. Here we present an approach in which theoretically driven models of task response are used to define constraints based on which the final model is derived computationally using the actual data. Specifically, we developed 'autohrf' — a package for the 'R' programming language that allows for data-driven estimation of HRF models. The package uses genetic algorithms to efficiently search for models that fit the underlying data well. The package uses automated parameter search to find the onset and duration of task predictors which result in the highest fitness of the resulting GLM based on the fMRI signal under predefined restrictions. We evaluate the usefulness of the 'autohrf' package on publicly available datasets of task-related fMRI activity. Our results suggest that by using 'autohrf' users can find better task related brain activity models in a quick and efficient manner.

Jure Dem<c5><a1>ar

autohrf

Automated Generation of Data-Informed GLM Models in Task-Based
fMRI Data Analysis

Jure Demšar

Nina Purg

Grega Repovš

fit_to_constraints function

<dl><dt>model_id</dt>
<dd>ID of the model.</dd>
<dt>d</dt>
<dd>A dataframe with the signal data: roi, t and y. ROI is the name of
the region, t is the timestamp and y the value of the signal.</dd>
<dt>model_constraints</dt>
<dd>A list of model specifications to use for fitting.
Each specification is represented as a data frame containing information
about it (event, start_time, end_time, min_duration and max_duration).</dd>
<dt>tr</dt>
<dd>MRI's repetition time.</dd>
<dt>roi_weights</dt>
<dd>A data frame with ROI weights: roi, weight. ROI is the
name of the region, weight a number that defines the importance of that roi,
the default weight for a ROI is 1. If set to 2 for a particular ROI that ROI
will be twice as important.</dd>
<dt>allow_overlap</dt>
<dd>Whether to allow overlap between events.</dd>
<dt>population</dt>
<dd>The size of the population in the genetic algorithm.</dd>
<dt>iter</dt>
<dd>Number of iterations in the genetic algorithm.</dd>
<dt>mutation_rate</dt>
<dd>The mutation rate in the genetic algorithm.</dd>
<dt>mutation_factor</dt>
<dd>The mutation factor in the genetic algorithm.</dd>
<dt>elitism</dt>
<dd>The degree of elitism (promote a percentage of the best
solutions) in the genetic algorithm.</dd>
<dt>hrf</dt>
<dd>Method to use for HRF generation.</dd>
<dt>t</dt>
<dd>The t parameter for Boynton or SPM HRF generation.</dd>
<dt>p_boynton</dt>
<dd>Parameters for the Boynton's HRF.</dd>
<dt>p_spm</dt>
<dd>Parameters for the SPM HRF.</dd>
<dt>f</dt>
<dd>Upsampling factor.</dd>
<dt>autohrf</dt>
<dd>Results of a previous autohrf run to continue.</dd>
<dt>verbose</dt>
<dd>Whether to print progress of the fitting process.</dd></dl>

Arguments

fit_to_constraints — fit_to_constraints

<dl>

<dt>model_id</dt>
<dd>ID of the model.</dd>


<dt>d</dt>
<dd>A dataframe with the signal data: roi, t and y. ROI is the name of
the region, t is the timestamp and y the value of the signal.</dd>


<dt>model_constraints</dt>
<dd>A list of model specifications to use for fitting.
Each specification is represented as a data frame containing information
about it (event, start_time, end_time, min_duration and max_duration).</dd>


<dt>tr</dt>
<dd>MRI's repetition time.</dd>


<dt>roi_weights</dt>
<dd>A data frame with ROI weights: roi, weight. ROI is the
name of the region, weight a number that defines the importance of that roi,
the default weight for a ROI is 1. If set to 2 for a particular ROI that ROI
will be twice as important.</dd>


<dt>allow_overlap</dt>
<dd>Whether to allow overlap between events.</dd>


<dt>population</dt>
<dd>The size of the population in the genetic algorithm.</dd>


<dt>iter</dt>
<dd>Number of iterations in the genetic algorithm.</dd>


<dt>mutation_rate</dt>
<dd>The mutation rate in the genetic algorithm.</dd>


<dt>mutation_factor</dt>
<dd>The mutation factor in the genetic algorithm.</dd>


<dt>elitism</dt>
<dd>The degree of elitism (promote a percentage of the best
solutions) in the genetic algorithm.</dd>


<dt>hrf</dt>
<dd>Method to use for HRF generation.</dd>


<dt>t</dt>
<dd>The t parameter for Boynton or SPM HRF generation.</dd>


<dt>p_boynton</dt>
<dd>Parameters for the Boynton's HRF.</dd>


<dt>p_spm</dt>
<dd>Parameters for the SPM HRF.</dd>


<dt>f</dt>
<dd>Upsampling factor.</dd>


<dt>autohrf</dt>
<dd>Results of a previous autohrf run to continue.</dd>


<dt>verbose</dt>
<dd>Whether to print progress of the fitting process.</dd>

</dl>

fit_to_constraints: fit_to_constraints

Description

Usage

Value

Arguments