Estimating the parameters of stochastic volatility models using option price data

This article describes a maximum likelihood method for estimating the parameters of the standard square-root stochastic volatility model and a variant of the model that includes jumps in equity prices. The model is fitted to data on the S&P 500 Index and the prices of vanilla options written on the index, for the period 1990 to 2011. The method is able to estimate both the parameters of the physical measure (associated with the index) and the parameters of the risk-neutral measure (associated with the options), including the volatility and jump risk premia. The estimation is implemented using a particle filter whose efficacy is demonstrated under simulation. The computational load of this estimation method, which previously has been prohibitive, is managed by the effective use of parallel computing using graphics processing units (GPUs). The empirical results indicate that the parameters of the models are reliably estimated and consistent with values reported in previous work. In particular, both the volatility risk premium and the jump risk premium are found to be significant.

Population health planning and health promotion: Opportunities and challenges to improve social determinants

Complex social factors and health issues challenge equitable health outcomes for many people, in particular those living in marginalised communities. Primary health care promises solutions through population health and health promotion approaches to improve social conditions (determinants) affecting health with emphasis on change at systems levels. Yet short-term efficiency focus policy decisions without long-term planning can undermine the effectiveness of primary health care. The workshop goal is to explore opportunities and share ideas about population health planning in Primary Health Networks and other community health care settings, so as to draw out opportunities, challenges and forward thinking health planning and health promotion strategies.

Stochastic growth models for analyzing crustacean data

The contemporary methodology for growth models of organisms is based on continuous trajectories and thus it hinders us from modelling stepwise growth in crustacean populations. Growth models for fish are normally assumed to follow a continuous function, but a different type of model is needed for crustacean growth. Crustaceans must moult in order for them to grow. The growth of crustaceans is a discontinuous process due to the periodical shedding of the exoskeleton in moulting. The stepwise growth of crustaceans through the moulting process makes the growth estimation more complex. Stochastic approaches can be used to model discontinuous growth or what are commonly known as "jumps" (Figure 1). However, in stochastic growth model we need to ensure that the stochastic growth model results in only positive jumps. In view of this, we will introduce a subordinator that is a special case of a Levy process. A subordinator is a non-decreasing Levy process, that will assist in modelling crustacean growth for better understanding of the individual variability and stochasticity in moulting periods and increments. We develop the estimation methods for parameter estimation and illustrate them with the help of a dataset from laboratory experiments. The motivational dataset is from the ornate rock lobster, Panulirus ornatus, which can be found between Australia and Papua New Guinea. Due to the presence of sex effects on the growth (Munday et al., 2004), we estimate the growth parameters separately for each sex. Since all hard parts are shed too often, the exact age determination of a lobster can be challenging. However, the growth parameters for the aforementioned moult processes from tank data being able to estimate through: (i) inter-moult periods, and (ii) moult increment. We will attempt to derive a joint density, which is made up of two functions: one for moult increments and the other for time intervals between moults. We claim these functions are conditionally independent given pre-moult length and the inter-moult periods. The variables moult increments and inter-moult periods are said to be independent because of the Markov property or conditional probability. Hence, the parameters in each function can be estimated separately. Subsequently, we integrate both of the functions through a Monte Carlo method. We can therefore obtain a population mean for crustacean growth (e. g. red curve in Figure 1). [GRAPHICS]

Study of Pinna nobilis growth from inner record: How biased are posterior adductor muscle scars estimates?

Previous studies have shown that the external growth records of the posterior adductor muscle scar (PAMS) of the bivalve Pinna nobilis are incomplete and do not produce accurate age estimations. We have developed a new methodology to study age and growth using the inner record of the PAMS, which avoids the necessity of costly in situ shell measurements or isotopic studies. Using the inner record we identified the positions of PAMS previously obscured by nacre and estimated the number of missing records in adult specimens with strong abrasion of the calcite layer in the anterior portion of the shell. The study of the PAMS and inner record of two shells that were 6 years old when collected showed that only 2 and 3 PAMS were observed, while 6 inner records could be counted, thus confirming our working methodology. Growth parameters of a P. nobilis population located in Moraira, Spain (western Mediterranean) were estimated with the new methodology and compared to those obtained using PAMS data and in situ measurements. For the comparisons, we applied different models considering the data alternatively as length-at-age (LA) and tag-recapture (TR). Among every method we tested to fit the Von Bertalanffy growth model, we observed that LA data from inner record fitted to the model using non-linear mixed effects and the estimation of missing records using the calcite width was the most appropriate. The equation obtained with this method, L = 573*(1 - e(-0.16(t-0.02))), is very similar to that calculated previously from in situ measurements for the same population.

