On the relationship between connections and the asymptotic properties of predictive distributions

In a recent paper, Komaki studied the second-order asymptotic properties of predictive distributions, using the Kullback-Leibler divergence as a loss function. He showed that estimative distributions with asymptotically efficient estimators can be improved by predictive distributions that do not belong to the model. The model is assumed to be a multidimensional curved exponential family. In this paper we generalize the result assuming as a loss function any f divergence. A relationship arises between alpha connections and optimal predictive distributions. In particular, using an alpha divergence to measure the goodness of a predictive distribution, the optimal shift of the estimate distribution is related to alpha-covariant derivatives. The expression that we obtain for the asymptotic risk is also useful to study the higher-order asymptotic properties of an estimator, in the mentioned class of loss functions.

An efficient algorithm to perform multiple testing in epistasis screening

Background: Research in epistasis or gene-gene interaction detection for human complex traits has grown over the last few years. It has been marked by promising methodological developments, improved translation efforts of statistical epistasis to biological epistasis and attempts to integrate different omics information sources into the epistasis screening to enhance power. The quest for gene-gene interactions poses severe multiple-testing problems. In this context, the maxT algorithm is one technique to control the false-positive rate. However, the memory needed by this algorithm rises linearly with the amount of hypothesis tests. Gene-gene interaction studies will require a memory proportional to the squared number of SNPs. A genome-wide epistasis search would therefore require terabytes of memory. Hence, cache problems are likely to occur, increasing the computation time. In this work we present a new version of maxT, requiring an amount of memory independent from the number of genetic effects to be investigated. This algorithm was implemented in C++ in our epistasis screening software MBMDR-3.0.3. We evaluate the new implementation in terms of memory efficiency and speed using simulated data. The software is illustrated on real-life data for Crohn’s disease. Results: In the case of a binary (affected/unaffected) trait, the parallel workflow of MBMDR-3.0.3 analyzes all gene-gene interactions with a dataset of 100,000 SNPs typed on 1000 individuals within 4 days and 9 hours, using 999 permutations of the trait to assess statistical significance, on a cluster composed of 10 blades, containing each four Quad-Core AMD Opteron(tm) Processor 2352 2.1 GHz. In the case of a continuous trait, a similar run takes 9 days. Our program found 14 SNP-SNP interactions with a multiple-testing corrected p-value of less than 0.05 on real-life Crohn’s disease (CD) data. Conclusions: Our software is the first implementation of the MB-MDR methodology able to solve large-scale SNP-SNP interactions problems within a few days, without using much memory, while adequately controlling the type I error rates. A new implementation to reach genome-wide epistasis screening is under construction. In the context of Crohn’s disease, MBMDR-3.0.3 could identify epistasis involving regions that are well known in the field and could be explained from a biological point of view. This demonstrates the power of our software to find relevant phenotype-genotype higher-order associations.

The novel proteins Rng8 and Rng9 regulate the myosin-V Myo51 during fission yeast cytokinesis.

The myosin-V family of molecular motors is known to be under sophisticated regulation, but our knowledge of the roles and regulation of myosin-Vs in cytokinesis is limited. Here, we report that the myosin-V Myo51 affects contractile ring assembly and stability during fission yeast cytokinesis, and is regulated by two novel coiled-coil proteins, Rng8 and Rng9. Both rng8Δ and rng9Δ cells display similar defects as myo51Δ in cytokinesis. Rng8 and Rng9 are required for Myo51's localizations to cytoplasmic puncta, actin cables, and the contractile ring. Myo51 puncta contain multiple Myo51 molecules and walk continuously on actin filaments in rng8(+) cells, whereas Myo51 forms speckles containing only one dimer and does not move efficiently on actin tracks in rng8Δ. Consistently, Myo51 transports artificial cargos efficiently in vivo, and this activity is regulated by Rng8. Purified Rng8 and Rng9 form stable higher-order complexes. Collectively, we propose that Rng8 and Rng9 form oligomers and cluster multiple Myo51 dimers to regulate Myo51 localization and functions.

Why should we care about soil fauna?

