On the Behaviour of Marginal and Conditional Akaike Information Criteria in Linear Mixed Models

In linear mixed models, model selection frequently includes the selection of random effects. Two versions of the Akaike information criterion (AIC) have been used, based either on the marginal or on the conditional distribution. We show that the marginal AIC is no longer an asymptotically unbiased estimator of the Akaike information, and in fact favours smaller models without random effects. For the conditional AIC, we show that ignoring estimation uncertainty in the random effects covariance matrix, as is common practice, induces a bias that leads to the selection of any random effect not predicted to be exactly zero. We derive an analytic representation of a corrected version of the conditional AIC, which avoids the high computational cost and imprecision of available numerical approximations. An implementation in an R package is provided. All theoretical results are illustrated in simulation studies, and their impact in practice is investigated in an analysis of childhood malnutrition in Zambia.

Applying Least Absolute Shrinkage Selection Operator and Akaike Information Criterion Analysis to Find the Best Multiple Linear Regression Models between Climate Indices and Components of Cow’s Milk

This study focuses on multiple linear regression models relating six climate indices (temperature humidity THI, environmental stress ESI, equivalent temperature index ETI, heat load HLI, modified HLI (HLI new), and respiratory rate predictor RRP) with three main components of cow’s milk (yield, fat, and protein) for cows in Iran. The least absolute shrinkage selection operator (LASSO) and the Akaike information criterion (AIC) techniques are applied to select the best model for milk predictands with the smallest number of climate predictors. Uncertainty estimation is employed by applying bootstrapping through resampling. Cross validation is used to avoid over-fitting. Climatic parameters are calculated from the NASA-MERRA global atmospheric reanalysis. Milk data for the months from April to September, 2002 to 2010 are used. The best linear regression models are found in spring between milk yield as the predictand and THI, ESI, ETI, HLI, and RRP as predictors with p-value < 0.001 and R2 (0.50, 0.49) respectively. In summer, milk yield with independent variables of THI, ETI, and ESI show the highest relation (p-value < 0.001) with R2 (0.69). For fat and protein the results are only marginal. This method is suggested for the impact studies of climate variability/change on agriculture and food science fields when short-time series or data with large uncertainty are available.

Inclusão da epistasia em modelo de avaliação genética de bovinos de corte compostos

Dados de bovinos compostos foram analisados para avaliar o efeito da epistasia nos modelos de avaliação genética. As características analisadas foram os pesos aos 205 (P205) e 390 dias (P390) e perímetro escrotal aos 390 dias (PE390). As análises foram realizadas pela metodologia de máxima verossimilhança considerando-se dois modelos: o modelo 1 incluiu como covariáveis os efeitos aditivos diretos e maternos, e os não aditivos das heterozigoses para os efeitos diretos e para o materno total, e o modelo 2 considerou também o efeito direto de epistasia. Para comparação dos modelos, foram utilizados o critério de informação de Akaike (AIC) e o critério de informação Bayesiano de Schwartz (BIC), e o teste de razão de verossimilhança. A inclusão da epistasia no modelo de avaliação genética pouco alterou as estimativas de componentes de (co)variâncias genéticas aditivas e, consequentemente, as herdabilidades. O teste de verossimilhança e o critério de Akaike sugeriram que o modelo 2, que inclui a epistasia, apresentou maior aderência aos dados para todas as características analisadas. O critério BIC indicou este modelo como o melhor apenas para P205. Para análise genética dessa população, o modelo que considerou o efeito de epistasia foi o mais adequado.

Random regression models to estimate test-day milk yield genetic parameters Holstein cows in Southeastern Brazil

