Synchronizability of EEG-Based Functional Networks in Early Alzheimer's Disease.

Recently graph theory and complex networks have been widely used as a mean to model functionality of the brain. Among different neuroimaging techniques available for constructing the brain functional networks, electroencephalography (EEG) with its high temporal resolution is a useful instrument of the analysis of functional interdependencies between different brain regions. Alzheimer's disease (AD) is a neurodegenerative disease, which leads to substantial cognitive decline, and eventually, dementia in aged people. To achieve a deeper insight into the behavior of functional cerebral networks in AD, here we study their synchronizability in 17 newly diagnosed AD patients compared to 17 healthy control subjects at no-task, eyes-closed condition. The cross-correlation of artifact-free EEGs was used to construct brain functional networks. The extracted networks were then tested for their synchronization properties by calculating the eigenratio of the Laplacian matrix of the connection graph, i.e., the largest eigenvalue divided by the second smallest one. In AD patients, we found an increase in the eigenratio, i.e., a decrease in the synchronizability of brain networks across delta, alpha, beta, and gamma EEG frequencies within the wide range of network costs. The finding indicates the destruction of functional brain networks in early AD.

Selecting statistical models to study the relationship between soybean yield and soil physical properties

Statistical models allow the representation of data sets and the estimation and/or prediction of the behavior of a given variable through its interaction with the other variables involved in a phenomenon. Among other different statistical models, are the autoregressive state-space models (ARSS) and the linear regression models (LR), which allow the quantification of the relationships among soil-plant-atmosphere system variables. To compare the quality of the ARSS and LR models for the modeling of the relationships between soybean yield and soil physical properties, Akaike's Information Criterion, which provides a coefficient for the selection of the best model, was used in this study. The data sets were sampled in a Rhodic Acrudox soil, along a spatial transect with 84 points spaced 3 m apart. At each sampling point, soybean samples were collected for yield quantification. At the same site, soil penetration resistance was also measured and soil samples were collected to measure soil bulk density in the 0-0.10 m and 0.10-0.20 m layers. Results showed autocorrelation and a cross correlation structure of soybean yield and soil penetration resistance data. Soil bulk density data, however, were only autocorrelated in the 0-0.10 m layer and not cross correlated with soybean yield. The results showed the higher efficiency of the autoregressive space-state models in relation to the equivalent simple and multiple linear regression models using Akaike's Information Criterion. The resulting values were comparatively lower than the values obtained by the regression models, for all combinations of explanatory variables.

Model-Free Reconstruction of Excitatory Neuronal Connectivity from Calcium Imaging Signals

A systematic assessment of global neural network connectivity through direct electrophysiological assays has remained technically infeasible, even in simpler systems like dissociated neuronal cultures. We introduce an improved algorithmic approach based on Transfer Entropy to reconstruct structural connectivity from network activity monitored through calcium imaging. We focus in this study on the inference of excitatory synaptic links. Based on information theory, our method requires no prior assumptions on the statistics of neuronal firing and neuronal connections. The performance of our algorithm is benchmarked on surrogate time series of calcium fluorescence generated by the simulated dynamics of a network with known ground-truth topology. We find that the functional network topology revealed by Transfer Entropy depends qualitatively on the time-dependent dynamic state of the network (bursting or non-bursting). Thus by conditioning with respect to the global mean activity, we improve the performance of our method. This allows us to focus the analysis to specific dynamical regimes of the network in which the inferred functional connectivity is shaped by monosynaptic excitatory connections, rather than by collective synchrony. Our method can discriminate between actual causal influences between neurons and spurious non-causal correlations due to light scattering artifacts, which inherently affect the quality of fluorescence imaging. Compared to other reconstruction strategies such as cross-correlation or Granger Causality methods, our method based on improved Transfer Entropy is remarkably more accurate. In particular, it provides a good estimation of the excitatory network clustering coefficient, allowing for discrimination between weakly and strongly clustered topologies. Finally, we demonstrate the applicability of our method to analyses of real recordings of in vitro disinhibited cortical cultures where we suggest that excitatory connections are characterized by an elevated level of clustering compared to a random graph (although not extreme) and can be markedly non-local.

Comparison of the power spectral changes of the voluntary surface electromyogram and M wave during intermittent maximal voluntary contractions.

