Fourth order pseudo maximum likelihood methods

We extend PML theory to account for information on the conditional moments up to order four, but without assuming a parametric model, to avoid a risk of misspecification of the conditional distribution. The key statistical tool is the quartic exponential family, which allows us to generalize the PML2 and QGPML1 methods proposed in Gourieroux et al. (1984) to PML4 and QGPML2 methods, respectively. An asymptotic theory is developed. The key numerical tool that we use is the Gauss-Freud integration scheme that solves a computational problem that has previously been raised in several fields. Simulation exercises demonstrate the feasibility and robustness of the methods [Authors]

Effects of the peroxisome proliferator-activated receptor (PPAR)-gamma agonist pioglitazone on renal and hormonal responses to salt in diabetic and hypertensive individuals.

Aims/Hypothesis: Glitazones are powerful insulin sensitisers prescribed for the treatment of type 2 diabetes. Their use is, however, associated with fluid retention and an increased risk of congestive heart failure. We previously demonstrated that pioglitazone increases proximal sodium reabsorption in healthy volunteers. This study examines the effects of pioglitazone on renal sodium handling in individuals prone to insulin resistance, i.e. those with diabetes and/or hypertension. Methods: In this double-blind randomised placebo-controlled four-way crossover study, we examined the effects of pioglitazone (45 mg daily during 6 weeks) or placebo on renal, systemic and hormonal responses to changes in sodium intake in 16 individuals, eight with type 2 diabetes and eight with hypertension. Results: Pioglitazone was associated with a rapid increase in body weight and an increase in diurnal proximal sodium reabsorption, without any change in renal haemodynamics or in the modulation of the renin-angiotensin aldosterone system to changes in salt intake. A compensatory increase in brain natriuretic peptide levels was observed. In spite of sodium retention, pioglitazone dissociated the blood-pressure response to salt and abolished salt sensitivity in salt-sensitive individuals. Conclusions/Interpretation: Pioglitazone increases diurnal proximal sodium retention in diabetic and hypertensive individuals. These effects cause fluid retention and may contribute to the increased incidence of congestive heart failure with glitazones.

Cytostatic lung perfusion results in heterogeneous spatial regional blood flow and drug distribution: evaluation of different cytostatic lung perfusion techniques in a porcine model.

OBJECTIVES: Comparison of doxorubicin uptake, leakage and spatial regional blood flow, and drug distribution was made for antegrade, retrograde, combined antegrade and retrograde isolated lung perfusion, and pulmonary artery infusion by endovascular inflow occlusion (blood flow occlusion), as opposed to intravenous administration in a porcine model. METHODS: White pigs underwent single-pass lung perfusion with doxorubicin (320 mug/mL), labeled 99mTc-microspheres, and Indian ink. Visual assessment of the ink distribution and perfusion scintigraphy of the perfused lung was performed. 99mTc activity and doxorubicin levels were measured by gamma counting and high-performance liquid chromatography on 15 tissue samples from each perfused lung at predetermined localizations. RESULTS: Overall doxorubicin uptake in the perfused lung was significantly higher (P = .001) and the plasma concentration was significantly lower (P < .0001) after all isolated lung perfusion techniques, compared with intravenous administration, without differences between them. Pulmonary artery infusion (blood flow occlusion) showed an equally high doxorubicin uptake in the perfused lung but a higher systemic leakage than surgical isolated lung perfusion (P < .0001). The geometric coefficients of variation of the doxorubicin lung tissue levels were 175%, 279%, 226%, and 151% for antegrade, retrograde, combined antegrade and retrograde isolated lung perfusion, and pulmonary artery infusion by endovascular inflow occlusion (blood flow occlusion), respectively, compared with 51% for intravenous administration (P = .09). 99mTc activity measurements of the samples paralleled the doxorubicin level measurements, indicating a trend to a more heterogeneous spatial regional blood flow and drug distribution after isolated lung perfusion and blood flow occlusion compared with intravenous administration. CONCLUSIONS: Cytostatic lung perfusion results in a high overall doxorubicin uptake, which is, however, heterogeneously distributed within the perfused lung.

Patterns of pain freedom response in 497 cases with classical trigeminal neuralgia and at least one year of follow-up treated by gamma knife radiosurgery

