Functional analysis: evaluation of response intensities--tailoring ANOVA for lists of expression subsets.

Background: Microarray data is frequently used to characterize the expression profile of a whole genome and to compare the characteristics of that genome under several conditions. Geneset analysis methods have been described previously to analyze the expression values of several genes related by known biological criteria (metabolic pathway, pathology signature, co-regulation by a common factor, etc.) at the same time and the cost of these methods allows for the use of more values to help discover the underlying biological mechanisms. Results: As several methods assume different null hypotheses, we propose to reformulate the main question that biologists seek to answer. To determine which genesets are associated with expression values that differ between two experiments, we focused on three ad hoc criteria: expression levels, the direction of individual gene expression changes (up or down regulation), and correlations between genes. We introduce the FAERI methodology, tailored from a two-way ANOVA to examine these criteria. The significance of the results was evaluated according to the self-contained null hypothesis, using label sampling or by inferring the null distribution from normally distributed random data. Evaluations performed on simulated data revealed that FAERI outperforms currently available methods for each type of set tested. We then applied the FAERI method to analyze three real-world datasets on hypoxia response. FAERI was able to detect more genesets than other methodologies, and the genesets selected were coherent with current knowledge of cellular response to hypoxia. Moreover, the genesets selected by FAERI were confirmed when the analysis was repeated on two additional related datasets. Conclusions: The expression values of genesets are associated with several biological effects. The underlying mathematical structure of the genesets allows for analysis of data from several genes at the same time. Focusing on expression levels, the direction of the expression changes, and correlations, we showed that two-step data reduction allowed us to significantly improve the performance of geneset analysis using a modified two-way ANOVA procedure, and to detect genesets that current methods fail to detect.

A methodology for the adaptation of a PMP at the determination of a PMF

This paper presents a new type of very fine grid hydrological model based on the spatiotemporal repartition of a PMP (Probable Maximum Precipitation) and on the topography. The goal is to estimate the influence of this rain on a PMF (Probable Maximum Flood) on a catchment area in Switzerland. The spatiotemporal distribution of the PMP was realized using six clouds modeled by the advection-diffusion equation. The equation shows the movement of the clouds over the terrain and also gives the evolution of the rain intensity in time. This hydrological modeling is followed by a hydraulic modeling of the surface and subterranean flow, done considering the factors that contribute to the hydrological cycle, such as the infiltration, the resurgence and the snowmelt. These added factors make the developed model closer to reality and also offer flexibility in the initial condition that is added to the factors concerning the PMP, such as the duration of the rain, the speed and direction of the wind. All these initial conditions taken together offer a complete image of the PMF.

The early Miocene rotation of the Corso-Sardinian block. New paleomagnetic constraints for the end of the motion

