Applying Least Absolute Shrinkage Selection Operator and Akaike Information Criterion Analysis to Find the Best Multiple Linear Regression Models between Climate Indices and Components of Cow’s Milk

This study focuses on multiple linear regression models relating six climate indices (temperature humidity THI, environmental stress ESI, equivalent temperature index ETI, heat load HLI, modified HLI (HLI new), and respiratory rate predictor RRP) with three main components of cow’s milk (yield, fat, and protein) for cows in Iran. The least absolute shrinkage selection operator (LASSO) and the Akaike information criterion (AIC) techniques are applied to select the best model for milk predictands with the smallest number of climate predictors. Uncertainty estimation is employed by applying bootstrapping through resampling. Cross validation is used to avoid over-fitting. Climatic parameters are calculated from the NASA-MERRA global atmospheric reanalysis. Milk data for the months from April to September, 2002 to 2010 are used. The best linear regression models are found in spring between milk yield as the predictand and THI, ESI, ETI, HLI, and RRP as predictors with p-value < 0.001 and R2 (0.50, 0.49) respectively. In summer, milk yield with independent variables of THI, ETI, and ESI show the highest relation (p-value < 0.001) with R2 (0.69). For fat and protein the results are only marginal. This method is suggested for the impact studies of climate variability/change on agriculture and food science fields when short-time series or data with large uncertainty are available.

On the effect of confounding in linear regression models: an approach based on the theory of quadratic forms

The main topic of this thesis is confounding in linear regression models. It arises when a relationship between an observed process, the covariate, and an outcome process, the response, is influenced by an unmeasured process, the confounder, associated with both. Consequently, the estimators for the regression coefficients of the measured covariates might be severely biased, less efficient and characterized by misleading interpretations. Confounding is an issue when the primary target of the work is the estimation of the regression parameters. The central point of the dissertation is the evaluation of the sampling properties of parameter estimators. This work aims to extend the spatial confounding framework to general structured settings and to understand the behaviour of confounding as a function of the data generating process structure parameters in several scenarios focusing on the joint covariate-confounder structure. In line with the spatial statistics literature, our purpose is to quantify the sampling properties of the regression coefficient estimators and, in turn, to identify the most prominent quantities depending on the generative mechanism impacting confounding. Once the sampling properties of the estimator conditionally on the covariate process are derived as ratios of dependent quadratic forms in Gaussian random variables, we provide an analytic expression of the marginal sampling properties of the estimator using Carlson’s R function. Additionally, we propose a representative quantity for the magnitude of confounding as a proxy of the bias, its first-order Laplace approximation. To conclude, we work under several frameworks considering spatial and temporal data with specific assumptions regarding the covariance and cross-covariance functions used to generate the processes involved. This study allows us to claim that the variability of the confounder-covariate interaction and of the covariate plays the most relevant role in determining the principal marker of the magnitude of confounding.

Analysis Of The Fluorescence Of Body Fluids On Different Surfaces And Times.

The use of screening techniques, such as an alternative light source (ALS), is important for finding biological evidence at a crime scene. The objective of this study was to evaluate whether biological fluid (blood, semen, saliva, and urine) deposited on different surfaces changes as a function of the age of the sample. Stains were illuminated with a Megamaxx™ ALS System and photographed with a Canon EOS Utility™ camera. Adobe Photoshop™ was utilized to prepare photographs for analysis, and then ImageJ™ was used to record the brightness values of pixels in the images. Data were submitted to analysis of variance using a generalized linear mixed model with two fixed effects (surface and fluid). Time was treated as a random effect (through repeated measures) with a first-order autoregressive covariance structure. Means of significant effects were compared by the Tukey test. The fluorescence of the analyzed biological material varied depending on the age of the sample. Fluorescence was lower when the samples were moist. Fluorescence remained constant when the sample was dry, up to the maximum period analyzed (60 days), independent of the substrate on which the fluid was deposited, showing the novelty of this study. Therefore, the forensic expert can detect biological fluids at the crime scene using an ALS even several days after a crime has occurred.

Modelling the Distribution of Ischaemic Stroke-Specific Survival Time Using an EM-based Mixture Approach with Random Effects Adjustment

