Robust response transformations based on optimal prediction

Nonlinear regression problems can often be reduced to linearity by transforming the response variable (e.g., using the Box-Cox family of transformations). The classic estimates of the parameter defining the transformation as well as of the regression coefficients are based on the maximum likelihood criterion, assuming homoscedastic normal errors for the transformed response. These estimates are nonrobust in the presence of outliers and can be inconsistent when the errors are nonnormal or heteroscedastic. This article proposes new robust estimates that are consistent and asymptotically normal for any unimodal and homoscedastic error distribution. For this purpose, a robust version of conditional expectation is introduced for which the prediction mean squared error is replaced with an M scale. This concept is then used to develop a nonparametric criterion to estimate the transformation parameter as well as the regression coefficients. A finite sample estimate of this criterion based on a robust version of smearing is also proposed. Monte Carlo experiments show that the new estimates compare favorably with respect to the available competitors.

O Comércio Intra-sectorial entre Portugal e os Países da União Europeia

A presente investigação tem como objectivo analisar a evolução do comércio intra-sectorial nas relações bilaterais entre Portugal e os países da União Europeia durante o período de 1995-2008 para os sectores: 31 e 34 e o impacto sobre das principais variáveis explicativas: diferença do rendimento per capita, diferença das dotações relativas, valor mínimo do rendimento per capita, valor máximo do rendimento per capita, BORDER (adjacência/fronteira), distância geográfica e dimensão do mercado (das economias) sobre o comércio intra-sectorial. A intensidade e a contribuição do comércio intra-sectorial é medida em termos absolutos e relativos através do índice proposto por Grubel e Lloyd (1975); a desagregação do comércio intra-sectorial horizontal e vertical é efectuada utilizando o critério de Abd-EL-Rahman (1991) e Grenaway et al. (1994) aplicado ao índice de Grubel e Lloyd. Os resultados obtidos mostraram que o comércio intra-sectorial vertical representa maior peso no total do comércio intra-sectorial nas relações bilaterais entre Portugal e a União Europeia, para ambos os sectores tanto em termos absolutos como em termos relativos. Em termos dos resultados econométricos verificou-se que a variável, proxy diferença do consumo eléctrico per capita utilizada para medir as dotações relativas possui um efeito negativo sobre o comércio intra-sectorial. As variáveis dimensão do mercado e valor mínimo do rendimento per capita apresentam um impacto positivo sobre o comércio intra-sectorial.

Accounting information and the prediction of farm viability

Until recently farm management made little use of accounting and agriculture has been largely excluded from the scope of accounting standards. This article examines the current use of accounting in agriculture and points theneed to establish accounting standards for agriculture. Empirical evidence shows that accounting can make a significant contribution to agricultural management and farm viability and could also be important for other agents involved in agricultural decision making. Existing literature on failureprediction models and farm viability prediction studies provide the starting point for our research, in which two dichotomous logit models were applied to subsamples of viable and unviable farms in Catalonia, Spain. The firstmodel considered only non-financial variables, while the other also considered financial ones. When accounting variables were added to the model, a significant reduction in deviance was observed.

Intravitreal ranibizumab for age-related macular degeneration : optical coherence tomography guided retreatment intervals and their intra- and inter-individual variability

Purpose: To determine whether the need for retreatment after an initial phase of 3 monthly intravitreal injections of ranibizumab shows an intra-individual regular rhythm and to what degree it varies between different patients. Methods: Prospective study with 42 patients with exudative AMD, treatment naïve. Loading dose of 3 monthly doses of ranibizumab (0,5 mg), followed by a 12 months pro re nata (PRN) regimen according to early exudative signs on HD-OCT Cirrus, Zeiss. The follow-up visits were intensified (week 4, 5, 6, 7, 8, 10, 12, 14, 16, 20, etc after each injection) in order to detect recurrences early, and injection followed within 3 days in cases of subretinal fluid, cysts, or central thickness increase of>50microns. Intervals were calculated between injections for the 12 month follow-up with PRN treatment. Variability was expressed as standard deviation (SD). Results: Visual acuity (VA) improved from a mean ETDRS score of 61.6 (SD 10.8) at baseline to 68.0 (SD 10.2) at month 3 and to 74.7(SD 9.0) at month 12. The 15 patients who have already completed the study showed maintenance of the VA improvement. Central foveal thickness improved from a mean value of 366 microns (baseline) to 253 microns (month 3), well maintained thereafter. Mean number of injections was 8.8 (SD 3.5,range 0-12) per 12 months of follow-up (after 3 doses), with mean individual treatment-recurrence (TR) intervals ranging from 28->365 days (mean 58). Intraindividual variability of TR intervals (SD) was 7.1 days as a mean value (range 1.7¡V22.6). It ranged within 20% of the mean intra-individual interval for 30 (91%) and within 15% for 21 patients (64%). The first interval was within 1 week of the mean intra-individual interval in 64% and within 2 weeks in 89% of patients. Conclusions: The majority of AMD patients showed a relatively stable rhythm for PRN injections of ranibizumab after initial loading phase, associated with excellent functional/anatomical results. The initial interval last loading dose-first recurrence may have a predictive value for further need of treatment, potentially facilitating follow-up and patient care.

