Statistical analysis of maximum likelihood estimator images of human brain FDG PET studies

The work presented evaluates the statistical characteristics of regional bias and expected error in reconstructions of real positron emission tomography (PET) data of human brain fluoro-deoxiglucose (FDG) studies carried out by the maximum likelihood estimator (MLE) method with a robust stopping rule, and compares them with the results of filtered backprojection (FBP) reconstructions and with the method of sieves. The task of evaluating radioisotope uptake in regions-of-interest (ROIs) is investigated. An assessment of bias and variance in uptake measurements is carried out with simulated data. Then, by using three different transition matrices with different degrees of accuracy and a components of variance model for statistical analysis, it is shown that the characteristics obtained from real human FDG brain data are consistent with the results of the simulation studies.

Treatment of penile carcinoma: to cut or not to cut?

PURPOSE: The aim of this study was to assess the outcome in patients with penile cancer. METHODS AND MATERIALS: A total of 60 patients with penile carcinoma were included. Of the patients, 45 (n = 27) underwent surgery, and 51 underwent definitive (n = 29) or postoperative (n = 22) radiotherapy (RT). Median follow-up was 62 months. RESULTS: Median time to locoregional relapse was 14 months. Local failure was observed in 3 of 23 patients (13%) treated with surgery with or without postoperative RT vs. in 19 of 33 patients (56%) given organ-sparing treatment (p = 0.0008). Of 22 local failures, 16 (73%) were salvaged with surgery. Of the 33 patients treated with definitive RT (n = 29) and the 4 patients refusing RT after excisional biopsy, local control was obtained with organ preservation in 13 (39%). In the remaining 20, 4 patients with local failure underwent salvage conservatively, resulting in an ultimate penis preservation rate of 17 of 33 (52%) patients treated with definitive RT. The 5-year and 10-year probability of surviving with an intact penis was 43% and 26%, respectively. There was no survival difference between the patients treated with definitive RT and primary surgery (56% vs. 53%; p = 0.16). In multivariate analysis, independent factors influencing survival were N-classification and pathologic grade. Surgery was the only independent predictor for better local control. CONCLUSION: Based on our study findings, in patients with penile cancer, local control is superior with surgery. However, there is no difference in survival between patients treated with surgery and those treated with definitive RT, with 52% organ preservation.

Quality of care and survival of haemodialysed patients in western Switzerland.

BACKGROUND: Many factors affect survival in haemodialysis (HD) patients. Our aim was to study whether quality of clinical care may affect survival in this population, when adjusted for demographic characteristics and co-morbidities. METHODS: We studied survival in 553 patients treated by chronic HD during March 2001 in 21 dialysis facilities in western Switzerland. Indicators of quality of care were established for anaemia control, calcium and phosphate product, serum albumin, pre-dialysis blood pressure (BP), type of vascular access and dialysis adequacy (spKt/V) and their baseline values were related to 3-year survival. The modified Charlson co-morbidity index (including age) and transplantation status were also considered as a predictor of survival. RESULTS: Three-year survival was obtained for 96% of the patients; 39% (211/541) of these patients had died. The 3-year survival was 50, 62 and 69%, respectively, in patients who had 0-2, 3 and >or=4 fulfilled indicators of quality of care (test for linear trend, P < 0.001). In a Cox multivariate analysis model, the absence of transplantation, a higher modified Charlson's score, decreased fulfilment of indicators of good clinical care and low pre-dialysis systolic BP were independent predictors of death. CONCLUSION: Good clinical care improves survival in HD patients, even after adjustment for availability of transplantation and co-morbidities.

Predictors of poor response to methotrexate in polyarticular-course juvenile idiopathic arthritis: analysis of the PRINTO methotrexate trial.

OBJECTIVES: To determine whether baseline demographic, clinical, articular and laboratory variables predict methotrexate (MTX) poor response in polyarticular-course juvenile idiopathic arthritis. METHODS: Patients newly treated for 6 months with MTX enrolled in the Paediatric Rheumatology International Trials Organization (PRINTO) MTX trial. Bivariate and logistic regression analyses were used to identify baseline predictors of poor response according to the American College of Rheumatology pediatric (ACR-ped) 30 and 70 criteria. RESULTS: In all, 405/563 (71.9%) of patients were women; median age at onset and disease duration were 4.3 and 1.4 years, respectively, with anti-nuclear antibody (ANA) detected in 259/537 (48.2%) patients. With multivariate logistic regression analysis, the most important determinants of ACR-ped 70 non-responders were: disease duration > 1.3 years (OR 1.93), ANA negativity (OR 1.77), Childhood Health Assessment Questionnaire (CHAQ) disability index > 1.125 (OR 1.65) and the presence of right and left wrist activity (OR 1.55). Predictors of ACR-ped 30 non-responders were: ANA negativity (OR 1.92), CHAQ disability index > 1.14 (OR 2.18) and a parent's evaluation of child's overall well-being < or = 4.69 (OR 2.2). CONCLUSION: The subgroup of patients with longer disease duration, ANA negativity, higher disability and presence of wrist activity were significantly associated with a poorer response to a 6-month MTX course.

