Reproducibility of infant fNIRS studies: a meta-analytic approach

Jessica Gemignani; Irene de la Cruz-Pavía; Anna Martinez; Caroline Nallet; Alessia Pasquini; Gaia Lucarini; Francesca Cavicchiolo; Judit Gervain

doi:10.1117/1.NPh.10.2.023518

10 March 2023 Reproducibility of infant fNIRS studies: a meta-analytic approach

Jessica Gemignani, Irene de la Cruz-Pavía, Anna Martinez, Caroline Nallet, Alessia Pasquini, Gaia Lucarini, Francesca Cavicchiolo, Judit Gervain

Author Affiliations +

Neurophotonics, Vol. 10, Issue 2, 023518 (March 2023). https://doi.org/10.1117/1.NPh.10.2.023518

Abstract

Significance

Concerns about the reproducibility of experimental findings have recently emerged in many disciplines, from psychology to medicine and neuroscience. As NIRS is a relatively recent brain imaging technique, the question of reproducibility has not yet been systematically addressed.

Aim

The current study seeks to test the replicability of effects observed in NIRS experiments assessing young infants’ rule-learning ability.

Approach

We conducted meta-analyses and mixed-effects modeling-based inferential statistics to determine whether effect sizes were replicable and comparable in a sample of 23 NIRS studies investigating infants’ abilities to process repetition- and diversity-based regularities in linguistic and nonlinguistic auditory and visual sequences. Additionally, we tested whether effect sizes were modulated by different factors such as the age of participants or the laboratory. We obtained NIRS data from 12 published and 11 unpublished studies. The 23 studies involved a total of 487 infants, aged between 0 and 9 months, tested in four different countries (Canada, France, Italy, and USA).

Results

Our most important finding is that study and laboratory were never significant moderators of variation in effect sizes, indicating that results replicated reliably across the different studies and labs included in the sample. We observed small-to-moderate effect sizes, similar to effect sizes found with other neuroimaging and behavioral techniques in the developmental literature. In line with existing findings, effect sizes were modulated by the participants’ age and differed across the different regularities tested, with repetition-based regularities giving rise to the strongest effects; in particular, the overall magnitude of this effect in the left temporal region was 0.27 when analyzing the entire dataset.

Conclusions

Meta-analysis is a useful tool for assessing replicability and cross-study variability. Here, we have shown that infant NIRS studies in the language domain replicate robustly across various NIRS machines, testing sites, and developmental populations.

1. Introduction

Since the introduction of near-infrared spectroscopy, its use has grown dramatically in the neurosciences, particularly in developmental cognitive neuroscience,¹^,² due to its many practical advantages: ease of use, portability, noninvasiveness, and low cost.

As the volume of fNIRS research in developmental neuroscience has increased, so have the concerns about the reproducibility of experimental findings, a core principle of scientific progress. Issues about reproducibility have emerged in many related areas, from psychology to medicine, as some key findings were found not to replicate reliably.³^–⁶ Reasons for the replication crisis include undocumented diversification in research methodology and data analysis techniques, nontransparent data inclusion and exclusion criteria, difficulty publishing null results, HARKing (hypothesizing after the results are known), p-hacking, and other dubious analysis and reporting practices. Some of these issues are further exacerbated in developmental research, as young children are challenging participants, which may lead to large amounts of missing data, unclear data rejection criteria, short study durations and thus small numbers of trials per participant, small sample sizes, and relatedly, low statistical power. Developmental NIRS research is no exception, especially because NIRS is a relatively new technique, and standardization of research practices has only started recently,⁷ e.g., with systematic comparisons of handware performance,⁸^,⁹ preprocessing methods,¹⁰^–¹² and statistical analyses.¹³^,¹⁴

NIRS by now has become a very commonly used neuroimaging technique in developmental neuroscience. It is thus timely to address the question of reproducibility in a systematic way to support good theory-building and to identify issues that weaken replicability.

One approach to assessing replicability, quantifying cross-lab variability and identifying potential moderating factors underlying variability is to conduct a meta-analysis of existing (published and/or unpublished) studies that address the same research question. A meta-analysis is a quantitative method for aggregating across experimental studies¹⁵ to reveal the average effect size of a specific phenomenon. In addition to quantifying replicability and variability across studies, meta-analyses have the advantage of pooling data over a larger sample size than typically possible in single studies, licensing more robust or more general conclusions.¹⁶ Meta-analyses have been conducted over various types of data in many different domains of psychology, neuroscience, and medicine¹⁷^–²² and are readily applicable to developmental data.²³

The goal of the current study is, therefore, to conduct a meta-analysis of infant NIRS studies in an attempt to test for reproducibility, quantify variability across studies, and identify moderators explaining the observed variability. We chose repetition-based rule learning as the phenomenon under investigation. The ability to extract and learn rules from speech is foundational for language development.²⁴ This question has thus received considerable attention in the language acquisition literature. In particular, a large number of behavioral and NIRS studies (for reviews, Refs. 25 and 26, respectively) asked whether infants are able to recognize and represent different sequences that include a repetition (e.g., AAB versus ABB) or distinguish them from diversity-based sequences, i.e., sequences in which all syllable are different (e.g., ABB versus ABC). These studies typically use artificial grammars that generate sequences with the relevant structure (e.g., ABB: “mubaba,” “penana,” etc., ABC: “mubage,” “penaku,” etc.).

We chose this research question because a relatively large number of published and unpublished NIRS studies addressed it using similar stimuli and experimental designs. Such a sample of studies allows for a good estimation of lab-based variation, e.g., differences in machine type, experimenter characteristics, NIRS-relevant population characteristics, such as hair density and color, etc. Many of these factors are not under the researchers’ direct control and often remain unreported in publications. This is the type of variation that meta-analyses are particularly well suited to assess.

Ours is the first meta-analytic study of infant NIRS data, and its aim is primarily methodological. We seek to provide a first quantitative assessment of variability in effect sizes across NIRS studies that differ along specific dimensions to address the issue of cross-lab/cross-study variability. In other words, we are less concerned here with the theoretical issues related to infants’ ability to learn linguistic rules, despite its relevance for developmental research. For these questions, we refer the reader to Refs. 25 and 26.

2. Methods

2.1.

Data

2.1.1.

Studies