INTRODUCTION: To compare the power spectral changes of the voluntary surface electromyogram (sEMG) and of the compound action potential (M wave) in the vastus medialis and vastus lateralis muscles during fatiguing contractions. METHODS: Interference sEMG and force were recorded during 48 intermittent 3-s isometric maximal voluntary contractions (MVC) from 13 young, healthy subjects. M waves and twitches were evoked using supramaximal femoral nerve stimulation between the successive MVCs. Mean frequency (F mean), and median frequency were calculated from the sEMG and M waves. Muscle fiber conduction velocity (MFCV) was computed by cross-correlation. RESULTS: The power spectral shift to lower frequencies was significantly greater for the voluntary sEMG than for the M waves (P < 0.05). Over the fatiguing protocol, the overall average decrease in MFCV (~25 %) was comparable to that of sEMG F mean (~22 %), but significantly greater than that of M-wave F mean (~9 %) (P < 0.001). The mean decline in MFCV was highly correlated with the mean decreases in both sEMG and M-wave F mean. CONCLUSIONS: The present findings indicated that, as fatigue progressed, central mechanisms could enhance the relative weight of the low-frequency components of the voluntary sEMG power spectrum, and/or the end-of-fiber (non-propagating) components could reduce the sensitivity of the M-wave spectrum to changes in conduction velocity.

EEG-based functional networks in schizophrenia.

Schizophrenia is often considered as a dysconnection syndrome in which, abnormal interactions between large-scale functional brain networks result in cognitive and perceptual deficits. In this article we apply the graph theoretic measures to brain functional networks based on the resting EEGs of fourteen schizophrenic patients in comparison with those of fourteen matched control subjects. The networks were extracted from common-average-referenced EEG time-series through partial and unpartial cross-correlation methods. Unpartial correlation detects functional connectivity based on direct and/or indirect links, while partial correlation allows one to ignore indirect links. We quantified the network properties with the graph metrics, including mall-worldness, vulnerability, modularity, assortativity, and synchronizability. The schizophrenic patients showed method-specific and frequency-specific changes especially pronounced for modularity, assortativity, and synchronizability measures. However, the differences between schizophrenia patients and normal controls in terms of graph theory metrics were stronger for the unpartial correlation method.

A two-dimensional doppler-radiometer for earth observation

Compared to synthetic aperture radars (SARs), the angular resolution of microwave radiometers is quite poor. Traditionally, it has been limited by the physical size of the antenna. However, the angular resolution can be improved by means of aperture synthesis interferometric techniques. A narrow beam is synthesized during the image formation processing of the cross-correlations measured at zero-lag between pairs of signals collected by an array of antennas. The angular resolution is then determined by the maximum antenna spacing normalized to the wavelength (baseline). The next step in improving the angular resolution is the Doppler-Radiometer, somehow related to the super-synthesis radiometers and the Radiometer-SAR. This paper presents the concept of a three-antenna Doppler-Radiometer for 2D imaging. The performance of this instrument is evaluated in terms of angular/spatial resolution and radiometric sensitivity, and an L-band illustrative example is presented.

Effect of radar rainfall time resolution on the predictive capability of a distributed hydrologic model

The performance of a hydrologic model depends on the rainfall input data, both spatially and temporally. As the spatial distribution of rainfall exerts a great influence on both runoff volumes and peak flows, the use of a distributed hydrologic model can improve the results in the case of convective rainfall in a basin where the storm area is smaller than the basin area. The aim of this study was to perform a sensitivity analysis of the rainfall time resolution on the results of a distributed hydrologic model in a flash-flood prone basin. Within such a catchment, floods are produced by heavy rainfall events with a large convective component. A second objective of the current paper is the proposal of a methodology that improves the radar rainfall estimation at a higher spatial and temporal resolution. Composite radar data from a network of three C-band radars with 6-min temporal and 2 × 2 km2 spatial resolution were used to feed the RIBS distributed hydrological model. A modification of the Window Probability Matching Method (gauge-adjustment method) was applied to four cases of heavy rainfall to improve the observed rainfall sub-estimation by computing new Z/R relationships for both convective and stratiform reflectivities. An advection correction technique based on the cross-correlation between two consecutive images was introduced to obtain several time resolutions from 1 min to 30 min. The RIBS hydrologic model was calibrated using a probabilistic approach based on a multiobjective methodology for each time resolution. A sensitivity analysis of rainfall time resolution was conducted to find the resolution that best represents the hydrological basin behaviour.

Magnetoestratigrafía de las sucesiones del Mioceno medio y superior del Vallés occidental (Depresión del Vallés-Penedés, NE España): Implicaciones biocronológicas y cronoestratigráficas.