The reasons why we care about soil fauna are related to their intrinsic, utilitarian and functional values. The intrinsic values embrace aesthetic or moral reasons for conserving below-ground biodiversity. Unfortunately, the protection of soil invertebrates has rarely been a criterion for avoiding changes in land use and management. Utilitarian, or direct use values, have been investigated more extensively for fungi, bacteria and marine invertebrates than for soil fauna. However, some traditional remedies, novel enzymes and pharmaceutical compounds have been derived from earthworms, termites and other groups, and gut symbionts may provide microbial strains with interesting properties for biotechnology. The functional importance of soil invertebrates in ecosystem processes has been a major focus of research in recent decades. It is suggested herein that it is rarely possible to identify the role of soil invertebrates as rate determinants of soil processes at plot and ecosystem scales of hectares and above because other biophysical controls override their effects. There are situations, however, where the activities of functional groups of soil animals, even of species, are synchronised in space or time by plant events, resource inputs, seasonality or other perturbations to the system, and their emergent effects are detectable as higher order controls.

Simulation of surface waves in porous media

We present a novel numerical algorithm for the simulation of seismic wave propagation in porous media, which is particularly suitable for the accurate modelling of surface wave-type phenomena. The differential equations of motion are based on Biot's theory of poro-elasticity and solved with a pseudospectral approach using Fourier and Chebyshev methods to compute the spatial derivatives along the horizontal and vertical directions, respectively. The time solver is a splitting algorithm that accounts for the stiffness of the differential equations. Due to the Chebyshev operator the grid spacing in the vertical direction is non-uniform and characterized by a denser spatial sampling in the vicinity of interfaces, which allows for a numerically stable and accurate evaluation of higher order surface wave modes. We stretch the grid in the vertical direction to increase the minimum grid spacing and reduce the computational cost. The free-surface boundary conditions are implemented with a characteristics approach, where the characteristic variables are evaluated at zero viscosity. The same procedure is used to model seismic wave propagation at the interface between a fluid and porous medium. In this case, each medium is represented by a different grid and the two grids are combined through a domain-decomposition method. This wavefield decomposition method accounts for the discontinuity of variables and is crucial for an accurate interface treatment. We simulate seismic wave propagation with open-pore and sealed-pore boundary conditions and verify the validity and accuracy of the algorithm by comparing the numerical simulations to analytical solutions based on zero viscosity obtained with the Cagniard-de Hoop method. Finally, we illustrate the suitability of our algorithm for more complex models of porous media involving viscous pore fluids and strongly heterogeneous distributions of the elastic and hydraulic material properties.

Impulse response recovery of linear systems through weighted cumulant slice

Identifiability of the so-called ω-slice algorithm is proven for ARMA linear systems. Although proofs were developed in the past for the simpler cases of MA and AR models, they were not extendible to general exponential linear systems. The results presented in this paper demonstrate a unique feature of the ω-slice method, which is unbiasedness and consistency when order is overdetermined, regardless of the IIR or FIR nature of the underlying system, and numerical robustness.

Designing online interaction to address disciplinary competencies: A cross-country comparison of faculty perspectives

This study was conducted at colleges in three countries (United States, Venezuela, and Spain) and across three academic disciplines (engineering, education, and business), to examine how experienced faculty define competencies for their discipline, and design instructional interaction for online courses. A qualitative research design employing in-depth interviews was selected. Results show that disciplinary knowledge takes precedence when faculty members select competencies to be developed in online courses for their respective professions. In all three disciplines, the design of interaction to correspond with disciplinary competencies was often influenced by contextual factors that modify faculty intention. Therefore, instructional design will vary across countries in the same discipline to address the local context, such as the needs and expectations of the learners, faculty perspectives, beliefs and values, and the needs of the institution, the community, and country. The three disciplines from the three countries agreed on the importance of the following competencies: knowledge of the field, higher order cognitive processes such as critical thinking, analysis, problem solving, transfer of knowledge, oral and written communication skills, team work, decision making, leadership and management skills, indicating far more similarities in competencies than differences between the three different applied disciplines. We found a lack of correspondence between faculty¿s intent to develop collaborative learning skills and the actual development of them. Contextual factors such as faculty prior experience in design, student reluctance to engage in collaborative learning, and institutional assessment systems that focus on individual performance were some of these reasons.