We aggregated 23 published and unpublished fNIRS studies, conducted in four different laboratories, testing brain responses of typically developing infants to two different types of linguistic regularities: repetition-based regularities (R) and diversity-, i.e., nonrepetition-based regularities (N). The studies were identified by searching through PubMed and Google Scholar, using the search strings “repetition-based regularity,” “rule learning,” “fNIRS,” and “infants.” Exclusion criteria included (i) testing atypical populations or (ii) using methods other than NIRS. Papers including more than one study were considered separate studies. Of the 43 hits, those that met either of the exclusion criteria were discarded, leaving 12 published studies.²⁷^–³¹ Additionally, 11 studies from the last author’s laboratory were added. These studies were not published in peer-reviewed articles, although some of them are available online in PhD dissertations (Table 1). We know of no other unpublished studies.

Table 1

List of studies included in each of the three meta-analyses (R versus 0, N versus 0, and R versus N). As described in Sec. 2.1.1, different sets of studies contributed to the three comparisons of interest, indicated by the color of the last three columns (bold = included, NA = not included); the columns also report which regularity was considered for each study.

The studies used similar methodology, e.g., similar stimuli and experimental designs, but it addressed slightly different theoretical questions and, as a result, varied along a few dimensions (e.g., the auditory versus visual nature of the stimuli, see Table 1). We used these factors as moderators in the meta-analysis. Furthermore, studies varied in whether they tested repetition-based regularities, diversity-based regularities, or both (Table 1). Consequently, we conducted three separate meta-analyses evaluating the effect sizes of (i) the comparison between brain responses to repetitions versus a zero baseline (“R versus 0”); (ii) the comparison between brain responses to diversity (nonrepetition) versus a zero baseline (“N versus 0”); and (iii) the comparison between brain responses to repetitions versus diversity (“R versus N”). A specific study may have contributed to just one or several of the comparisons (see the last column of Table 1). As a result, 23 studies were included in the final analysis for the R versus 0 comparison, 19 in the analysis of the R versus N comparison, and 17 in the analysis of the N versus 0 comparison. Details can be found in Table 1. Additionally, Table 2 reports the most relevant technical details of NIRS data acquisition.

Table 2

Technical characteristics of the NIRS machines employed in each study.

ID	Study	Lab	CW-NIRS	Fs (Hz)	Wavelengths (nm)
1	0m-Speech-A_A-A_C	Trieste	Hitachi ETG-4000	10	695, 830
2	0m-Speech-AA-AB	Trieste	Hitachi ETG-4000	10	695, 830
3	6m-VisualNonLinguistic-AA-AB	Paris	NIRx NIRScout 8x16	15.625	760, 850
4	6m-VisualLinguistic-AA-AB	Paris	NIRx NIRScout 8x16	15.625	760, 850
5	6m-Speech-AA-AB	Paris	NIRx NIRScout 8x16	15.625	760, 850
6	0m-Speech-AAB-ABC	Vancouver	Hitachi ETG-4000	10	695, 830
7	0m-Speech-AAB-ABB-AltNonAlt-R=AAB	Vancouver	Hitachi ETG-4000	10	695, 830
8	0m-Speech-AAB-ABB-AltNonAlt-R=ABB	Vancouver	Hitachi ETG-4000	10	695, 830
9	0m-Speech-AAB-ABB-SimpleBlocks-R=AAB	Vancouver	Hitachi ETG-4000	10	695, 830
10	0m-Speech-AAB-ABB-SimpleBlocks-R=ABB	Vancouver	Hitachi ETG-4000	10	695, 830
11	0m-Speech-AAB-ABC-ReplicationNIRx	Paris	NIRx NIRScout 8x16	10	760, 850
12	0m-Tones-AAB-ABC	Paris	Hitachi ETG-4000	10	695, 830
13	6m-Speech-AAB-ABC-Entropy	Paris	NIRx NIRScout 8x16	15.625	760, 850
14	0m-Speech-ABA-ABC	Trieste	Hitachi ETG 4000	10	695, 830
15	0m-Speech-ABB-ABC	Trieste	Hitachi ETG 4000	10	695, 830
16	0m-Speech-ABB-ABC-RepeatedContext-R=C	Paris	NIRx NIRScout 10x8	10.417	760, 850
17	0m-Speech-ABB-ABC-RepeatedContext-R=V	Paris	NIRx NIRScout 10x8	10.417	760, 850
18	0m-Speech-ABB-ABC-VariableContext-R=C	Paris	NIRx NIRScout 10x8	10.417	760, 850
19	0m-Speech-ABB-ABC-VariableContext-R=V	Paris	NIRx NIRScout 10x8	10.417	760, 850
20	0m-Speech-ABB-ABC-CV	Paris	NIRx NIRScout 10x8	10.417	760, 850
21	6m-Speech-ABB-ABC-AltNonAlt	Paris	NIRx NIRScout 8x16	15.625	760, 850
22	7m-Speech-ABB-ABC	Boston	Hitachi ETG-4000	10	695, 830
23	9m-Speech-ABB-ABC	Boston	Hitachi ETG 4000	10	695, 830

The study comprised data from a total of 487 infants, aged between 0 and 9 months, tested in four different countries (Canada, France, Italy, and USA). Table 1 provides information about each study’s individual sample size.

2.1.2.

Materials

All included studies used two artificial grammars in the form of bisyllablic or trisyllabic sequences. Table 1 reports the specific structures employed in each study. Other stimulus characteristics can be found in the respective publications and were highly similar across studies.

2.1.3.

Procedure

In all studies, infants were tested with an NIRS device (the brand, wavelengths, and sampling frequencies are listed in Table 2), and sound stimuli were administered through two loud speakers.

Eight to ten sources and eight detectors were placed on the infants’ heads bilaterally, with a 2.5 to 3 cm source–detector distance, forming 10 to 12 channels per hemisphere (Fig. 1). The anatomical localization of the resulting array is described in detail in Ref. 32. On the basis of the original study²⁸ that first tested newborns’ abilities to detect repetition-based sequences, and given the well-documented relevance of this area for speech and language processing, we focused on the left temporal lobe, which we operationalized as the cluster of channels 3 and 6 (Fig. 1). Choosing the most relevant region of interest (ROI) is a research question in and of itself, with both anatomical and functional localization methods now available in the literature.⁷ Although this issue is relevant for studies comparing different datasets and headgears, it falls outside the scope of the current study. Here, we have chosen to use a predefined ROI, the left temporal region. This area is activated in the large majority of the included studies. Further research will need to address other, study-specific methods for determining the ROI for a meta-analysis.

Fig. 1

Optode arrangement employed in the studies included 8 or 10 sources and 8 detectors, forming a total of 20 or 24 channels. The region of interest that we focused on is the left temporal lobe, formed by channels 3 and 6; it is highlighted in red.

2.2.

Data Analysis

2.2.1.

fNIRS preprocessing

