Diagnosis of variability and trends in a global precipitation dataset using a physically motivated statistical model

A physically motivated statistical model is used to diagnose variability and trends in wintertime ( October - March) Global Precipitation Climatology Project (GPCP) pentad (5-day mean) precipitation. Quasi-geostrophic theory suggests that extratropical precipitation amounts should depend multiplicatively on the pressure gradient, saturation specific humidity, and the meridional temperature gradient. This physical insight has been used to guide the development of a suitable statistical model for precipitation using a mixture of generalized linear models: a logistic model for the binary occurrence of precipitation and a Gamma distribution model for the wet day precipitation amount. The statistical model allows for the investigation of the role of each factor in determining variations and long-term trends. Saturation specific humidity q(s) has a generally negative effect on global precipitation occurrence and with the tropical wet pentad precipitation amount, but has a positive relationship with the pentad precipitation amount at mid- and high latitudes. The North Atlantic Oscillation, a proxy for the meridional temperature gradient, is also found to have a statistically significant positive effect on precipitation over much of the Atlantic region. Residual time trends in wet pentad precipitation are extremely sensitive to the choice of the wet pentad threshold because of increasing trends in low-amplitude precipitation pentads; too low a choice of threshold can lead to a spurious decreasing trend in wet pentad precipitation amounts. However, for not too small thresholds, it is found that the meridional temperature gradient is an important factor for explaining part of the long-term trend in Atlantic precipitation.

Diagnosis of variability and trends in a global precipitation dataset using a physically motivated statistical model

A physically motivated statistical model is used to diagnose variability and trends in wintertime ( October - March) Global Precipitation Climatology Project (GPCP) pentad (5-day mean) precipitation. Quasi-geostrophic theory suggests that extratropical precipitation amounts should depend multiplicatively on the pressure gradient, saturation specific humidity, and the meridional temperature gradient. This physical insight has been used to guide the development of a suitable statistical model for precipitation using a mixture of generalized linear models: a logistic model for the binary occurrence of precipitation and a Gamma distribution model for the wet day precipitation amount. The statistical model allows for the investigation of the role of each factor in determining variations and long-term trends. Saturation specific humidity q(s) has a generally negative effect on global precipitation occurrence and with the tropical wet pentad precipitation amount, but has a positive relationship with the pentad precipitation amount at mid- and high latitudes. The North Atlantic Oscillation, a proxy for the meridional temperature gradient, is also found to have a statistically significant positive effect on precipitation over much of the Atlantic region. Residual time trends in wet pentad precipitation are extremely sensitive to the choice of the wet pentad threshold because of increasing trends in low-amplitude precipitation pentads; too low a choice of threshold can lead to a spurious decreasing trend in wet pentad precipitation amounts. However, for not too small thresholds, it is found that the meridional temperature gradient is an important factor for explaining part of the long-term trend in Atlantic precipitation.

Detection of frailty in Weibull lifetime data using Outlier tests

Heterogeneity in lifetime data may be modelled by multiplying an individual's hazard by an unobserved frailty. We test for the presence of frailty of this kind in univariate and bivariate data with Weibull distributed lifetimes, using statistics based on the ordered Cox-Snell residuals from the null model of no frailty. The form of the statistics is suggested by outlier testing in the gamma distribution. We find through simulation that the sum of the k largest or k smallest order statistics, for suitably chosen k , provides a powerful test when the frailty distribution is assumed to be gamma or positive stable, respectively. We provide recommended values of k for sample sizes up to 100 and simple formulae for estimated critical values for tests at the 5% level.

Modeling the variability of single-cell lag times for Listeria innocua populations after sublethal and lethal heat treatments

