Metabifurcation analysis of a mean field model of the cortex

Mean field models (MFMs) of cortical tissue incorporate salient, average features of neural masses in order to model activity at the population level, thereby linking microscopic physiology to macroscopic observations, e.g., with the electroencephalogram (EEG). One of the common aspects of MFM descriptions is the presence of a high-dimensional parameter space capturing neurobiological attributes deemed relevant to the brain dynamics of interest. We study the physiological parameter space of a MFM of electrocortical activity and discover robust correlations between physiological attributes of the model cortex and its dynamical features. These correlations are revealed by the study of bifurcation plots, which show that the model responses to changes in inhibition belong to two archetypal categories or “families”. After investigating and characterizing them in depth, we discuss their essential differences in terms of four important aspects: power responses with respect to the modeled action of anesthetics, reaction to exogenous stimuli such as thalamic input, and distributions of model parameters and oscillatory repertoires when inhibition is enhanced. Furthermore, while the complexity of sustained periodic orbits differs significantly between families, we are able to show how metamorphoses between the families can be brought about by exogenous stimuli. We here unveil links between measurable physiological attributes of the brain and dynamical patterns that are not accessible by linear methods. They instead emerge when the nonlinear structure of parameter space is partitioned according to bifurcation responses. We call this general method “metabifurcation analysis”. The partitioning cannot be achieved by the investigation of only a small number of parameter sets and is instead the result of an automated bifurcation analysis of a representative sample of 73,454 physiologically admissible parameter sets. Our approach generalizes straightforwardly and is well suited to probing the dynamics of other models with large and complex parameter spaces.

Scenario and modelling uncertainty in global mean temperature change derived from emission-driven global climate models

We compare future changes in global mean temperature in response to different future scenarios which, for the first time, arise from emission-driven rather than concentration-driven perturbed parameter ensemble of a global climate model (GCM). These new GCM simulations sample uncertainties in atmospheric feedbacks, land carbon cycle, ocean physics and aerosol sulphur cycle processes. We find broader ranges of projected temperature responses arising when considering emission rather than concentration-driven simulations (with 10–90th percentile ranges of 1.7 K for the aggressive mitigation scenario, up to 3.9 K for the high-end, business as usual scenario). A small minority of simulations resulting from combinations of strong atmospheric feedbacks and carbon cycle responses show temperature increases in excess of 9 K (RCP8.5) and even under aggressive mitigation (RCP2.6) temperatures in excess of 4 K. While the simulations point to much larger temperature ranges for emission-driven experiments, they do not change existing expectations (based on previous concentration-driven experiments) on the timescales over which different sources of uncertainty are important. The new simulations sample a range of future atmospheric concentrations for each emission scenario. Both in the case of SRES A1B and the Representative Concentration Pathways (RCPs), the concentration scenarios used to drive GCM ensembles, lies towards the lower end of our simulated distribution. This design decision (a legacy of previous assessments) is likely to lead concentration-driven experiments to under-sample strong feedback responses in future projections. Our ensemble of emission-driven simulations span the global temperature response of the CMIP5 emission-driven simulations, except at the low end. Combinations of low climate sensitivity and low carbon cycle feedbacks lead to a number of CMIP5 responses to lie below our ensemble range. The ensemble simulates a number of high-end responses which lie above the CMIP5 carbon cycle range. These high-end simulations can be linked to sampling a number of stronger carbon cycle feedbacks and to sampling climate sensitivities above 4.5 K. This latter aspect highlights the priority in identifying real-world climate-sensitivity constraints which, if achieved, would lead to reductions on the upper bound of projected global mean temperature change. The ensembles of simulations presented here provides a framework to explore relationships between present-day observables and future changes, while the large spread of future-projected changes highlights the ongoing need for such work.

The numbers of sample design for evaluation of the soil fertility

The determination of the amount of sample units that will compose the sample express the optimization of the workforce, and reduce errors inherent in the report of recommendation and evaluation of soil fertility. This study aimed to determine in three systems use and soil management, the numbers of units samples design, needed to form the composed sample, for evaluation of soil fertility. It was concluded that the number of sample units needed to compose the composed sample to determination the attributes of organic matter, pH, P, K, Ca, Mg, Al and H+Al and base saturation of soil vary by use and soil management and error acceptable to the mean estimate. For the same depth of collected, increasing the number of sample units, reduced the percentage error in estimating the average, allowing the recommendation of 14, 14 and 11 sample in management with native vegetation, pasture cultivation and corn, respectively, for a error 20% on the mean estimate.