Purpose: To identify patterns of initially pain freedom response in patients with classical trigeminal neuralgia (CTN) with Gamma Knife surgery (GKS) and to compare their associated hypoesthesia and recurrence rates. Methods: In this study we analysed only 497 patients treated between July 1992 and November 2010, with a follow-up longer than 1 year, after excluding megadolichobasilar artery and multiple sclerosis related trigeminal neuralgia, as well as the second GKS treatments so as to have only cases with CTN and single radiosurgical treatment. GKS using a Gamma Knife (model B or C or Perfexion) was performed, based on both MRI and computer tomography (CT) targeting. A single 4-mm isocenter was positioned in the cisternal portion of the trigeminal nerve at a median distance of 7.8 mm (range 4.5-14) anteriorly to the emergence of the nerve. A median maximum dose of 85 Gy (range 70-90) was delivered. After empiric methods but also by using a chart with clear cut-off periods of pain free distribution, we were able to divide the initially pain free patients into 3 separate groups: within the first 48 hours, after 48 hours till 30 days and after 30 days, respectively. Results: The median follow- up period was 43.75 months (range 12-174.41). Four hundred and fifty-four patients (91.75%) were initially pain free in a median time of 10 days (range 1-459): 169 (37.2%) became pain free within the first 48 hours (pf<=48 h), compared to 194 (42.8%) between the 3-rd day and the day 30 (pf (>48 h, <=30 d)), inclusively and 91 (20%) patients after 30 days (pf>30d). Differences in terms of postoperative hypoesthesia were found with a p value of 0.014 as follows: the group pf<=48 h had 18 (13.7%) compared to pf (>48 h, <=30 d) with 30 (19%) and pf>30d with 22 (30.5%) patients developing a postoperative GKS hypoesthesia. One hundred and fifty seven (34.4%) patients initially pain free experienced a recurrence with a median delay of 24 months (range 0.62-150.06). There were no statistically significant differences between the three groups concerning recurrence (p value of 0.515). Conclusions: An important number of patients (169 cases, 37.2%) became pain free in the first 48 hours. Hypoesthesia rate was higher within the group becoming pain free after 30 days with a statistically significant difference between the three populations (p= 0.014). Further analysis will eventually help to elucidate the differences observed among groups.

Alternative ways to estimate change points in multinomial sequences. An application to an authorship attribution problem

The statistical analysis of literary style is the part of stylometry that compares measurable characteristicsin a text that are rarely controlled by the author, with those in other texts. When thegoal is to settle authorship questions, these characteristics should relate to the author’s style andnot to the genre, epoch or editor, and they should be such that their variation between authors islarger than the variation within comparable texts from the same author.For an overview of the literature on stylometry and some of the techniques involved, see for exampleMosteller and Wallace (1964, 82), Herdan (1964), Morton (1978), Holmes (1985), Oakes (1998) orLebart, Salem and Berry (1998).Tirant lo Blanc, a chivalry book, is the main work in catalan literature and it was hailed to be“the best book of its kind in the world” by Cervantes in Don Quixote. Considered by writterslike Vargas Llosa or Damaso Alonso to be the first modern novel in Europe, it has been translatedseveral times into Spanish, Italian and French, with modern English translations by Rosenthal(1996) and La Fontaine (1993). The main body of this book was written between 1460 and 1465,but it was not printed until 1490.There is an intense and long lasting debate around its authorship sprouting from its first edition,where its introduction states that the whole book is the work of Martorell (1413?-1468), while atthe end it is stated that the last one fourth of the book is by Galba (?-1490), after the death ofMartorell. Some of the authors that support the theory of single authorship are Riquer (1990),Chiner (1993) and Badia (1993), while some of those supporting the double authorship are Riquer(1947), Coromines (1956) and Ferrando (1995). For an overview of this debate, see Riquer (1990).Neither of the two candidate authors left any text comparable to the one under study, and thereforediscriminant analysis can not be used to help classify chapters by author. By using sample textsencompassing about ten percent of the book, and looking at word length and at the use of 44conjunctions, prepositions and articles, Ginebra and Cabos (1998) detect heterogeneities that mightindicate the existence of two authors. By analyzing the diversity of the vocabulary, Riba andGinebra (2000) estimates that stylistic boundary to be near chapter 383.Following the lead of the extensive literature, this paper looks into word length, the use of the mostfrequent words and into the use of vowels in each chapter of the book. Given that the featuresselected are categorical, that leads to three contingency tables of ordered rows and therefore tothree sequences of multinomial observations.Section 2 explores these sequences graphically, observing a clear shift in their distribution. Section 3describes the problem of the estimation of a suden change-point in those sequences, in the followingsections we propose various ways to estimate change-points in multinomial sequences; the methodin section 4 involves fitting models for polytomous data, the one in Section 5 fits gamma modelsonto the sequence of Chi-square distances between each row profiles and the average profile, theone in Section 6 fits models onto the sequence of values taken by the first component of thecorrespondence analysis as well as onto sequences of other summary measures like the averageword length. In Section 7 we fit models onto the marginal binomial sequences to identify thefeatures that distinguish the chapters before and after that boundary. Most methods rely heavilyon the use of generalized linear models