Optical density measurements were used to estimate the effect of heat treatments on the single-cell lag times of Listeria innocua fitted to a shifted gamma distribution. The single-cell lag time was subdivided into repair time ( the shift of the distribution assumed to be uniform for all cells) and adjustment time (varying randomly from cell to cell). After heat treatments in which all of the cells recovered (sublethal), the repair time and the mean and the variance of the single-cell adjustment time increased with the severity of the treatment. When the heat treatments resulted in a loss of viability (lethal), the repair time of the survivors increased with the decimal reduction of the cell numbers independently of the temperature, while the mean and variance of the single-cell adjustment times remained the same irrespective of the heat treatment. Based on these observations and modeling of the effect of time and temperature of the heat treatment, we propose that the severity of a heat treatment can be characterized by the repair time of the cells whether the heat treatment is lethal or not, an extension of the F value concept for sublethal heat treatments. In addition, the repair time could be interpreted as the extent or degree of injury with a multiple-hit lethality model. Another implication of these results is that the distribution of the time for cells to reach unacceptable numbers in food is not affected by the time-temperature combination resulting in a given decimal reduction.

Secular trends in daily precipitation characteristics: greenhouse gas simulation with a coupled AOGCM

Secular trends of daily precipitation characteristics are considered in the transient climate change experiment with a coupled atmosphere-ocean general circulation model ECHAM4/OPYC3 for 1900-2099. The climate forcing is due to increasing concentrations of the greenhouse gases in the atmosphere. Mean daily precipitation, precipitation intensity, probability of wet days and parameters of the gamma distribution are analyzed. Particular attention is paid to the changes of heavy precipitation, Analysis of the annual mean precipitation trends for 1900-1999 revealed general agreement with observations with significant positive trends in mean precipitation over continental areas. In the 2000-2099 period precipitation trend patterns followed the tendency obtained for 1900-1999 but with significantly increased magnitudes. Unlike the annual mean precipitation trends for which negative values were found for some continental areas, the mean precipitation intensity and scale parameter of the fitted gamma distribution increased over all land territories . Negative trends in the number of wet days were found over most of the land areas except high latitudes in the Northern Hemisphere. The shape parameter of the gamma distribution in general revealed a slight negative trend in the areas of the precipitation increase. Investigation of daily precipitation revealed an unproportional increase of heavy precipitation events for the land areas including local maxima in Europe and the eastern United States.

A new family of generalized distributions

Kumaraswamy [Generalized probability density-function for double-bounded random-processes, J. Hydrol. 462 (1980), pp. 79-88] introduced a distribution for double-bounded random processes with hydrological applications. For the first time, based on this distribution, we describe a new family of generalized distributions (denoted with the prefix `Kw`) to extend the normal, Weibull, gamma, Gumbel, inverse Gaussian distributions, among several well-known distributions. Some special distributions in the new family such as the Kw-normal, Kw-Weibull, Kw-gamma, Kw-Gumbel and Kw-inverse Gaussian distribution are discussed. We express the ordinary moments of any Kw generalized distribution as linear functions of probability weighted moments (PWMs) of the parent distribution. We also obtain the ordinary moments of order statistics as functions of PWMs of the baseline distribution. We use the method of maximum likelihood to fit the distributions in the new class and illustrate the potentiality of the new model with an application to real data.

Younger siblings, C-Reactive protein, and risk of age-related macular degeneration

In this study, we examined the relationship between exposure to siblings and 1) the risk of age-related macular degeneration (AMD) and 2) C-reactive protein levels. We retrospectively analyzed pooled cross-sectional data from 2 studies: the Cardiovascular Health and Age-Related Maculopathy Study (2001–2002) and the Age-Related Maculopathy Statin Study (2004–2006). Associations between number of siblings and AMD were assessed by using multinomial logistic regression. Associations between number of siblings and C-reactive protein levels were examined by using a generalized linear model for γ distribution. A higher number of younger siblings was associated with significantly lower odds of early AMD in those with a family history of AMD (odds ratio = 0.2, 95% confidence interval: 0.1, 0.8) (P = 0.022) but was unrelated to AMD for those who had no family history of the disease (odds ratio = 1.0, 95% confidence interval: 0.9, 1.2) (P = 0.874). A higher number of younger siblings correlated with lower C-reactive protein levels (β = −0.19, 95% confidence interval: −0.38, −0.01) (P = 0.036). This supports the theory that immune modulation contributes to AMD pathogenesis and suggests that exposure to younger siblings might be protective when there is a family history of AMD.