For studies 1 to 21 (Table 1), NIRS data was preprocessed in the same way as in the original publications (or sources, e.g., unpublished PhD dissertations), which was similar across most studies. Briefly, light intensities were converted to optical densities and to hemoglobin concentration changes using the modified Beer–Lambert Law with the absorption coefficients $μ_{a}$ , ${mm}^{- 1} \times {mM}^{- 1}$ : $μ_{a} (HbO, 695 nm) = 0.0955$ , $μ_{a} (HbO, 760 nm) = 0.1496$ , $μ_{a} (HbO, 830 nm) = 0.2320$ , $μ_{a} (HbO, 850 nm) = 0.2526$ ; $μ_{a} (HbR, 695 nm) = 0.4513$ , $μ_{a} (HbR, 760 nm) = 0.3865$ , $μ_{a} (HbR, 830 nm) = 0.1792$ ; and $μ_{a} (HbR, 850 nm) = 0.1798$ . The product of the optical pathlength and the differential pathlength factor was set to 1, resulting in concentration changes being expressed in $mM \times mm$ . A bandpass filter between 0.01 and 0.7 Hz was applied to concentration changes using an $f f t$ digital filter. Then, as illustrated in Ref. 10, blocks of single-trial data were rejected if they contained motion artifacts or if the light intensity reached the saturation value, with motion artifacts defined as signal changes larger than $0.1 mM \times mm$ over 0.2 s. The artifact detection and trial rejection procedures were performed independently for each channel, and channels with less than at least two valid blocks were discarded from the analysis. The trial inclusion rate for each study ranged between 52% and 100% (M: 65.1%, SD: 12.8%). Finally, for the nonrejected blocks, a baseline was linearly fit between the mean of the 5 s preceding the onset of the block and the mean of the 5 s preceding the onset of the next one. Blocks were then averaged within each infant to obtain channel-wise block averages for each condition as well as across infants to obtain grand averages. This preprocessing routine has been shown to yield an accurate recovery of the infant hemodynamic response.¹⁰ Study-level grand averages were employed to compute study-level effect sizes, whereas individual trial averages were employed to compute infant-level effect sizes (Sec. 2.2.2).

For studies 22 and 23 (Table 1), we obtained each subject’s channel-wise average activation to each experimental condition, i.e., we obtained preprocessed data from the authors and had no access to the raw data. For these studies, we could, therefore, only compute the study-level effect size but not the individual infant-level effect size. Data processing for these datasets is described in the original publication.²⁷

2.2.2.

Calculation of effect sizes

We performed two analyses using two different and complementary approaches: a meta-analytical approach that analyzes study-level effect sizes and a mixed-effects modeling approach that analyzes infant-level effect sizes. The meta-analytic framework estimates the variability of effect sizes across studies. It can be conducted even when only group-level averages, but not individual participant data, are available and when procedures or data types are not standardized.²³ When trial-level data for each participant is available, it is possible to compute individual effect sizes and perform a mixed-effects model, yielding a more sensitive measure of within-study variability. In the current study, participant-level NIRS data was available for 21 of the 23 included studies, and it accounted for the entirety of two age groups (newborns and 6-month-olds).

We conducted both meta-analyses and mixed-effects models for the 21 studies for which individual trial-level data was available, whereas the two studies for which we only had individual averages were only entered into the meta-analysis. We did not use NIRS data as the dependent variable in either the meta-analyses or the mixed-effects models because relative concentrations, which are the NIRS measures obtained in the included studies, are not necessarily comparable across different participants and different machines. We used effect sizes as the dependent variable instead.

For both analyses, we first averaged the time series of the hemodynamic response within a block over a time window starting at the onset of the stimulus and lasting up to 15 s after the end of the stimulation block. Effect sizes were then computed for three comparisons of interest: (i) R versus 0, i.e., the comparison between the repetition condition and the zero baseline; (ii) N versus 0, i.e., the comparison between the nonrepetition condition and the zero baseline; (iii) R versus N, i.e., a comparison between repetition and nonrepetition responses.

The meta-analytic effect sizes ( $d_{study}$ ) were computed by averaging a participant’s responses in all trials of a given condition (Fig. 2). These individual means were then averaged and divided by their standard deviation. The sampling variances of the meta-analytic effect sizes were computed as $V_{d} = 2 / n + d_{study}^{2} / 4 n$ , with $n$ being the number of participants¹⁵; i.e., effect sizes were weighted by the number of participants in a study.³³

Fig. 2

Schematic illustration of infant-level and study-level effect sizes calculations. Activation refers to the R versus 0, N versus 0, and R versus N contrasts computed as the average of the HRF along its time course. Magenta and cyan indicate repetition trials (HbO and HbR, respectively), and red and blue indicate nonrepetition trials (HbO and HbR, respectively).

The individual effect sizes ( $d_{infant}$ ) were computed by dividing each infant’s average activation in a given condition by its standard deviation across trials (Fig. 2).

Both individual and meta-analytic effect sizes were calculated for each channel and hemoglobin component independently. Then, we averaged effect sizes across channels within our predefined ROI, i.e., the left temporal area.

2.2.3.

Statistical analysis

Meta-analysis

Study-level effect sizes were analyzed by employing meta-analytic random-effects models with the metafor package in R.³⁴ Models were fit using restricted maximum likelihood (REML).

In addition to this overall meta-analysis, we also conducted moderated meta-analyses to test the effects of specific factors on cross-study effect sizes and variability. Specifically, we explored whether the effect was significantly moderated by the lab (Boston/Trieste/Paris/Vancouver), as this allowed us to directly address our main question of cross-lab replicability. We also tested the moderator age with four levels (newborns/6-/7-/9-months). Infants’ ability to learn repetition-based regularities has been reported to improve between 5 and 11 months³⁵^,³⁶ at the behavioral level, and by 6 months of age, infants show much stronger brain responses to diversity-based regularities than at birth.²⁶ The ability to process these structural regularities may thus undergo developmental change, which we sought to model by adding age as a moderator to the meta-analyses.

Finally, we analyzed whether the effect size was moderated by repetition position within trisyllabic sequences (initial/final). Existing results suggest that final repetitions may be easier to process than initial ones, possibly due to recency effects in memory, although this advantage may be relatively weak.³⁷ We thus included repetition position for the comparisons in which it was relevant, namely R versus 0 and R versus N.

Importantly, these three analyses explored each moderating variable independently, and their joint impact on the effect was analyzed through mixed-effects modeling (see below).

