A note on the use of the generalized odds ratio in meta-analysis of association studies involving bi- and tri-allelic polymorphisms

Abstract Background The generalized odds ratio (GOR) was recently suggested as a genetic model-free measure for association studies. However, its properties were not extensively investigated. We used Monte Carlo simulations to investigate type-I error rates, power and bias in both effect size and between-study variance estimates of meta-analyses using the GOR as a summary effect, and compared these results to those obtained by usual approaches of model specification. We further applied the GOR in a real meta-analysis of three genome-wide association studies in Alzheimer's disease. Findings For bi-allelic polymorphisms, the GOR performs virtually identical to a standard multiplicative model of analysis (e.g. per-allele odds ratio) for variants acting multiplicatively, but augments slightly the power to detect variants with a dominant mode of action, while reducing the probability to detect recessive variants. Although there were differences among the GOR and usual approaches in terms of bias and type-I error rates, both simulation- and real data-based results provided little indication that these differences will be substantial in practice for meta-analyses involving bi-allelic polymorphisms. However, the use of the GOR may be slightly more powerful for the synthesis of data from tri-allelic variants, particularly when susceptibility alleles are less common in the populations (≤10%). This gain in power may depend on knowledge of the direction of the effects. Conclusions For the synthesis of data from bi-allelic variants, the GOR may be regarded as a multiplicative-like model of analysis. The use of the GOR may be slightly more powerful in the tri-allelic case, particularly when susceptibility alleles are less common in the populations.

Relation of Uric Acid to Serum Levels of High-Sensitivity C-Reactive Protein, Triglycerides, and High-Density Lipoprotein Cholesterol and to Hepatic Steatosis

Increased uric acid (UA) is strongly linked to cardiovascular disease. However, the independent role of UA is still debated because it is associated with several cardiovascular risk factors including obesity and metabolic syndrome. This study assessed the association of UA with increased high-sensitivity C-reactive protein (hs-CRP), increased ratio of triglyceride to high-density lipoprotein cholesterol (TG/HDL), sonographically detected hepatic steatosis, and their clustering in the presence and absence of obesity and metabolic syndrome. We evaluated 3,518 employed subjects without clinical cardiovascular disease from November 2008 through July 2010. Prevalence of tis-CRP >= 3 mg/L was 19%, that of TG/HDL >= 3 was 44%, and that of hepatic steatosis was 43%. In multivariable logistic regression after adjusting for traditional cardiovascular risk factors and confounders, highest versus lowest UA quartile was associated with hs-CRP >= 3 mg/L (odds ratio [OR] 1.52, 95% confidence interval [CI] 1.01 to 2.28, p = 0.04), TG/HDL >= 3 (OR 3.29, 95% CI 2.36 to 4.60, p <0.001), and hepatic steatosis (OR 3.10, 95% CI 2.22 to 4.32, p <0.001) independently of obesity and metabolic syndrome. Association of UA with hs-CRP >= 3 mg/L became nonsignificant in analyses stratified by obesity. Ascending UA quartiles compared to the lowest UA quartile demonstrated a graded increase in the odds of having 2 or 3 of these risk conditions and a successive decrease in the odds of having none. In conclusion, high UA levels were associated with increased TG/HDL and hepatic steatosis independently of metabolic syndrome and obesity and with increased hs-CRP independently of metabolic syndrome. (C) 2012 Elsevier Inc. All rights reserved. (Am J Cardiol 2012;110:1787-1792)

Parameter recovery for a skew-normal IRT model under a Bayesian approach: hierarchical framework, prior and kernel sensitivity and sample size

