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

Development and Validation of Single Nucleotide Polymorphisms (SNPs) Markers from Two Transcriptome 454-Runs of Turbot (Scophthalmus maximus) Using High-Throughput Genotyping

The turbot (Scophthalmus maximus) is a commercially valuable flatfish and one of the most promising aquaculture species in Europe. Two transcriptome 454-pyrosequencing runs were used in order to detect Single Nucleotide Polymorphisms (SNPs) in genesrelated to immune response and gonad differentiation. A total of 866 true SNPs were detected in 140 different contigs representing 262,093 bp as a whole. Only one true SNP was analyzed in each contig. One hundred and thirteen SNPs out of the 140 analyzed were feasible (genotyped), while Ш were polymorphic in a wild population. Transition/transversion ratio (1.354) was similar to that observed in other fish studies. Unbiased gene diversity (He) estimates ranged from 0.060 to 0.510 (mean = 0.351), minimum allele frequency (MAF) from 0.030 to 0.500 (mean = 0.259) and all loci were in Hardy-Weinberg equilibrium after Bonferroni correction. A large number of SNPs (49) were located in the coding region, 33 representing synonymous and 16 non-synonymous changes. Most SNP-containing genes were related to immune response and gonad differentiation processes, and could be candidates for functional changes leading to phenotypic changes. These markers will be useful for population screening to look for adaptive variation in wild and domestic turbot

Characterization of thirty two microsatellite loci for three Atlanto-Mediterranean echinoderm species

Thirty two microsatellites were optimized from 454 pyrosequencing libraries for three Atlanto-Mediterranean echinoderms: Coscinasterias tenuispina, Echinaster sepositus and Arbacia lixula. We observed different frequency of microsatellite types (di-, tri-, tetra- and pentanucleotide) throughout the genome of the species, but no significant differences were observed in allele richness among different microsatellite repeats. No loci showed linkage disequilibrium. Heterozygosity deficit and departure from Hardy Weinberg equilibrium were observed for some loci, in two species, probably due to high levels of inbreeding. Heterozygosity excess observed in C. tenuispina could be explained by selection against homozygotes and/or outcrossing.

Inversion polymorphism in two Serbian natural populations of Drosophila subobscura: Analysis of long-term changes

To study whether inversions (or arrangements) by themselves or karyotypes are the main global warming adaptive target of natural selection, two Drosophila subobscura Serbian populations (Apatin and Petnica) were re analyzed using different statistical approaches. Both populations were sampled in an approximately 15 years period: Apatin in 1994 and 2008 + 2009 and Petnica in 1995 and 2010. For all chromosomes, the four collections studied were in Hardy-Weinberg equilibrium. Thus, it seems that inversions (or arrangements) combined at random to constitute populations" karyotypes. However, there were differences in karyotypic fre quencies along the years, although they were significant only for Apatin population. It is possible to conclude that inversions (or arrangements) are likely the target of natural selection, because they presented long term changes, but combine at random to generate the corresponding karyotypic combinations. As a consequence, the frequencies of karyotypes also change along time.

Characterization of thirty two microsatellite loci for three Atlanto-Mediterranean echinoderm species

Thirty two microsatellites were optimized from 454 pyrosequencing libraries for three Atlanto-Mediterranean echinoderms: Coscinasterias tenuispina, Echinaster sepositus and Arbacia lixula. We observed different frequency of microsatellite types (di-, tri-, tetra- and pentanucleotide) throughout the genome of the species, but no significant differences were observed in allele richness among different microsatellite repeats. No loci showed linkage disequilibrium. Heterozygosity deficit and departure from Hardy Weinberg equilibrium were observed for some loci, in two species, probably due to high levels of inbreeding. Heterozygosity excess observed in C. tenuispina could be explained by selection against homozygotes and/or outcrossing.

Quantifying rock fabrics: a test of autocorrelation of the spatial distribution of cristals