Although our dataset included both published and unpublished studies, for the sake of completeness, we also generated funnel plots and the fail-safe-N as estimates of publication bias.³⁴ Funnel plots can identify certain forms of publication bias. Estimates are expected to be randomly sampled around the mean if there is no bias in the published or selected studies, whereas an asymmetry in the distribution of the effects may indicate either true heterogeneity in the phenomenon under study or that studies with nonsignificant findings remained unpublished.¹⁶ The fail-safe-N is another useful measure to detect bias in the literature as it quantifies how many unpublished studies with null results would have to exist for the overall effect size to be zero. We computed the funnel plots and the fail-safe-N using the metafor R package,³⁴ following recommendations in the literature.³⁸

Linear mixed-effects modeling

Linear mixed-effects models were computed over infant-level effect sizes. The planned random effects structure consisted of random intercepts for study, lab, and age. In the case of failure to converge, it was iteratively simplified by pruning first the random intercept for age, then for lab, and finally for study. Candidate fixed effects were incrementally included in the fixed effects structure, and the resulting models were compared. The best-fitting model was chosen based on the AIC (Akaike Information Criterion) value.

Application of analyses to study sets

Analyses were performed over different sets of studies as a function of the factors to which a given study contributed (Table 1). Specifically, three sets were created on the basis of sample sizes (Fig. 3): (i) the entire dataset (Set 1; R versus 0: 23 studies, N versus 0: 17 studies, R versus N: 19 studies); (ii) a set of studies using speech stimuli (Set 2; R versus 0: 20 studies, N versus 0: 14 studies, R versus N: 16 studies); and (iii) a set of studies using speech stimuli with adjacent repetitions (Set 3; R versus 0: 16 studies, R versus N: 13 studies). The meta-analytic approach was applied to all studies, and mixed-effects modeling was employed only for studies for which raw data was available (studies 1 to 21).

Fig. 3

Organization of the dataset into the three sets of studies. The numbers overlaid on the bars indicate the number of participants in each category. Age is color coded (green is newborns, and violet is 6- to 9-month-olds).

3. Results

3.1.

Meta-Analysis

3.1.1.

Unmoderated models

R versus 0

When analyzing the entire dataset (set 1), the overall magnitude of the effect was 0.271 (95% $CI = [0.144, 0.398]$ , $z = 4.20$ , $p < 0.001$ ). Breaking it down by age, the effect size was 0.273 for newborns, 0.289 for 6-month-olds, and 0.266 for 7- to 9-month-olds. When considering studies employing speech stimuli (Set 2), the meta-analytic effect size was 0.282 ([0.146, 0.418], $z = 4.07$ , $p < 0.001$ , 0.282 for newborns, 0.327 for 6-month-olds, and 0.28 for 7- to 9-month-olds). For studies with adjacent trisyllabic repetitions in speech (set 3), the overall magnitude of the effect was 0.288 ([0.136, 0.439], $z = 3.73$ , $p < 0.001$ ), with an effect size of 0.289 for newborns, 0.348 for 6-month-olds, and 0.283 for 7- to 9-month-olds. Corresponding forest plots are shown in the left panel of Fig. 4.

Fig. 4

Forest plots of the meta-analytic effect sizes for brain activation in the left temporal area to repetition-based (R versus 0) and diversity-based (N versus 0) regularities, as well as for the difference in activation between them (R versus N) for HbO. Each study’s estimate is indicated by the corresponding square. Error bars indicate the 95% confidence interval. Diamonds show the summary estimates of each set of studies, with the center of the diamond corresponding to the estimate and the left and right edges indicating the confidence interval limits. Not all studies contributed to all comparisons, as described in Sec. 2.1.1 and Table 1. Corresponding forest plots obtained on HbR time traces are reported in the Supplementary Material (Fig. S1).

N versus 0

The overall magnitude of the effect, computed on the entire dataset, was 0.18 (95% CI = [0.03, 0.33], $z = 2.35$ , $p < 0.05$ ). Estimates were 0.14 for newborns, 0.20 for 6-month-olds, and 0.25 for 7- to 9-month-olds. For Set 2, the estimated meta-analytic effect size was 0.17 (95% CI = [0.01, 0.34], $z = 2.074$ , $p < 0.05$ ), with 0.14 for newborns, 0.18 for 6-month-olds, and 0.26 for 7- to 9-month-olds (Fig. 4).

R versus N

The overall magnitude of the effect was 0.08 (95% $CI = [- 0.06, 0.22]$ , $z = 1.14$ , ns), with a larger effect for newborns than for older infants (newborns: 0.12, 6-month-olds: 0.04, 7- to 9-month-olds: 0). For set 2, we found an overall meta-analytic effect size of 0.10 (95% $CI = [- 0.05, 0.25]$ , $z = 1.28$ , ns), with estimates of 0.133 for newborns, 0.072 for 6-month-olds, and 0.01 for 7- to 9-month-olds (Fig. 4). For set 3, the overall effect size was 0.14 (95% $CI = [- 0.04, 0.32]$ , $z = 1.56$ , ns), with subgroup estimates that were again larger in newborns (0.19) than for older infants (6-month-olds: 0.12, 7- to 9-month-olds: 0.02).

3.1.2.

Moderated analyses

R versus 0

The lab-moderated models yielded a nonsignificant effect of lab in the three sets of data, as did the age-moderated models for age. Crucially, lab and age were highly collinear because some ages were only tested in some labs. Therefore, these two moderating effects are better interpreted when analyzed jointly through mixed-effects modeling (Sec. 3.2). In set 3, the only set in which this moderator was analyzed, the model yielded no significant effect of repetition position.

N versus 0

The lab-moderated models yielded a nonsignificant effect of lab in all sets of data, as did the age-moderated models for age.

R versus N

Lab-moderated models yielded nonsignificant effects of lab in all sets of data, as did the age-moderated models for age in sets 1 and 2. For set 3, age showed a trend toward significance, with the intercept being 0.19 for newborns ( $[- 0.03, 0.40]$ , $z = 1.72$ , $p < 0.1$ ), and the decreases in effect magnitude were 0.07 for 6-month-olds, 0.36 for 7-month-olds, and 0.29 for 9-month-olds (Table S6 in the Supplementary Material). The model including repetition position as a moderator for set 3 yielded no significant effect.

3.1.3.

Estimates of bias

Funnel plots

Figure 5 shows funnel plots for the three comparisons of interest over the entire dataset. In addition to visual inspection, we also carried out a rank correlation test³⁹ to detect asymmetries in the funnel plots. Asymmetries were not detected in any comparison (HbO, R versus 0: Kendall’s $τ = 0.05$ , ns; $N$ versus 0: Kendall’s $τ = 0.22$ , ns; R versus N: Kendall’s $τ = 0.04$ , ns; HbR, R versus 0: Kendall’s $τ = 0.11$ , ns; N versus 0: Kendall’s τ = 0.03, ns; R versus N: Kendall’s $τ = 0.11$ , ns).