Item response theory (IRT) comprises a set of statistical models which are useful in many fields, especially when there is an interest in studying latent variables (or latent traits). Usually such latent traits are assumed to be random variables and a convenient distribution is assigned to them. A very common choice for such a distribution has been the standard normal. Recently, Azevedo et al. [Bayesian inference for a skew-normal IRT model under the centred parameterization, Comput. Stat. Data Anal. 55 (2011), pp. 353-365] proposed a skew-normal distribution under the centred parameterization (SNCP) as had been studied in [R. B. Arellano-Valle and A. Azzalini, The centred parametrization for the multivariate skew-normal distribution, J. Multivariate Anal. 99(7) (2008), pp. 1362-1382], to model the latent trait distribution. This approach allows one to represent any asymmetric behaviour concerning the latent trait distribution. Also, they developed a Metropolis-Hastings within the Gibbs sampling (MHWGS) algorithm based on the density of the SNCP. They showed that the algorithm recovers all parameters properly. Their results indicated that, in the presence of asymmetry, the proposed model and the estimation algorithm perform better than the usual model and estimation methods. Our main goal in this paper is to propose another type of MHWGS algorithm based on a stochastic representation (hierarchical structure) of the SNCP studied in [N. Henze, A probabilistic representation of the skew-normal distribution, Scand. J. Statist. 13 (1986), pp. 271-275]. Our algorithm has only one Metropolis-Hastings step, in opposition to the algorithm developed by Azevedo et al., which has two such steps. This not only makes the implementation easier but also reduces the number of proposal densities to be used, which can be a problem in the implementation of MHWGS algorithms, as can be seen in [R.J. Patz and B.W. Junker, A straightforward approach to Markov Chain Monte Carlo methods for item response models, J. Educ. Behav. Stat. 24(2) (1999), pp. 146-178; R. J. Patz and B. W. Junker, The applications and extensions of MCMC in IRT: Multiple item types, missing data, and rated responses, J. Educ. Behav. Stat. 24(4) (1999), pp. 342-366; A. Gelman, G.O. Roberts, and W.R. Gilks, Efficient Metropolis jumping rules, Bayesian Stat. 5 (1996), pp. 599-607]. Moreover, we consider a modified beta prior (which generalizes the one considered in [3]) and a Jeffreys prior for the asymmetry parameter. Furthermore, we study the sensitivity of such priors as well as the use of different kernel densities for this parameter. Finally, we assess the impact of the number of examinees, number of items and the asymmetry level on the parameter recovery. Results of the simulation study indicated that our approach performed equally as well as that in [3], in terms of parameter recovery, mainly using the Jeffreys prior. Also, they indicated that the asymmetry level has the highest impact on parameter recovery, even though it is relatively small. A real data analysis is considered jointly with the development of model fitting assessment tools. The results are compared with the ones obtained by Azevedo et al. The results indicate that using the hierarchical approach allows us to implement MCMC algorithms more easily, it facilitates diagnosis of the convergence and also it can be very useful to fit more complex skew IRT models.

Low Bone Mineral Density in Middle-Aged Breast Cancer Survivors: Prevalence and Associated Factors

Background: The aim of this study was to investigate the prevalence of low bone mineral density (BMD) and associated factors in middle-aged breast cancer survivors (BCS). Patients and Methods: A cross-sectional study was conducted with 70 BCS of 45-65 years of age undergoing complete oncology treatment. Logistic regression models were used to identify factors associated with low BMD (osteopenia and osteoporosis taken together as a single group). Results: The mean age of participants was 53.2 +/- 5.9 years. BMD was low at the femoral neck in 28.6% of patients and at the lumbar spine in 45.7%. Body mass index <= 30 kg/m(2) (adjusted odds ratio (OR) 3.43; 95% confidence interval (CI) 1.0-11.3) and postmenopausal status (OR adjusted 20.42; 95% CI 2.0-201.2) were associated with low BMD at the lumbar spine. Femoral neck measurements, age > 50 years (OR 3.41; 95% CI 1.0-11.6), and time since diagnosis > 50 months (OR adjusted 3.34; 95% CI 1.0-11.3) increased the likelihood of low BMD. Conclusion: These findings show that low BMD is common in middle-aged BCS. Factors were identified that may affect BMD in BCS and should be considered when implementing strategies to minimize bone loss in middle-aged women with breast cancer.