Different concentrations of Mg++ ions affect nuclear matrix protein distribution during thermal stabilization of isolated nuclei.

The nuclear matrix, a proteinaceous network believed to be a scaffolding structure determining higher-order organization of chromatin, is usually prepared from intact nuclei by a series of extraction steps. In most cell types investigated the nuclear matrix does not spontaneously resist these treatments but must be stabilized before the application of extracting agents. Incubation of isolated nuclei at 37C or 42C in buffers containing Mg++ has been widely employed as stabilizing agent. We have previously demonstrated that heat treatment induces changes in the distribution of three nuclear scaffold proteins in nuclei prepared in the absence of Mg++ ions. We studied whether different concentrations of Mg++ (2.0-5 mM) affect the spatial distribution of nuclear matrix proteins in nuclei isolated from K562 erythroleukemia cells and stabilized by heat at either 37C or 42C. Five proteins were studied, two of which were RNA metabolism-related proteins (a 105-kD component of splicing complexes and an RNP component), one a 126-kD constituent of a class of nuclear bodies, and two were components of the inner matrix network. The localization of proteins was determined by immunofluorescent staining and confocal scanning laser microscope. Mg++ induced significant changes of antigen distribution even at the lowest concentration employed, and these modifications were enhanced in parallel with increase in the concentration of the divalent cation. The different sensitivity to heat stabilization and Mg++ of these nuclear proteins might reflect a different degree of association with the nuclear scaffold and can be closely related to their functional or structural role.

Phenomenological models of holographic superconductors and hall currents

We study general models of holographic superconductivity parametrized by four arbitrary functions of a neutral scalar field of the bulk theory. The models can accommodate several features of real superconductors, like arbitrary critical temperatures and critical exponents in a certain range, and perhaps impurities or boundary or thickness effects. We find analytical expressions for the critical exponents of the general model and show that they satisfy the Rushbrooke identity. An important subclass of models exhibit second order phase transitions. A study of the specific heat shows that general models can also describe holographic superconductors undergoing first, second and third (or higher) order phase transitions. We discuss how small deformations of the HHH model can lead to the appearance of resonance peaks in the conductivity, which increase in number and become narrower as the temperature is gradually decreased, without the need for tuning mass of the scalar to be close to the Breitenlohner-Freedman bound. Finally, we investigate the inclusion of a generalized ¿theta term¿ producing Hall effect without magnetic field.

The role of intramolecular barriers on the glass transition of polymers: computer simulations versus mode coupling theory

We present computer simulations of a simple bead-spring model for polymer melts with intramolecular barriers. By systematically tuning the strength of the barriers, we investigate their role on the glass transition. Dynamic observables are analyzed within the framework of the mode coupling theory (MCT). Critical nonergodicity parameters, critical temperatures, and dynamic exponents are obtained from consistent fits of simulation data to MCT asymptotic laws. The so-obtained MCT λ-exponent increases from standard values for fully flexible chains to values close to the upper limit for stiff chains. In analogy with systems exhibiting higher-order MCT transitions, we suggest that the observed large λ-values arise form the interplay between two distinct mechanisms for dynamic arrest: general packing effects and polymer-specific intramolecular barriers. We compare simulation results with numerical solutions of the MCT equations for polymer systems, within the polymer reference interaction site model (PRISM) for static correlations. We verify that the approximations introduced by the PRISM are fulfilled by simulations, with the same quality for all the range of investigated barrier strength. The numerical solutions reproduce the qualitative trends of simulations for the dependence of the nonergodicity parameters and critical temperatures on the barrier strength. In particular, the increase in the barrier strength at fixed density increases the localization length and the critical temperature. However the qualitative agreement between theory and simulation breaks in the limit of stiff chains. We discuss the possible origin of this feature.

Reaction time asymmetries between expansion and contraction