Genetic association tests: a method for the joint analysis of family and case-control data.

With the trend in molecular epidemiology towards both genome-wide association studies and complex modelling, the need for large sample sizes to detect small effects and to allow for the estimation of many parameters within a model continues to increase. Unfortunately, most methods of association analysis have been restricted to either a family-based or a case-control design, resulting in the lack of synthesis of data from multiple studies. Transmission disequilibrium-type methods for detecting linkage disequilibrium from family data were developed as an effective way of preventing the detection of association due to population stratification. Because these methods condition on parental genotype, however, they have precluded the joint analysis of family and case-control data, although methods for case-control data may not protect against population stratification and do not allow for familial correlations. We present here an extension of a family-based association analysis method for continuous traits that will simultaneously test for, and if necessary control for, population stratification. We further extend this method to analyse binary traits (and therefore family and case-control data together) and accurately to estimate genetic effects in the population, even when using an ascertained family sample. Finally, we present the power of this binary extension for both family-only and joint family and case-control data, and demonstrate the accuracy of the association parameter and variance components in an ascertained family sample.

Relationships between microbial activity and soil physical and chemical properties in native and reforested Araucaria angustifolia forests in the state of São Paulo, Brazil

Araucaria angustifolia (Bert.) O. Kuntze is the main component of the Mixed Ombrophilous forest and, in the State of São Paulo, it is associated with a high diversity of soil organisms, essential for the maintenance of soil quality, making the conservation of this ecosystem a major and pressing challenge. The objective of this study was to identify the physical and chemical properties that are most closely correlated with dehydrogenase enzyme activity, basal respiration and microbial biomass under native (NF) and replanted (RF) Araucaria angustifolia forests in three regions of the state of São Paulo, in winter and summer. The main differentiating factors between the areas were also determined. Each forest was represented by three true replications; at each site, from around the araucaria trees, 15 soil samples (0-20 cm) were collected to evaluate the soil physical, chemical and microbiological properties. At the same points, forest litter was sampled to assess mass and chemical properties. The following microbiological properties were evaluated: microbial biomass carbon (MBC), basal respiration (CO2-C), metabolic quotient (Q: CO2), dehydrogenase enzyme activity (DHA) as well as the physical properties (moisture, bulk density, macroporosity and total porosity), soil chemical properties [pH, organic carbon (org-C), P, Ca, K, Mg, Al, H+Al], litter dry mass, and C, N and S contents. The data were subjected to analysis of variance (TWO-WAY: ANOVA). A Canonical Discriminant Analysis (CDA) and a Canonical Correlation Analysis (CCA) were also performed. In the soil under NF, the values of K, P, soil macroporosity, and litter dry mass were higher and Q: CO2 and DHA lower, regardless of the sampling period, and DHA was lower in winter. In the RF areas, the levels of moisture, porosity and Q: CO2 were higher in both sampling periods, and DHA was higher in winter. The MBC was only higher under NF in the summer, while the litter contents of C, N and S were greater in winter. In winter, CCA showed a high correlation of DHA with CO2-C, pH and H+Al, while in the summer org-C, moisture, Mg, pH and litter C were more associated with DHA and CO2-C. The CDA indicated H+Al, available P, total porosity, litter S content, and soil moisture as the most discriminating variables between NF and RF, but moisture was the most relevant, in both seasons and CO2-C only in winter. The combined analysis of CCA and CDA showed that the contribution of the microbiological variables to a differentiation of the areas was small at both samplings, which may indicate that the period after reforestation was long enough to allow an almost complete recovery of the microbial activity.

Safety and functional outcome of thrombolysis in dissection-related ischemic stroke: a meta-analysis of individual patient data.