A simple method for estimating growth parameters from multiple length-frequency data in presence of continuous recruitment

The extended recruitment season for short-lived species such as prawns biases the estimation of growth parameters from length-frequency data when conventional methods are used. We propose a simple method for overcoming this bias given a time series of length-frequency data. The difficulties arising from extended recruitment are eliminated by predicting the growth of the succeeding samples and the length increments of the recruits in previous samples. This method requires that some maximum size at recruitment can be specified. The advantages of this multiple length-frequency method are: it is simple to use; it requires only three parameters; no specific distributions need to be assumed; and the actual seasonal recruitment pattern does not have to be specified. We illustrate the new method with length-frequency data on the tiger prawn Penaeus esculentus from the north-western Gulf of Carpentaria, Australia.

An extravariation model for improving confidence intervals of population size estimates from removal data

We propose a new model for estimating the size of a population from successive catches taken during a removal experiment. The data from these experiments often have excessive variation, known as overdispersion, as compared with that predicted by the multinomial model. The new model allows catchability to vary randomly among samplings, which accounts for overdispersion. When the catchability is assumed to have a beta distribution, the likelihood function, which is refered to as beta-multinomial, is derived, and hence the maximum likelihood estimates can be evaluated. Simulations show that in the presence of extravariation in the data, the confidence intervals have been substantially underestimated in previous models (Leslie-DeLury, Moran) and that the new model provides more reliable confidence intervals. The performance of these methods was also demonstrated using two real data sets: one with overdispersion, from smallmouth bass (Micropterus dolomieu), and the other without overdispersion, from rat (Rattus rattus).

Growth curves with time-dependent explanatory variables

A class of growth models incorporating time-dependent factors and stochastic perturbations are introduced. The proposed model includes the existing growth models used in fisheries as special cases. Particular attention is given to growth of a population (in average weight or length) from which observations are taken randomly each time and the analysis of tag-recapture data. Two real data sets are used for illustration: (a) to estimate the seasonal effect and population density effect on growth of farmed prawn (Penaeus monodon) from weight data and (b) to assess the effect of tagging on growth of barramundi (Lates calcarifer)

Maximum likelihood estimation of mortality and growth with individual variability from multiple length-frequency data

We consider estimation of mortality rates and growth parameters from length-frequency data of a fish stock and derive the underlying length distribution of the population and the catch when there is individual variability in the von Bertalanffy growth parameter L-infinity. The model is flexible enough to accommodate 1) any recruitment pattern as a function of both time and length, 2) length-specific selectivity, and 3) varying fishing effort over time. The maximum likelihood method gives consistent estimates, provided the underlying distribution for individual variation in growth is correctly specified. Simulation results indicate that our method is reasonably robust to violations in the assumptions. The method is applied to tiger prawn data (Penaeus semisulcatus) to obtain estimates of natural and fishing mortality.

Population structure, mortality and growth of Pinna nobilis Linnaeus, 1758 (Mollusca, Bivalvia) at different depths in Moraira bay (Alicante, Western Mediterranean)

An investigation to characterize the causes of Pinna nobilis population structure in Moraira bay (Western Mediterranean) was developed. Individuals of two areas of the same Posidonia meadow, located at different depths (A1, -13 and A2, -6 m), were inventoried, tagged, their positions accurately recorded and monitored from July 1997 to July 2002. On each area, different aspects of population demography were studied (i.e. spatial distribution, size structure, displacement evidences, mortality, growth and shell orientation). A comparison between both groups of individuals was carried out, finding important differences between them. In A1, the individuals were more aggregated and mean and maximum size were higher (A1, 10.3 and A2, 6 individuals/100 m(2); A1, x = 47.2 +/- 9.9; A2, x = 29.8 +/- 7.4 cm, P < 0.001, respectively). In A2, growth rate and mortality were higher, the latter concentrated on the largest individuals, in contrast to A1, where the smallest individuals had the higher mortality rate [A1, L = 56.03(1 - e(-0.17t)); A2, L = 37.59(1 - e(-0.40t)), P < 0.001; mean annual mortality A1: 32 dead individuals out of 135, 23.7% and A2: 16 dead individuals out of 36, 44.4%, and total mortality coefficients (z), z(A1(-30)) = 0.28, z(A1(31-45)) = 0.05, z(A1(46-)) = 0.08; z(A2(-30)) = 0.15, z(A2(31-45)) = 0.25]. A common shell orientation N-S, coincident with the maximum shore exposure, was observed in A2. Spatial distribution in both areas showed not enough evidence to discard a random distribution of the individuals, despite the greater aggregation on the deeper area (A1) (A1, chi(2) = 0.41, df = 3, P > 0.5, A2, chi(2)= 0.98, df = 2 and 0.3 < P < 0.5). The obtained results have demonstrated that the depth-related size segregation usually shown by P. nobilis is mainly caused by differences in mortality and growth among individuals located at different depths, rather than by the active displacement of individuals previously reported in the literature. Furthermore, dwarf individuals are observed in shallower levels and as a consequence, the relationship between size and age are not comparable even among groups of individuals inhabiting the same meadow at different depths. The final causes of the differences on mortality and growth are also discussed.

