Indirect likelihood inference

Given a sample from a fully specified parametric model, let Zn be a given finite-dimensional statistic - for example, an initial estimator or a set of sample moments. We propose to (re-)estimate the parameters of the model by maximizing the likelihood of Zn. We call this the maximum indirect likelihood (MIL) estimator. We also propose a computationally tractable Bayesian version of the estimator which we refer to as a Bayesian Indirect Likelihood (BIL) estimator. In most cases, the density of the statistic will be of unknown form, and we develop simulated versions of the MIL and BIL estimators. We show that the indirect likelihood estimators are consistent and asymptotically normally distributed, with the same asymptotic variance as that of the corresponding efficient two-step GMM estimator based on the same statistic. However, our likelihood-based estimators, by taking into account the full finite-sample distribution of the statistic, are higher order efficient relative to GMM-type estimators. Furthermore, in many cases they enjoy a bias reduction property similar to that of the indirect inference estimator. Monte Carlo results for a number of applications including dynamic and nonlinear panel data models, a structural auction model and two DSGE models show that the proposed estimators indeed have attractive finite sample properties.

Cerebellar transcranial direct current stimulation modulates verbal working memory.

BACKGROUND: Neuroimaging studies show cerebellar activations in a wide range of cognitive tasks and patients with cerebellar lesions often present cognitive deficits suggesting a cerebellar role in higher-order cognition. OBJECTIVE: We used cathodal transcranial direct current stimulation (tDCS), known to inhibit neuronal excitability, over the cerebellum to investigate if cathodal tDCS impairs verbal working memory, an important higher-order cognitive faculty. METHOD: We tested verbal working memory as measured by forward and backward digit spans in 40 healthy young participants before and after applying cathodal tDCS (2 mA, stimulation duration 25 min) to the right cerebellum using a randomized, sham-controlled, double-blind, cross-over design. In addition, we tested the effect of cerebellar tDCS on word reading, finger tapping and a visually cued sensorimotor task. RESULTS: In line with lower digit spans in patients with cerebellar lesions, cerebellar tDCS reduced forward digit spans and blocked the practice dependent increase in backward digit spans. No effects of tDCS on word reading, finger tapping or the visually cued sensorimotor task were found. CONCLUSION: Our results support the view that the cerebellum contributes to verbal working memory as measured by forward and backward digit spans. Moreover, the induction of reversible "virtual cerebellar lesions" in healthy individuals by means of tDCS may improve our understanding of the mechanistic basis of verbal working memory deficits in patients with cerebellar lesions.

Zuckerman's reivsed alternative five-factor model: Validation of the Zuckerman-Kuhlman-Aluja Personality Questionnaire in four French-speaking countries

The aim of this study was to analyze the replicability of Zuckerman's revised Alternative Five-factor model in a French-speaking context by validating the Zuckerman-Kuhlman-Aluja Personality Questionnaire (ZKA-PQ) simultaneously in 4 French-speaking countries. The total sample was made up of 1,497 subjects from Belgium, Canada, France, and Switzerland. The internal consistencies for all countries were generally similar to those found for the normative U.S. and Spanish samples. A factor analysis confirmed that the normative structure replicated well and was stable within this French-speaking context. Moreover, multigroup confirmatory factor analyses have shown that the ZKA-PQ reaches scalar invariance across these 4 countries. Mean scores were slightly different for women and men, with women scoring higher on Neuroticism but lower on Sensation Seeking. Globally, mean score differences across countries were small. Overall, the ZKA-PQ seems an interesting alternative to assess both lower and higher order personality traits for applied or research purposes.

A randomized controlled trial of the effectiveness of a Computer-Assisted Cognitive Remediation (CACR) program in adolescents with psychosis or at high risk of psychosis: Short term and long term outcomes