The distribution pattern of Lutzomyia longipalpis (Diptera: Psychodidae) in the peridomiciles of a sector with canine and human visceral leishmaniasis transmission in the municipality of Dracena, São Paulo, Brazil

The specimen distribution pattern of a species can be used to characterise a population of interest and also provides area-specific guidance for pest management and control. In the municipality of Dracena, in the state of São Paulo, we analysed 5,889 Lutzomyia longipalpis specimens collected from the peridomiciles of 14 houses in a sector where American visceral leishmaniasis (AVL) is transmitted to humans and dogs. The goal was to analyse the dispersion and a theoretical fitting of the species occurrence probability. From January-December 2005, samples were collected once per week using CDC light traps that operated for 12-h periods. Each collection was considered a sub-sample and was evaluated monthly. The standardised Morisita index was used as a measure of dispersion. Adherence tests were performed for the log-series distribution. The number of traps was used to adjust the octave plots. The quantity of Lu. longipalpis in the sector was highly aggregated for each month of the year, adhering to a log-series distribution for 11 of the 12 months analysed. A sex-stratified analysis demonstrated a pattern of aggregated dispersion adjusted for each month of the year. The classes and frequencies of the traps in octaves can be employed as indicators for entomological surveillance and AVL control.

Use of Two-Part Regression Calibration Model to Correct for Measurement Error in Episodically Consumed Foods in a Single-Replicate Study Design: EPIC Case Study.

In epidemiologic studies, measurement error in dietary variables often attenuates association between dietary intake and disease occurrence. To adjust for the attenuation caused by error in dietary intake, regression calibration is commonly used. To apply regression calibration, unbiased reference measurements are required. Short-term reference measurements for foods that are not consumed daily contain excess zeroes that pose challenges in the calibration model. We adapted two-part regression calibration model, initially developed for multiple replicates of reference measurements per individual to a single-replicate setting. We showed how to handle excess zero reference measurements by two-step modeling approach, how to explore heteroscedasticity in the consumed amount with variance-mean graph, how to explore nonlinearity with the generalized additive modeling (GAM) and the empirical logit approaches, and how to select covariates in the calibration model. The performance of two-part calibration model was compared with the one-part counterpart. We used vegetable intake and mortality data from European Prospective Investigation on Cancer and Nutrition (EPIC) study. In the EPIC, reference measurements were taken with 24-hour recalls. For each of the three vegetable subgroups assessed separately, correcting for error with an appropriately specified two-part calibration model resulted in about three fold increase in the strength of association with all-cause mortality, as measured by the log hazard ratio. Further found is that the standard way of including covariates in the calibration model can lead to over fitting the two-part calibration model. Moreover, the extent of adjusting for error is influenced by the number and forms of covariates in the calibration model. For episodically consumed foods, we advise researchers to pay special attention to response distribution, nonlinearity, and covariate inclusion in specifying the calibration model.

Generalized food-deceptive orchid species flower earlier and occur at lower altitudes than rewarding ones