Different asymmetries between expansion and contraction (radial motions) have been reported in the literature. Often these patterns have been regarded as implying different channels for each type of radial direction (outward versus inwards) operating at a higher level of visual motion processing. In two experiments (detection and discrimination tasks) we report reaction time asymmetries between expansion and contraction. Power functions were fitted to the data. While an exponent of 0.5 accounted for the expansion data better, a value of unity yielded the best fit for the contraction data. Instead of interpreting these differences as corresponding to different higher order motion detectors, we regard these findings as reflecting the fact that expansion and contraction tap two distinct psychophysical input channels underlying the processing of fast and slow velocities respectively.

Intrahemispheric cortico-cortical connections of the human auditory cortex.

The human auditory cortex comprises the supratemporal plane and large parts of the temporal and parietal convexities. We have investigated the relevant intrahemispheric cortico-cortical connections using in vivo DSI tractography combined with landmark-based registration, automatic cortical parcellation and whole-brain structural connection matrices in 20 right-handed male subjects. On the supratemporal plane, the pattern of connectivity was related to the architectonically defined early-stage auditory areas. It revealed a three-tier architecture characterized by a cascade of connections from the primary auditory cortex to six adjacent non-primary areas and from there to the superior temporal gyrus. Graph theory-driven analysis confirmed the cascade-like connectivity pattern and demonstrated a strong degree of segregation and hierarchy within early-stage auditory areas. Putative higher-order areas on the temporal and parietal convexities had more widely spread local connectivity and long-range connections with the prefrontal cortex; analysis of optimal community structure revealed five distinct modules in each hemisphere. The pattern of temporo-parieto-frontal connectivity was partially asymmetrical. In conclusion, the human early-stage auditory cortical connectivity, as revealed by in vivo DSI tractography, has strong similarities with that of non-human primates. The modular architecture and hemispheric asymmetry in higher-order regions is compatible with segregated processing streams and lateralization of cognitive functions.

A thematic analysis of delusion with religious content in schizophrenia: open, closed, and mixed dynamics

The aim of the present study was to elicit how patients with delusions with religious contents conceptualized or experienced their spirituality and religiousness. Sixty-two patients with present or past religious delusions went through semistructured interviews, which were analyzed using the three coding steps described in the grounded theory. Three major themes were found in religious delusions: ''spiritual identity,'' ''meaning of illness,'' and ''spiritual figures.'' One higher-order concept was found: ''structure of beliefs.'' We identified dynamics that put these personal beliefs into a constant reconstruction through interaction with the world and others (i.e., open dynamics) and conversely structural dynamics that created a complete rupture with the surrounding world and others (i.e., closed structural dynamics); those dynamics may coexist. These analyses may help to identify psychological functions of delusions with religious content and, therefore, to better conceptualize interventions when dealing with it in psychotherapy.

25 anos de quimiometria no Brasil

Chemometric activities in Brazil are described according to three phases: before the existence of microcomputers in the 1970s, through the initial stages of microcomputer use in the 1980s and during the years of extensive microcomputer applications of the ´90s and into this century. Pioneering activities in both the university and industry are emphasized. Active research areas in chemometrics are cited including experimental design, pattern recognition and classification, curve resolution for complex systems and multivariate calibration. New trends in chemometrics, especially higher order methods for treating data, are emphasized.

Indirect likelihood inference (revised)

Standard indirect Inference (II) estimators take a given finite-dimensional statistic, Z_{n} , and then estimate the parameters by matching the sample statistic with the model-implied population moment. We here propose a novel estimation method that utilizes all available information contained in the distribution of Z_{n} , not just its first moment. This is done by computing the likelihood of Z_{n}, and then estimating the parameters by either maximizing the likelihood or computing the posterior mean for a given prior of the parameters. These are referred to as the maximum indirect likelihood (MIL) and Bayesian Indirect Likelihood (BIL) estimators, respectively. We show that the IL estimators are first-order equivalent to the corresponding moment-based II estimator that employs the optimal weighting matrix. However, due to higher-order features of Z_{n} , the IL estimators are higher order efficient relative to the standard II estimator. The likelihood of Z_{n} will in general be unknown and so simulated versions of IL estimators are developed. Monte Carlo results for a structural auction model and a DSGE model show that the proposed estimators indeed have attractive finite sample properties.