The present study investigates the short- and long-term outcomes of a computer-assisted cognitive remediation (CACR) program in adolescents with psychosis or at high risk. 32 adolescents participated in a blinded 8-week randomized controlled trial of CACR treatment compared to computer games (CG). Clinical and neuropsychological evaluations were undertaken at baseline, at the end of the program and at 6-month. At the end of the program (n = 28), results indicated that visuospatial abilities (Repeatable Battery for the Assessment of Neuropsychological Status, RBANS; P = .005) improved signifi cantly more in the CACR group compared to the CG group. Furthermore, other cognitive functions (RBANS), psychotic symptoms (Positive and Negative Symptom Scale) and psychosocial functioning (Social and Occupational Functioning Assessment Scale) improved signifi cantly, but at similar rates, in the two groups. At long term (n = 22), cognitive abilities did not demonstrated any amelioration in the control group while, in the CACR group, signifi cant long-term improvements in inhibition (Stroop; P = .040) and reasoning (Block Design Test; P = .005) were observed. In addition, symptom severity (Clinical Global Improvement) decreased signifi cantly in the control group (P = .046) and marginally in the CACR group (P = .088). To sum up, CACR can be successfully administered in this population. CACR proved to be effective over and above CG for the most intensively trained cognitive ability. Finally, on the long-term, enhanced reasoning and inhibition abilities, which are necessary to execute higher-order goals or to adapt behavior to the ever-changing environment, were observed in adolescents benefi ting from a CACR.

Musical training intensity yields opposite effects on grey matter density in cognitive versus sensorimotor networks.

Using optimized voxel-based morphometry, we performed grey matter density analyses on 59 age-, sex- and intelligence-matched young adults with three distinct, progressive levels of musical training intensity or expertise. Structural brain adaptations in musicians have been repeatedly demonstrated in areas involved in auditory perception and motor skills. However, musical activities are not confined to auditory perception and motor performance, but are entangled with higher-order cognitive processes. In consequence, neuronal systems involved in such higher-order processing may also be shaped by experience-driven plasticity. We modelled expertise as a three-level regressor to study possible linear relationships of expertise with grey matter density. The key finding of this study resides in a functional dissimilarity between areas exhibiting increase versus decrease of grey matter as a function of musical expertise. Grey matter density increased with expertise in areas known for their involvement in higher-order cognitive processing: right fusiform gyrus (visual pattern recognition), right mid orbital gyrus (tonal sensitivity), left inferior frontal gyrus (syntactic processing, executive function, working memory), left intraparietal sulcus (visuo-motor coordination) and bilateral posterior cerebellar Crus II (executive function, working memory) and in auditory processing: left Heschl's gyrus. Conversely, grey matter density decreased with expertise in bilateral perirolandic and striatal areas that are related to sensorimotor function, possibly reflecting high automation of motor skills. Moreover, a multiple regression analysis evidenced that grey matter density in the right mid orbital area and the inferior frontal gyrus predicted accuracy in detecting fine-grained incongruities in tonal music.

Time-Delayed Theory of the Neolithic Transition in Europe

The classical wave-of-advance model of the neolithic transition (i.e., the shift from hunter-gatherer to agricultural economies) is based on Fisher's reaction-diffusion equation. Here we present an extension of Einstein's approach to Fickian diffusion, incorporating reaction terms. On this basis we show that second-order terms in the reaction-diffusion equation, which have been neglected up to now, are not in fact negligible but can lead to important corrections. The resulting time-delayed model agrees quite well with observations

Noninvasive functional and structural exploration of human subcortical structures with high resolution

Projecte de recerca elaborat a partir d’una estada al Max Planck Institute for Human Cognitive and Brain Sciences, Alemanya, entre 2010 i 2012. El principal objectiu d’aquest projecte era estudiar en detall les estructures subcorticals, en concret, el rol dels ganglis basals en control cognitiu durant processament lingüístic i no-lingüístic. Per tal d’assolir una diferenciació minuciosa en els diferents nuclis dels ganglis basals s’utilitzà ressonància magnètica d’ultra-alt camp i alta resolució (7T-MRI). El còrtex prefrontal lateral i els ganglis basals treballant conjuntament per a mitjançar memòria de treball i la regulació “top-down” de la cognició. Aquest circuit regula l’equilibri entre respostes automàtiques i d’alt-ordre cognitiu. Es crearen tres condicions experimentals principals: frases/seqüències noambigües, no-gramatical i ambigües. Les frases/seqüències no-ambigües haurien de provocar una resposta automàtica, mentre les frases/seqüències ambigües i no-gramaticals produïren un conflicte amb la resposta automàtica, i per tant, requeririen una resposta de d’alt-ordre cognitiu. Dins del domini de la resposta de control, la ambigüitat i no-gramaticalitat representen dues dimensions diferents de la resolució de conflicte, mentre per una frase/seqüència temporalment ambigua existeix una interpretació correcte, aquest no és el cas per a les frases/seqüències no-gramaticals. A més, el disseny experimental incloïa una manipulació lingüística i nolingüística, la qual posà a prova la hipòtesi que els efectes són de domini-general; així com una manipulació semàntica i sintàctica que avaluà les diferències entre el processament d’ambigüitat/error “intrínseca” vs. “estructural”. Els resultats del primer experiment (sintax-lingüístic) mostraren un gradient rostroventralcaudodorsal de control cognitiu dins del nucli caudat, això és, les regions més rostrals sostenint els nivells més alts de processament cognitiu

