Unconstrained optimization with MINTOOLKIT for GNU octave

The paper documents MINTOOLKIT for GNU Octave. MINTOOLKIT provides functions for minimization and numeric differentiation. The main algorithms are BFGS, LBFGS, and simulated annealing. Examples are given.

Estimation of dynamic latent variable models using simulated nonparametric moments

Given a model that can be simulated, conditional moments at a trial parameter value can be calculated with high accuracy by applying kernel smoothing methods to a long simulation. With such conditional moments in hand, standard method of moments techniques can be used to estimate the parameter. Since conditional moments are calculated using kernel smoothing rather than simple averaging, it is not necessary that the model be simulable subject to the conditioning information that is used to define the moment conditions. For this reason, the proposed estimator is applicable to general dynamic latent variable models. Monte Carlo results show that the estimator performs well in comparison to other estimators that have been proposed for estimation of general DLV models.

Improving the risk concept: a revision of Arrow-Pratt theory in the context of controlled dynamic stochastic environments

In the literature on risk, one generally assume that uncertainty is uniformly distributed over the entire working horizon, when the absolute risk-aversion index is negative and constant. From this perspective, the risk is totally exogenous, and thus independent of endogenous risks. The classic procedure is "myopic" with regard to potential changes in the future behavior of the agent due to inherent random fluctuations of the system. The agent's attitude to risk is rigid. Although often criticized, the most widely used hypothesis for the analysis of economic behavior is risk-neutrality. This borderline case must be envisaged with prudence in a dynamic stochastic context. The traditional measures of risk-aversion are generally too weak for making comparisons between risky situations, given the dynamic �complexity of the environment. This can be highlighted in concrete problems in finance and insurance, context for which the Arrow-Pratt measures (in the small) give ambiguous.

Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

Improving the effort concept: a revision of the traditional approach in the context of controlled Dynamic Stochastic Environments

The objective of this paper is to re-evaluate the attitude to effort of a risk-averse decision-maker in an evolving environment. In the classic analysis, the space of efforts is generally discretized. More realistic, this new approach emploies a continuum of effort levels. The presence of multiple possible efforts and performance levels provides a better basis for explaining real economic phenomena. The traditional approach (see, Laffont, J. J. & Tirole, J., 1993, Salanie, B., 1997, Laffont, J.J. and Martimort, D, 2002, among others) does not take into account the potential effect of the system dynamics on the agent's behavior to effort over time. In the context of a Principal-agent relationship, not only the incentives of the Principal can determine the private agent to allocate a good effort, but also the evolution of the dynamic system. The incentives can be ineffective when the environment does not incite the agent to invest a good effort. This explains why, some effici

Static and dynamic shoulder stabilizer adaptations in javelin throwers: A preliminary study

BACKGROUND: Adaptations to Internal (IR) and external (ER) rotator shoulder muscles improving overhead throwing kinematics could lead to muscular strength imbalances and be considered an intrinsic risk factor for shoulder injury, as well as modified shoulder range of motion (RoM). OBJECTIVE: To establish profiles of internal and external rotation RoM and isokinetic IR and ER strength in adolescent- and national-level javelin throwers. METHODS: Fourteen healthy subjects were included in this preliminary cross-sectional study, 7 javelin throwers (JTG) and 7 nonathletes (CG). Passive internal and external rotation RoM were measured at 90 degrees of shoulder abduction. Isokinetic strength of dominant and non-dominant IR and ER was evaluated during concentric (60, 120 and 240 degrees/s) and eccentric (60 degrees/s) contractions by Con-Trex (R) dynamometer with the subject in a seated position with 45 degrees of shoulder abduction in the scapular plane. RESULTS: We reported significantly lower internal rotation and significantly higher external rotation RoM in JTG than in CG. Concentric and eccentric IR and ER strength were significantly higher for the dominant shoulder side in JTG (P < 0.05), without significant differences in ER/IR ratios. CONCLUSIONS: The main finding of this preliminary study confirmed static and dynamic shoulder stabilizer adaptations due to javelin throw practice in a population of adolescent- and national-level javelin throwers.

RADIC II: a fault tolerant architecture with flexible dynamic redundancy