Background and Purpose-The safety and efficacy of thrombolysis in cervical artery dissection (CAD) are controversial. The aim of this meta-analysis was to pool all individual patient data and provide a valid estimate of safety and outcome of thrombolysis in CAD.Methods-We performed a systematic literature search on intravenous and intra-arterial thrombolysis in CAD. We calculated the rates of pooled symptomatic intracranial hemorrhage and mortality and indirectly compared them with matched controls from the Safe Implementation of Thrombolysis in Stroke-International Stroke Thrombolysis Register. We applied multivariate regression models to identify predictors of excellent (modified Rankin Scale=0 to 1) and favorable (modified Rankin Scale=0 to 2) outcome.Results-We obtained individual patient data of 180 patients from 14 retrospective series and 22 case reports. Patients were predominantly female (68%), with a mean +/- SD age of 46 +/- 11 years. Most patients presented with severe stroke (median National Institutes of Health Stroke Scale score=16). Treatment was intravenous thrombolysis in 67% and intra-arterial thrombolysis in 33%. Median follow-up was 3 months. The pooled symptomatic intracranial hemorrhage rate was 3.1% (95% CI, 1.3 to 7.2). Overall mortality was 8.1% (95% CI, 4.9 to 13.2), and 41.0% (95% CI, 31.4 to 51.4) had an excellent outcome. Stroke severity was a strong predictor of outcome. Overlapping confidence intervals of end points indicated no relevant differences with matched controls from the Safe Implementation of Thrombolysis in Stroke-International Stroke Thrombolysis Register.Conclusions-Safety and outcome of thrombolysis in patients with CAD-related stroke appear similar to those for stroke from all causes. Based on our findings, thrombolysis should not be withheld in patients with CAD. (Stroke. 2011;42:2515-2520.)

Pharmacogenetics-based population pharmacokinetic analysis of etravirine in HIV-1 infected individuals.

OBJECTIVES: Etravirine (ETV) is metabolized by cytochrome P450 (CYP) 3A, 2C9, and 2C19. Metabolites are glucuronidated by uridine diphosphate glucuronosyltransferases (UGT). To identify the potential impact of genetic and non-genetic factors involved in ETV metabolism, we carried out a two-step pharmacogenetics-based population pharmacokinetic study in HIV-1 infected individuals. MATERIALS AND METHODS: The study population included 144 individuals contributing 289 ETV plasma concentrations and four individuals contributing 23 ETV plasma concentrations collected in a rich sampling design. Genetic variants [n=125 single-nucleotide polymorphisms (SNPs)] in 34 genes with a predicted role in ETV metabolism were selected. A first step population pharmacokinetic model included non-genetic and known genetic factors (seven SNPs in CYP2C, one SNP in CYP3A5) as covariates. Post-hoc individual ETV clearance (CL) was used in a second (discovery) step, in which the effect of the remaining 98 SNPs in CYP3A, P450 cytochrome oxidoreductase (POR), nuclear receptor genes, and UGTs was investigated. RESULTS: A one-compartment model with zero-order absorption best characterized ETV pharmacokinetics. The average ETV CL was 41 (l/h) (CV 51.1%), the volume of distribution was 1325 l, and the mean absorption time was 1.2 h. The administration of darunavir/ritonavir or tenofovir was the only non-genetic covariate influencing ETV CL significantly, resulting in a 40% [95% confidence interval (CI): 13-69%] and a 42% (95% CI: 17-68%) increase in ETV CL, respectively. Carriers of rs4244285 (CYP2C19*2) had 23% (8-38%) lower ETV CL. Co-administered antiretroviral agents and genetic factors explained 16% of the variance in ETV concentrations. None of the SNPs in the discovery step influenced ETV CL. CONCLUSION: ETV concentrations are highly variable, and co-administered antiretroviral agents and genetic factors explained only a modest part of the interindividual variability in ETV elimination. Opposing effects of interacting drugs effectively abrogate genetic influences on ETV CL, and vice-versa.

Une méta-analyse GWA de la glycémie après 2h d'HGPO révèle que GIPR est associé avec la sécrétion de l'insuline en réponse au glucose et que l'ADCY5 est un nouveau gène de suseptibilité au diabète de type 2. [A meta-analysis of GWA blood glucose after 2 h of OGTT revealed that GIPR is associated with the secretion of insulin based on glucose and ADCY5 is a new susceptibility gene for type 2 diabetes.]