Fig. 5

Meta-analytic effect sizes by age in sets 1 to 3. Gray bands indicate 95% confidence intervals, computed from the standard errors of each random-effects model. The top panel reports effect sizes obtained from HbO time traces, and the bottom panel shows the corresponding plots for HbR.

Fail-safe-N

The estimated fail-safe-N for the R versus 0 comparison over set 1 was 126, indicating that 126 studies with a null effect would have to exist in the “file drawer” for the overall effect size to become zero. The fact that this number is much greater than the number of studies in the analysis suggests that the effect of repetition-based regularities on brain responses is robust.

3.2.

Mixed-Effects Modeling

Distributions of individual effect sizes are shown in Fig. 6. Grand averages are shown in Fig. 8.

Fig. 6

Funnel plots of effect sizes against standard errors, for each of the three comparisons of interest, obtained on HbO (top panel) and HbR time traces (bottom panel); the white region shows the 95% confidence interval around the estimates, and the gray region shows the 95% to 99% interval; colors indicate the lab where the study was carried out, and shapes indicate the age group of each study.

3.2.1.

Oxygenated hemoglobin (HbO)

For the R versus 0 comparison over set 1, the best fitting model included a random intercept for StudyID and fixed effects for lab and age, and it yielded a significant main effect of age (F(1, 410) = 3.96, $p < 0.05$ ), carried by a marginally larger effect for 6-month-olds than for newborns (estimate: $- 0.166$ , $t (19) = - 1.97$ , $p = 0.06$ ). Over set 2, the best-fitting model included a random intercept for StudyID and fixed effects for lab and age. This model also yielded a main effect of age ( $F (1, 343) = 5.48$ , $p < 0.05$ ), carried by a larger effect in 6-month-olds than in newborns (estimate $- 0.27$ , $t (26) = - 2.33$ , $p = 0.02$ ). Over set 3, the best-fitting model included a random intercept for StudyID and fixed effects for lab and age. The model yielded a significant main effect of age [ $F (1,266) = 7.02$ , $p < 0.01$ ], carried by a larger effect in 6-month-olds than in newborns [estimate $- 0.36$ , $- (20) = - 2.63$ , $p = 0.01$ ].

For the N versus 0 comparison, for both sets 1 and 2, the best fitting models included a random intercept for StudyID and fixed effects for lab and age, but in both cases, no significant effects were found.

For the R versus N comparison, the best fitting models for sets 1 to 3 included a random intercept for StudyID and fixed effects for lab and age, but no significant effects were found.

3.2.2.

Deoxygenated hemoglobin (HbR)