Emotionally negative stimuli can overcome attentional deficits in patients with visuo-spatial hemineglect.

Left unilateral spatial neglect resulting from right brain damage is characterized by loss of awareness for stimuli in the contralesional side of space, despite intact visual pathways. We examined using fMRI whether patients with neglect are more likely to consciously detect in the neglected hemifield, emotionally negative complex scenes rather than visually similar neutral pictures and if so, what neural mechanisms mediate this effect. Photographs of emotional and neutral scenes taken from the IAPS were presented in a divided visual field paradigm. As expected, the detection rate for emotional stimuli presented in the neglected field was higher than for neutral ones. Successful detection of emotional scenes as opposed to neutral stimuli in the left visual field (LVF) produced activations in the parahippocampal and anterior cingulate areas in the right hemisphere. Detection of emotional stimuli presented in the intact right visual field (RVF) activated a distributed network of structures in the left hemisphere, including anterior and posterior cingulate cortex, insula, as well as visual striate and extrastriate areas. LVF-RVF contrasts for emotional stimuli revealed activations in right and left attention related prefrontal areas whereas RVF-LVF comparison showed activations in the posterior cingulate and extrastriate visual cortex in the left hemisphere. An additional analysis contrasting detected vs. undetected emotional LVF stimuli showed involvement of left anterior cingulate, right frontal and extrastriate areas. We hypothesize that beneficial role of emotion in overcoming neglect is achieved by activation of frontal and limbic lobe networks, which provide a privileged access of emotional stimuli to attention by top-down modulation of processing in the higher-order extrastriate visual areas. Our results point to the importance of top-down regulatory role of the frontal attentional systems, which might enhance visual activations and lead to greater salience of emotional stimuli for perceptual awareness.

Regional child care availability and fertility decisions in Spain

In this paper I explore two hypotheses: (1) Formal child care availability for children under three has a positive effect across contexts, according to the degree of adaptation of social institutions to changes in gender roles. Event history models with regional fixed effects are applied to data from the European Community Household Panel (1994-2001). The results show a significant and positive effect of regional day care availability on both, first and higher order births, while results are consistent with the second hypothesis only for second or higher order births.

Towards robust network based complex systems : from evolutionary cellular automata to biological models

Abstract Sitting between your past and your future doesn't mean you are in the present. Dakota Skye Complex systems science is an interdisciplinary field grouping under the same umbrella dynamical phenomena from social, natural or mathematical sciences. The emergence of a higher order organization or behavior, transcending that expected of the linear addition of the parts, is a key factor shared by all these systems. Most complex systems can be modeled as networks that represent the interactions amongst the system's components. In addition to the actual nature of the part's interactions, the intrinsic topological structure of underlying network is believed to play a crucial role in the remarkable emergent behaviors exhibited by the systems. Moreover, the topology is also a key a factor to explain the extraordinary flexibility and resilience to perturbations when applied to transmission and diffusion phenomena. In this work, we study the effect of different network structures on the performance and on the fault tolerance of systems in two different contexts. In the first part, we study cellular automata, which are a simple paradigm for distributed computation. Cellular automata are made of basic Boolean computational units, the cells; relying on simple rules and information from- the surrounding cells to perform a global task. The limited visibility of the cells can be modeled as a network, where interactions amongst cells are governed by an underlying structure, usually a regular one. In order to increase the performance of cellular automata, we chose to change its topology. We applied computational principles inspired by Darwinian evolution, called evolutionary algorithms, to alter the system's topological structure starting from either a regular or a random one. The outcome is remarkable, as the resulting topologies find themselves sharing properties of both regular and random network, and display similitudes Watts-Strogtz's small-world network found in social systems. Moreover, the performance and tolerance to probabilistic faults of our small-world like cellular automata surpasses that of regular ones. In the second part, we use the context of biological genetic regulatory networks and, in particular, Kauffman's random Boolean networks model. In some ways, this model is close to cellular automata, although is not expected to perform any task. Instead, it simulates the time-evolution of genetic regulation within living organisms under strict conditions. The original model, though very attractive by it's simplicity, suffered from important shortcomings unveiled by the recent advances in genetics and biology. We propose to use these new discoveries to improve the original model. Firstly, we have used artificial topologies believed to be closer to that of gene regulatory networks. We have also studied actual biological organisms, and used parts of their genetic regulatory networks in our models. Secondly, we have addressed the improbable full synchronicity of the event taking place on. Boolean networks and proposed a more biologically plausible cascading scheme. Finally, we tackled the actual Boolean functions of the model, i.e. the specifics of how genes activate according to the activity of upstream genes, and presented a new update function that takes into account the actual promoting and repressing effects of one gene on another. Our improved models demonstrate the expected, biologically sound, behavior of previous GRN model, yet with superior resistance to perturbations. We believe they are one step closer to the biological reality.