The paleomagnetic investigations carried out in the 70's on Oligo-Miocene volcanics of Sardinia have demonstrated that the island was turned by 35-30 degrees clockwise from 33 Ma up to 3-1-20.5 Ma and rotated counterclockwise in a few million years [De Jong et al., 1969, 1973; Bobier et Coulon, 1970; Coulon et al., 1974; Manzoni, 1974, 1975; Bellon rr nl.. 1977: Edel et Lortscher, 1977; Edel, 1979, 1980]. Since then, the end of the rotation fixed at 19 Ma by Montigny er al. [1981] was the subject of discussions and several studies associating paleomagnetism and radiometric dating were undertaken [Assorgia er al., 1994: Vigliotti et Langenheim, 1995: Deino et al., 1997; Gattacceca rt Deino, 1999]. This is a contribution to this debate that is hampered by thr important secular variation recorded in the volcanics. The only way to get our of this problem is to sample series of successive flows as completely as possible, and to reduce the effect of secular variation by the calculation of means. Sampling was performed north of Bonorva in 5 pyroclastic flows that belong to the upper ignimbritic series SI2 according to Coulon rr nl. [1974] or LBLS according to Assorgia et al. [1997] (fig. I). Ar-40/Ar-39 dating of biotites from the debris flow (MDF) has yielded an age or 18.35 +/- 0.03 Ma [Dubois, 2000]. Five of the investigated sites are located beneath the debris flow ITV, TVB, TVD, SPM85, SPM86), one site was cured in the matrix of the debris flow (MDF) and one in 4 metric blocks included in the flow (DFC). Another site was sampled in the upper ash flow (PDM) that marks the end of the pyroclastic activity, just before the marine transgression. According to micropaleontological and radiometric dating this transgression has occurred between 18.35 and 17.6 Ma [Dubois, 2000]. After removal of a soft viscous component, the thermal demagnetization generally shows a univectorial behaviour of the remanent magnetization (fig. 2a). The maximum unblocking temperatures of 580-620 degrees (tab. I) and a rapid saturation below 100 mT (fig. 3) indicate that the carrier of the characteristic magnetization is magnetite. The exception comes: from the upper site PDM in which were found two characteristic components, one with a normal polarity and low unblocking temperatures up to 350 degreesC and one with a reversed polarity and maximum unblocking temperatures at 580-600 degreesC of magnetite. After calculation of a mean direction for each flow, the mean << Al >> direction 4 degrees /57 degrees (alpha (95) = 13 degrees) computed with the mean directions for the 5 flows may be considered as weakly affected by secular variation. But the results require a more careful examination. The declinations are N to NNW beneath the debris flow. NNW in the debris flow. and NNE (or SSW) above the debris flow, The elongated distribution of the directions obtained at sites TVB and TVD. scattered from the mean direction of TV to the mean direction of MDF is interpreted as due to partial overprinting during the debris How volcanic episode, The low temperature component PDMa is likely related to the alteration seen on thin sections and is also viewed as an overprint. As NNE/SSW directions occur as well below (mean direction << B >> : 5 degrees /58 degrees) as above the debris flow (PDMb : 200 degrees/-58 degrees). the NNW directions (<< C >> : 337 degrees /64 degrees) associated with the debris flow volcanism may be interpreted as resulting from a magnetic field excursion. According to the polarity scale of Cande and Kent [1992, 1995] and the radiometric age of MDF, the directions with normal polarity (TV, TVB, TVD, SPM85. SPM86a. MDF. DFC) may represent the period 5En. while the directions with reversed polarity PDMb and SPM86b were likely acquired during the period 5Dr. Using the mean << Al >> direction, the mean << B >>, or the PDM direction (tab. I). the deviation in declination with the direction of stable Europe 6.4 degrees /58.7 degrees (alpha (95) = 8 degrees) for a selection of 4 middle Tertiary poles by Besse et Courtillot [1991] or 7 degrees /56 degrees (alpha (95) = 3 degrees) for 19 poles listed by Edel [1980] can be considered as negligible. Using the results from the uppermost ignimbritic layer of Anglona also emplaced around 18.3 Ma [Odin rt al.. 1994]. the mean direction << E >> (3 degrees /51.5 degrees) leads to the same conclusion. On the contrary, when taking into account all dated results available for the period 5En (mean direction << D >> 353 degrees /56 degrees for 45 sites) (tab. II). the deviation 13 degrees is much more significant. As the rotation of Sardinia started around 21-20.5 Ma. the assumption of a constant velocity of rotation and the deviations of the Sardinia directions with respect to the stable Europe direction locate the end of the motion between 18.3 and 17.2 or 16.7 Ma (fig. 4). During the interval 18.35-17.5 Ma, the marine transgression took place. At the same period a NE-SW shortening interpreted as resulting from the collision of Sardinia with Apulia affected different parts of the island [Letouzey et al., 1982]. Consequently, the new paleomagnetic results and the tectono-sedimentary evolution are in favour of an end of the rotation at 17.5-18 Ma.

The truncated mean of an asymmetric distribution

This paper investigates a simple procedure to estimate robustly the mean of an asymmetric distribution. The procedure removes the observations which are larger or smaller than certain limits and takes the arithmetic mean of the remaining observations, the limits being determined with the help of a parametric model, e.g., the Gamma, the Weibull or the Lognormal distribution. The breakdown point, the influence function, the (asymptotic) variance, and the contamination bias of this estimator are explored and compared numerically with those of competing estimates.

Asymmetry of Functional Connectivity in Schizophrenia at Different Spatial Scales