The demand for computational power has been leading the improvement of the High Performance Computing (HPC) area, generally represented by the use of distributed systems like clusters of computers running parallel applications. In this area, fault tolerance plays an important role in order to provide high availability isolating the application from the faults effects. Performance and availability form an undissociable binomial for some kind of applications. Therefore, the fault tolerant solutions must take into consideration these two constraints when it has been designed. In this dissertation, we present a few side-effects that some fault tolerant solutions may presents when recovering a failed process. These effects may causes degradation of the system, affecting mainly the overall performance and availability. We introduce RADIC-II, a fault tolerant architecture for message passing based on RADIC (Redundant Array of Distributed Independent Fault Tolerance Controllers) architecture. RADIC-II keeps as maximum as possible the RADIC features of transparency, decentralization, flexibility and scalability, incorporating a flexible dynamic redundancy feature, allowing to mitigate or to avoid some recovery side-effects.

Three essays on the psychology of investment and financial markets

Préface My thesis consists of three essays where I consider equilibrium asset prices and investment strategies when the market is likely to experience crashes and possibly sharp windfalls. Although each part is written as an independent and self contained article, the papers share a common behavioral approach in representing investors preferences regarding to extremal returns. Investors utility is defined over their relative performance rather than over their final wealth position, a method first proposed by Markowitz (1952b) and by Kahneman and Tversky (1979), that I extend to incorporate preferences over extremal outcomes. With the failure of the traditional expected utility models in reproducing the observed stylized features of financial markets, the Prospect theory of Kahneman and Tversky (1979) offered the first significant alternative to the expected utility paradigm by considering that people focus on gains and losses rather than on final positions. Under this setting, Barberis, Huang, and Santos (2000) and McQueen and Vorkink (2004) were able to build a representative agent optimization model which solution reproduced some of the observed risk premium and excess volatility. The research in behavioral finance is relatively new and its potential still to explore. The three essays composing my thesis propose to use and extend this setting to study investors behavior and investment strategies in a market where crashes and sharp windfalls are likely to occur. In the first paper, the preferences of a representative agent, relative to time varying positive and negative extremal thresholds are modelled and estimated. A new utility function that conciliates between expected utility maximization and tail-related performance measures is proposed. The model estimation shows that the representative agent preferences reveals a significant level of crash aversion and lottery-pursuit. Assuming a single risky asset economy the proposed specification is able to reproduce some of the distributional features exhibited by financial return series. The second part proposes and illustrates a preference-based asset allocation model taking into account investors crash aversion. Using the skewed t distribution, optimal allocations are characterized as a resulting tradeoff between the distribution four moments. The specification highlights the preference for odd moments and the aversion for even moments. Qualitatively, optimal portfolios are analyzed in terms of firm characteristics and in a setting that reflects real-time asset allocation, a systematic over-performance is obtained compared to the aggregate stock market. Finally, in my third article, dynamic option-based investment strategies are derived and illustrated for investors presenting downside loss aversion. The problem is solved in closed form when the stock market exhibits stochastic volatility and jumps. The specification of downside loss averse utility functions allows corresponding terminal wealth profiles to be expressed as options on the stochastic discount factor contingent on the loss aversion level. Therefore dynamic strategies reduce to the replicating portfolio using exchange traded and well selected options, and the risky stock.

Cell-matrix interactions in Schistosomal portal fibrosis: a dynamic event

In recent years, one of the most significant progress in the understanding of liver diseases was the demonstration that liver fibrosis is a dynamic process resulting from a balance between synthesis and degradation of several matrix components, collagen in particular. Thus, fibrosis has been found to be a very early event during liver diseases, be it of toxic, viral or parasitic origin, and to be spontaneously reversible, either partially or totally. In liver fibrosis cell matrix interactions are dependent on the existence of the many factors (sometimes acting in combination) which produce the same events at the cellular and molecular levels. These events are: (i) the recruitment of fiber-producing cells, (ii) their proliferation, (iii) the secretion of matrix constituents of the extracellular matrix, and (iv) the remodeling and degradation of the newly formed matrix. All these events represent, at least in principle, a target for a therapeutic intervention aimed at influencing the experimentally induced hepatic fibrosis. In this context, hepatosplenic schistosomiasis is of particular interest, being an immune cell-mediated granulomatous disease and a model of liver fibrosis allowing extensive studies in human and animals as well as providing original in vitro models.