Introduction: Les études GVvA (Genome-wide association ,-studies) ont identifié et confirmé plus de 20 gènes de susceptibilité au DT2 et ont contribué à mieux comprendre sa physiopathologie. L'hyperglycémie à jeun (GJ), et 2 heures après une HGPO (G2h) sont les deux mesures cliniques du diagnostic du DT2. Nous avons identifié récemment la G6P du pancréas (G6PC2) comme déterminant de la variabilité physiologique de la GJ puis Ie récepteur à la mélatonine (MTNRIB) qui de plus lie la régulation du rythme circadien au DT2. Dans ce travail nous avons étudié la génétique de la G2h à l'aide de l'approche GWA. Résultats: Nous avons réalisé une méta-analyse GWA dans le cadre de MAGIC (Meta-Analysis of Glucose and Insulin related traits Consortium) qui a inclus 9 études GWA (N=15'234). La réplication de 29 loci (N=6958-30 121, P < 10-5 ) a confirmé 5 nouveaux loci; 2 étant connus comme associés avec Ie DT2 (TCF7L2, P = 1,6 X 10-10 ) et la GJ (GCKR, p = 5,6 X 10-10 ); alors que GIPR (p= 5,2 X 10-12), VSP13C (p= 3,9 X 10-8) et ADCY5 (p = 1,11 X 10-15 ) sont inédits. GIPR code Ie récepteur au GIP (gastric inhibitory polypeptide) qui est sécrété par les ceIlules intestinales pour stimuler la sécrétion de l'insuline en réponse au glucose (l'effet incrétine). Les porteurs du variant GIPR qui augmente la G2h ont également un indice insulinogénique plus bas, (p= 1,0 X 10-17) mais ils ne présentent aucune modification de leur glycémie suite à une hyperglycémie provoquée par voie veineuse (p= 0,21). Ces résultats soutiennent un effet incrétine du locus GIPR qui expliquerait ~9,6 % de la variance total de ce trait. La biologie de ADCY5 et VPS13C et son lien avec l'homéostasie du glucose restent à élucider. GIPR n'est pas associé avec le risque de DT2 indiquant qu'il influence la variabilité physiologique de la G2h alors que le locus ADCY5 est associé avec le DT2 (OR = 1,11, P = 1,5 X 10-15). Conclusion: Notre étude démontre que l'étude de la G2h est une approche efficace d'une part pour la compréhension de la base génétique de la physiologie de ce trait clinique important et d'autre part pour identifier de nouveaux gènes de susceptibilité au DT2.

Density, climate and varying return points: an analysis of long-term population fluctuations in the threatened European tree frog.

Experimental research has identified many putative agents of amphibian decline, yet the population-level consequences of these agents remain unknown, owing to lack of information on compensatory density dependence in natural populations. Here, we investigate the relative importance of intrinsic (density-dependent) and extrinsic (climatic) factors impacting the dynamics of a tree frog (Hyla arborea) population over 22 years. A combination of log-linear density dependence and rainfall (with a 2-year time lag corresponding to development time) explain 75% of the variance in the rate of increase. Such fluctuations around a variable return point might be responsible for the seemingly erratic demography and disequilibrium dynamics of many amphibian populations.

Multivariate analyses of genotype x environment interaction of popcorn

The objectives of this work were to evaluate the genotype x environment (GxE) interaction for popcorn and to compare two multivariate analyses methods. Nine popcorn cultivars were sown on four dates one month apart during each of the agricultural years 1998/1999 and 1999/2000. The experiments were carried out using randomized block designs, with four replicates. The cv. Zélia contributed the least to the GxE interaction. The cv. Viçosa performed similarly to cv. Rosa-claro. Optimization of GxE was obtained for cv. CMS 42 for a favorable mega-environment, and for cv. CMS 43 for an unfavorable environment. Multivariate analysis supported the results from the method of Eberhart & Russell. The graphic analysis of the Additive Main effects and Multiplicative Interaction (AMMI) model was simple, allowing conclusions to be made about stability, genotypic performance, genetic divergence between cultivars, and the environments that optimize cultivar performance. The graphic analysis of the Genotype main effects and Genotype x Environment interaction (GGE) method added to AMMI information on environmental stratification, defining mega-environments and the cultivars that optimized performance in those mega-environments. Both methods are adequate to explain the genotype x environment interactions.

Analysis of sensitivity with respect to a compositional parameter : A comment on 'Local sensitivity analysis for compositional data with application to soil texture in hydrologic modelling' by L. Loosvelt and co-authors

