Linear least squares estimation of the first order moving average parameter

We propose an iterative procedure to minimize the sum of squares function which avoids the nonlinear nature of estimating the first order moving average parameter and provides a closed form of the estimator. The asymptotic properties of the method are discussed and the consistency of the linear least squares estimator is proved for the invertible case. We perform various Monte Carlo experiments in order to compare the sample properties of the linear least squares estimator with its nonlinear counterpart for the conditional and unconditional cases. Some examples are also discussed

Fecal pellets collection as a method for assessing egesta of the marine cave-dwelling mysid Hemimysis speluncula

Egesta of a cave-dwelling mysid (Hemimysis speluncola Ledoyer, 1963) was studied in a submarine cave of Medes Islands, NW Mediterranean by in situ fecal pellet collecting. Fecal pellet production and gut fullness of mysids during incubation experiments are used to estimate mysid egestion rates. Intrinsic factors related with the natural history of this species such as population structure, density of mysids, daily rhythms and pellet decomposition rates are tested for their influence on the egestion rate. The effects of methodological artifacts, such as the stress induced by both incubation and preservation procedures, are also studied. An average mysid egests about 2.5 pellets per day into the cave. The time of day is the main factor affecting egestion. The highest deposition rate is between 2 to 4 hours after sunrise when about 38 % of the total daily pellet production becomes egested. Fecal pellet morphology changes with mysid demographic classes: immature mysids produce slender and thick pellets, whereas mature mysids produce only thick pellets. Immature classes show higher percentages of full guts than mature ones. Mysid density in the incubators does not affect the results on gut fullness, but it causes a decrease in the number of pellets collected after incubation. Coprorhexia seems to be the only plausible process to explain this paradox. The incubation procedure does not increase deposition rate significantly. Time of incubation is critical because the half-life of fecal pellets is about 2.5 hours. Fixation with liquid nitrogen decreases gut fullness and also deposition rates. Higher values are obtained with 70 % ethanol and 5 % formalin solutions which show very similar results for both gut fullness and pellet deposition rates. Nevertheless, ethanol is not suitable as fixative because it enhances the opacity of the body. Several suggestions are given in order to optimize the reliability of further in situ experiments for evaluation of egesta of Hemimysis speluncola in submarine caves.

Geomorphological method in the elaboration of hazard maps for flash-floods in the municipality of Jucuarán (El Salvador)

This work deals with the elaboration of flood hazard maps. These maps reflect the areas prone to floods based on the effects of Hurricane Mitch in the Municipality of Jucuarán of El Salvador. Stream channels located in the coastal range in the SE of El Salvador flow into the Pacific Ocean and generate alluvial fans. Communities often inhabit these fans can be affected by floods. The geomorphology of these stream basins is associated with small areas, steep slopes, well developed regolite and extensive deforestation. These features play a key role in the generation of flash-floods. This zone lacks comprehensive rainfall data and gauging stations. The most detailed topographic maps are on a scale of 1:25 000. Given that the scale was not sufficiently detailed, we used aerial photographs enlarged to the scale of 1:8000. The effects of Hurricane Mitch mapped on these photographs were regarded as the reference event. Flood maps have a dual purpose (1) community emergency plans, (2) regional land use planning carried out by local authorities. The geomorphological method is based on mapping the geomorphological evidence (alluvial fans, preferential stream channels, erosion and sedimentation, man-made terraces). Following the interpretation of the photographs this information was validated on the field and complemented by eyewitness reports such as the height of water and flow typology. In addition, community workshops were organized to obtain information about the evolution and the impact of the phenomena. The superimposition of this information enables us to obtain a comprehensive geomorphological map. Another aim of the study was the calculation of the peak discharge using the Manning and the paleohydraulic methods and estimates based on geomorphologic criterion. The results were compared with those obtained using the rational method. Significant differences in the order of magnitude of the calculated discharges were noted. The rational method underestimated the results owing to short and discontinuous periods of rainfall data with the result that probabilistic equations cannot be applied. The Manning method yields a wide range of results because of its dependence on the roughness coefficient. The paleohydraulic method yielded higher values than the rational and Manning methods. However, it should be pointed out that it is possible that bigger boulders could have been moved had they existed. These discharge values are lower than those obtained by the geomorphological estimates, i.e. much closer to reality. The flood hazard maps were derived from the comprehensive geomorphological map. Three categories of hazard were established (very high, high and moderate) using flood energy, water height and velocity flow deduced from geomorphological and eyewitness reports.

Rapid Biodiversity Assessment and Monitoring Method for Highly Diverse Benthic Communities: A Case Study of Mediterranean Coralligenous Outcrops