A general threshold stress hybrid hazard model for lifetime data

In this paper we propose a hybrid hazard regression model with threshold stress which includes the proportional hazards and the accelerated failure time models as particular cases. To express the behavior of lifetimes the generalized-gamma distribution is assumed and an inverse power law model with a threshold stress is considered. For parameter estimation we develop a sampling-based posterior inference procedure based on Markov Chain Monte Carlo techniques. We assume proper but vague priors for the parameters of interest. A simulation study investigates the frequentist properties of the proposed estimators obtained under the assumption of vague priors. Further, some discussions on model selection criteria are given. The methodology is illustrated on simulated and real lifetime data set.

Weibull and generalised exponential overdispersion models with an application to ozone air pollution

We consider the problem of estimating the mean and variance of the time between occurrences of an event of interest (inter-occurrences times) where some forms of dependence between two consecutive time intervals are allowed. Two basic density functions are taken into account. They are the Weibull and the generalised exponential density functions. In order to capture the dependence between two consecutive inter-occurrences times, we assume that either the shape and/or the scale parameters of the two density functions are given by auto-regressive models. The expressions for the mean and variance of the inter-occurrences times are presented. The models are applied to the ozone data from two regions of Mexico City. The estimation of the parameters is performed using a Bayesian point of view via Markov chain Monte Carlo (MCMC) methods.

IRF6 is a Risk Factor for Nonsyndromic Cleft Lip in the Brazilian Population

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is a complex disorder with a worldwide incidence estimated at 1:700. Among the putative susceptibility loci, the IRF6 gene and a region at 8q24.21 have been corroborated in different populations. To test the role of IRF6 in NSCL/P predisposition in the Brazilian population, we conducted a structured association study with the SNPs rs642961 and rs590223, respectively, located at 5' and 3' of the IRF6 gene and not in strong linkage disequilibrium (LD), in patients from five different Brazilian locations. We also evaluated the effect of these SNPs in IRF6 expression in mesenchymal stem cells (MSC). We observed association between rs642961 and cleft lip only (CLO) (P = 0.009; odds ratio (OR) for AA genotype = 1.83 [95% Confidence interval (CI), 0.64-5.31]; OR for AG genotype = 1.72 [95% CI, 1.03-2.84]). This association seems to be driven by the affected patients from Barbalha, a location which presents the highest heritability estimate (H-2 = 0.85), and the A allele at rs642961 is acting through a dominant model. No association was detected for the SNP rs590223. We did not find any correlation between expression levels and genotypes of the two loci, and it is possible that these SNPs have a functional role in some specific period of embryogenesis. (C) 2012 Wiley Periodicals, Inc.

VIP expression in vessels of the vaginal wall: relationship with estrogen status

Objective Vasoactive intestinal peptide (VIP) is a neuropeptide with elevated expression in regions that control urogenital functions. Estrogen appears to modulate VIP expression in various organs, but this effect has not been demonstrated in the vaginal wall. The aim of this study was to evaluate the influence of estrogen status on VIP expression in vessels of the vaginal wall. Methods Surgical specimens were removed from the vaginal walls of 18 premenopausal women and 12 postmenopausal women who were given surgery for genital prolapse grade I or II. Vaginal specimens were stained with estrogen receptor-alpha (ER-alpha) and VIP antibodies. Levels of follicle stimulating hormone (FSH), estradiol, prolactin, fasting glucose and serum thyroxine stimulating hormone were also measured. Estrogen status was assessed on the basis of FSH and ER-alpha scores. Results The vaginal walls of premenopausal women had significantly higher ER-alpha scores than those of menopausal women (premenopausal group, 3.6 +/- 2.2; menopausal group, 1.4 +/- 1.8; p = 0.01). Premenopausal women also had significantly higher levels of VIP in the vaginal wall than menopausal women (p = 0.02). Increasing age was associated with lower level of VIP staining (odds ratio 0.88; 95% confidence interval 0.78-0.99). Conclusion Levels of ER-alpha and VIP expression in the posterior vaginal wall were higher in premenopausal than in menopausal women, but VIP expression was not associated with estrogen status. Age was an independent predictor of VIP staining in vaginal wall biopsies.