Aims Food-deceptive pollination, in which plants do not offer any food reward to their pollinators, is common within the Orchidaceae. As food-deceptive orchids are poorer competitors for pollinator visitation than rewarding orchids, their occurrence in a given habitat may be more constrained than that of rewarding orchids. In particular, the success of deceptive orchids strongly relies on several biotic factors such as interactions with co-flowering rewarding species and pollinators, which may vary with altitude and over time. Our study compares generalized food-deceptive (i.e. excluding sexually deceptive) and rewarding orchids to test whether (i) deceptive orchids flower earlier compared to their rewarding counterparts and whether (ii) the relative occurrence of deceptive orchids decreases with increasing altitude. Methods To compare the flowering phenology of rewarding and deceptive orchids, we analysed data compiled from the literature at the species level over the occidental Palaearctic area. Since flowering phenology can be constrained by the latitudinal distribution of the species and by their phylogenetic relationships, we accounted for these factors in our analysis. To compare the altitudinal distribution of rewarding and deceptive orchids, we used field observations made over the entire Swiss territory and over two Swiss mountain ranges. Important Findings We found that deceptive orchid species start flowering earlier than rewarding orchids do, which is in accordance with the hypotheses of exploitation of naive pollinators and/or avoidance of competition with rewarding co-occurring species. Also, the relative frequency of deceptive orchids decreases with altitude, suggesting that deception may be less profitable at high compared to low altitude.

A stratified approach for modeling the distribution of a threatened ant species in the Swiss National Park

We present models predicting the potential distribution of a threatened ant species, Formica exsecta Nyl., in the Swiss National Park ( SNP). Data to fit the models have been collected according to a random-stratified design with an equal number of replicates per stratum. The basic aim of such a sampling strategy is to allow the formal testing of biological hypotheses about those factors most likely to account for the distribution of the modeled species. The stratifying factors used in this study were: vegetation, slope angle and slope aspect, the latter two being used as surrogates of solar radiation, considered one of the basic requirements of F. exsecta. Results show that, although the basic stratifying predictors account for more than 50% of the deviance, the incorporation of additional non-spatially explicit predictors into the model, as measured in the field, allows for an increased model performance (up to nearly 75%). However, this was not corroborated by permutation tests. Implementation on a national scale was made for one model only, due to the difficulty of obtaining similar predictors on this scale. The resulting map on the national scale suggests that the species might once have had a broader distribution in Switzerland. Reasons for its particular abundance within the SNP might possibly be related to habitat fragmentation and vegetation transformation outside the SNP boundaries.

Determining the best population-level alcohol consumption model and its impact on estimates of alcohol-attributable harms.

BACKGROUND: The goals of our study are to determine the most appropriate model for alcohol consumption as an exposure for burden of disease, to analyze the effect of the chosen alcohol consumption distribution on the estimation of the alcohol Population- Attributable Fractions (PAFs), and to characterize the chosen alcohol consumption distribution by exploring if there is a global relationship within the distribution. METHODS: To identify the best model, the Log-Normal, Gamma, and Weibull prevalence distributions were examined using data from 41 surveys from Gender, Alcohol and Culture: An International Study (GENACIS) and from the European Comparative Alcohol Study. To assess the effect of these distributions on the estimated alcohol PAFs, we calculated the alcohol PAF for diabetes, breast cancer, and pancreatitis using the three above-named distributions and using the more traditional approach based on categories. The relationship between the mean and the standard deviation from the Gamma distribution was estimated using data from 851 datasets for 66 countries from GENACIS and from the STEPwise approach to Surveillance from the World Health Organization. RESULTS: The Log-Normal distribution provided a poor fit for the survey data, with Gamma and Weibull distributions providing better fits. Additionally, our analyses showed that there were no marked differences for the alcohol PAF estimates based on the Gamma or Weibull distributions compared to PAFs based on categorical alcohol consumption estimates. The standard deviation of the alcohol distribution was highly dependent on the mean, with a unit increase in alcohol consumption associated with a unit increase in the mean of 1.258 (95% CI: 1.223 to 1.293) (R2 = 0.9207) for women and 1.171 (95% CI: 1.144 to 1.197) (R2 = 0. 9474) for men. CONCLUSIONS: Although the Gamma distribution and the Weibull distribution provided similar results, the Gamma distribution is recommended to model alcohol consumption from population surveys due to its fit, flexibility, and the ease with which it can be modified. The results showed that a large degree of variance of the standard deviation of the alcohol consumption Gamma distribution was explained by the mean alcohol consumption, allowing for alcohol consumption to be modeled through a Gamma distribution using only average consumption.

The truncated mean of an asymmetric distribution

This paper investigates a simple procedure to estimate robustly the mean of an asymmetric distribution. The procedure removes the observations which are larger or smaller than certain limits and takes the arithmetic mean of the remaining observations, the limits being determined with the help of a parametric model, e.g., the Gamma, the Weibull or the Lognormal distribution. The breakdown point, the influence function, the (asymptotic) variance, and the contamination bias of this estimator are explored and compared numerically with those of competing estimates.