A total of 152,145 weekly test-day milk yield records from 7317 first lactations of Holstein cows distributed in 93 herds in southeastern Brazil were analyzed. Test-day milk yields were classified into 44 weekly classes of DIM. The contemporary groups were defined as herd-year-week of test-day. The model included direct additive genetic, permanent environmental and residual effects as random and fixed effects of contemporary group and age of cow at calving as covariable, linear and quadratic effects. Mean trends were modeled by a cubic regression on orthogonal polynomials of DIM. Additive genetic and permanent environmental random effects were estimated by random regression on orthogonal Legendre polynomials. Residual variances were modeled using third to seventh-order variance functions or a step function with 1, 6,13,17 and 44 variance classes. Results from Akaike`s and Schwarz`s Bayesian information criterion suggested that a model considering a 7th-order Legendre polynomial for additive effect, a 12th-order polynomial for permanent environment effect and a step function with 6 classes for residual variances, fitted best. However, a parsimonious model, with a 6th-order Legendre polynomial for additive effects and a 7th-order polynomial for permanent environmental effects, yielded very similar genetic parameter estimates. (C) 2008 Elsevier B.V. All rights reserved.

Comparison of different nonlinear functions to describe Nelore cattle growth

This work aims to compare different nonlinear functions for describing the growth curves of Nelore females. The growth curve parameters, their (co) variance components, and environmental and genetic effects were estimated jointly through a Bayesian hierarchical model. In the first stage of the hierarchy, 4 nonlinear functions were compared: Brody, Von Bertalanffy, Gompertz, and logistic. The analyses were carried out using 3 different data sets to check goodness of fit while having animals with few records. Three different assumptions about SD of fitting errors were considered: constancy throughout the trajectory, linear increasing until 3 yr of age and constancy thereafter, and variation following the nonlinear function applied in the first stage of the hierarchy. Comparisons of the overall goodness of fit were based on Akaike information criterion, the Bayesian information criterion, and the deviance information criterion. Goodness of fit at different points of the growth curve was compared applying the Gelfand`s check function. The posterior means of adult BW ranged from 531.78 to 586.89 kg. Greater estimates of adult BW were observed when the fitting error variance was considered constant along the trajectory. The models were not suitable to describe the SD of fitting errors at the beginning of the growth curve. All functions provided less accurate predictions at the beginning of growth, and predictions were more accurate after 48 mo of age. The prediction of adult BW using nonlinear functions can be accurate when growth curve parameters and their (co) variance components are estimated jointly. The hierarchical model used in the present study can be applied to the prediction of mature BW in herds in which a portion of the animals are culled before adult age. Gompertz, Von Bertalanffy, and Brody functions were adequate to establish mean growth patterns and to predict the adult BW of Nelore females. The Brody model was more accurate in predicting the birth weight of these animals and presented the best overall goodness of fit.

Analysis of residence patterns of Sperm whales (Physeter macrocephalus) in Azores Islands using opportunistic data

27th Annual Conference of the European Cetacean Society. Setúbal, Portugal, 8-10 April 2013.

Gestão da cadeia de abastecimento de bens de consumo baseada em modelos de previsão lineares

Dissertação apresentada ao Instituto Politécnico do Porto para obtenção do Grau de Mestre em Logística

Performance of state space and ARIMA models for consumer retail sales forecasting

Forecasting future sales is one of the most important issues that is beyond all strategic and planning decisions in effective operations of retail businesses. For profitable retail businesses, accurate demand forecasting is crucial in organizing and planning production, purchasing, transportation and labor force. Retail sales series belong to a special type of time series that typically contain trend and seasonal patterns, presenting challenges in developing effective forecasting models. This work compares the forecasting performance of state space models and ARIMA models. The forecasting performance is demonstrated through a case study of retail sales of five different categories of women footwear: Boots, Booties, Flats, Sandals and Shoes. On both methodologies the model with the minimum value of Akaike's Information Criteria for the in-sample period was selected from all admissible models for further evaluation in the out-of-sample. Both one-step and multiple-step forecasts were produced. The results show that when an automatic algorithm the overall out-of-sample forecasting performance of state space and ARIMA models evaluated via RMSE, MAE and MAPE is quite similar on both one-step and multi-step forecasts. We also conclude that state space and ARIMA produce coverage probabilities that are close to the nominal rates for both one-step and multi-step forecasts.

Spatial analysis and investigation of fire events occurrences in the Valencian Community, Spain

Dissertation submitted in partial fulfillment of the requirements for the Degree of Master of Science in Geospatial Technologies.

Spectrum sensing through software defined radio

Dissertação apresentada para obtenção do Grau de Mestre em Engenharia Electrotécnica e de Computadores, pela Universidade Nova de Lisboa, Faculdade de Ciências e Tecnologia

Ireland’s understudied flatfish: reproduction, age and growth of the dab Limanda limanda (L.) in Irish coastal waters

In Ireland, although flatfish form a valuable fishery, little is known about the smallest, the dab Limanda limanda. In this study, a variety of parameters of reproductive development, including ovarian phase description, gonadosomatic index (GSI), hepatosomatic index (HSI), relative condition (Kn) and oocyte size were analysed to provide information on the dab’s reproductive cycle and spawning periods. Sampling were collected monthly over an 18-month period using bottom trawls of the Irish coastline. A six phase macroscopic guide was developed for both sexes of dab, and verified using histology. In comparisons of macroscopic and microscopic phases, there was high agreement in the proposed female guide (86%), with males demonstratively lower (62%). No significant bias was observed between the the two reproductive methods. When the male macroscopic guide was examined, misclassification was high in phase 5 and phase 5 (41%), with 96% of misclassification occurring in adjacent phases. The sampled population was primarily composed of females, with ratios of females to males 1:0.6, although the predominance of females was less noticeable during the reproductive season. Oocyte growth in dab follows asynchronous development, and spawn over a protracted period indicating a batch spawning strategy. Spawning occurred mainly in early spring, with total regeneration of gonads by May. The length at which 50% of the population was reproductively mature was identified as 14cm and 17cm, for male and female dab, respectively. Precision and bias in age determinations using whole otoliths to age dab was investigated using six age readers from various institutions. Low levels of precision were obtained (CV: 10-23%) inferring the need for an alternative methodology. Precision and bias was influence by the level of experience of the reader, with ageing error attributed to interpretative differences and difficulty in edge determination. Sectioned otolith age determinations were subsequently compared to whole otolith age determinations using two age readers experienced in dab ageing. Although increased precision was observed in whole otoliths from previous estimates (CV=0%, 0% APE), sectioned otoliths were used for growth models. This was based on multinominal logistic regression on age length keys developed using both ageing methods. Biological data (length and age) for both sexes was applied to four growth models, where the Akaike criterion and Multi model Inference indicated the logistic model as having the best fit to the collected data. In general, female dab attained a longer length then males, with growth rates significantly different between the two sexes. Length weight relationships between the two sexes were also significantly different.

Global literature of fiddler crabs, genus Uca (Decapoda, Ocypodidae): trends and future directions

Crabs of the genus Uca Leach, 1814 are characterized by having strong sexual dimorphism and a global distribution. Currently, 97 species have been described and analyzed under several aspects, including population ecology, physiology and ethology. However, there is no general summary of the information from the various literatures. The aim of this study is to perform a scientometric analysis of fiddler crab studies. For this we searched papers available in the Thomson ISI database that contained the words "Uca" OR "fiddler* crab*" between the years 1991 and 2007. For each paper, we researched and recorded the following characteristics: publication year; journal of publication; the first author's nationality; the country where the study was conducted; study type; species studied; and the work area. Our results indicated that there was no increase in the number of articles through the years considered. The Journal of Experimental Marine Biology and Ecology published most of the articles on Uca, indicating the importance of this group as a model for testing ecological hypotheses using experimental approaches. Our results also showed that United States had the highest number of authors and published studies on Uca, following the overall trend in dominance on scientific research. Furthermore, using models with three variables (per capita income, number of species of Uca and extent of coastal countries) we observed that, according to the Akaike Information Criterion, the per capita income was the most important correlate for the number of articles per country (both the author's country and country of study). Additionally, our results show that the species U. pugilator (distributed on the East Coast of the North American continent) was the species most singularly referenced in the papers considered. Moreover, our results indicate that most studies on Uca use a descriptive and local scale. The majority of papers in our literature search reflect studies in population biology, followed by behavioral and physiological characteristics.

Multimodel inference and multimodel averaging in empirical modeling of occupational exposure levels.

Empirical modeling of exposure levels has been popular for identifying exposure determinants in occupational hygiene. Traditional data-driven methods used to choose a model on which to base inferences have typically not accounted for the uncertainty linked to the process of selecting the final model. Several new approaches propose making statistical inferences from a set of plausible models rather than from a single model regarded as 'best'. This paper introduces the multimodel averaging approach described in the monograph by Burnham and Anderson. In their approach, a set of plausible models are defined a priori by taking into account the sample size and previous knowledge of variables influent on exposure levels. The Akaike information criterion is then calculated to evaluate the relative support of the data for each model, expressed as Akaike weight, to be interpreted as the probability of the model being the best approximating model given the model set. The model weights can then be used to rank models, quantify the evidence favoring one over another, perform multimodel prediction, estimate the relative influence of the potential predictors and estimate multimodel-averaged effects of determinants. The whole approach is illustrated with the analysis of a data set of 1500 volatile organic compound exposure levels collected by the Institute for work and health (Lausanne, Switzerland) over 20 years, each concentration having been divided by the relevant Swiss occupational exposure limit and log-transformed before analysis. Multimodel inference represents a promising procedure for modeling exposure levels that incorporates the notion that several models can be supported by the data and permits to evaluate to a certain extent model selection uncertainty, which is seldom mentioned in current practice.

She’s a femme fatale: low-density larval development produces good disease vectors

Two hypotheses for how conditions for larval mosquitoes affect vectorial capacity make opposite predictions about the relationship of adult size and frequency of infection with vector-borne pathogens. Competition among larvae produces small adult females. The competition-susceptibility hypothesis postulates that small females are more susceptible to infection and predicts frequency of infection should decrease with size. The competition-longevity hypothesis postulates that small females have lower longevity and lower probability of becoming competent to transmit the pathogen and thus predicts frequency of infection should increase with size. We tested these hypotheses for Aedes aegypti in Rio de Janeiro, Brazil, during a dengue outbreak. In the laboratory, longevity increases with size, then decreases at the largest sizes. For field-collected females, generalised linear mixed model comparisons showed that a model with a linear increase of frequency of dengue with size produced the best Akaike’s information criterion with a correction for small sample sizes (AICc). Consensus prediction of three competing models indicated that frequency of infection increases monotonically with female size, consistent with the competition-longevity hypothesis. Site frequency of infection was not significantly related to site mean size of females. Thus, our data indicate that uncrowded, low competition conditions for larvae produce the females that are most likely to be important vectors of dengue. More generally, ecological conditions, particularly crowding and intraspecific competition among larvae, are likely to affect vector-borne pathogen transmission in nature, in this case via effects on longevity of resulting adults. Heterogeneity among individual vectors in likelihood of infection is a generally important outcome of ecological conditions impacting vectors as larvae.

U.S. subprime financial crisis contagion on BRIC and European Union stock markets

ABSTRACTThe Copula Theory was used to analyze contagion among the BRIC (Brazil, Russia, India and China) and European Union stock markets with the U.S. Equity Market. The market indexes used for the period between January 01, 2005 and February 27, 2010 are: MXBRIC (BRIC), MXEU (European Union) and MXUS (United States). This article evaluated the adequacy of the main copulas found in the financial literature using log-likelihood, Akaike information and Bayesian information criteria. This article provides a groundbreaking study in the area of contagion due to the use of conditional copulas, allowing to calculate the correlation increase between indexes with non-parametric approach. The conditional Symmetrized Joe-Clayton copula was the one that fitted better to the considered pairs of returns. Results indicate evidence of contagion effect in both markets, European Union and BRIC members, with a 5% significance level. Furthermore, there is also evidence that the contagion of U.S. financial crisis was more pronounced in the European Union than in the BRIC markets, with a 5% significance level. Therefore, stock portfolios formed by equities from the BRIC countries were able to offer greater protection during the subprime crisis. The results are aligned with recent papers that present an increase in correlation between stock markets, especially in bear markets.