Cointegration: Bayesian Significance Test

To estimate causal relationships, time series econometricians must be aware of spurious correlation, a problem first mentioned by Yule (1926). To deal with this problem, one can work either with differenced series or multivariate models: VAR (VEC or VECM) models. These models usually include at least one cointegration relation. Although the Bayesian literature on VAR/VEC is quite advanced, Bauwens et al. (1999) highlighted that "the topic of selecting the cointegrating rank has not yet given very useful and convincing results". The present article applies the Full Bayesian Significance Test (FBST), especially designed to deal with sharp hypotheses, to cointegration rank selection tests in VECM time series models. It shows the FBST implementation using both simulated and available (in the literature) data sets. As illustration, standard non informative priors are used.

Income inequality and elderly self-rated health in Sao Paulo, Brazil

Purpose: To test the association between income inequality and elderly self-rated health and to propose a pathway to explain the relationship. Methods: We analyzed a sample of 2143 older individuals (60 years of age and over) from 49 distritos of the Municipality of Sao Paulo, Brazil. Bayesian multilevel logistic models were performed with poor self-rated health as the outcome variable. Results: Income inequality (measured by the Gini coefficient) was found to be associated with poor self-rated health after controlling for age, sex, income and education (odds ratio, 1.19; 95% credible interval, 1.01-1.38). When the practice of physical exercise and homicide rate were added to the model, the Gini coefficient lost its statistical significance (P>.05). We fitted a structural equation model in which income inequality affects elderly health by a pathway mediated by violence and practice of physical exercise. Conclusions: The health of older individuals may be highly susceptible to the socioeconomic environment of residence, specifically to the local distribution of income. We propose that this association may be mediated by fear of violence and lack of physical activity. (C) 2012 Elsevier Inc. All rights reserved.

Bayesian analysis of random regression models using B-splines to model test-day milk yield of Holstein cattle in Brazil

The objective of this paper is to model variations in test-day milk yields of first lactations of Holstein cows by RR using B-spline functions and Bayesian inference in order to fit adequate and parsimonious models for the estimation of genetic parameters. They used 152,145 test day milk yield records from 7317 first lactations of Holstein cows. The model established in this study was additive, permanent environmental and residual random effects. In addition, contemporary group and linear and quadratic effects of the age of cow at calving were included as fixed effects. Authors modeled the average lactation curve of the population with a fourth-order orthogonal Legendre polynomial. They concluded that a cubic B-spline with seven random regression coefficients for both the additive genetic and permanent environment effects was to be the best according to residual mean square and residual variance estimates. Moreover they urged a lower order model (quadratic B-spline with seven random regression coefficients for both random effects) could be adopted because it yielded practically the same genetic parameter estimates with parsimony. (C) 2012 Elsevier B.V. All rights reserved.

A Bayesian destructive weighted Poisson cure rate model and an application to a cutaneous melanoma data