A comment about the article “Local sensitivity analysis for compositional data with application to soil texture in hydrologic modelling” writen by L. Loosvelt and co-authors. The present comment is centered in three specific points. The first one is related to the fact that the authors avoid the use of ilr-coordinates. The second one refers to some generalization of sensitivity analysis when input parameters are compositional. The third tries to show that the role of the Dirichlet distribution in the sensitivity analysis is irrelevant

Characterization and multivariate classification of grapes and wines of two Cabernet Sauvignon clones

The objective of this work was to assess and characterize two clones, 169 and 685, of Cabernet Sauvignon grapes and to evaluate the wine produced from these grapes. The experiment was carried out in São Joaquim, SC, Brazil, during the 2009 harvest season. During grape ripening, the evolution of physical-chemical properties, phenolic compounds, organic acids, and anthocyanins was evaluated. During grape harvest, yield components were determined for each clone. Individual and total phenolics, individual and total anthocyanins, and antioxidant activity were evaluated for wine. The clones were also assessed regarding the duration of their phenological cycle. During ripening, the evolution of phenolic compounds and of physical-chemical parameters was similar for both clones; however, during harvest, significant differences were observed regarding yield, number of bunches per plant and berries per bunch, leaf area, and organic acid, polyphenol, and anthocyanin content. The wines produced from these clones showed significant differences regarding chemical composition. The clones showed similar phenological cycle and responses to bioclimatic parameters. Principal component analysis shows that clone 685 is strongly correlated with color characteristics, mainly monomeric anthocyanins, while clone 169 is correlated with individual phenolic compounds.

Factors predisposing to dislocation after primary total hip arthroplasty: a multivariate analysis.

We conducted this study to determine the relative influence of various mechanical and patient-related factors on the incidence of dislocation after primary total hip asthroplasty (THA). Of 2,023 THAs, 21 patients who had at least 1 dislocation were compared with a control group of 21 patients without dislocation, matched for age, gender, pathology, and year of surgery. Implant positioning, seniority of the surgeon, American Society of Anesthesiologists (ASA) score, and diminished motor coordination were recorded. Data analysis included univariate and multivariate methods. The dislocation risk was 6.9 times higher if total anteversion was not between 40 degrees and 60 degrees and 10 times higher in patients with high ASA scores. Surgeons should pay attention to total anteversion (cup and stem) of THA. The ASA score should be part of the preoperative assessment of the dislocation risk.

Predictors of mortality from pump failure and sudden cardiac death in patients with systolic heart failure and left ventricular dyssynchrony: results of the CARE-HF trial.

BACKGROUND: Determining a specific death cause may facilitate individualized therapy in patients with heart failure (HF). Cardiac resynchronization therapy (CRT) decreased mortality in the Cardiac Resynchronization in Heart Failure trial by reducing pump failure and sudden cardiac death (SCD). This study analyzes predictors of specific causes of death. METHODS AND RESULTS: Univariate and multivariate analyses used 8 baseline and 3-month post-randomization variables to predict pump failure and SCD (categorized as "definite," "probable," and "possible"). Of 255 deaths, 197 were cardiovascular. There were 71 SCDs with a risk reduction by CRT of 0.47 (95% confidence interval 0.29-0.76; P = .002) with similar reductions in SCD classified as definite, probable, and possible. Univariate SCD predictors were 3-month HF status (mitral regurgitation [MR] severity, plasma brain natriuretic peptide [BNP], end-diastolic volume, and systolic blood pressure), whereas randomization to CRT decreased risk. Multivariate SCD predictors were randomization to CRT 0.56 (0.53-0.96, P = .035) and 3-month MR severity 1.82 (1.77-2.60, P = .0012). Univariate pump failure death predictors related to baseline HF state (quality of life score, interventricular mechanical delay, end-diastolic volume, plasma BNP, MR severity, and systolic pressure), whereas randomization to CRT and nonischemic cardiomyopathy decreased risk; multivariate predictors of pump failure death were baseline plasma BNP and systolic pressure and randomization to CRT. CONCLUSION: CRT decreased SCD in patients with systolic HF and ventricular dyssynchrony. SCD risk was increased with increased severity of MR (including the 3-month value for MR as a time-dependent covariate) and reduced by randomization to CRT. HF death was increased related to the level of systolic blood pressure, log BNP, and randomization to CRT. These results emphasize the importance and interdependence of HF severity to mortality from pump failure and SCD.