Generalized linear and generalized additive models in studies of species distributions: setting the scene

An important statistical development of the last 30 years has been the advance in regression analysis provided by generalized linear models (GLMs) and generalized additive models (GAMs). Here we introduce a series of papers prepared within the framework of an international workshop entitled: Advances in GLMs/GAMs modeling: from species distribution to environmental management, held in Riederalp, Switzerland, 6-11 August 2001.We first discuss some general uses of statistical models in ecology, as well as provide a short review of several key examples of the use of GLMs and GAMs in ecological modeling efforts. We next present an overview of GLMs and GAMs, and discuss some of their related statistics used for predictor selection, model diagnostics, and evaluation. Included is a discussion of several new approaches applicable to GLMs and GAMs, such as ridge regression, an alternative to stepwise selection of predictors, and methods for the identification of interactions by a combined use of regression trees and several other approaches. We close with an overview of the papers and how we feel they advance our understanding of their application to ecological modeling.

Aitchison Geometry for Probability and Likelihood as a new approach to mathematical statistics

The Aitchison vector space structure for the simplex is generalized to a Hilbert space structure A2(P) for distributions and likelihoods on arbitrary spaces. Centralnotations of statistics, such as Information or Likelihood, can be identified in the algebraical structure of A2(P) and their corresponding notions in compositional data analysis, such as Aitchison distance or centered log ratio transform.In this way very elaborated aspects of mathematical statistics can be understoodeasily in the light of a simple vector space structure and of compositional data analysis. E.g. combination of statistical information such as Bayesian updating,combination of likelihood and robust M-estimation functions are simple additions/perturbations in A2(Pprior). Weighting observations corresponds to a weightedaddition of the corresponding evidence.Likelihood based statistics for general exponential families turns out to have aparticularly easy interpretation in terms of A2(P). Regular exponential families formfinite dimensional linear subspaces of A2(P) and they correspond to finite dimensionalsubspaces formed by their posterior in the dual information space A2(Pprior).The Aitchison norm can identified with mean Fisher information. The closing constant itself is identified with a generalization of the cummulant function and shown to be Kullback Leiblers directed information. Fisher information is the local geometry of the manifold induced by the A2(P) derivative of the Kullback Leibler information and the space A2(P) can therefore be seen as the tangential geometry of statistical inference at the distribution P.The discussion of A2(P) valued random variables, such as estimation functionsor likelihoods, give a further interpretation of Fisher information as the expected squared norm of evidence and a scale free understanding of unbiased reasoning

Cellular distribution of glucose and monocarboxylate transporters in human brain white matter and multiple sclerosis lesions.

To ensure efficient energy supply to the high demanding brain, nutrients are transported into brain cells via specific glucose (GLUT) and monocarboxylate transporters (MCT). Mitochondrial dysfunction and altered glucose metabolism are thought to play an important role in the progression of neurodegenerative diseases, including multiple sclerosis (MS). Here, we investigated the cellular localization of key GLUT and MCT proteins in human brain tissue of non-neurological controls and MS patients. We show that in control brain tissue GLUT and MCT proteins were abundantly expressed in a variety of central nervous system cells, particularly in microglia and endothelial cells. In active MS lesions, GLUTs and MCTs were highly expressed in infiltrating leukocytes and reactive astrocytes. Astrocytes manifest increased MCT1 staining and maintain GLUT expression in inactive lesions, whereas demyelinated axons exhibit significantly reduced GLUT3 and MCT2 immunoreactivity in inactive lesions. Finally, we demonstrated that the co-transcription factor peroxisome proliferator-activated receptor gamma co-activator 1-alpha (PGC-1α), an important protein involved in energy metabolism, is highly expressed in reactive astrocytes in active MS lesions. Overexpression of PGC-1α in astrocyte-like cells resulted in increased production of several GLUT and MCT proteins. In conclusion, we provide for the first time a comprehensive overview of key nutrient transporters in white matter brain samples. Moreover, our data demonstrate an altered expression of these nutrient transporters in MS brain tissue, including a marked reduction of axonal GLUT3 and MCT2 expression in chronic lesions, which may impede efficient nutrient supply to the hypoxic demyelinated axons thereby contributing to the ongoing neurodegeneration in MS. GLIA 2014;62:1125-1141.

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.