The magnetostratigraphic analysis of the middle to late Miocene continental deposits from the Valles-Penedes basin, combined with its well-documented fossil mammal record, provides a well-resoluted chronology for the upper basin infill. It is based on the biostratigraphic and magnetostratigraphic cross-correlation of 18 sections throughout the alluvial and transitional/shallow marine sequences in the Western Valles area. The biostratigraphic framework consists of 24 mammal localities of upper Aragonian and Vallesian age. Correlation of the studied sections to the geomagnetic polarity time scale (GPTS) is based on the distinctive pattern of local magnetozones, as well as the radiometric age of the late Vallesian fauna from the Bicorp Basin (9.6 + 0.3 Ma) and the known relationship of the late Vallesian assemblages with marine beds belonging to the planktonic forarninifera N16 zone. It has led to an absolute dating of the fauna1 events and a precise chronostratigraphy of the Vallesian marnrnal stage in its type area. The Hipparion First Appearance Datum (FAD) defines the lower Vallesian boundary and is dated at 11.1 Ma, at the base of chron C5r. ln. It is in good agreement with radiometric ages from the early Hipparion bearing sites in the Vienna Basin (1 1.1 * 0.5 Ma) and the classic Howenegg locality in Germany (10.8 * 0.3 Ma). It also agrees with the age of the turkish localities of Yailacilar (1 1.6 + 0.25 Ma) and Yenieskihisar-2 (1 1.1 * 0.2 Ma) with absence of Hipparion. Al1 these support the isochrony of the dispersa1 of Hipparion throughout the Mediterranean region. A possible isochrony at a larger geographical scale (Old World, Mesogea) must await more reliable ages of the Hipparion FAD in Asia and Africa. The Cricetulodon FAD that defines the MN9a/MN9b boundary occurs at the middle part of C5n. Assuming an on average constant sedimentation rate, this datum has an age of approximately 10.4 Ma. The earlyllate Vallesian boundary is marked by one of the most distinct fauna1 events of the late Neogene: the dispersa1 of the muridae Progonomys into Europe and North Africa, which coincides with an important macromarnmal turnover. The first extensive appearance of Progonomys in Europe (MN9íMN10 boundary) is dated at 9.7 Ma (C4Ar3r), showing a remarkable diachrony with the Himalayan region. F9i d lly, the FAD of Rotundomys bressnnus occurs in the upper part of C4Ar.ln (9.2-9.3 Ma). The Vallesian spans 2.4 Myr, from 11.1 Ma (CSr.ln) to 8.7 Ma (C4An), and correlates to the early Tortonian.

Radio continuum and near-infrared study of the MGRO J2019+37 region

MGRO J2019+37 is an unidentified extended source of very high energy gamma-rays originally reported by the Milagro Collaboration as the brightest TeV source in the Cygnus region. Its extended emission could be powered by either a single or several sources. The GeV pulsar AGL J2020.5+3653 , discovered by AGILE and associated with PSR J2021+3651 , could contribute to the emission from MGRO J2019+37 . Aims. Our aim is to identify radio and near-infrared sources in the field of the extended TeV source MGRO J2019+37 , and study potential counterparts to explain its emission. Methods. We surveyed a region of about 6 square degrees with the Giant Metrewave Radio Telescope (GMRT) at the frequency 610 MHz. We also observed the central square degree of this survey in the near-infrared -band using the 3.5 m telescope in Calar Alto. Archival X-ray observations of some specific fields are included. VLBI observations of an interesting radio source were performed. We explored possible scenarios to produce the multi-TeV emission from MGRO J2019+37 and studied which of the sources could be the main particle accelerator. Results. We present a catalogue of 362 radio sources detected with the GMRT in the field of MGRO J2019+37 , and the results of a cross-correlation of this catalog with one obtained at near-infrared wavelengths, which contains ~3105 sources, as well as with available X-ray observations of the region. Some peculiar sources inside the ~1° uncertainty region of the TeV emission from MGRO J2019+37 are discussed in detail, including the pulsar PSR J2021+3651 and its pulsar wind nebula PWN G75.2+0.1 , two new radio-jet sources, the H II region Sh 2-104 containing two star clusters, and the radio source NVSS J202032+363158 . We also find that the hadronic scenario is the most likely in case of a single accelerator, and discuss the possible contribution from the sources mentioned above. Conclusions. Although the radio and GeV pulsar PSR J2021+3651 / AGL J2020.5+3653 and its associated pulsar wind nebula PWN G75.2+0.1 can contribute to the emission from MGRO J2019+37 , extrapolation of the GeV spectrum does not explain the detected multi-TeV flux. Other sources discussed here could contribute to the emission of the Milagro source.