Explaining High Economic Growth in Small Tourism Countries with a Dynamic General Equilibrium Model

This paper shows that tourism specialisation can help to explain the observed high growth rates of small countries. For this purpose, two models of growth and trade are constructed to represent the trade relations between two countries. One of the countries is large, rich, has an own source of sustained growth and produces a tradable capital good. The other is a small poor economy, which does not have an own engine of growth and produces tradable tourism services. The poor country exports tourism services to and imports capital goods from the rich economy. In one model tourism is a luxury good, while in the other the expenditure elasticity of tourism imports is unitary. Two main results are obtained. In the long run, the tourism country overcomes decreasing returns and permanently grows because its terms of trade continuously improve. Since the tourism sector is relatively less productive than the capital good sector, tourism services become relatively scarcer and hence more expensive than the capital good. Moreover, along the transition the growth rate of the tourism economy holds well above the one of the rich country for a long time. The growth rate differential between countries is particularly high when tourism is a luxury good. In this case, there is a faster increase in the tourism demand. As a result, investment of the small economy is boosted and its terms of trade highly improve.

A musculoskeletal shoulder model based on pseudo-inverse and null-space optimization.

The goal of the present work was assess the feasibility of using a pseudo-inverse and null-space optimization approach in the modeling of the shoulder biomechanics. The method was applied to a simplified musculoskeletal shoulder model. The mechanical system consisted in the arm, and the external forces were the arm weight, 6 scapulo-humeral muscles and the reaction at the glenohumeral joint, which was considered as a spherical joint. The muscle wrapping was considered around the humeral head assumed spherical. The dynamical equations were solved in a Lagrangian approach. The mathematical redundancy of the mechanical system was solved in two steps: a pseudo-inverse optimization to minimize the square of the muscle stress and a null-space optimization to restrict the muscle force to physiological limits. Several movements were simulated. The mathematical and numerical aspects of the constrained redundancy problem were efficiently solved by the proposed method. The prediction of muscle moment arms was consistent with cadaveric measurements and the joint reaction force was consistent with in vivo measurements. This preliminary work demonstrated that the developed algorithm has a great potential for more complex musculoskeletal modeling of the shoulder joint. In particular it could be further applied to a non-spherical joint model, allowing for the natural translation of the humeral head in the glenoid fossa.

Dynamic Option Adjusted Spread and the Value of Mortgage Backed Securities (Draft 1)

We extend a reduced form model for pricing pass-through mortgage backed securities (MBS) and provide a novel hedging tool for investors in this market. To calculate the price of an MBS, traders use what is known as option-adjusted spread (OAS). The resulting OAS value represents the required basis points adjustment to reference curve discounting rates needed to match an observed market price. The OAS suffers from some drawbacks. For example, it remains constant until the maturity of the bond (thirty years in mortgage-backed securities), and does not incorporate interest rate volatility. We suggest instead what we call dynamic option adjusted spread (DOAS). The latter allows investors in the mortgage market to account for both prepayment risk and changes of the yield curve.

Dynamic Option Adjusted Spread and the Value of Mortgage Backed Securities (Draft 2)

We extend a reduced form model for pricing pass-through mortgage backed securities (MBS) and provide a novel hedging tool for investors in this market. To calculate the price of an MBS, traders use what is known as option-adjusted spread (OAS). The resulting OAS value represents the required basis points adjustment to reference curve discounting rates needed to match an observed market price. The OAS suffers from some drawbacks. For example, it remains constant until the maturity of the bond (thirty years in mortgage-backed securities), and does not incorporate interest rate volatility. We suggest instead what we call dynamic option adjusted spread (DOAS), which allows investors in the mortgage market to account for both prepayment risk and changes of the yield curve.

IInflation and Globalisation: A Dynamic Factor Model with Stochastic Volatility

National inflation rates reflect domestic and international (regional and global) influences. The relative importance of these components remains a controversial empirical issue. We extend the literature on inflation co-movement by utilising a dynamic factor model with stochastic volatility to account for shifts in the variance of inflation and endogenously determined regional groupings. We find that most of inflation variability is explained by the country specific disturbance term. Nevertheless, the contribution of the global component in explaining industrialised countries’ inflation rates has increased over time.