A two-component survival mixture model is proposed to analyse a set of ischaemic stroke-specific mortality data. The survival experience of stroke patients after index stroke may be described by a subpopulation of patients in the acute condition and another subpopulation of patients in the chronic phase. To adjust for the inherent correlation of observations due to random hospital effects, a mixture model of two survival functions with random effects is formulated. Assuming a Weibull hazard in both components, an EM algorithm is developed for the estimation of fixed effect parameters and variance components. A simulation study is conducted to assess the performance of the two-component survival mixture model estimators. Simulation results confirm the applicability of the proposed model in a small sample setting. Copyright (C) 2004 John Wiley Sons, Ltd.

Marginal adaptation in vitro and clinical outcome of Class V restorations

OBJECTIVES: We examined the correlation between the quantitative margin analysis of two laboratory test methods (Berlin, Zurich) and the clinical outcome in Class V restorations. METHODS: Prospective clinical studies with an observation period of at least 18 months were searched in the literature, for which laboratory data were also available. The clinical outcome variables were retention loss, marginal discoloration, detectable margins and secondary caries. Forty-four clinical studies matched the inclusion criteria, including 34 adhesive systems for which laboratory data were also present. For both laboratory test methods and the clinical studies, an index was formulated to better compare the in vitro and in vivo results. Linear mixed models which included a random study effect were calculated. As most clinical data were available for 12 and 24 months, the main analysis was restricted to these recall intervals. RESULTS: The comparative analysis revealed a weak correlation between the clinical index and both in vitro indices. The correlation was statistically significant for the Berlin method but not for the Zurich method and only present if studies were compared which used the same composite in the in vitro and in vivo study. When defining specific cut-off values, the prognosis for the good clinical performance of an adhesive system based on in vitro results was 78% (Berlin) or 100% (Zurich). For poor performance it was 67% and 60%, respectively. No correlation was found between both in vitro methods. SIGNIFICANCE: The surrogate parameter "marginal adaptation" of restorations placed in extracted teeth has a mediocre value to predict the clinical performance of an adhesive system in cervical cavities. The composite is an important factor for a successful prediction. The comparison between in vitro/in vivo is sometimes hampered by the great variability of clinical results on the same adhesive system.

Genomics of the divergence continuum in an African plant biodiversity hotspot, I: drivers of population divergence in Restio capensis (Restionaceae).

Understanding the drivers of population divergence, speciation and species persistence is of great interest to molecular ecology, especially for species-rich radiations inhabiting the world's biodiversity hotspots. The toolbox of population genomics holds great promise for addressing these key issues, especially if genomic data are analysed within a spatially and ecologically explicit context. We have studied the earliest stages of the divergence continuum in the Restionaceae, a species-rich and ecologically important plant family of the Cape Floristic Region (CFR) of South Africa, using the widespread CFR endemic Restio capensis (L.) H.P. Linder & C.R. Hardy as an example. We studied diverging populations of this morphotaxon for plastid DNA sequences and >14 400 nuclear DNA polymorphisms from Restriction site Associated DNA (RAD) sequencing and analysed the results jointly with spatial, climatic and phytogeographic data, using a Bayesian generalized linear mixed modelling (GLMM) approach. The results indicate that population divergence across the extreme environmental mosaic of the CFR is mostly driven by isolation by environment (IBE) rather than isolation by distance (IBD) for both neutral and non-neutral markers, consistent with genome hitchhiking or coupling effects during early stages of divergence. Mixed modelling of plastid DNA and single divergent outlier loci from a Bayesian genome scan confirmed the predominant role of climate and pointed to additional drivers of divergence, such as drift and ecological agents of selection captured by phytogeographic zones. Our study demonstrates the usefulness of population genomics for disentangling the effects of IBD and IBE along the divergence continuum often found in species radiations across heterogeneous ecological landscapes.

Clinical performance of cervical restorations: a meta-analysis.