Introduction: The interhemispheric asymmetries that originate from connectivity-related structuring of the cerebral cortex are compromised in schizophrenia (SZ). Recently, we have revealed the whole-head topography of EEG synchronization in SZ (Jalili et al. 2007; Knyazeva et al. 2008). Here we extended the analysis to assess the abnormality in the asymmetry of synchronization, which is further motivated by the evidence that the interhemispheric asymmetries suspected to be abnormal in SZ originate from the connectivity-related structuring of the cortex. Methods: Thirteen right-handed SZ patients and thirteen matched controls, participated in this study and the multichannel (128) EEGs were recorded for 3-5 minutes at rest. Then, Laplacian EEG (LEEG) were calculated using a 2-D spline. The LEEGs were analysis through calculating the power spectral density using Welch's average periodogram method. Furthermore, using a state-space based multivariate synchronization measure, S-estimator, we analyzed the correlate of the functional cortico-cortical connectivity in SZ patients compared to the controls. The values of S-estimator were obtained at three different special scales: first-order neighbors for each sensor location, second-order neighbors, and the whole hemisphere. The synchronization measures based on LEEG of alpha and beta bands were applied and tuned to various spatial scales including local, intraregional, and long-distance levels. To assess the between-group differences, we used a permutation version of Hotelling's T2 test. For correlation analysis, Spearman Rank Correlation was calculated. Results: Compared to the controls, who had rightward asymmetry at a local level (LEEG power), rightward anterior and leftward posterior asymmetries at an intraregional level (first- and second-order S-estimator), and rightward global asymmetry (hemispheric S-estimator), SZ patients showed generally attenuated asymmetry, the effect being strongest for intraregional synchronization. This deviation in asymmetry across the anterior-to-posterior axis is consistent with the cerebral form of the so-called Yakovlevian or anticlockwise cerebral torque. Moreover, the negative occipital and positive frontal asymmetry values suggest higher regional synchronization among the left occipital and the right frontal locations relative to their symmetrical counterparts. Correlation analysis linked the posterior intraregional and hemispheric abnormalities to the negative SZ symptoms, whereas the asymmetry of LEEG power appeared to be weakly coupled to clinical ratings. The posterior intraregional abnormalities of asymmetry were shown to increase with the duration of the disease. The tentative links between these findings and gross anatomical asymmetries, including the cerebral torque and gyrification pattern in normal subjects and SZ patients, are discussed. Conclusions: Overall, our findings reveal the abnormalities in the synchronization asymmetry in SZ patients and heavy involvement of the right hemisphere in these abnormalities. These results indicate that anomalous asymmetry of cortico-cortical connections in schizophrenia is amenable to electrophysiological analysis.

The Sociology of educational mismatch

This paper studies the theoretical relationships between core research lines of sociology such as intergenerational mobility, class structure, cultural capital and educational mismatches. By educational mismatch we mean two things. Firstly an individual can be horizontally mismatched whereby their field of study is inadequate for the job. Another direction of educational mismatch is the so called vertical mismatch where worker possesses more/less education than the job requires resulting in over-/under-education. While analyzing the educational mismatches I keep present the conclusions of Rational Action Theory on individuals’ rational choices in their educational careers. I arrive to conclusions where the influences between educational mismatches and social classes are bidirectional and one can establish fairly clear theoretical links between class of origins and likelihood of being educationally mismatched.

Geological transect across the Northwestern Himalaya in eastern Ladakh and Lahul (A model for the continental collision of India and Asia)

