Hardy-Weinberg Equilibrium and the Ternary Plot

The Hardy-Weinberg law, formulated about 100 years ago, states that under certainassumptions, the three genotypes AA, AB and BB at a bi-allelic locus are expected to occur inthe proportions p2, 2pq, and q2 respectively, where p is the allele frequency of A, and q = 1-p.There are many statistical tests being used to check whether empirical marker data obeys theHardy-Weinberg principle. Among these are the classical xi-square test (with or withoutcontinuity correction), the likelihood ratio test, Fisher's Exact test, and exact tests in combinationwith Monte Carlo and Markov Chain algorithms. Tests for Hardy-Weinberg equilibrium (HWE)are numerical in nature, requiring the computation of a test statistic and a p-value.There is however, ample space for the use of graphics in HWE tests, in particular for the ternaryplot. Nowadays, many genetical studies are using genetical markers known as SingleNucleotide Polymorphisms (SNPs). SNP data comes in the form of counts, but from the countsone typically computes genotype frequencies and allele frequencies. These frequencies satisfythe unit-sum constraint, and their analysis therefore falls within the realm of compositional dataanalysis (Aitchison, 1986). SNPs are usually bi-allelic, which implies that the genotypefrequencies can be adequately represented in a ternary plot. Compositions that are in exactHWE describe a parabola in the ternary plot. Compositions for which HWE cannot be rejected ina statistical test are typically “close" to the parabola, whereas compositions that differsignificantly from HWE are “far". By rewriting the statistics used to test for HWE in terms ofheterozygote frequencies, acceptance regions for HWE can be obtained that can be depicted inthe ternary plot. This way, compositions can be tested for HWE purely on the basis of theirposition in the ternary plot (Graffelman & Morales, 2008). This leads to nice graphicalrepresentations where large numbers of SNPs can be tested for HWE in a single graph. Severalexamples of graphical tests for HWE (implemented in R software), will be shown, using SNPdata from different human populations

Intrinsic test of independence in contingency tables

A condition needed for testing nested hypotheses from a Bayesianviewpoint is that the prior for the alternative model concentratesmass around the small, or null, model. For testing independencein contingency tables, the intrinsic priors satisfy this requirement.Further, the degree of concentration of the priors is controlled bya discrete parameter m, the training sample size, which plays animportant role in the resulting answer regardless of the samplesize.In this paper we study robustness of the tests of independencein contingency tables with respect to the intrinsic priors withdifferent degree of concentration around the null, and comparewith other “robust” results by Good and Crook. Consistency ofthe intrinsic Bayesian tests is established.We also discuss conditioning issues and sampling schemes,and argue that conditioning should be on either one margin orthe table total, but not on both margins.Examples using real are simulated data are given

Pattern of recurrence of early breast cancer is different according to intrinsic subtype and proliferation index.

INTRODUCTION Recurrence risk in breast cancer varies throughout the follow-up time. We examined if these changes are related to the level of expression of the proliferation pathway and intrinsic subtypes. METHODS Expression of estrogen and progesterone receptor, Ki-67, human epidermal growth factor receptor 2 (HER2), epidermal growth factor receptor (EGFR) and cytokeratin 5/6 (CK 5/6) was performed on tissue-microarrays constructed from a large and uniformly managed series of early breast cancer patients (N = 1,249). Subtype definitions by four biomarkers were as follows: luminal A (ER + and/or PR+, HER2-, Ki-67 <14), luminal B (ER + and/or PR+, HER2-, Ki-67 ≥14), HER2-enriched (any ER, any PR, HER2+, any Ki-67), triple-negative (ER-, PR-, HER2-, any Ki-67). Subtype definitions by six biomarkers were as follows: luminal A (ER + and/or PR+, HER2-, Ki-67 <14, any CK 5/6, any EGFR), luminal B (ER + and/or PR+, HER2-, Ki-67 ≥14, any CK 5/6, any EGFR), HER2-enriched (ER-, PR-, HER2+, any Ki-67, any CK 5/6, any EGFR), Luminal-HER2 (ER + and/or PR+, HER2+, any Ki-67, any CK 5/6, any EGFR), Basal-like (ER-, PR-, HER2-, any Ki-67, CK5/6+ and/or EGFR+), triple-negative nonbasal (ER-, PR-, HER2-, any Ki-67, CK 5/6-, EGFR-). Each four- or six-marker defined intrinsic subtype was divided in two groups, with Ki-67 <14% or with Ki-67 ≥14%. Recurrence hazard rate function was determined for each intrinsic subtype as a whole and according to Ki-67 value. RESULTS Luminal A displayed a slow risk increase, reaching its maximum after three years and then remained steady. Luminal B presented most of its relapses during the first five years. HER2-enriched tumors show a peak of recurrence nearly twenty months post-surgery, with a greater risk in Ki-67 ≥14%. However a second peak occurred at 72 months but the risk magnitude was greater in Ki-67 <14%. Triple negative tumors with low proliferation rate display a smooth risk curve, but with Ki-67 ≥14% show sharp peak at nearly 18 months. CONCLUSIONS Each intrinsic subtype has a particular pattern of relapses over time which change depending on the level of activation of the proliferation pathway assessed by Ki-67. These findings could have clinical implications both on adjuvant treatment trial design and on the recommendations concerning the surveillance of patients.