OBJECTIVES: To carry out a meta-analysis in order to assess the influencing factors on retention loss and marginal discoloration of cervical restorations made of composites and glass ionomer (derivates). METHODS: The literature was searched for prospective clinical studies on cervical restorations with an observation period of at least 18 months. RESULTS: Fifty clinical studies involving 40 adhesive systems matched the inclusion criteria. On average, 10% of the cervical fillings were lost and 24% exhibited marginal discoloration after 3 years. The variability ranged from 0% to 50% for retention loss and from 0% to 74% for marginal discoloration. Hardly any secondary caries was detected. When linear mixed models with a study and experiment effect were used, the analysis revealed that the adhesive/restorative class had the most significant influence, with 2-step self-etching adhesive systems performing best and 1-step self-etching adhesive systems performing worst; 3-step etch-and-rinse systems, glass ionomers/resin-modified glass ionomers, 2-step etch-and-rinse systems and polyacid-modified resin composites were ranked in between. Restorations placed in teeth whose dentin/enamel had been prepared/roughened showed a statistically significant higher retention rate than those placed in teeth with unprepared dentin (p<0.05). Beveling of the enamel and the type of isolation used (rubberdam/cotton rolls) had no significant influence. SIGNIFICANCE: The clinical performance of cervical restorations is significantly influenced by the type of adhesive system used and/or the adhesive class to which the system belonged and whether the dentin/enamel is prepared or not. 2-Step self-etching- and 3-step etch&rinse systems shall be chosen over 1-step self-etching systems and glass ionomer derivates. The dentin (and enamel) surface shall be roughened before placement of the restoration.

Macronutrient composition of the diet and prospective weight change in participants of the EPIC-PANACEA study.

BACKGROUND The effect of the macronutrient composition of the usual diet on long term weight maintenance remains controversial. METHODS 373,803 subjects aged 25-70 years were recruited in 10 European countries (1992-2000) in the PANACEA project of the EPIC cohort. Diet was assessed at baseline using country-specific validated questionnaires and weight and height were measured at baseline and self-reported at follow-up in most centers. The association between weight change after 5 years of follow-up and the iso-energetic replacement of 5% of energy from one macronutrient by 5% of energy from another macronutrient was assessed using multivariate linear mixed-models. The risk of becoming overweight or obese after 5 years was investigated using multivariate Poisson regressions stratified according to initial Body Mass Index. RESULTS A higher proportion of energy from fat at the expense of carbohydrates was not significantly associated with weight change after 5 years. However, a higher proportion of energy from protein at the expense of fat was positively associated with weight gain. A higher proportion of energy from protein at the expense of carbohydrates was also positively associated with weight gain, especially when carbohydrates were rich in fibre. The association between percentage of energy from protein and weight change was slightly stronger in overweight participants, former smokers, participants ≥60 years old, participants underreporting their energy intake and participants with a prudent dietary pattern. Compared to diets with no more than 14% of energy from protein, diets with more than 22% of energy from protein were associated with a 23-24% higher risk of becoming overweight or obese in normal weight and overweight subjects at baseline. CONCLUSION Our results show that participants consuming an amount of protein above the protein intake recommended by the American Diabetes Association may experience a higher risk of becoming overweight or obese during adult life.

Characterisation of Culex quinquefasciatus (Diptera: Culicidae) larval habitats at ground level and temporal fluctuations of larval abundance in Córdoba, Argentina

The aims of this study were to characterise the ground-level larval habitats of the mosquito Culex quinquefasciatus, to determine the relationships between habitat characteristics and larval abundance and to examine seasonal larval-stage variations in Córdoba city. Every two weeks for two years, 15 larval habitats (natural and artificial water bodies, including shallow wells, drains, retention ponds, canals and ditches) were visited and sampled for larval mosquitoes. Data regarding the water depth, temperature and pH, permanence, the presence of aquatic vegetation and the density of collected mosquito larvae were recorded. Data on the average air temperatures and accumulated precipitation during the 15 days prior to each sampling date were also obtained. Cx. quinquefasciatus larvae were collected throughout the study period and were generally most abundant in the summer season. Generalised linear mixed models indicated the average air temperature and presence of dicotyledonous aquatic vegetation as variables that served as important predictors of larval densities. Additionally, permanent breeding sites supported high larval densities. In Córdoba city and possibly in other highly populated cities at the same latitude with the same environmental conditions, control programs should focus on permanent larval habitats with aquatic vegetation during the early spring, when the Cx. quinquefasciatus population begins to increase.

Association of alcohol consumption and HIV surrogate markers in participants of the swiss HIV cohort study.