In this article, we propose a new Bayesian flexible cure rate survival model, which generalises the stochastic model of Klebanov et al. [Klebanov LB, Rachev ST and Yakovlev AY. A stochastic-model of radiation carcinogenesis - latent time distributions and their properties. Math Biosci 1993; 113: 51-75], and has much in common with the destructive model formulated by Rodrigues et al. [Rodrigues J, de Castro M, Balakrishnan N and Cancho VG. Destructive weighted Poisson cure rate models. Technical Report, Universidade Federal de Sao Carlos, Sao Carlos-SP. Brazil, 2009 (accepted in Lifetime Data Analysis)]. In our approach, the accumulated number of lesions or altered cells follows a compound weighted Poisson distribution. This model is more flexible than the promotion time cure model in terms of dispersion. Moreover, it possesses an interesting and realistic interpretation of the biological mechanism of the occurrence of the event of interest as it includes a destructive process of tumour cells after an initial treatment or the capacity of an individual exposed to irradiation to repair altered cells that results in cancer induction. In other words, what is recorded is only the damaged portion of the original number of altered cells not eliminated by the treatment or repaired by the repair system of an individual. Markov Chain Monte Carlo (MCMC) methods are then used to develop Bayesian inference for the proposed model. Also, some discussions on the model selection and an illustration with a cutaneous melanoma data set analysed by Rodrigues et al. [Rodrigues J, de Castro M, Balakrishnan N and Cancho VG. Destructive weighted Poisson cure rate models. Technical Report, Universidade Federal de Sao Carlos, Sao Carlos-SP. Brazil, 2009 (accepted in Lifetime Data Analysis)] are presented.

Low sclerostin levels: a predictive marker of persistent inflammation in ankylosing spondylitis during anti-tumor necrosis factor therapy?

Abstract Introduction Sclerostin levels have been reported to be low in ankylosing spondylitis (AS), but there is no data regarding the possible role of this Wnt inhibitor during anti-tumor necrosis factor (TNF) therapy. The present study longitudinally evaluated sclerostin levels, inflammatory markers and bone mineral density (BMD) in AS patients under anti-TNF therapy. Methods Thirty active AS patients were assessed at baseline, 6 and 12 months after anti-TNF therapy regarding clinical parameters, inflammatory markers, BMD and baseline radiographic damage (mSASSS). Thirty age- and sex-matched healthy individuals comprised the control group. Patients' sclerostin levels, sclerostin binding low-density lipoprotein receptor-related protein 6 (LRP6) and BMD were evaluated at the same time points and compared to controls. Results At baseline, AS patients had lower sclerostin levels (60.5 ± 32.7 vs. 96.7 ± 52.9 pmol/L, P = 0.002) and comparable sclerostin binding to LRP6 (P = 0.387) than controls. Improvement of Bath Ankylosing Spondylitis Disease Activity Index (BASDAI), Bath Ankylosing Spondylitis Functional Index (BASFI), Bath Ankylosing Spondylitis Metrology Index (BASMI), Ankylosing Spondylitis quality of life (ASQoL) was observed at baseline vs. 6 vs. 12 months (P < 0.01). Concomitantly, a gradual increase in spine BMD (P < 0.001) and a positive correlation between baseline mSASSS and spine BMD was found (r = 0.468, P < 0.01). Inflammatory parameters reduction was observed comparing baseline vs. 6 vs. 12 months (P <0.01). Sclerostin levels progressively increased [baseline (60.5 ± 32.7) vs. 6 months (67.1 ± 31.9) vs. 12 months (72.7 ± 32.3) pmol/L, P <0.001]. At 12 months, the sclerostin levels remained significantly lower in patients compared to controls (72.7 ± 32.3 vs. 96.70 ± 52.85 pmol/L, P = 0.038). Moreover, sclerostin serum levels at 12 months were lower in the 10 patients with high C reactive protein (CRP) (≥ 5 mg/l) compared to the other 20 patients with normal CRP (P = 0.004). Of note, these 10 patients with persistent inflammation also had lower sclerostin serum levels at baseline compared to the other patients (P = 0.023). Univariate logistic regression analysis demonstrated that AS patients with lower sclerostin serum levels had an increased risk to have high CRP at 12 months (odds ratio = 7.43, 95% CI 1.23 to 45.01, P = 0.020) than those with higher sclerostin values. Conclusions Persistent low sclerostin levels may underlie continuous inflammation in AS patients under anti-TNF therapy.