Distribution of neuronal and metabolic markers in the human auditory cortex

SUMMARY The human auditory cortex, located on the supratemporal plane of the temporal lobe, is divided in a primary auditory area and several non-primary areas surrounding it. These different areas show anatomical and functional differences. Many studies have focussed on auditory areas in non-human primates, using investigation techniques such as electrophysiological recordings, tracing of neural connections, or immunohistochemical and histochemical staining. Some of these studies have suggested parallel and hierarchical organization of the cortical auditory areas as well as subcortical auditory relays. In humans, only few studies have investigated these regions immunohistochemically, but activation and lesion studies speak in favour of parallel and hierarchical organization, very similar to that of non-human primates. Calcium-binding proteins and metabolic markers were used to investigate possible correlates of hierarchical and parallel organization in man. Calcium-binding proteins, parvalbumin, calretinin and calbindin, modulate the concentration of intracellular free calcium ions and were found in distinct subpopulations of GABAergic neurons in non-human primates species. In our study, their distribution showed several differences between auditory areas: the primary auditory area was darkly stained for both parvalbumin and calbindin, and their expression rapidly decreased while moving away from the primary area. This staining pattern suggests a hierarchical organization of the areas, in which the more darkly stained areas could correspond to an earlier integration level and the areas showing light staining may correspond to higher level integration areas. Parallel organization of primary and non-primary auditory areas was suggested by the complementarity, within a given area, between parvalbumin and calbindin expression across layers. To investigate the possible differences in the energetic metabolism of the cortical auditory areas, several metabolic markers were used: cytochrome oxidase and LDH1 were used as oxidative metabolism markers and LDH5 was used as glycolytic metabolism marker. The results obtained show a difference in the expression of enzymes involved in oxidative metabolism between areas. In the primary auditory area the oxidative metabolism markers were maximally expressed in layer IV. In contrast, higher order areas showed maximal staining in supragranular layers. The expression of LDH5 varied in patches, but did not differ between the different hierarchical auditory areas. The distribution of the two LDH enzymes isoforms also provides information about cellular aspects of metabolic organization, since neurons expressed the LDH1 isoform whereas astrocytes express primarily LDH5, but some astrocytes also contained the LDH1 isoform. This cellular distribution pattern supports the hypothesis of the existence of an astrocyte-neuron lactate shuttle, previously suggested in rodent studies, and in particular of lactate transfer from astrocytes, which produce lactate from the glucose obtained from the circulation, to neurons that use lactate as energy substrate. In conclusion, the hypothesis of parallel and hierarchical organization of the auditory areas can be supported by CaBPs, cytochrome oxidase and LDH1 distribution. Moreover, the two LDHs cellular distribution pattern support the hypothesis of an astrocyte-neuron lactate shuttle in human cortex.

Neocortical maturation during adolescence: change in neuronal soma dimension.

During adolescence, cognitive abilities increase robustly. To search for possible related structural alterations of the cerebral cortex, we measured neuronal soma dimension (NSD = width times height), cortical thickness and neuronal densities in different types of neocortex in post-mortem brains of five 12-16 and five 17-24 year-olds (each 2F, 3M). Using a generalized mixed model analysis, mean normalized NSD comparing the age groups shows layer-specific change for layer 2 (p < .0001) and age-related differences between categorized type of cortex: primary/primary association cortex (BA 1, 3, 4, and 44) shows a generalized increase; higher-order regions (BA 9, 21, 39, and 45) also show increase in layers 2 and 5 but decrease in layers 3, 4, and 6 while limbic/orbital cortex (BA 23, 24, and 47) undergoes minor decrease (BA 1, 3, 4, and 44 vs. BA 9, 21, 39, and 45: p = .036 and BA 1, 3, 4, and 44 vs. BA 23, 24, and 47: p = .004). These data imply the operation of cortical layer- and type-specific processes of growth and regression adding new evidence that the human brain matures during adolescence not only functionally but also structurally.