BACKGROUND: Alcohol consumption may affect the course of HIV infection and/or antiretroviral therapy (ART). The authors investigated the association between self-reported alcohol consumption and HIV surrogate markers in both treated and untreated individuals. DESIGN: Prospective cohort study. METHODS: Over a 7-year period, the authors analyzed 2 groups of individuals in the Swiss HIV Cohort Study: (1) ART-naïve individuals remaining off ART and (2) individuals initiating first ART. For individuals initiating first ART, time-dependent Cox proportional hazards models were used to assess the association between alcohol consumption, virological failure, and ART interruption. For both groups, trajectories of log-transformed CD4 cell counts were analyzed using linear mixed models with repeated measures. RESULTS: The authors included 2982 individuals initiating first ART and 2085 ART naives. In individuals initiating first ART, 241 (8%) experienced virological failure. Alcohol consumption was not associated with virological failure. ART interruption was noted in 449 (15%) individuals and was more prevalent in severe compared with none/light health risk drinkers [hazard ratio: 2.24, 95% confidence interval: 1.42 to 3.52]. The association remained significant even after adjusting for nonadherence. The authors did not find an association between alcohol consumption and change in CD4 cell count over time in either group. CONCLUSIONS: No effect of alcohol consumption on either virological failure or CD4 cell count in both groups of ART-initiating and ART-naive individuals was found. However, severe drinkers were more likely to interrupt ART. Efforts on ART continuation should be especially implemented in individuals reporting high alcohol consumption.

The effect of antiretroviral treatment of different durations in primary HIV infection.