A finite sample improvement of the fixed effects estimator : applied to technical infficiency

The FE ('fixed effects') estimator of technical inefficiency performs poorly when N ('number of firms') is large and T ('number of time observations') is small. We propose estimators of both the firm effects and the inefficiencies, which have small sample gains compared to the traditional FE estimator. The estimators are based on nonparametric kernel regression of unordered variables, which includes the FE estimator as a special case. In terms of global conditional MSE ('mean square error') criterions, it is proved that there are kernel estimators which are efficient to the FE estimators of firm effects and inefficiencies, in finite samples. Monte Carlo simulations supports our theoretical findings and in an empirical example it is shown how the traditional FE estimator and the proposed kernel FE estimator lead to very different conclusions about inefficiency of Indonesian rice farmers.

A new chart based on sample variances for monitoring the covariance matrix of multivariate processes

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Control charts for monitoring the mean vector and the covariance matrix of bivariate processes

In this article, we propose new control charts for monitoring the mean vector and the covariance matrix of bivariate processes. The traditional tools used for this purpose are the T (2) and the |S| charts. However, these charts have two drawbacks: (1) the T (2) and the |S| statistics are not easy to compute, and (2) after a signal, they do not distinguish the variable affected by the assignable cause. As an alternative to (1), we propose the MVMAX chart, which only requires the computation of sample means and sample variances. As an alternative to (2), we propose the joint use of two charts based on the non-central chi-square statistic (NCS statistic), named as the NCS charts. Once the NCS charts signal, the user can immediately identify the out-of-control variable. In general, the synthetic MVMAX chart is faster than the NCS charts and the joint T (2) and |S| charts in signaling processes disturbances.

Monitoring a wandering mean with an np chart

Este artigo considera um gráfico np x proposto por Wu et al. (2009) para controle de média de processo como uma alternativa ao uso do gráfico de. O que distingue do gráfico de controle np x é o fato das unidades amostrais serem classificadas como unidades de primeiro ou de segunda classe de acordo com seus limites discriminantes. O gráfico tradicional np é um caso particular do gráfico np x quando os limites discriminantes coincidem com os limites de especificação e unidade de primeira (segunda) classe é um item conforme (não conforme). Estendendo o trabalho de Reynolds Junior, Arnold e Baik (1996), consideramos que a média de processo oscila mesmo na ausência de alguma causa especial. As propriedades de Cadeia de Markov foram adotadas para avaliar o desempenho do gráfico np x no monitoramento de média de processos que oscila. de modo geral, o gráfico np x requer amostras duas vezes maior para superar desempenho do gráfico (enquanto que o gráfico tradicional np necessita tamanho de amostras cinco ou seis vezes maior).

Joint (X)over-bar and R charts with variable sample sizes and sampling intervals

Recent studies have shown that the (X) over bar chart with variable sampling intervals (VSI) and/or with variable sample sizes (VSS) detects process shifts faster than the traditional (X) over bar chart. This article extends these studies for processes that are monitored by both the (X) over bar and R charts. A Markov chain model is used to determine the properties of the joint (X) over bar and R charts with variable sample sizes and sampling intervals (VSSI). The VSSI scheme improves the joint (X) over bar and R control chart performance in terms of the speed with which shifts in the process mean and/or variance are detected.

Treatment effects produced by the Bionator appliance. Comparison with an untreated Class II sample

The purpose of this retrospective investigation was to evaluate the dentoalveolar and skeletal cephalometric changes of the Bionator appliance on individuals with a Class II division 1 malocclusion. Lateral cephalograms of 44 patients were divided into two equal groups. The control group comprised 22 untreated Class II children (11 males, 11 females), with an initial mean age of 8 years 7 months who were followed without treatment for a period of 13 months. The Bionator group (111 males, 11 females) had an initial mean age of 10 years 8 months, and were treated for a mean period of 16 months. Lateral cephalometric headfilms were obtained of each patient and control at the beginning and end of treatment.The results showed that there were no changes in forward growth of the maxilla in the experimental group compared with the control group. However, the Bionator treatment produced a statistically significant increase in mandibular protrusion, and in total mandibular and body lengths. There were no statistically significant differences in craniofacial growth direction between the Bionator group and the control group, although the treated patients demonstrated a greater increase in posterior face height. The Bionator appliance produced labial tipping of the lower incisors and lingual inclination of the upper incisors, as well as a significant increase (P < 0.01) in mandibular posterior dentoalveolar height. The major effects of the Bionator appliance were dentoalveolar, with a smaller significant skeletal effect. The results indicate that the correction of a Class II division 1 malocclusion with the Bionator appliance is achieved not only by a combination of mandibular skeletal effects, but also by significant dentoalveolar changes.