Further insights on the French WISC-IV factor structure through Bayesian structural equation modeling

The interpretation of the Wechsler Intelligence Scale for Children-Fourth Edition (WISC-IV) is based on a 4-factor model, which is only partially compatible with the mainstream Cattell-Horn-Carroll (CHC) model of intelligence measurement. The structure of cognitive batteries is frequently analyzed via exploratory factor analysis and/or confirmatory factor analysis. With classical confirmatory factor analysis, almost all crossloadings between latent variables and measures are fixed to zero in order to allow the model to be identified. However, inappropriate zero cross-loadings can contribute to poor model fit, distorted factors, and biased factor correlations; most important, they do not necessarily faithfully reflect theory. To deal with these methodological and theoretical limitations, we used a new statistical approach, Bayesian structural equation modeling (BSEM), among a sample of 249 French-speaking Swiss children (8-12 years). With BSEM, zero-fixed cross-loadings between latent variables and measures are replaced by approximate zeros, based on informative, small-variance priors. Results indicated that a direct hierarchical CHC-based model with 5 factors plus a general intelligence factor better represented the structure of the WISC-IV than did the 4-factor structure and the higher order models. Because a direct hierarchical CHC model was more adequate, it was concluded that the general factor should be considered as a breadth rather than a superordinate factor. Because it was possible for us to estimate the influence of each of the latent variables on the 15 subtest scores, BSEM allowed improvement of the understanding of the structure of intelligence tests and the clinical interpretation of the subtest scores.

Solving heterogeneous-agent models by projection and perturbation

The paper proposes a numerical solution method for general equilibrium models with a continuum of heterogeneous agents, which combines elements of projection and of perturbation methods. The basic idea is to solve first for the stationary solutionof the model, without aggregate shocks but with fully specified idiosyncratic shocks. Afterwards one computes a first-order perturbation of the solution in the aggregate shocks. This approach allows to include a high-dimensional representation of the cross-sectional distribution in the state vector. The method is applied to a model of household saving with uninsurable income risk and liquidity constraints. The model includes not only productivity shocks, but also shocks to redistributive taxation, which cause substantial short-run variation in the cross-sectional distribution of wealth. If those shocks are operative, it is shown that a solution method based on very few statistics of the distribution is not suitable, while the proposed method can solve the model with high accuracy, at least for the case of small aggregate shocks. Techniques are discussed to reduce the dimension of the state space such that higher order perturbations are feasible.Matlab programs to solve the model can be downloaded.

Tax reforms and inequality: theoretical and empirical implications

In this paper we examine the effect of tax policy on the relationship between inequality and growth in a two-sector non-scale model. With non-scale models, the longrun equilibrium growth rate is determined by technological parameters and it is independent of macroeconomic policy instruments. However, this fact does not imply that fiscal policy is unimportant for long-run economic performance. It indeed has important effects on the different levels of key economic variables such as per capita stock of capital and output. Hence, although the economy grows at the same rate across steady states, the bases for economic growth may be different.The model has three essential features. First, we explicitly model skill accumulation, second, we introduce government finance into the production function, and we introduce an income tax to mirror the fiscal events of the 1980¿s and 1990¿s in the US. The fact that the non-scale model is associated with higher order dynamics enables it to replicate the distinctly non-linear nature of inequality in the US with relative ease. The results derived in this paper attract attention to the fact that the non-scale growth model does not only fit the US data well for the long-run (Jones, 1995b) but also that it possesses unique abilities in explaining short term fluctuations of the economy. It is shown that during transition the response of the relative simulated wage to changes in the tax code is rather non-monotonic, quite in accordance to the US inequality pattern in the 1980¿s and early 1990¿s.More specifically, we have analyzed in detail the dynamics following the simulation of an isolated tax decrease and an isolated tax increase. So, after a tax decrease the skill premium follows a lower trajectory than the one it would follow without a tax decrease. Hence we are able to reduce inequality for several periods after the fiscal shock. On the contrary, following a tax increase, the evolution of the skill premium remains above the trajectory carried on by the skill premium under a situation with no tax increase. Consequently, a tax increase would imply a higher level of inequality in the economy