OBJECTIVES: To compare immunological, virological and clinical outcomes in persons initiating combination antiretroviral therapy (cART of different durations within 6 months of seroconversion (early treated) with those who deferred therapy (deferred group). DESIGN: CD4 cell and HIV-RNA measurements for 'early treated' individuals following treatment cessation were compared with the corresponding ART-free period for the 'deferred' group using piecewise linear mixed models. Individuals identified during primary HIV infection were included if they seroconverted from 1st January 1996 and were at least 15 years of age at seroconversion. Those with at least 2 CD4 less than 350 cells/microl or AIDS within the first 6 months following seroconversion were excluded. RESULTS: Of 348 'early treated' patients, 147 stopped cART following treatment for at least 6 (n = 38), more than 6-12 (n = 40) or more than 12 months (n = 69). CD4 cell loss was steeper for the first 6 months following cART cessation, but subsequent loss rate was similar to the 'deferred' group (n = 675, P = 0.26). Although those treated for more than 12 months appeared to maintain higher CD4 cell counts following cART cessation, those treated for 12 months or less had CD4 cell counts 6 months after cessation comparable to those in the 'deferred' group. There was no difference in HIV-RNA set points between the 'early' and 'deferred' groups (P = 0.57). AIDS rates were similar but death rates, mainly due to non-AIDS causes, were higher in the 'deferred' group (P = 0.05). CONCLUSION: Transient cART, initiated within 6 months of seroconversion, seems to have no effect on viral load set point and limited beneficial effect on CD4 cell levels in individuals treated for more than 12 months. Its long-term effects remain inconclusive and need further investigation.

Wear of two artificial tooth materials in vivo: a 12-month pilot study.

STATEMENT OF PROBLEM: Wear of methacrylate artificial teeth resulting in vertical loss is a problem for both dentists and patients. PURPOSE: The purpose of this study was to quantify wear of artificial teeth in vivo and to relate it to subject and tooth variables. MATERIAL AND METHODS: Twenty-eight subjects treated with complete dentures received 2 artificial tooth materials (polymethyl methacrylate (PMMA)/double-cross linked PMMA fillers; 35%/59% (SR Antaris DCL, SR Postaris DCL); experimental 48%/46%). At baseline and after 12 months, impressions of the dentures were poured with improved stone. After laser scanning, the casts were superimposed and matched. Maximal vertical loss (mm) and volumetric loss (mm(3)) were calculated for each tooth and log-transformed to reduce variability. Volumetric loss was related to the occlusally active surface area. Linear mixed models were used to study the influence of the factors jaw, tooth, and material on adjusted (residual) wear values (alpha=.05). RESULTS: Due to drop outs (n=5) and unmatchable casts (n=3), 69% of all teeth were analyzed. Volumetric loss had a strong linear relationship to surface area (P<.001); this was less pronounced for vertical loss (P=.004). The factor showing the highest influence was the subject. Wear was tooth dependent (increasing from incisors to molars). However, these differences diminished once the wear rates were adjusted for occlusal area, and only a few remained significant (anterior versus posterior maxillary teeth). Another influencing factor was the age of the subject. CONCLUSIONS: Clinical wear of artificial teeth is higher than previously measured or expected. The presented method of analyzing wear of artificial teeth using a laser-scanning device seemed suitable.

Patient- and therapy-related factors on the wear of denture teeth : results of a clinical trial

OBJECTIVE: When we examined a previously published prospective multi-center clinical trial in which complete denture-wearers were followed over a period of 2 years, we found that about 30% of the variability in the clinical wear data of denture teeth was due to unknown characteristics of the subjects. In the second part of the study, we try to identify which patient- and therapy-related factors may explain some of this variability. METHODS: The clinical wear data of denture teeth at different recall times (6, 12, 18, 24 months) in 89 subjects (at baseline) were correlated with the following parameters, which may all have an influence on the wear of denture teeth: age, gender, bruxism as reported by the subjects, number of prostheses used so far, time since last extraction, smoking, fit of dentures as judged by the subject and the clinician, average denture wearing time and wearing of denture during the night. To evaluate the influence of the different patient- and therapy-related variables, both a univariate analysis (one extra factor to the model) and a multivariate analysis were carried out using linear mixed models with the variable Log mean as the outcome. RESULTS: None of the patient- and therapy-related parameters showed a statistically significant effect on the wear of denture teeth. There was, however, a trend for women to show less wear compared to men and a trend of decreasing wear with increasing age. SIGNIFICANCE: Further research is required to identify the factors which are responsible for the high variability observed between the subjects regarding clinical wear data.

Risk of cardiovascular events and blood pressure control in hypertensive HIV-infected patients: Swiss HIV Cohort Study (SHCS).

BACKGROUND: Prevalence of hypertension in HIV infection is high, and information on blood pressure control in HIV-infected individuals is insufficient. We modeled blood pressure over time and the risk of cardiovascular events in hypertensive HIV-infected individuals. METHODS: All patients from the Swiss HIV Cohort Study with confirmed hypertension (systolic or diastolic blood pressure above 139 or 89 mm Hg on 2 consecutive visits and presence of at least 1 additional cardiovascular risk factor) between April 1, 2000 and March 31, 2011 were included. Patients with previous cardiovascular events, already on antihypertensive drugs, and pregnant women were excluded. Change in blood pressure over time was modeled using linear mixed models with repeated measurement. RESULTS: Hypertension was diagnosed in 2595 of 10,361 eligible patients. Of those, 869 initiated antihypertensive treatment. For patients treated for hypertension, we found a mean (95% confidence interval) decrease in systolic and diastolic blood pressure of -0.82 (-1.06 to -0.58) mm Hg and -0.89 (-1.05 to -0.73) mm Hg/yr, respectively. Factors associated with a decline in systolic blood pressure were baseline blood pressure, presence of chronic kidney disease, cardiovascular events, and the typical risk factors for cardiovascular disease. In patients with hypertension, increase in systolic blood pressure [(hazard ratio 1.18 (1.06 to 1.32) per 10 mm Hg increase], total cholesterol, smoking, age, and cumulative exposure to protease inhibitor-based and triple nucleoside regimens were associated with cardiovascular events. CONCLUSIONS: Insufficient control of hypertension was associated with increased risk of cardiovascular events indicating the need for improved management of hypertension in HIV-infected individuals.

Regression analysis of spatial data.

Many of the most interesting questions ecologists ask lead to analyses of spatial data. Yet, perhaps confused by the large number of statistical models and fitting methods available, many ecologists seem to believe this is best left to specialists. Here, we describe the issues that need consideration when analysing spatial data and illustrate these using simulation studies. Our comparative analysis involves using methods including generalized least squares, spatial filters, wavelet revised models, conditional autoregressive models and generalized additive mixed models to estimate regression coefficients from synthetic but realistic data sets, including some which violate standard regression assumptions. We assess the performance of each method using two measures and using statistical error rates for model selection. Methods that performed well included generalized least squares family of models and a Bayesian implementation of the conditional auto-regressive model. Ordinary least squares also performed adequately in the absence of model selection, but had poorly controlled Type I error rates and so did not show the improvements in performance under model selection when using the above methods. Removing large-scale spatial trends in the response led to poor performance. These are empirical results; hence extrapolation of these findings to other situations should be performed cautiously. Nevertheless, our simulation-based approach provides much stronger evidence for comparative analysis than assessments based on single or small numbers of data sets, and should be considered a necessary foundation for statements of this type in future.