The detailed geological mapping and structural study of a complete transect across the northwestern Himalaya allow to describe the tectonic evolution of the north Indian continental margin during the Tethys ocean opening and the Himalayan Orogeny. The Late Paleozoic Tethys rifting is associated with several tectonomagmatic events. In Upper Lahul and SE Zanskar, this extensional phase is recorded by Lower Carboniferous synsedimentary transtensional faults, a Lower Permian stratigraphic unconformity, a Lower Permian granitic intrusion and middle Permian basaltic extrusions (Panjal Traps). In eastern Ladakh, a Permian listric normal fault is also related to this phase. The scarcity of synsedimentary faults and the gradual increase of the Permian syn-rift sediment thickness towards the NE suggest a flexural type margin. The collision of India and Asia is characterized by a succession of contrasting orogenic phases. South of the Suture Zone, the initiation of the SW vergent Nyimaling-Tsarap Nappe corresponds to an early phase of continental underthrusting. To the S, in Lahul, an opposite underthrusting within the Indian plate is recorded by the NE vergent Tandi Syncline. This structure is associated with the newly defined Shikar Beh Nappe, now partly eroded, which is responsible for the high grade (amphibolite facies) regional metamorphism of South Lahul. The main thrusting of the Nyimaling-Tsarap Nappe followed the formation of the Shikar Beh Nappe. The Nyimaling-Tsarap Nappe developed by ductile shear of the upper part of the subducted Indian continental margin and is responsible for the progressive regional metamorphism of SE Zanskar, reaching amphibolite facies below the frontal part of the nappe, near Sarchu. In Upper Lahul, the frontal parts of the Nyimaling-Tsarap and Shikar Beh nappes are separated by a zone of low grade metamorphic rocks (pumpellyite-actinolite facies to lower greenschist facies). At high structural level, the Nyimaling-Tsarap Nappe is characterized by imbricate structures, which grade into a large ductile shear zone with depth. The related crustal shortening is about 87 km. The root zone and the frontal part of this nappe have been subsequently affected by two zones of dextral transpression and underthrusting: the Nyimaling Shear Zone and the Sarchu Shear Zone. These shear zones are interpreted as consequences of the counterclockwise rotation of the continental underthrusting direction of India relative to Asia, which occurred some 45 and 36 Ma ago, according to plate tectonic models. Later, a phase of NE vergent `'backfolding'' developed on these two zones of dextral transpression, creating isoclinal folds in SE Zanskar and more open folds in the Nyimaling Dome and in the Indus Molasse sediments. During a late stage of the Himalayan Orogeny, the frontal part of the Nyimaling-Tsarap Nappe underwent an extension of about 15 km. This phase is represented by two types of structures, responsible for the tectonic unroofing of the amphibolite facies rocks of the Sarchu area: the Sarchu high angle Normal Fault, cutting a first set of low angle normal faults, which have been created by reactivation of older thrust planes related to the Nyimaling-Tsarap Nappe.

An Analysis of the Link between Ethanol, Energy, and Crop Markets, November 2006

This study analyzes the impact of price shocks in three input and output markets critical to ethanol: gasoline, corn, and sugar. We investigate the impact of these shocks on ethanol and related agricultural markets in the United States and Brazil. We find that the composition of a country’s vehicle fleet determines the direction of the response of ethanol consumption to changes in the gasoline price. We also find that a change in feedstock costs affects the profitability of ethanol producers and the domestic ethanol price. In Brazil, where two commodities compete for sugarcane, changes in the sugar market affect the competing ethanol market.

Dispersal strategies, few dominating or many coexisting: the effect of environmental spatial structure and multiple sources of mortality.

Interspecific competition, life history traits, environmental heterogeneity and spatial structure as well as disturbance are known to impact the successful dispersal strategies in metacommunities. However, studies on the direction of impact of those factors on dispersal have yielded contradictory results and often considered only few competing dispersal strategies at the same time. We used a unifying modeling approach to contrast the combined effects of species traits (adult survival, specialization), environmental heterogeneity and structure (spatial autocorrelation, habitat availability) and disturbance on the selected, maintained and coexisting dispersal strategies in heterogeneous metacommunities. Using a negative exponential dispersal kernel, we allowed for variation of both species dispersal distance and dispersal rate. We showed that strong disturbance promotes species with high dispersal abilities, while low local adult survival and habitat availability select against them. Spatial autocorrelation favors species with higher dispersal ability when adult survival and disturbance rate are low, and selects against them in the opposite situation. Interestingly, several dispersal strategies coexist when disturbance and adult survival act in opposition, as for example when strong disturbance regime favors species with high dispersal abilities while low adult survival selects species with low dispersal. Our results unify apparently contradictory previous results and demonstrate that spatial structure, disturbance and adult survival determine the success and diversity of coexisting dispersal strategies in competing metacommunities.

Equilibrium asset prices and savings of heterogeneous agents in the presence of incomplete markets and portfolio constraints

We study the quantitative properties of a dynamic general equilibrium model in which agents face both idiosyncratic and aggregate income risk, state-dependent borrowing constraints that bind in some but not all periods and markets are incomplete. Optimal individual consumption-savings plans and equilibrium asset prices are computed under various assumptions about income uncertainty. Then we investigate whether our general equilibrium model with incomplete markets replicates two empirical observations: the high correlation between individual consumption and individual income, and the equity premium puzzle. We find that, when the driving processes are calibrated according to the data from wage income in different sectors of the US economy, the results move in the direction of explaining these observations, but the model falls short of explaining the observed correlations quantitatively. If the incomes of agents are assumed independent of each other, the observations can be explained quantitatively.

The genetical theory of social behaviour.

We survey the population genetic basis of social evolution, using a logically consistent set of arguments to cover a wide range of biological scenarios. We start by reconsidering Hamilton's (Hamilton 1964 J. Theoret. Biol. 7, 1-16 (doi:10.1016/0022-5193(64)90038-4)) results for selection on a social trait under the assumptions of additive gene action, weak selection and constant environment and demography. This yields a prediction for the direction of allele frequency change in terms of phenotypic costs and benefits and genealogical concepts of relatedness, which holds for any frequency of the trait in the population, and provides the foundation for further developments and extensions. We then allow for any type of gene interaction within and between individuals, strong selection and fluctuating environments and demography, which may depend on the evolving trait itself. We reach three conclusions pertaining to selection on social behaviours under broad conditions. (i) Selection can be understood by focusing on a one-generation change in mean allele frequency, a computation which underpins the utility of reproductive value weights; (ii) in large populations under the assumptions of additive gene action and weak selection, this change is of constant sign for any allele frequency and is predicted by a phenotypic selection gradient; (iii) under the assumptions of trait substitution sequences, such phenotypic selection gradients suffice to characterize long-term multi-dimensional stochastic evolution, with almost no knowledge about the genetic details underlying the coevolving traits. Having such simple results about the effect of selection regardless of population structure and type of social interactions can help to delineate the common features of distinct biological processes. Finally, we clarify some persistent divergences within social evolution theory, with respect to exactness, synergies, maximization, dynamic sufficiency and the role of genetic arguments.

A generalization of histogram type estimators

We introduce simple nonparametric density estimators that generalize theclassical histogram and frequency polygon. The new estimators are expressed as linear combination of density functions that are piecewisepolynomials, where the coefficients are optimally chosen in order to minimize the integrated square error of the estimator. We establish the asymptotic behaviour of the proposed estimators, and study theirperformance in a simulation study.

The effect of intergroup competition on group coordination: An experimental study

We report an experiment on the effect of intergroup competition on group coordination in the minimal-effort game (Van Huyck et al., 1990). The competition was between two 7-person groups. Each player in each group independently chose an integer from 1 to 7. The group with the higher minimum won the competition and each of its members was paid according to the game s original payoff matrix. Members of the losing group were paid nothing. In case of a tie, each player was paid half the payoff in the original matrix. This treatment was contrasted with two control treatments where each of the two groups played an independent coordination game, either with or without information about the minimum chosen by the outgroup. Although the intergroup competition does not change the set of strict equilibria, we found that it improved collective rationality by moving group members in the direction of higher-payoff equilibria. Merely providing group members with information about the minimal-effort level in the other group was not sufficient to generate this effect.

Technological change and the wealth of nations

We discuss a unified theory of directed technological change and technologyadoption that can shed light on the causes of persistent productivity differencesacross countries. In our model, new technologies are designed in advanced countries and diffuse endogenously to less developed countries. Our framework is richenough to highlight three broad reasons for productivity differences: inappropriatetechnologies, policy-induced barriers to technology adoption, and within-countrymisallocations across sectors due to policy distortions. We also discuss the effectsof two aspects of globalization, trade in goods and migration, on the wealth ofnations through their impact on the direction of technical progress. By doing so,we illustrate some of the equalizing and unequalizing forces of globalization.

Minimum distance estimation of stationary and non-stationary ARFIMA processes

A new parametric minimum distance time-domain estimator for ARFIMA processes is introduced in this paper. The proposed estimator minimizes the sum of squared correlations of residuals obtained after filtering a series through ARFIMA parameters. The estimator iseasy to compute and is consistent and asymptotically normally distributed for fractionallyintegrated (FI) processes with an integration order d strictly greater than -0.75. Therefore, it can be applied to both stationary and non-stationary processes. Deterministic components are also allowed in the DGP. Furthermore, as a by-product, the estimation procedure provides an immediate check on the adequacy of the specified model. This is so because the criterion function, when evaluated at the estimated values, coincides with the Box-Pierce goodness of fit statistic. Empirical applications and Monte-Carlo simulations supporting the analytical results and showing the good performance of the estimator in finite samples are also provided.