Improving generalized regression analysis for the spatial prediction of forest communities

Aim This study used data from temperate forest communities to assess: (1) five different stepwise selection methods with generalized additive models, (2) the effect of weighting absences to ensure a prevalence of 0.5, (3) the effect of limiting absences beyond the environmental envelope defined by presences, (4) four different methods for incorporating spatial autocorrelation, and (5) the effect of integrating an interaction factor defined by a regression tree on the residuals of an initial environmental model. Location State of Vaud, western Switzerland. Methods Generalized additive models (GAMs) were fitted using the grasp package (generalized regression analysis and spatial predictions, http://www.cscf.ch/grasp). Results Model selection based on cross-validation appeared to be the best compromise between model stability and performance (parsimony) among the five methods tested. Weighting absences returned models that perform better than models fitted with the original sample prevalence. This appeared to be mainly due to the impact of very low prevalence values on evaluation statistics. Removing zeroes beyond the range of presences on main environmental gradients changed the set of selected predictors, and potentially their response curve shape. Moreover, removing zeroes slightly improved model performance and stability when compared with the baseline model on the same data set. Incorporating a spatial trend predictor improved model performance and stability significantly. Even better models were obtained when including local spatial autocorrelation. A novel approach to include interactions proved to be an efficient way to account for interactions between all predictors at once. Main conclusions Models and spatial predictions of 18 forest communities were significantly improved by using either: (1) cross-validation as a model selection method, (2) weighted absences, (3) limited absences, (4) predictors accounting for spatial autocorrelation, or (5) a factor variable accounting for interactions between all predictors. The final choice of model strategy should depend on the nature of the available data and the specific study aims. Statistical evaluation is useful in searching for the best modelling practice. However, one should not neglect to consider the shapes and interpretability of response curves, as well as the resulting spatial predictions in the final assessment.

Prediction of postnatal outcome in fetuses with pulmonary hypoplasia by fetal-MRI: Comparison between lung volumetry and pulmonary signal intensity values

Purpose: Pulmonary hypoplasia is a determinant parameter for extra-uterine life. In the last years, MRI appears as a complement to US in order to evaluate the degree of pulmonary hypoplasia in foetuses with congenital anomalies, by using different methods - fetal lung volumetry (FLV), lung-to-liver signal intensity ratio (LLSIR)-. But until now, information about the correlation between the MRI prediction and the real postnatal outcome is limited. Methods and materials: We retrospectively reviewed the fetal MRI performed at our Institution in the last 8 years and selected the cases with suspicion of fetal pulmonary hypoplasia (n = 30). The pulmonary volumetry data of these foetuses were collected and the lung-to-liver signal intensity ratio (LLSIR) measures performed. These data were compared with those obtained from a control group of 25 foetuses considered as normal at MRI. The data of the study group were also correlated with the autopsy records or the post-natal clinical information of the patients. Results: As expected, the control group showed higher FLV and LLSIR values than the problem group at all gestational ages. Higher values of FLV and LLSIR were associated with a better post-natal outcome. Sensitivity, specificity, positive and negative predictive values and accuracy for the relative LLSIR and the relative FLV showed no significant differences. Conclusion: Our data show that not only the FLV but also the relative LLSIR inform about the degree of fetal lung development. This information may help to predict the fetal outcome and to evaluate the need for neonatal intensive care.

Influência da competição intra-específica entre fêmeas e da ausência de hospedeiro no parasitismo de Telenomus remus Nixon (Hymenoptera, Scelionidae) sobre ovos de Spodoptera frugiperda (J. E. Smith) (Lepidoptera, Noctuidae)

A competição intraespecífica e a ausência de hospedeiros podem comprometer o sucesso de programas de controle biológico. Assim, este trabalho objetivou avaliar o efeito de diferentes densidades de fêmeas (1, 3, 6 e 9 fêmeas/ 100 ovos) e da ausência de hospedeiro no parasitismo de Telenomus remus Nixon (Hymenoptera, Scelionidae) em ovos de Spodoptera frugiperda (J. E. Smith) (Lepidoptera, Noctuidae). Observou-se aumento no parasitismo de ovos de S. frugiperda por T. remus a partir da densidade 3 fêmeas/ 100 ovos tanto em condições de laboratório (18%) quanto em campo (32%), embora o parasitismo causado por densidades maiores que 1 fêmea/100 ovos não foi significativamente diferente. Todavia, isso não causou superparasitismo. Por outro lado, o número de fêmeas produzidas na progênie diminuiu (39%) com a maior densidade de fêmeas sob condições de laboratório. Contudo, não verificou-se diferença significativa sob condições de campo. Mesmo após 9 dias de ausência de hospedeiro não houve redução nos índices de parasitismo, permanecendo em torno de 95 ovos parasitados/fêmea/dia. Todavia, no décimo dia pode-se observar redução (40%) no número médio de ovos parasitados. Desse modo, a liberação de três fêmeas por 100 ovos é adequada para se obter elevados níveis de parasitismo.

The intra-generational redistributive effects of social security

This paper provides a quantitative evaluation of the intra--cohortredistributive elements of the United States social security system in thecontext of a computable general equilibrium model. I determine how thewell--being of individuals that differ across {\sl gender, race} and {\sl education}is affected by government social security policy. I find that females, whitesand non--college graduates stand less to gain (lose) from reductions(increases) in the size of social security than males, non--whites andcollege graduates, respectively. Differences in mortality risk and laborproductivity translate into differences in the magnitudes of capitalaccumulation and labor supply distortions, that are responsible for theobserved welfare difference between types. Results imply that the currentprogram is lifetime progressive across gender and education, yet lifetimeregressive across race.

Genome-Wide Prediction Methods in Highly Diverse and Heterozygous Species: Proof-of-Concept through Simulation in Grapevine.

Nowadays, genome-wide association studies (GWAS) and genomic selection (GS) methods which use genome-wide marker data for phenotype prediction are of much potential interest in plant breeding. However, to our knowledge, no studies have been performed yet on the predictive ability of these methods for structured traits when using training populations with high levels of genetic diversity. Such an example of a highly heterozygous, perennial species is grapevine. The present study compares the accuracy of models based on GWAS or GS alone, or in combination, for predicting simple or complex traits, linked or not with population structure. In order to explore the relevance of these methods in this context, we performed simulations using approx 90,000 SNPs on a population of 3,000 individuals structured into three groups and corresponding to published diversity grapevine data. To estimate the parameters of the prediction models, we defined four training populations of 1,000 individuals, corresponding to these three groups and a core collection. Finally, to estimate the accuracy of the models, we also simulated four breeding populations of 200 individuals. Although prediction accuracy was low when breeding populations were too distant from the training populations, high accuracy levels were obtained using the sole core-collection as training population. The highest prediction accuracy was obtained (up to 0.9) using the combined GWAS-GS model. We thus recommend using the combined prediction model and a core-collection as training population for grapevine breeding or for other important economic crops with the same characteristics.

A simple randomized algorithm for consistent sequential prediction of ergodic time series

We present a simple randomized procedure for the prediction of a binary sequence. The algorithm uses ideas from recent developments of the theory of the prediction of individual sequences. We show that if thesequence is a realization of a stationary and ergodic random process then the average number of mistakes converges, almost surely, to that of the optimum, given by the Bayes predictor.

Intra-individual stability of neuromotor tasks from 6 to 18 years: a longitudinal study.

This study investigates the intra-individual stability of the speed of several motor tasks and the intensity of associated movements in 256 children (131 girls, 125 boys) from the Zurich generational study using the Zurich neuromotor assessment battery (ZNA) over a 12-year period from the age of 6 to 18 years. The stability was assessed by correlograms of standard deviation scores calculated from age- and gender-adjusted normative values and compared with standing height and full scale intelligence quotient (IQ). While motor tasks of hand, finger and foot (HFT) and contralateral associated movements (CAM) exhibited a moderate stability (summary measure as correlation coefficients between two measurements made 4 years apart: .61 and .60), other tasks (dynamic balance, static balance and pegboard) were only weakly stable (.46, .47 and .49). IQ and height were more stable than neuromotor components (.72 and .86). We conclude that the moderately stable HFT and CAM may reflect "motor traits", while the stability of the pegboard and balance tasks is weaker because these skills are more experience related and state-dependent.

Inter- and intra-individual variation of oxytocin secretion in oestrous mares exposed to stallions

The neuro-peptide hormone oxytocin regulates several reproductive mechanisms in mammals, such as uterine contractions during parturition and milk ejection in the lactating mammary gland. Oxytocin may also influence behavior and behavioral strategies, e.g. pair bonding, social recognition, maternal behavior, trust building, or anxiety. Teasing oestrous mares by a stallion provokes the release of oxytocin. We therefore tested whether such elevated oxytocin levels reveal possible mate preferences as determined in typical preference tests.

Intra-puparial development of the females of Chrysomya albiceps (Wiedemann) (Diptera, Calliphoridae)

Intra-puparial development of the females of Chrysomya albiceps (Wiedemann) (Diptera, Calliphoridae). The chronology and morphological changes that take place during intra-puparial development of Chrysomya albiceps is described based on 254 specimens reared in the laboratory. Larvae were obtained from the eggs laid by a single female. The pre-pupae were separated according to the reduction of larval length and the degree of pigmentation and sclerotization of the cuticle. After pupation, 10 individuals were fixed in Carnoy's solution and preserved in 70% ethanol, 10 individuals were fixed every 3 hours up to complete the first 24 hours (n = 80), the remaining individuals were fixed every six hours up to the 90th hour (n = 110) when 54 females emerged. The pupae were immersed in 5% formic acid for 48 hours and maintained in 70% ethanol, and then dissected and analyzed. C. albiceps shows four intra-puparial stages, each of which were described and compared with those described for Musca domestica, Calliphora erythrocephala, Sarcophaga bullata, Cuterebra tenebrosa, Oestrus ovis and Dermatobia hominis. Four developmental stages may be described: (1) the larva-pupa apolysis, after three hours; (2) the criptocephalic pupa, after six hours; (3) the phanerocephalic pupa, after nine hours; (4) the pharate pupa, after nine hours. The pharate adult is completely formed after 81 hours.

Prediction of functional outcome 18 months after a first psychotic episode: a proton magnetic resonance spectroscopy study.

CONTEXT: Recent magnetic resonance imaging studies have attempted to relate volumetric brain measurements in early schizophrenia to clinical and functional outcome some years later. These studies have generally been negative, perhaps because gray and white matter volumes inaccurately assess the underlying dysfunction that might be predictive of outcome. OBJECTIVE: To investigate the predictive value of frontal and temporal spectroscopy measures for outcome in patients with first-episode psychoses. DESIGN: Left prefrontal cortex and left mediotemporal lobe voxels were assessed using proton magnetic resonance spectroscopy to provide the ratio of N-acetylaspartate (NAA) and choline-containing compounds to creatine and phosphocreatine (Cr) (NAA/Cr ratio). These data were used to predict outcome at 18 months after admission, as assessed by a systematic medical record audit. SETTING: Early psychosis clinic. PARTICIPANTS: Forty-six patients with first-episode psychosis. MAIN OUTCOME MEASURES: We used regression models that included age at imaging and duration of untreated psychosis to predict outcome scores on the Global Assessment of Functioning Scale, Clinical Global Impression scales, and Social and Occupational Functional Assessment Scale, as well as the number of admissions during the treatment period. We then further considered the contributions of premorbid function and baseline level of negative symptoms. RESULTS: The only spectroscopic predictor of outcome was the NAA/Cr ratio in the prefrontal cortex. Low scores on this variable were related to poorer outcome on all measures. In addition, the frontal NAA/Cr ratio explained 17% to 30% of the variance in outcome. CONCLUSIONS: Prefrontal neuronal dysfunction is an inconsistent feature of early psychosis; rather, it is an early marker of poor prognosis across the first years of illness. The extent to which this can be used to guide treatment and whether it predicts outcome some years after first presentation are questions for further research.

Egg laying site selection by a host plant specialist leaf miner moth at two intra-plant levels in the northern Chilean Atacama Desert

Egg laying site selection by a host plant specialist leaf miner moth at two intra-plant levels in the northern Chilean Atacama Desert. The spatial distribution of the immature stages of the leaf miner Angelabella tecomae Vargas & Parra, 2005 was determined at two intra-plant levels (shoot and leaflet) on the shrub Tecoma fulva fulva (Cav.) D. Don (Bignoniaceae) in the Azapa valley, northern Chilean Atacama Desert. An aggregated spatial pattern was detected for all the immature stages along the shoot, with an age dependent relative position: eggs and first instar larvae were clumped at apex; second, third and fourth instar larvae were mostly found at intermediate positions; meanwhile the spinning larva and pupa were clumped at basis. This pattern suggests that the females select new, actively growing leaflets for egg laying. At the leaflet level, the immature stages were found more frequently at underside. Furthermore, survivorship was higher for larvae from underside mines. All these results highlight the importance of an accurate selection of egg laying site in the life history of this highly specialized leaf miner. By contrast, eventual wrong choices in the egg laying site selection may be associated with diminished larval survivorship. The importance of the continuous availability of new plant tissue in this highly human modified arid environment is discussed in relation with the observed patterns.