Estimation of growth parameters from multiple-recapture data

This article develops a method for analysis of growth data with multiple recaptures when the initial ages for all individuals are unknown. The existing approaches either impute the initial ages or model them as random effects. Assumptions about the initial age are not verifiable because all the initial ages are unknown. We present an alternative approach that treats all the lengths including the length at first capture as correlated repeated measures for each individual. Optimal estimating equations are developed using the generalized estimating equations approach that only requires the first two moment assumptions. Explicit expressions for estimation of both mean growth parameters and variance components are given to minimize the computational complexity. Simulation studies indicate that the proposed method works well. Two real data sets are analyzed for illustration, one from whelks (Dicathais aegaota) and the other from southern rock lobster (Jasus edwardsii) in South Australia.

Estimating equations for parameters in stochastic growth models from tag-recapture data

James (1991, Biometrics 47, 1519-1530) constructed unbiased estimating functions for estimating the two parameters in the von Bertalanffy growth curve from tag-recapture data. This paper provides unbiased estimating functions for a class of growth models that incorporate stochastic components and explanatory variables. a simulation study using seasonal growth models indicates that the proposed method works well while the least-squares methods that are commonly used in the literature may produce substantially biased estimates. The proposed model and method are also applied to real data from tagged rack lobsters to assess the possible seasonal effect on growth.

An improved Fabens method for estimation of growth parameters in the vonBertalanffy model with individual asymptotes

In the analysis of tagging data, it has been found that the least-squares method, based on the increment function known as the Fabens method, produces biased estimates because individual variability in growth is not allowed for. This paper modifies the Fabens method to account for individual variability in the length asymptote. Significance tests using t-statistics or log-likelihood ratio statistics may be applied to show the level of individual variability. Simulation results indicate that the modified method reduces the biases in the estimates to negligible proportions. Tagging data from tiger prawns (Penaeus esculentus and Penaeus semisulcatus) and rock lobster (Panulirus ornatus) are analysed as an illustration.

Variable Size Population NSGA-II: VPNSGA-II

Multi-objective optimization is an active field of research with broad applicability in aeronautics. This report details a variant of the original NSGA-II software aimed to improve the performances of such a widely used Genetic Algorithm in finding the optimal Pareto-front of a Multi-Objective optimization problem for the use of UAV and aircraft design and optimsaiton. Original NSGA-II works on a population of predetermined constant size and its computational cost to evaluate one generation is O(mn^2 ), being m the number of objective functions and n the population size. The basic idea encouraging this work is that of reduce the computational cost of the NSGA-II algorithm by making it work on a population of variable size, in order to obtain better convergence towards the Pareto-front in less time. In this work some test functions will be tested with both original NSGA-II and VPNSGA-II algorithms; each test will be timed in order to get a measure of the computational cost of each trial and the results will be compared.

Direct shear testing on unsaturated silty soils to investigate the effects of drying and wetting on shear strength parameters at low suction

A modified conventional direct shear device was used to measure unsaturated shear strength of two silty soils at low suction values (0 ~ 50 kPa) that were achieved by following drying and wetting paths of soil water characteristic curves (SWCCs). The results revealed that the internal friction angle of the soils was not significantly affected by either the suction or the drying wetting SWCCs. The apparent cohesion of soil increased with a decreasing rate as suction increased. Shear stress-shear displacement curves obtained from soil specimens subjected to the same net normal stress and different suction values showed a higher initial stiffness and a greater peak stress as suction increased. A soil in wetting exhibited slightly higher peak shear stress and more contractive volume change behavior than that of soil in drying at the same net normal stress and suction.