Synthetic control chart for monitoring the pprocess mean and variance

In this article, we consider the synthetic control chart with two-stage sampling (SyTS chart) to control the process mean and variance. During the first stage, one item of the sample is inspected; if its value X, is close to the target value of the process mean, then the sampling is interrupted. Otherwise, the sampling goes on to the second stage, where the remaining items are inspected and the statistic T = Sigma [x(i) - mu(0) + xi sigma(0)](2) is computed taking into account all items of the sample. The design parameter is function of X-1. When the statistic T is larger than a specified value, the sample is classified as nonconforming. According to the synthetic procedure, the signal is based on Conforming Run Length (CRL). The CRL is the number of samples taken from the process since the previous nonconforming sample until the occurrence of the next nonconforming sample. If the CRL is sufficiently small, then a signal is generated. A comparative study shows that the SyTS chart and the joint X and S charts with double sampling are very similar in performance. However, from the practical viewpoint, the SyTS chart is more convenient to administer than the joint charts.

X charts with supplementary samples to control the mean and variance

A standard X chart for controlling a process takes regular individual observations, for instance every half hour. This article proposes a modification of the X chart that allows one to take supplementary samples. The supplementary sample is taken (and the (X) over bar and R values computed) when the current value of X falls outside the control limits. With the supplementary sample, the signal of out-of-control is given by an (X) over bar value outside the (X) over bar chart's control limits or an R value outside the R chart's control limit. The proposed chart is designed to hold the supplementary sample frequency, during the in-control period, as low as 5% or less. In this context, the practitioner might prefer to verify an out-of-control condition by simply comparing the (X) over bar and R values with the control limits. In other words, without plotting the (X) over bar and R points. The X chart with supplementary samples has two major advantages when compared with the standard (X) over bar and A charts: (a) the user will be plotting X values instead of (X) over bar and R values; (b) the shifts in the process mean and/or changes in the process variance are detected faster.

Monitoring the process mean and variance using a synthetic control chart with two-stage testing

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

The Variable Sample Size (X)over-bar Chart with Estimated Parameters

The VSS X chart, dedicated to the detection of small to moderate mean shifts in the process, has been investigated by several researchers under the assumption of known process parameters. In practice, the process parameters are rarely known and are usually estimated from an in-control Phase I data set. In this paper, we evaluate the (run length) performances of the VSS chart when the process parameters are estimated, we compare them in the case where the process parameters are assumed known and we propose specific optimal control chart parameters taking the number of Phase I samples into account.

Joint X̄ and R Charts with Variable Sample Sizes and Sampling Intervals

Recent studies have shown that the X̄ chart with variable sampling intervals (VSI) and/or with variable sample sizes (VSS) detects process shifts faster than the traditional X̄ chart. This article extends these studies for processes that are monitored by both the X̄ and R charts. A Markov chain model is used to determine the properties of the joint X and R charts with variable sample sizes and sampling intervals (VSSI). The VSSI scheme improves the joint X̄ and R control chart performance in terms of the speed with which shifts in the process mean and/or variance are detected.

A single chart with supplementary runs rules for monitoring the mean vector and the covariance matrix of multivariate processes

The MRMAX chart is a single chart based on the standardized sample means and sample ranges for monitoring the mean vector and the covariance matrix of multivariate processes. User's familiarity with the computation of these statistics is a point in favor of the MRMAX chart. As a single chart, the recently proposed MRMAX chart is very appropriate for supplementary runs rules. In this article, we compare the supplemented MRMAX chart and the synthetic MRMAX chart with the standard MRMAX chart. The supplementary and the synthetic runs rules enhance the performance of the MRMAX chart. © 2013 Elsevier Ltd.