For the R versus 0 comparison over set 1, the best fitting model included a random intercept for StudyID and fixed effects for lab and age, and it yielded a significant main effect of age ( $F (1,414) = 6.34$ , $p < 0.05$ ), carried by a stronger activation in 6-month-olds than in newborns (estimate 0.49, $t (19) = 2.49$ , $p < 0.05$ ). A similar result was found for Set 2 ( $F (1,348 = 8.45, p < 0.01)$ , carried again by a larger effect size in 6-month-olds than in newborns (estimate 0.78, $t (23) = 2.89$ , $p < 0.01$ ), as well as for set 3 ( $F (1,270) = 11.95$ , $p < 0.001$ , estimate 0 to 6 month = 1.2, $t (17) = 3.43$ , $p < 0.01$ ).

The best-fitting model comparing N versus 0 included a random intercept for StudyID and fixed effects for lab and age, but it yielded no significant effects for either set 1 or set 2.

For the R versus N comparison, the best fitting model similarly included a random intercept for StudyID and fixed effects for lab and age, but it did not yield significant effects for any of the sets of data.

4. Discussion

We have conducted meta-analyses and mixed-effects modeling-based inferential statistics to determine whether effect sizes were replicable in a sample of 23 NIRS studies investigating 0- to 9-month-old infants’ abilities to process repetition- and diversity-based regularities in linguistic and nonlinguistic auditory and visual sequences. In addition to quantifying effect sizes and their variation across studies, we have investigated whether these are modulated by different factors such as the age of participants or the laboratory, with the latter standing as a proxy for a set of various dimensions along which laboratories differed, e.g., NIRS machine, population characteristics, etc. We obtained NIRS data from 12 published and 11 unpublished studies. The 23 studies involved a total of 487 infants, aged between 0 and 9 months, tested in four different countries (Canada, France, Italy, and USA). The purpose of the study was primarily methodological, i.e., to test the replicability of experimental effects in infant NIRS research, an issue that has not been addressed in the developmental neuroscience literature, but which is essential for establishing NIRS as a reliable tool for infant brain imaging.

We tested three comparisons—infants’ responses to repetition-based regularities with respect to baseline (R versus 0), infants’ responses to diversity-based regularities with respect to baseline (N versus 0), and infants’ responses to repetition- versus diversity-based regularities (R versus N)—as all three are relevant to various aspects of language development.²⁶

We used two statistical approaches: a meta-analytic one, which addresses variability and replicability at the study level, and linear mixed-effects modeling, which tests the significance of the factors of interest (study, lab, and age) over individual effect sizes, while taking into account the nested nature of sampling (study nested in lab) by the random effects structure. Importantly, in both cases, we used effect sizes (study-level effect sizes for meta-analyses and individual-level effect sizes for the mixed-effects models) and not hemoglobin concentrations as dependent variables. In other words, our analyses address not only whether infants respond to a certain structural regularity but also, and more importantly, whether and how the responses vary across studies and laboratories. If these factors are found to modulate effect sizes, then there is considerable methodological variation across studies, which may undermine replicability and raise concerns about the reliability of infant NIRS methodology.

4.1.

Factors of Variability

4.1.1.

Cross-study variability and replicability

Our first and most important result is, therefore, that the factors/moderators study and laboratory never showed any significant effects for any of the comparisons in any of the analyses, meta-analytic or inferential. Effect sizes were statistically indistinguishable across different studies and laboratories. Importantly, this was true despite the fact that, as reported in the publications, studies and labs differed considerably in several factors known to impact NIRS data as well as along many undocumented factors that have less well known impacts on NIRS data (e.g., experimenter characteristics, features of the testing space, season/time of day of testing, etc.).

Estimates of bias, of which we computed the fail-safe-N and funnel plots, also confirm that the effects are robust and show no particular biases related to the selective publication of positive results or otherwise.

Our results, therefore, suggest that infants’ ability to process structural regularities in linguistic and nonlinguistic sequences can be tested in comparable and replicable ways across studies conducted with different NIRS machines, near-infrared light wavelengths, infants with different hair quality and color, etc. These results show, for the first time, that NIRS studies replicate robustly, even with the youngest infants.

4.1.2.

Theoretically relevant sources of variability

In addition to the above-discussed methodological factors, for which we expected no significant variation if infant NIRS studies were to be reliably replicable, we also included two additional factors in our analyses, age and repetition position, which have been suggested to modulate infants’ responses to structural regularities in theoretically relevant ways.

We found that age indeed modulated effect sizes. In the meta-analyses for the R versus N comparison, age showed a trend toward significance, with infants exhibiting a decreasing difference in their responses to repetition- versus diversity-based regularities with increasing age. These results are in accordance with the literature²⁶ and derive from the response to repetition-based regularities remaining stable across developmental time with the response to diversity-based regularities increasing (Figs. 7 and 8). Although this tendency was numerically present in all sets of the data, it only approached significance in Set 3, i.e., studies using trisyllabic speech sequences with adjacent repetitions because, in this restricted set, variability attributable to other factors was reduced; thus the effect of age could surface more clearly. In the mixed-effects models, we also observed a significant main effect of age. However, in this case, we saw an age-related increase in the responses to repetition-based sequences compared with the baseline, with 6-month-olds’ effect sizes being greater than those of newborns. The effect was marginal over the entire set of studies, but became significant for sets 2 and 3, which were more homogeneous and thus had less variance. Although it may seem that the results of the meta-analyses and the mixed-effects models diverge, it needs to be remembered that the mixed-effects models did not include data from 7- and 9-month-olds, so this is only apparent given the lack of trial-by-trial data for those studies. As a result, the meta-analyses have detected a larger developmental trend, which the mixed-effects models could not find in the absence of relevant data. Relevantly, however, the increase in effect sizes between 0 and 6 months in the R versus 0 comparison is also present numerically in the meta-analyses, and indeed its magnitude increases as the sets of studies get more homogeneous.

Fig. 7

Distributions of individual effect sizes for the three comparisons of interest (R versus 0, N versus 0, and R versus N) for HbO (upper panel) and HbR (lower panel).

Fig. 8

Grand average hemodynamic responses across studies in the left temporal area for the three comparisons of interest (R versus 0, N versus 0, and R versus N).

We thus observed two developmental trends: a decrease in the differential response to repetition- versus diversity-based regularities, especially starting at 7 months, and an increase in the response to repetitions. Both were particularly strong for the least variable, most homogenous set of studies, which were those that tested adjacent repetitions in trisyllabic speech sequences. Whether similar trends may also be observed for other stimuli, e.g., for visual sequences, could not be determined as the number of studies using other types of stimuli was insufficient. The two observed developmental trends converge with existing findings as infants’ ability to detect repetitions has been shown to improve during the second half of the first year of life at the behavioral level,³⁵^,³⁶ and the differential response between repetition-based and diversity-based patterns has been suggested to decrease by 6 months of age.²⁶

By contrast, our results show no evidence for repetition position to modulate the effect sizes of infants’ responses in any of the analyses. This effect is indeed relatively weak, and has not systematically been found in all behavioral studies.²⁵^,²⁶^,³⁷ We may thus not have sufficient statistical power to detect it as a lot fewer studies included in our sample tested initial repetition compared with the final ones, or the effect may genuinely be absent at the neural level.

4.2.

Observed Effect Sizes

Our results have established that, at least for the linguistic tasks under investigation, effects obtained in infants using NIRS replicate across studies and labs, while they are modulated by developmental factors. It is thus now relevant to quantify the magnitude of these effects. Effect sizes are indeed becoming increasingly important statistical measures as experimental neuroscience moves away from the much criticized and highly restrictive interpretation of study results in terms of the strict dichotomy of significant and nonsignificant $p$ -values.⁴⁰^–⁴²

The effect size of infants’ responses to repetition-based regularities compared with the baseline was found to be close to 0.3, with only slight differences as a function of study sets and participants’ age, as discussed above. This is thus a moderate effect. Responses to diversity-based regularities tend to be weaker, varying between 0.14 and 0.20, i.e., they are small-to-moderate effects. The effects of the differential response between the two regularities show considerable variation and tend to be small, between 0.04 and 0.19.

Small-to-moderate effect sizes are quite common in the developmental literature, even in technically less challenging behavioral tasks.²³^,³³ That infant neuroimaging, often presenting greater challenges both for researchers and participants, produces similar effect sizes can, therefore, be considered an important finding, demonstrating the reliability of infant NIRS methodology.

4.3.

Limitations and Future Directions

Our study sought to quantify the sizes and variability of the effects found in a sample of NIRS studies with young infants to assess their replicability. Consequently, we needed to make some methodological choices to reduce the dimensionality of the data and to facilitate comparisons across studies. The most important of these was the decision to focus on an anatomically defined, predetermined region of interest (ROI), the left temporal area. This choice was motivated by various factors. First, this area has been shown to be strongly involved in speech and language processing from the earliest ages, as shown by converging evidence from various brain imaging modalities such as fMRI⁴³ and NIRS.⁴⁴ Second, this area showed the strongest activation in the first NIRS study testing the processing of repetition- and diversity-based regularities in newborns.²⁸ Third, given the headgears and optode configurations used in the current studies, channels overlaid on the temporal areas seemed to vary less with respect to the underlying brain anatomy than other channels, even despite age-related changes in head and brain size, as shown by localization analyses.³¹^,³² This region of interest is thus useful for obtaining comparable data across many different studies and so was well suited for our purposes. However, for a more theoretically oriented meta-analysis, i.e., if the goal is to better understand the neural mechanisms underlying infants’ rule learning abilities, it will be more suitable to derive functionally-based, data-driven ROIs specifically for each study as it may be the case that the strongest effects are found in different ROIs in different studies, e.g., as a function of age-related changes or changes in the nature or sensory modality of the stimuli. Using such an approach is beyond the scope of the current study, but future research may address this issue.

5. Conclusion

The replication crisis has raised serious questions about the reliability and robustness of empirical findings in many disciplines from sociology and economics to psychology and neuroscience. NIRS, being a relatively recent brain imaging technique, with little standardization in research and analysis practices across laboratories, is in need of a systematic assessment of study replicability. In a meta-analysis of 23 NIRS studies on infants’ rule-learning abilities, we have shown, for the first time, that effects are robust and replicable across different studies and laboratories, with small-to-moderate effect sizes, as is typical in developmental research. This approach also allowed us to identify moderating factors relevant to theory-building, such as developmental age.

Disclosures

No conflicts of interest, financial or otherwise, are declared by the authors

Acknowledgments

This study was funded by the ERC Consolidator Grant “BabyRhythm” nr. 773202 to Judit Gervain, the Marie Curie Individual Fellowship EF-ST “BabyMindReader” nr. 101031716 to Jessica Gemignani, and the Basque Foundation for Science Ikerbasque and the Spanish Ministry of Science and Innovation [grant nr. PID2019-105100RJ-I00] to Irene de la Cruz-Pavía.

References

1.

R. N. Aslin, “Questioning the questions that have been asked about the infant brain using near-infrared spectroscopy,” Cogn. Neuropsychol., 29 (1–2), 7 –33 https://doi.org/10.1080/02643294.2012.654773 COGNEP 1464-0627 (2012). Google Scholar

2.

J. Gervain et al., “Near-infrared spectroscopy: a report from the McDonnell infant methodology consortium,” Dev. Cognit. Neurosci., 1 (1), 22 –46 https://doi.org/10.1016/j.dcn.2010.07.004 (2011). Google Scholar

3.

B. A. Nosek, J. R. Spies and M. Motyl, “Scientific utopia: II. Restructuring incentives and practices to promote truth over publishability,” Perspect. Psychol. Sci., 7 (6), 615 –631 https://doi.org/10.1177/1745691612459058 (2012). Google Scholar

4.

M. R. Munafò et al., “A manifesto for reproducible science,” Nat. Hum. Behav., 1 0021 https://doi.org/10.1038/s41562-016-0021 (2017). Google Scholar

5.

Open Science Collaboration, “Estimating the reproducibility of psychological science,” Science, 349 (6251), aac4716 https://doi.org/10.1126/science.aac4716 SCIEAS 0036-8075 (2015). Google Scholar

6.

R. A. Klein et al., “Investigating variation in replicability,” Soc. Psychol., 45 142 –152 https://doi.org/10.1027/1864-9335/a000178 (2014). Google Scholar

7.

M. A. Yücel et al., “Best practices for fNIRS publications,” Neurophotonics, 8 (1), 012101 https://doi.org/10.1117/1.NPh.8.1.012101 (2021). Google Scholar

8.

A. Pifferi et al., “Mechanically switchable solid inhomogeneous phantom for performance tests in diffuse imaging and spectroscopy,” J. Biomed. Opt., 20 (12), 121304 https://doi.org/10.1117/1.JBO.20.12.121304 JBOPFO 1083-3668 (2015). Google Scholar

9.

P. Lanka et al., “The BITMAP exercise–a multi-laboratory performance assessment campaign of diffuse optical instrumentation,” in Eur. Conf. Biomed. Opt., 11074_44 (2019). Google Scholar

10.

J. Gemignani and J. Gervain, “Comparing different pre-processing routines for infant fNIRS data,” Dev. Cognit. Neurosci., 48 100943 https://doi.org/10.1016/j.dcn.2021.100943 (2021). Google Scholar

11.

R. Di Lorenzo et al., “Recommendations for motion correction of infant fNIRS data applicable to multiple data sets and acquisition systems,” NeuroImage, 200 511 –527 https://doi.org/10.1016/j.neuroimage.2019.06.056 NEIMEF 1053-8119 (2019). Google Scholar

12.

P. Pinti et al., “Current status and issues regarding pre-processing of fNIRS neuroimaging data: an investigation of diverse signal filtering methods within a general linear model framework,” Front. Hum. Neurosci., 12 505 https://doi.org/10.3389/fnhum.2018.00505 (2019). Google Scholar

13.

J. W. Barker, A. Aarabi and T. J. Huppert, “Autoregressive model based algorithm for correcting motion and serially correlated errors in fNIRS,” Biomed. Opt. Express, 4 (8), 1366 –1379 https://doi.org/10.1364/BOE.4.001366 BOEICL 2156-7085 (2013). Google Scholar

14.

T. J. Huppert, “Commentary on the statistical properties of noise and its implication on general linear models in functional near-infrared spectroscopy,” Neurophotonics, 3 (1), 010401 https://doi.org/10.1117/1.NPh.3.1.010401 (2016). Google Scholar

15.

M. Borenstein et al., Introduction to Meta-Analysis, John Wiley & Sons( (2021). Google Scholar

16.

M. Lewis et al., “A quantitative synthesis of early language acquisition using meta-analysis,” (2016). https://doi.org/10.31234/osf.io/htsjm Google Scholar

17.

S. Tsuji and A. Cristia, “Perceptual attunement in vowels: a meta-analysis,” Dev. Psychobiol., 56 (2), 179 –191 https://doi.org/10.1002/dev.21179 DEPBA5 0012-1630 (2014). Google Scholar

18.

M. A. Qureshi, “A meta-analysis: age and second language grammar acquisition,” System, 60 147 –160 https://doi.org/10.1016/j.system.2016.06.001 (2016). Google Scholar

19.

V. P. Kwok et al., “A meta-analytic study of the neural systems for auditory processing of lexical tones,” Front. Hum. Neurosci., 11 375 https://doi.org/10.3389/fnhum.2017.00375 (2017). Google Scholar

20.

E. C. Ferstl et al., “The extended language network: a meta-analysis of neuroimaging studies on text comprehension,” Hum. Brain Mapp., 29 (5), 581 –593 https://doi.org/10.1002/hbm.20422 HBRME7 1065-9471 (2008). Google Scholar

21.

C. Dunst, E. Gorman and D. Hamby, “Preference for infant-directed speech in preverbal young children,” Cent. Early Lit. Learn., 5 (1), 1 –13 (2012). Google Scholar

22.

C. Bergmann and A. Cristia, “Development of infants’ segmentation of words from native speech: a meta-analytic approach,” Dev. Sci., 19 (6), 901 –917 https://doi.org/10.1111/desc.12341 (2016). Google Scholar

23.

ManyBabies Consortium, “Quantifying sources of variability in infancy research using the infant-directed speech preference,” Adv. Methods Pract. Psychol. Sci., 3 (1), 24 –52 https://doi.org/10.1177/2515245919900809 (2020). Google Scholar

24.

G. F. Marcus et al., “Rule learning by seven-month-old infants,” Science, 283 (5398), 77 –80 https://doi.org/10.1126/science.283.5398.77 SCIEAS 0036-8075 (1999). Google Scholar

25.

H. Rabagliati, B. Ferguson and C. Lew-Williams, “The profile of abstract rule learning in infancy: Meta-analytic and experimental evidence,” Dev. Sci., 22 (1), e12704 https://doi.org/10.1111/desc.12704 (2019). Google Scholar

26.

I. de la Cruz-Pavía and J. Gervain, “Infants’ perception of repetition-based regularities in speech: a look from the perspective of the same/different distinction,” Curr. Opin. Behav. Sci., 37 125 –132 https://doi.org/10.1016/j.cobeha.2020.11.014 (2021). Google Scholar

27.

J. B. Wagner et al., “Neural processing of repetition and non-repetition grammars in 7- and 9-month-old infants,” Front. Psychol., 2 168 https://doi.org/10.3389/fpsyg.2011.00168 1664-1078 (2011). Google Scholar

28.

J. Gervain et al., “The neonate brain detects speech structure,” Proc. Natl. Acad. Sci. U. S. A., 105 (37), 14222 –14227 https://doi.org/10.1073/pnas.0806530105 (2008). Google Scholar

29.

J. Gervain, I. Berent and J. F. Werker, “Binding at birth: The newborn brain detects identity relations and sequential position in speech,” J. Cognit. Neurosci., 24 564 –574 https://doi.org/10.1162/jocn_a_00157 JCONEO 0898-929X (2012). Google Scholar

30.

C. Bouchon, T. Nazzi and J. Gervain, “Hemispheric asymmetries in repetition enhancement and suppression effects in the newborn brain,” PLoS One, 10 (10), e0140160 https://doi.org/10.1371/journal.pone.0140160 POLNCL 1932-6203 (2015). Google Scholar

31.

I. Berent et al., “Infants differentially extract rules from language,” Sci. Rep., 11 20001 https://doi.org/10.1038/s41598-021-99539-8 SRCEC3 2045-2322 (2021). Google Scholar

32.

N. Abboub, T. Nazzi and J. Gervain, “Prosodic grouping at birth,” Brain Lang., 162 46 –59 https://doi.org/10.1016/j.bandl.2016.08.002 (2016). Google Scholar

33.

C. Bergmann et al., “Promoting replicability in developmental research through meta-analyses: Insights from language acquisition research,” Child Dev., 89 (6), 1996 –2009 https://doi.org/10.1111/cdev.13079 CHDEAW 0009-3920 (2018). Google Scholar

34.

W. Viechtbauer, “Conducting meta-analyses in R with the metafor package,” J. Stat. Software, 36 (3), 1 –48 https://doi.org/10.18637/jss.v036.i03 (2010). Google Scholar

35.

M. C. Frank et al., “Information from multiple modalities helps five-month-olds learn abstract rules,” Dev. Sci., 12 (4), 504 –509 https://doi.org/10.1111/j.1467-7687.2008.00794.x (2009). Google Scholar

36.

S. P. Johnson et al., “Abstract rule learning for visual sequences in 8-and 11-month-olds,” Infancy, 14 (1), 2 –18 https://doi.org/10.1080/15250000802569611 (2009). Google Scholar

37.

A. D. Endress, M. Nespor and J. Mehler, “Perceptual and memory constraints on language acquisition,” Trends Cognit. Sci., 13 (8), 348 –353 https://doi.org/10.1016/j.tics.2009.05.005 (2009). Google Scholar

38.

R. Rosenthal, “The file drawer problem and tolerance for null results,” Psychol. Bull., 86 (3), 638 –641 https://doi.org/10.1037/0033-2909.86.3.638 PSBUAI 0033-2909 (1979). Google Scholar

39.

C. B. Begg and M. Mazumdar, “Operating characteristics of a rank correlation test for publication bias,” Biometrics, 50 1088 –1101 https://doi.org/10.2307/2533446 BIOMB6 0006-341X (1994). Google Scholar

40.

G. Cumming, “The new statistics: why and how,” Psychol. Sci., 25 (1), 7 –29 https://doi.org/10.1177/0956797613504966 1467-9280 (2014). Google Scholar

41.

G. Cumming, Understanding the New Statistics: Effect sizes, Confidence Intervals, and Meta-Analysis, Routledge( (2013). Google Scholar

42.

I. Deviatko, “New data, new statistics: from reproducibility crisis toward new requirements to data analysis and presentation in social sciences,” Sotsiologicheskie Issled, 12 30 –38 https://doi.org/10.31857/S013216250003163-2 (2019). Google Scholar

43.

G. Dehaene-Lambertz, S. Dehaene and L. Hertz-Pannier, “Functional neuroimaging of speech perception in infants,” Science, 298 (5600), 2013 –2015 https://doi.org/10.1126/science.1077066 SCIEAS 0036-8075 (2002). Google Scholar

44.

M. Pena et al., “Sounds and silence: an optical topography study of language recognition at birth,” Proc. Natl. Acad. Sci. U. S. A., 100 (20), 11702 –11705 https://doi.org/10.1073/pnas.1934290100 (2003). Google Scholar

Biography

Jessica Gemignani is a Marie Skodowska-Curie Postdoctoral Fellow at the University of Padua, Italy. She received her degree in biomedical engineering from the University of Pisa in 2014 and her PhD in neurotechnology from Technical University of Berlin in 2019. Her research focuses on the development and validation of methods for the analysis of multivariate neuroimaging data acquired on newborns and infants, particularly in the field of language acquisition.

Irene de la Cruz-Pavía is a research fellow at Basque Foundation for Science Ikerbasque and JIN researcher at the Department of Linguistics and Basque Studies, University of the Basque Country UPV/EHU, Spain. Her research investigates linguistic processing throughout the lifespan from early infancy to old age, in monolingual and bilingual populations with a strong focus on the acquisition and processing of grammar.

Judit Gervain is a full professor at Developmental Psychology, University of Padua, Italy, and a senior research scientist at CNRS, France. Her research focuses early speech perception and language acquisition in typically and atypically developing infants. Her work is published in leading journals, such as Science Advances, Nature Communications, PNAS, and Current Biology. She is an associate editor at Developmental Science, Cognition, and Neurophotonics. Her work is currently funded by an ERC Consolidator Grant.

Biographies of the other authors are not available.

CC BY: © The Authors. Published by SPIE under a Creative Commons Attribution 4.0 International License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.

Citation Download Citation

Jessica Gemignani, Irene de la Cruz-Pavía, Anna Martinez, Caroline Nallet, Alessia Pasquini, Gaia Lucarini, Francesca Cavicchiolo, and Judit Gervain "Reproducibility of infant fNIRS studies: a meta-analytic approach," Neurophotonics 10(2), 023518 (10 March 2023). https://doi.org/10.1117/1.NPh.10.2.023518

Received: 30 July 2022; Accepted: 10 February 2023; Published: 10 March 2023

Access the abstract

JOURNAL ARTICLE
21 PAGES

DOWNLOAD PAPER SAVE TO MY LIBRARY