Radio continuum and near-ingrared study of the MGRO J2019+37 region

Context. MGRO J2019+37 is an unidentified extended source of very high energy gamma-rays originally reported by the Milagro Collaboration as the brightest TeV source in the Cygnus region. Its extended emission could be powered by either a single or several sources. The GeV pulsar AGL J2020.5+3653, discovered by AGILE and associated with PSR J2021+3651, could contribute to the emission from MGRO J2019+37. Our aim is to identify radio and near-infrared sources in the field of the extended TeV source MGRO J2019+37, and study potential counterparts to explain its emission. Methods: We surveyed a region of about 6 square degrees with the Giant Metrewave Radio Telescope (GMRT) at the frequency 610 MHz. We also observed the central square degree of this survey in the near-infrared Ks-band using the 3.5 m telescope in Calar Alto. Archival X-ray observations of some specific fields are included. VLBI observations of an interesting radio source were performed. We explored possible scenarios to produce the multi-TeV emission from MGRO J2019+37 and studied which of the sources could be the main particle accelerator. Results: We present a catalogue of 362 radio sources detected with the GMRT in the field of MGRO J2019+37, and the results of a cross-correlation of this catalog with one obtained at near-infrared wavelengths, which contains ∼3 × 105 sources, as well as with available X-ray observations of the region. Some peculiar sources inside the ∼1◦ uncertainty region of the TeV emission from MGRO J2019+37 are discussed in detail, including the pulsar PSR J2021+3651 and its pulsar wind nebula PWN G75.2+0.1, two new radio-jet sources, the Hii region Sh 2-104 containing two star clusters, and the radio source NVSS J202032+363158. We also find that the hadronic scenario is the most likely in case of a single accelerator, and discuss the possible contribution from the sources mentioned above. Conclusions: Although the radio and GeV pulsar PSR J2021+3651 / AGL J2020.5+3653 and its associated pulsar wind nebula PWN G75.2+0.1 can contribute to the emission from MGRO J2019+37, extrapolation of the GeV spectrum does not explain the detected multi-TeV flux. Other sources discussed here could contribute to the emission of the Milagro source

Local noise analysis of a Schottky contact: combined thermionic-emissiondiffusion theory

A theoretical model for the noise properties of Schottky barrier diodes in the framework of the thermionic-emission¿diffusion theory is presented. The theory incorporates both the noise inducedby the diffusion of carriers through the semiconductor and the noise induced by the thermionicemission of carriers across the metal¿semiconductor interface. Closed analytical formulas arederived for the junction resistance, series resistance, and contributions to the net noise localized indifferent space regions of the diode, all valid in the whole range of applied biases. An additionalcontribution to the voltage-noise spectral density is identified, whose origin may be traced back tothe cross correlation between the voltage-noise sources associated with the junction resistance andthose for the series resistance. It is argued that an inclusion of the cross-correlation term as a newelement in the existing equivalent circuit models of Schottky diodes could explain the discrepanciesbetween these models and experimental measurements or Monte Carlo simulations.

Local noise analysis of a Schottky contact: combined thermionic-emissiondiffusion theory

A theoretical model for the noise properties of Schottky barrier diodes in the framework of the thermionic-emission¿diffusion theory is presented. The theory incorporates both the noise inducedby the diffusion of carriers through the semiconductor and the noise induced by the thermionicemission of carriers across the metal¿semiconductor interface. Closed analytical formulas arederived for the junction resistance, series resistance, and contributions to the net noise localized indifferent space regions of the diode, all valid in the whole range of applied biases. An additionalcontribution to the voltage-noise spectral density is identified, whose origin may be traced back tothe cross correlation between the voltage-noise sources associated with the junction resistance andthose for the series resistance. It is argued that an inclusion of the cross-correlation term as a newelement in the existing equivalent circuit models of Schottky diodes could explain the discrepanciesbetween these models and experimental measurements or Monte Carlo simulations.

Kinetic modeling of mechanisms of industrially important organic reactions in gas and liquid phase