MIinimal N-point diameters and f-best-packing constants in Rd[superindex]

Vegeu el resum a l'inici del document de l'arxiu adjunt

On concentrators and related approximation constants

Pippenger [Pi77] showed the existence of (6m,4m,3m,6)-concentrator for each positive integer m using a probabilistic method. We generalize his approach and prove existence of (6m,4m,3m,5.05)-concentrator (which is no longer regular, but has fewer edges). We apply this result to improve the constant of approximation of almost additive set functions by additive set functions from 44.5 (established by Kalton and Roberts in [KaRo83] to 39. We show a more direct connection of the latter problem to the Whitney type estimate for approximation of continuous functions on a cube in &b&R&/b&&sup&d&/sup& by linear functions, and improve the estimate of this Whitney constant from 802 (proved by Brudnyi and Kalton in [BrKa00] to 73.

Bayesian nash equilibrium in ''linear'' cournot models with private information about costs

Calculating explicit closed form solutions of Cournot models where firms have private information about their costs is, in general, very cumbersome. Most authors consider therefore linear demands and constant marginal costs. However, within this framework, the nonnegativity constraint on prices (and quantities) has been ignored or not properly dealt with and the correct calculation of all Bayesian Nash equilibria is more complicated than expected. Moreover, multiple symmetric and interior Bayesianf equilibria may exist for an open set of parameters. The reason for this is that linear demand is not really linear, since there is a kink at zero price: the general ''linear'' inverse demand function is P (Q) = max{a - bQ, 0} rather than P (Q) = a - bQ.

Social capital and the equilibrium number of entrepreneurs

Social capital is viewed either as a proprietary asset that serves private interests, including those of entrepreneurs, or as a collective asset that supports trust-based transactions saving on transaction costs both in markets and within the boundaries of firms, and benefiting society as a whole. This paper explains the relative specialization between entrepreneurs and market-governed exchanges as a result of the interaction between social capital that lowers transaction costs, and the scale economies of ability in managerial jobs (Lucas 1978). The main hypothesis formulated in the paper is that higher social capital will benefit the hierarchy relatively more than the market as a governance mechanism, and therefore in economies with higher social capital, the equilibrium number of entrepreneurs will be lower and their average span of control larger than in economies with lower social capital. The empirical evidence, with data from the Spanish Autonomous Communities, is consistent with this prediction.

Mean intrinsic tensile properties of stone groundwood fibers from softwood

Stone groundwood (SGW) is a fibrous matter commonly prepared in a high yield process, and mainly used for papermaking applications. In this work, the use of SGW fibers is explored as reinforcing element of polypropylene (PP) composites. Due to its chemical and superficial features, the use of coupling agents is needed for a good adhesion and stress transfer across the fiber-matrix interface. The intrinsic strength of the reinforcement is a key parameter to predict the mechanical properties of the composite and to perform an interface analysis. The main objective of the present work was the determination of the intrinsic tensile strength of stone groundwood fibers. Coupled and non-coupled PP composites from stone groundwood fibers were prepared. The influence of the surface morphology and the quality at interface on the final properties of the composite was analyzed and compared to that of fiberglass PP composites. The intrinsic tensile properties of stone groundwood fibers, as well as the fiber orientation factor and the interfacial shear strength of the current composites were determined

Equilibrium of drops on inclined hydrophilic surfaces

Why does not gravity make drops slip down the inclined surfaces, e.g., plant leaves? The current explanation is based on the existence of surface inhomogeneities, which cause a sustaining force that pins the contact line. Following this theory, the drop remains in equilibrium until a critical value of the sustaining force is reached. We propose an alternative analysis, from the point of view of energy balance, for the particular case in which the drop leaves a liquid film behind. The critical angle of the inclined surface at which the drop slips down is predicted. This result does not depend explicitly on surface inhomogeneities, but only on the drop size and surface tensions. There is good agreement with experiments for contact angles below 90° where the formation of the film is expected, whereas for greater contact angles great discrepancies arise

Disseny, implementació i optimització d'un algorisme per al càlcul de les constants d'isocronia

S'ha desenvolupat una eina de programari capaç de calcular les constants d'isocronia. S'obté també per pantalla dades referents al cost temporal de l'algorisme de càlcul. Les constants d'isocronia obtingudes són recuperables pel Mathematica per a la seva posterior utilització. El càlcul de les constants d'isocronia és molt més òptim temporalment que el càlcul realitzat pel Mathematica. S'ha tingut en compte també l'optimització en la utilització d'espai de memòria.

Several neuronal and axonal types form long intrinsic connections in the cat primary auditory cortical field (AI).

Intrinsic connections in the cat primary auditory field (AI) as revealed by injections of Phaseolus vulgaris leucoagglutinin (PHA-L) or biocytin, had an anisotropic and patchy distribution. Neurons, labelled retrogradely with PHA-L were concentrated along a dorsoventral stripe through the injection site and rostral to it; the spread of rostrally located neurons was greater after injections into regions of low rather than high characteristic frequencies. The intensity of retrograde labelling varied from weak and granular to very strong and Golgi-like. Out of 313 Golgi like retrogradely labelled neurons 79.6% were pyramidal, 17.2% multipolar, 2.6% bipolar, and 0.6% bitufted; 13.4% were putatively inhibitory, i.e. aspiny or sparsely spiny multipolar, or bitufted. Individual anterogradely labelled intrinsic axons were reconstructed for distances of 2 to 7 mm. Five main types were distinguished on the basis of the branching pattern and the location of synaptic specialisations. Type 1 axons travelled horizontally within layers II to VI and sent collaterals at regular intervals; boutons were only present in the terminal arborizations of these collaterals. Type 2 axons also travelled horizontally within layers II to VI and had rather short and thin collateral branches; boutons or spine-like protrusions occurred in most parts of the axon. Type 3 axons travelled obliquely through the cortex and formed a single terminal arborization, the only site where boutons were found. Type 4 axons travelled for some distance in layer I; they formed a heterogeneous group as to their collaterals and synaptic specializations. Type 5 axons travelled at the interface between layer VI and the white matter; boutons en passant, spine-like protrusions, and thin short branches with boutons en passant were frequent all along their trajectory. Thus, only some axonal types sustain the patchy pattern of intrinsic connectivity, whereas others are involved in a more diffuse connectivity.

Kinetics of crystal growth and equilibrium domains in eclogite of the Sesia Zone, Western Alps

Kinetics of crystal growth and equilibrium domains in eclogite of the Sesia Zone, Western Alps Darbellay Bastien Institut de Minéralogie et Géochimie Résumé grand public Comme toute matière, la roche est sensible à son environnement et cherche à s'adapter pour acquérir un état stable (état d'équilibre). Les changements des conditions physiques (température et pression) vont ainsi impliquer des modifications dans la roche. Le métamorphisme est l'étude de ces changements. Les minéraux qui constituent la roche peuvent modifier, leur structure, leur chimie ou être remplacer par d'autres minéraux plus stables. Il est ainsi crucial de déterminer les processus responsables et limitant de la croissance minérale. Trois processus permettent la croissance ; (1) la dissolution des éléments du réactant, (2) le transport de ces éléments vers le site de croissance, (3) l'incorporation de ces éléments dans la nouvelle structure. Cette thèse se focalise sur les structures des minéraux de haute pression (forme, zonation chimique, structure interne) pour essayer de déterminer les facteurs importants à l'origine de leur état final. Les zones d'étude se situent dans la zone de Sésia. La première partie traite de la problématique liée à l'incorporation d'un élément dans une structure minérale. A l'image de la croissance humaine, les irrégularités minéralogiques permettent de mettre en lumière un dysfonctionnement de la croissance due à un excès ou à une carence d'un élément. Bien dosé, cet élément est cependant essentiel à la croissance. Les zoisites (épidotes) des métabasites de la région de Cima di Bonze montrent une zonation chimique en sablier. Dans cette zonation la teneur en fer excède la capacité maximum que peut contenir la structure orthorhombique de la zoisite. Des défauts de structure permettent l'accommodation de cet excès. La zoisite peut ainsi adapter sa structure pour permettre l'incorporation d'une relativement grande quantité de fer. Les études précédentes montraient, pour des conditions similaires, la formation de deux épidotes distinctes. La deuxième partie se penche sur la compétition entre le minéral qui fait sa croissance et les minéraux (réactants) qui l'entourent. Les métapélites de la région du Monte Mucrone contiennent des grenats atollaires. Des études détaillées de la texture et de la zonation chimique du grenat ainsi qu'une modélisation thermodynamique ont permis de mieux cerner les facteurs importants responsables de la forme atollaire. Cette structure est obtenue par un changement du comportement de la croissance du grenat le long d'un chemin P-T hercynien. Dans un premier stade, le grenat croît rapidement et consume peu le quartz de la matrice. La croissance se fait ainsi le long des jointures des grains de quartz. Dans un second temps, les changements de conditions PT donnent une croissance lente du grenat et une forte consommation du quartz. Le grenat peut ainsi développer sa forme dodécaédrale classique. La troisième partie s'intéresse aux distances de transport par diffusion d'un élément (ici l'argon) durant la haute pression. Pour ce faire, un profile d'âges 40Ar/39Ar sur biotite a été mesuré depuis un veine de haute pression riche en argon jusque dans son encaissant (granitoïd du Monte Mucrone). Le profile montre une répartition des âges suivant une courbe de diffusion. Le transport se fait sur une longueur de deux centimètre avec l'aide d'un fluide. Il est réduit à une échelle millimétrique quand la phase fluide disparaît. Cette étude montre ainsi les difficultés de transport des éléments durant la haute pression ne permettant pas un rééquilibrage de la roche à grande échelle. Kinetics of crystal growth and equilibrium domains in eclogite of the Sesia Zone, Western Alps Darbellay Bastien Institut de Minéralogie et Géochimie Résumé de thèse Les processus de croissance (diffusion des éléments et les réactions d'interface) et les conditions dans lesquelles les minéraux grandissent (température, pression, fluide, composition chimique de la roche), déterminent la texture ainsi que la zonation des minéraux. Cette thèse se focalise, par le biais de textures peu communes, sur trois différents processus impliqués dans la croissance minérale à haute pression (Zone de Sésia, Alpes de l'Ouest, Italie). L'incorporation d'un élément dans une structure minérale ne peut se faire que dans des sites en accord avec la taille et la charge ionique de l'élément. De plus, la balance de charge doit être maintenue dans le minéral. La régularité de la structure cristalline fixe ainsi une limite maximum de concentration d'un élément donné. Les zoisites provenant des métabasites de la région de Cima di Bonze montrent des zonations en sablier caractérisées par une concentration anormale en fer. La zonation se marque par une différente teinte de biréfringence et par un plus grand angle d'extinction que le reste de la zoisite. Une inter-croissance de clinozoisite à l'intérieur de la structure orthorhombique de la zoisite peut ainsi être suspectée. Les analyses XRD (diffraction des rayons x) ainsi que les analyses Raman ne confirment pas cette suspicion. Seules les analyses TEM (microscope à électrons transmis) montrent des défauts de structure pouvant être interprétés comme des modules de clinozoisite. Ils ne peuvent cependant pas être considérés comme une phase thermodynamique. Un nouveau trou d'immiscibilité entre deux zoisite (X ep= 0.1 and Xep = 0.15) a ainsi pu être établi. Dans les métapélites la région du Monte Mucrone, des grenats fortement zonés montrent une évolution texturale singulière. Ils présentent une forme initiale de `champignon' qui se développe pour former une structure atollaire finale. L'étude conjuguée de la structure 3D et des zonations, ainsi que l'établissement d'un model thermodynamique, indiquent que ces structures proviennent de deux épisodes de croissances : (1) La croissance du grenat durant un chemin prograde hercynien (de 525 °C et 6.2 kbar à 640 °C et 9 kbar) permet la formation des textures atollaires. Elles sont le résultat d'une croissance poecilitique initiale suivie d'une croissance idiomorphique du grenat. (2) La structure est rendue plus complexe par la cristallisation d'un grenat homogène tout autour ainsi qu'à l'intérieur du grenat hercynien durant la haute pression alpine (550 °C and 20 kbar). L'arrivée de l'eau durant la haute pression facilite le transport d'éléments et permet une cristallisation rapide du grenat. La diffusion peut être un facteur limitant de la croissance minéralogique. Elle a aussi une grande importance pour la géochronologie. Une veine de haute pression à l'intérieur du granitoïde du Monte Mucrone a été étudiée dans le but de déterminer la distance de diffusion de l'argon. Le profile d'âges 40Ar/39Ar sur biotites, établi de la veine vers le métagranitoïde, suit une courbe de diffusion. Les âges sont élevés proche de la veine (800 Ma) puis décroissent jusqu'à des âges homogènes (170-150 Ma) à deux centimètres de la veine. La présence de fluide, marqué par de hautes concentrations en chlore, permet une diffusion centimétrique. Cependant, la distance est réduite à une échelle millimétrique quand le fluide est absent. Les très faibles distances de diffusion préservent les âges pré-alpins et impliquent un événement géologique pour les âges de 170-150 Ma. Kinetics of crystal growth and equilibrium domains in eclogite of the Sesia Zone, Western Alps Darbellay Bastien Institut de Minéralogie et Géochimie Thesis abstract Rock textures and zonings are the consequence of growth processes (element diffusion and interface reaction) steered by the environment in which they grew (pressure, differential stress, temperature, fluid and rock composition). The thesis presented here focuses on three different topics, each of it dealing with aspects of mineral growth processes during subduction, in a high-pressure environment. All studies were conducted in the Sesia Zone of the Western European Alps, Italy. The first study addresses the crystallography and geochemistry of element incorporation in zoisite, one of the major hydrous minerals found in subduction zone rocks. Elements can be incorporated into a mineral structure only on crystallographic sites that offer enough space for the ion and the overall charge balance has to be maintained. Element concentrations are hence limited. Incorporation of some elements produces complex zoning, including hourglass like patterns, which are the focus of the first contribution. Zoisites from Cima di Bonze (Sesia Zone) show spectacular hourglass zoning defined by Fe-content variations. The hourglass zones have a distinct birefringence and a different extinction angle than the regular part of the zoisite. We show by detailed XRD (X-ray diffraction) and confocal Raman analyses that the high Fe-zones are nevertheless zoisite, and not clinozoisite as one might expect. High resolution TEM (transmission electron microscopy) analyses show planar defects on (100) that can be interpreted as small-scale clinozoisite modules. However, these clinozoisites cannot be interpreted as a distinctive thermodynamic phase and the entire mineral has to be considered as zoisite. The miscibility gap between two zoisites (Xep = 0.1 and Xep = 0.15) can be then definite at 550 ± 50°C and 14 to 20 Kbar. Strongly zoned garnets in quartz rich metapelite from the Monte Mucrone area (Sesia Zone) show evolution form 3D mushroom to atoll structure. The second contribution presents textural investigations, garnet zoning and thermodynamic modeling that demonstrate that atoll garnets are the result of two distinctive growth events. (1) Garnet atoll structure is already formed during a prograde Hercynian path from 525 °C and 6.2 kbar to 640 °C and 9 kbar. It results in an initial poikilitic growth followed by a final idiomorphic growth event. (2) Alpine HP garnet are homogenous (550 °C and 20 kbar) and grew around and also inside the Hercynian garnet. Lack of prograde Alpine garnet and fast growth of the HP garnet is explained by the absence of water during much of the prograde path. Water saturation was only observed towards the end, close towards the peak metamorphic conditions. Diffusion could be a limiting factor for crystal growth. It has also a great importance in geochronology. HP vein inside the metagranitoide of the Monte Mucrone (~300 Ma) was investigated to determine argon diffusion scales during high-pressure metamorphism. 40Ar/39Ar biotite ages profile from the vein toward the metagranodiorite show a diffusion curve: old ages (800 Ma) located close to the vein decrease until homogenous 170-150 Ma ages are obtained, two centimeter away from the vein. Centimeter-scale diffusion occurs with help of a fluid phase marked by high chlorine concentrations. Argon diffusion is reduced to a millimeter scale when free fluid is absent. Very short diffusion distance permits to preserve pre-Alpine ages. The 170-150 Ma ages are considered to be geologic meaningful, probably resulting from the extensional tectonics linked to opening of the Tethian ocean.

Feeding and the Equilibrium Feeder Animal Price-Weight Schedule, June 2005

The feeder animal price is a derivative in the sense that its value depends upon the price of animals for the consumption market. It also depends upon the biological growth technology and feed costs. Daily maintenance costs are of particular interest to the husbander because they can be avoided through accelerated feeding. In this paper, the optimal feeding path under equilibrium feeder animal prices is established. This analysis is used to gain a better understanding of feeding decisions, regulation in feedstuff markets, and the consequences of genetic innovations. It is shown that days on feed can increase or decrease with a genetic innovation or other improvement in feed conversion efficiency. The structure of comparative prices for feeder animals at different weights, the early slaughter decision, and equilibrium in feeder animal markets are also developed. Feeder animal prices can increase over a weight interval if biological feed efficiency parameters are low over the interval.