Predicting future purchases with the Poisson log-normal model

The negative binomial distribution (NBD) has been widely used in marketing for modeling purchase frequency counts, particularly in packaged goods contexts. A key managerially relevant use of this model is Conditional Trend Analysis (CTA)—a method of benchmarking future sales utilizing the NBD conditional expectation. CTA allows brand managers to identify whether the sales change in a second period is accounted for by previous non-, light, or heavy buyers of the brand. Although a useful tool, the conditional prediction of the NBD suffers from a bias: it under predicts what the period-one non-buyer class will do in period two and over predicts the sales contribution of existing buyers. In addition, the NBD's assumption of a gamma-distributed mean purchase rate lacks theoretical support—it is not possible to explain why a gamma distribution should hold. This paper therefore proposes an alternative model using a log-normal distribution in place of the gamma distribution, hence creating a Poisson log-normal (PLN) distribution. The PLN distribution has a stronger theoretical grounding than the NBD as it has a natural interpretation relying on the central limit theorem. Empirical analysis of brands in multiple categories shows that the PLN distribution gives better predictions than the NBD.

Avaliação genética da eficiência reprodutiva em vacas pardo-suíças por meio da análise de sobrevivência

Among the traits of economic importance to dairy cattle livestock those related to sexual precocity and longevity of the herd are essential to the success of the activity, because the stayability time of a cow in a herd is determined by their productive and reproductive lives. In Brazil, there are few studies about the reproductive efficiency of Swiss-Brown cows and no study was found using the methodology of survival analysis applied to this breed. Thus, in the first chapter of this study, the age at first calving from Swiss-Brown heifers was analyzed as the time until the event by the nonparametric method of Kaplan-Meier and the gamma shared frailty model, under the survival analysis methodology. Survival and hazard rate curves associated with this event were estimated and identified the influence of covariates on such time. The mean and median times at the first calving were 987.77 and 1,003 days, respectively, and significant covariates by the Log-Rank test, through Kaplan-Meier analysis, were birth season, calving year, sire (cow s father) and calving season. In the analysis by frailty model, the breeding values and the frailties of the sires (fathers) for the calving were predicted modeling the risk function of each cow as a function of the birth season as fixed covariate and sire as random covariate. The frailty followed the gamma distribution. Sires with high and positive breeding values possess high frailties, what means shorter survival time of their daughters to the event, i.e., reduction in the age at first calving of them. The second chapter aimed to evaluate the longevity of dairy cows using the nonparametric Kaplan-Meier and the Cox and Weibull proportional hazards models. It were simulated 10,000 records of the longevity trait from Brown-Swiss cows involving their respective times until the occurrence of five consecutive calvings (event), considered here as typical of a long-lived cow. The covariates considered in the database were age at first calving, herd and sire (cow s father). All covariates had influence on the longevity of cows by Log-Rank and Wilcoxon tests. The mean and median times to the occurrence of the event were 2,436.285 and 2,437 days, respectively. Sires that have higher breeding values also have a greater risk of that their daughters reach the five consecutive calvings until 84 months

Clinical correlates of social adjustment in patients with obsessive-compulsive disorder

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Modelagem estatística do frete da madeira no Uruguai

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Eficiência de estimadores, geradores e algoritmos na simulação de dados diários de precipitação pluviométrica utilizando a distribuição gama

Pós-graduação em Agronomia (Energia na Agricultura) - FCA

Comparação entre os estimadores de mínimos quadrados ordinários e mínimos desvios absolutos em modelos de regressão linear simples -: uma aplicação na energia na agricultura

Pós-graduação em Agronomia (Energia na Agricultura) - FCA

Comparação dos efeitos da hidrocinesioterapia e do treinamento físico realizado no solo em hipertensas

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Estudo da idade da vaca ao último parto para avaliar longevidade em rebanhos da raça nelore por análise de sobrevivência

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)