This dissertation is based on 5 articles which deal with reaction mechanisms of the following selected industrially important organic reactions: 1. dehydrocyclization of n-butylbenzene to produce naphthalene 2. dehydrocyclization of 1-(p-tolyl)-2-methylbutane (MB) to produce 2,6-dimethylnaphthalene 3. esterification of neopentyl glycol (NPG) with different carboxylic acids to produce monoesters 4. skeletal isomerization of 1-pentene to produce 2-methyl-1-butene and 2-methyl-2-butene The results of initial- and integral-rate experiments of n-butylbenzene dehydrocyclization over selfmade chromia/alumina catalyst were applied when investigating reaction 2. Reaction 2 was performed using commercial chromia/alumina of different acidity, platina on silica and vanadium/calcium/alumina as catalysts. On all catalysts used for the dehydrocyclization, major reactions were fragmentation of MB and 1-(p-tolyl)-2-methylbutenes (MBes), dehydrogenation of MB, double bond transfer, hydrogenation and 1,6-cyclization of MBes. Minor reactions were 1,5-cyclization of MBes and methyl group fragmentation of 1,6- cyclization products. Esterification reactions of NPG were performed using three different carboxylic acids: propionic, isobutyric and 2-ethylhexanoic acid. Commercial heterogeneous gellular (Dowex 50WX2), macroreticular (Amberlyst 15) type resins and homogeneous para-toluene sulfonic acid were used as catalysts. At first NPG reacted with carboxylic acids to form corresponding monoester and water. Then monoester esterified with carboxylic acid to form corresponding diester. In disproportionation reaction two monoester molecules formed NPG and corresponding diester. All these three reactions can attain equilibrium. Concerning esterification, water was removed from the reactor in order to prevent backward reaction. Skeletal isomerization experiments of 1-pentene were performed over HZSM-22 catalyst. Isomerization reactions of three different kind were detected: double bond, cis-trans and skeletal isomerization. Minor side reaction were dimerization and fragmentation. Monomolecular and bimolecular reaction mechanisms for skeletal isomerization explained experimental results almost equally well. Pseudohomogeneous kinetic parameters of reactions 1 and 2 were estimated by usual least squares fitting. Concerning reactions 3 and 4 kinetic parameters were estimated by the leastsquares method, but also the possible cross-correlation and identifiability of parameters were determined using Markov chain Monte Carlo (MCMC) method. Finally using MCMC method, the estimation of model parameters and predictions were performed according to the Bayesian paradigm. According to the fitting results suggested reaction mechanisms explained experimental results rather well. When the possible cross-correlation and identifiability of parameters (Reactions 3 and 4) were determined using MCMC method, the parameters identified well, and no pathological cross-correlation could be seen between any parameter pair.

Management zones definition using soil chemical and physical attributes in a soybean area

Several equipments and methodologies have been developed to make available precision agriculture, especially considering the high cost of its implantation and sampling. An interesting possibility is to define management zones aim at dividing producing areas in smaller management zones that could be treated differently, serving as a source of recommendation and analysis. Thus, this trial used physical and chemical properties of soil and yield aiming at the generation of management zones in order to identify whether they can be used as recommendation and analysis. Management zones were generated by the Fuzzy C-Means algorithm and their evaluation was performed by calculating the reduction of variance and performing means tests. The division of the area into two management zones was considered appropriate for the present distinct averages of most soil properties and yield. The used methodology allowed the generation of management zones that can serve as source of recommendation and soil analysis; despite the relative efficiency has shown a reduced variance for all attributes in divisions in the three sub-regions, the ANOVA did not show significative differences among the management zones.

Asymmetric Smiles, Leverage Effects and Structural Parameters.

In this paper, we characterize the asymmetries of the smile through multiple leverage effects in a stochastic dynamic asset pricing framework. The dependence between price movements and future volatility is introduced through a set of latent state variables. These latent variables can capture not only the volatility risk and the interest rate risk which potentially affect option prices, but also any kind of correlation risk and jump risk. The standard financial leverage effect is produced by a cross-correlation effect between the state variables which enter into the stochastic volatility process of the stock price and the stock price process itself. However, we provide a more general framework where asymmetric implied volatility curves result from any source of instantaneous correlation between the state variables and either the return on the stock or the stochastic discount factor. In order to draw the shapes of the implied volatility curves generated by a model with latent variables, we specify an equilibrium-based stochastic discount factor with time non-separable preferences. When we calibrate this model to empirically reasonable values of the parameters, we are able to reproduce the various types of implied volatility curves inferred from option market data.