A novel test of spatial independence of the distribution of crystals or phases in rocksbased on compositional statistics is introduced. It improves and generalizes the commonjoins-count statistics known from map analysis in geographic information systems.Assigning phases independently to objects in RD is modelled by a single-trial multinomialrandom function Z(x), where the probabilities of phases add to one and areexplicitly modelled as compositions in the K-part simplex SK. Thus, apparent inconsistenciesof the tests based on the conventional joins{count statistics and their possiblycontradictory interpretations are avoided. In practical applications we assume that theprobabilities of phases do not depend on the location but are identical everywhere inthe domain of de nition. Thus, the model involves the sum of r independent identicalmultinomial distributed 1-trial random variables which is an r-trial multinomialdistributed random variable. The probabilities of the distribution of the r counts canbe considered as a composition in the Q-part simplex SQ. They span the so calledHardy-Weinberg manifold H that is proved to be a K-1-affine subspace of SQ. This isa generalisation of the well-known Hardy-Weinberg law of genetics. If the assignmentof phases accounts for some kind of spatial dependence, then the r-trial probabilitiesdo not remain on H. This suggests the use of the Aitchison distance between observedprobabilities to H to test dependence. Moreover, when there is a spatial uctuation ofthe multinomial probabilities, the observed r-trial probabilities move on H. This shiftcan be used as to check for these uctuations. A practical procedure and an algorithmto perform the test have been developed. Some cases applied to simulated and realdata are presented.Key words: Spatial distribution of crystals in rocks, spatial distribution of phases,joins-count statistics, multinomial distribution, Hardy-Weinberg law, Hardy-Weinbergmanifold, Aitchison geometry

The structural approach of a natrex model on equilibrium exchange rates

Following a general macroeconomic approach, this paper sets a closed micro-founded structural model to determine the long run real exchange rate of a developed economy. In particular, the analysis follows the structure of a Natrex model. The main contribution of this research paper is the development of a solid theoretical framework that analyse in depth the basis of the real exchange rate and the details of the equilibrium dynamics after any shock influencing the steady state. In our case, the intertemporal factors derived from the stock-flow relationship will be particularly determinant. The main results of the paper can be summarised as follows. In first place, a complete well-integrated structural model for long-run real exchange rate determination is developed from first principles. Moreover, within the concrete dynamics of the model, it is found that some convergence restrictions will be necessary. On one hand, for the medium run convergence the sensitivity of the trade balance to changes in real exchange rate should be higher that the correspondent one to the investment decisions. On the other hand, and regarding long-run convergence, it is also necessary both that there exists a negative relationship between investment and capital stock accumulation and that the global saving of the economy depends positively on net foreign debt accumulation. In addition, there are also interesting conclusions about the effects that certain shocks over the exogenous variables of the model have on real exchange rates.

The effect of candidate quality on electoral equilibrium: an experimental study

When two candidates of different quality compete in a one dimensional policy space, the equilibrium outcomes are asymmetric and do not correspond to the median. There are three main effects. First, the better candidate adopts more centrist policies than the worse candidate. Second, the equilibrium is statistical, in the sense that it predicts a probability distribution of outcomes rather than a single degenerate outcome. Third, the equilibrium varies systematically with the level of uncertainty about the location of the median voter. We test these three predictions using laboratory experiments, and find strong support for all three. We also observe some biases and show that they canbe explained by quantal response equilibrium.

General equilibrium with asymmetric information: a dual approach

We study markets where the characteristics or decisions of certain agents are relevant but not known to their trading partners. Assuming exclusive transactions, the environment is described as a continuum economy with indivisible commodities. We characterize incentive efficient allocations as solutions to linear programming problems and appeal to duality theory to demonstrate the generic existence of external effects in these markets. Because under certain conditions such effects may generate non-convexities, randomization emerges as a theoretic possibility. In characterizing market equilibria we show that, consistently with the personalized nature of transactions, prices are generally non-linear in the underlying consumption. On the other hand, external effects may have critical implications for market efficiency. With adverse selection, in fact, cross-subsidization across agents with different private information may be necessary for optimality, and so, the market need not even achieve an incentive efficient allocation. In contrast, for the case of a single commodity, we find that when informational asymmetries arise after the trading period (e.g. moral hazard; ex post hidden types) external effects are fully internalized at a market equilibrium.

Consumption externalities, habit formation, and equilibrium efficiency

We analyze the welfare properties of the competitive equilibrium in a capital accumulation model where individual preferences are subjected to both habit formation and consumption spillovers. We also discuss how consumption externalities and habits interact to generate an inefficient dynamic equilibrium. Finally, we characterize optimal tax policies aimed to restore efficient decentralized paths.

Real equilibrium exchange rates. A panel data approach for advanced and emerging economies.

Based on an behavioral equilibrium exchange rate model, this paper examines the determinants of the real effective exchange rate and evaluates the degree of misalignment of a group of currencies since 1980. Within a panel cointegration setting, we estimate the relationship between exchange rate and a set of economic fundamentals, such as traded-nontraded productivity differentials and the stock of foreign assets. Having ascertained the variables are integrated and cointegrated, the long-run equilibrium value of the fundamentals are estimated and used to derive equilibrium exchange rates and misalignments. Although there is statistical homogeneity, some structural differences were found to exist between advanced and emerging economies.

Equilibrium play and best response to (Stated) beliefs in constant sum games

We report experimental results on one-shot two person 3x3 constant sum games played by non-economists without previous experience in the laboratory. Although strategically our games are very similar to previous experiments in which game theory predictions fail dramatically, 80% of actions taken in our experiment coincided with the prediction of the unique Nash equilibrium in pure strategies and 73% of actions were best responses to elicited beliefs. We argue how social preferences, presentation effects and belief elicitation procedures may influence how subjects play in simple but non trivial games and explain the diferences we observe with respect to previous work.

Weighted Hardy spaces for the unit disc: approximation properties

In this paper we study basic properties of the weighted Hardy space for the unit disc with the weight function satisfying Muckenhoupt's (Aq) condition, and study related approximation problems (expansion, moment and interpolation) with respect to two incomplete systems of holomorphic functions in this space.

A Computable general equilibrium approach to hypothetical extractions and missing links

Identifying key sectors or key locations in an interconnected economy is of paramount importance for improving policy planning and directing economic strategy. Hence the relevance of categorizing them and hence the corresponding need of evaluating their potential synergies in terms of their global economic thrust. We explain in this paper that standard measures based on gross outputs do not and cannot capture the relevant impact due to self- imposed modeling limitations. In fact, common gross output measures will be systematically downward biased. We argue that an economy wide Computable General Equilibrium (CGE) approach provides a modeling platform that overcomes these limitations since it provides (i) a more comprehensive measure of linkages and (ii) an alternate way of accounting for links' relevance that is in consonance with standard macromagnitudes in the National Income and Product Accounts.

Non-hyperbolic boundary equilibrium bifurcations in planar Filippov systems: a case study approach

Boundary equilibrium bifurcations in piecewise smooth discontinuous systems are characterized by the collision of an equilibrium point with the discontinuity surface. Generically, these bifurcations are of codimension one, but there are scenarios where the phenomenon can be of higher codimension. Here, the possible collision of a non-hyperbolic equilibrium with the boundary in a two-parameter framework and the nonlinear phenomena associated with such collision are considered. By dealing with planar discontinuous (Filippov) systems, some of such phenomena are pointed out through specific representative cases. A methodology for obtaining the corresponding bi-parametric bifurcation sets is developed.