Increasing anthropogenic pressures urge enhanced knowledge and understanding of the current state of marine biodiversity. This baseline information is pivotal to explore present trends, detect future modifications and propose adequate management actions for marine ecosystems. Coralligenous outcrops are a highly diverse and structurally complex deep-water habitat faced with major threats in the Mediterranean Sea. Despite its ecological, aesthetic and economic value, coralligenous biodiversity patterns are still poorly understood. There is currently no single sampling method that has been demonstrated to be sufficiently representative to ensure adequate community assessment and monitoring in this habitat. Therefore, we propose a rapid non-destructive protocol for biodiversity assessment and monitoring of coralligenous outcrops providing good estimates of its structure and species composition, based on photographic sampling and the determination of presence/absence of macrobenthic species. We used an extensive photographic survey, covering several spatial scales (100s of m to 100s of km) within the NW Mediterranean and including 2 different coralligenous assemblages: Paramuricea clavata (PCA) and Corallium rubrum assemblage (CRA). This approach allowed us to determine the minimal sampling area for each assemblage (5000 cm² for PCA and 2500 cm²for CRA). In addition, we conclude that 3 replicates provide an optimal sampling effort in order to maximize the species number and to assess the main biodiversity patterns of studied assemblages in variability studies requiring replicates. We contend that the proposed sampling approach provides a valuable tool for management and conservation planning, monitoring and research programs focused on coralligenous outcrops, potentially also applicable in other benthic ecosystems

Identifying quantitative operation principles in metabolic pathways: a systematic method for searching feasible enzyme activity patterns leading to cellular adaptive responses

Background: Optimization methods allow designing changes in a system so that specific goals are attained. These techniques are fundamental for metabolic engineering. However, they are not directly applicable for investigating the evolution of metabolic adaptation to environmental changes. Although biological systems have evolved by natural selection and result in well-adapted systems, we can hardly expect that actual metabolic processes are at the theoretical optimum that could result from an optimization analysis. More likely, natural systems are to be found in a feasible region compatible with global physiological requirements. Results: We first present a new method for globally optimizing nonlinear models of metabolic pathways that are based on the Generalized Mass Action (GMA) representation. The optimization task is posed as a nonconvex nonlinear programming (NLP) problem that is solved by an outer- approximation algorithm. This method relies on solving iteratively reduced NLP slave subproblems and mixed-integer linear programming (MILP) master problems that provide valid upper and lower bounds, respectively, on the global solution to the original NLP. The capabilities of this method are illustrated through its application to the anaerobic fermentation pathway in Saccharomyces cerevisiae. We next introduce a method to identify the feasibility parametric regions that allow a system to meet a set of physiological constraints that can be represented in mathematical terms through algebraic equations. This technique is based on applying the outer-approximation based algorithm iteratively over a reduced search space in order to identify regions that contain feasible solutions to the problem and discard others in which no feasible solution exists. As an example, we characterize the feasible enzyme activity changes that are compatible with an appropriate adaptive response of yeast Saccharomyces cerevisiae to heat shock Conclusion: Our results show the utility of the suggested approach for investigating the evolution of adaptive responses to environmental changes. The proposed method can be used in other important applications such as the evaluation of parameter changes that are compatible with health and disease states.

Hollywood Films in the Non-Western World: What Are the Criteria Followed by the Chinese Government When Choosing Hollywood Film Imports?

This dissertation argues that the Government of the People’s Republic of China, when it made the decision to import a quota of Hollywood films in 1994 to revive the failing domestic film industry, had different possible criteria in mind. This project has studied four of them: first, importing films that gave a negative image of the United States; second, importing films that featured Chinese talent or themes; third, importing films that were box office hits in the United States; and fourth, importing films with a strong technological innovation ingredient. In order to find out the most important criteria for the Chinese Government, this dissertation offers a dataset that analyzes a population of 262 Hollywood films released in the PRC between 1994 and 2010. For each unit, a method has been developed to compile data that will determine whether the film reflects any of the four criteria, and findings in the form of yearly percentages have been drawn. Results show that, out of the four studied criteria, the two predominant reasons for China to import films were technological innovation and box office hits. This tells us that, at this point, the Chinese Government shows more interest in obtaining big revenues and learning from technically innovative American films than in delivering underlying political messages. This dissertation contributes to the existing literature by analyzing the content of all the films imported by China between 1994 and 2010, while integrating in the analysis variables based on the existing knowledge.

New approximation to the third-order density : application to the calculation of correlated multicenter indices

In this work we present the formulas for the calculation of exact three-center electron sharing indices (3c-ESI) and introduce two new approximate expressions for correlated wave functions. The 3c-ESI uses the third-order density, the diagonal of the third-order reduced density matrix, but the approximations suggested in this work only involve natural orbitals and occupancies. In addition, the first calculations of 3c-ESI using Valdemoro's, Nakatsuji's and Mazziotti's approximation for the third-order reduced density matrix are also presented for comparison. Our results on a test set of molecules, including 32 3c-ESI values, prove that the new approximation based on the cubic root of natural occupancies performs the best, yielding absolute errors below 0.07 and an average absolute error of 0.015. Furthemore, this approximation seems to be rather insensitive to the amount of electron correlation present in the system. This newly developed methodology provides a computational inexpensive method to calculate 3c-ESI from correlated wave functions and opens new avenues to approximate high-order reduced density matrices in other contexts, such as the contracted Schrödinger equation and the anti-Hermitian contracted Schrödinger equation