Prediction with measurement errors in finite populations

We address the problem of selecting the best linear unbiased predictor (BLUP) of the latent value (e.g., serum glucose fasting level) of sample subjects with heteroskedastic measurement errors. Using a simple example, we compare the usual mixed model BLUP to a similar predictor based on a mixed model framed in a finite population (FPMM) setup with two sources of variability, the first of which corresponds to simple random sampling and the second, to heteroskedastic measurement errors. Under this last approach, we show that when measurement errors are subject-specific, the BLUP shrinkage constants are based on a pooled measurement error variance as opposed to the individual ones generally considered for the usual mixed model BLUP. In contrast, when the heteroskedastic measurement errors are measurement condition-specific, the FPMM BLUP involves different shrinkage constants. We also show that in this setup, when measurement errors are subject-specific, the usual mixed model predictor is biased but has a smaller mean squared error than the FPMM BLUP which points to some difficulties in the interpretation of such predictors. (C) 2011 Elsevier By. All rights reserved.

Efeitos gênicos de caracteres associados à eficiência no uso de nitrogênio em milho

Os objetivos deste trabalho foram determinar o controle genético da eficiência no uso do nitrogênio (EUN), identificar a importância das eficiências na absorção (EAN) e na utilização (EUtN) na sua composição,  e quantificar relação entre produção de matéria seca da parte aérea (MPS) e do sistema radicular com a EUN e  com seus componentes. Foram avaliadas 41 combinações híbridas em duas disponibilidades de N: baixa (BN)  e alta (AN). Utilizou-se o delineamento de blocos ao acaso com duas repetições, em arranjo fatorial simples  (combinação híbrida x disponibilidade de N). As análises estatísticas foram realizadas por meio das equações  de modelos mistos. Correlações de elevada magnitude foram detectadas entre EAN e EUN, bem como entre  essas eficiências e a MPS, tanto em BN como em AN. Em ambas as disponibilidades de N, efeitos genéticos  aditivos apresentaram maior importância para os caracteres associados à EUN. Dessa forma, a seleção baseada  no desempenho individual de linhagens quanto à MPS pode possibilitar a obtenção de genótipos com alta  EUN. Independentemente da disponibilidade de N, a EAN é o componente mais importante da EUN.

The effect of implant design on insertion torque and immediate micromotion

Objectives: To evaluate the effect of insertion torque on micromotion to a lateral force in three different implant designs. Material and methods: Thirty-six implants with identical thread design, but different cutting groove design were divided in three groups: (1) non-fluted (no cutting groove, solid screw-form); (2) fluted (901 cut at the apex, tap design); and (3) Blossomt (Patent pending) (non-fluted with engineered trimmed thread design). The implants were screwed into polyurethane foam blocks and the insertion torque was recorded after each turn of 901 by a digital torque gauge. Controlled lateral loads of 10N followed by increments of 5 up to 100N were sequentially applied by a digital force gauge on a titanium abutment. Statistical comparison was performed with two-way mixed model ANOVA that evaluated implant design group, linear effects of turns and displacement loads, and their interaction. Results: While insertion torque increased as a function of number of turns for each design, the slope and final values increased (Po0.001) progressively from the Blossomt to the fluted to the non-fluted design (M +/- standard deviation [SD] = 64.1 +/- 26.8, 139.4 +/- 17.2, and 205.23 +/- 24.3 Ncm, respectively). While a linear relationship between horizontal displacement and lateral force was observed for each design, the slope and maximal displacement increased (Po0.001) progressively from the Blossomt to the fluted to the non-fluted design (M +/- SD 530 +/- 57.7, 585.9 +/- 82.4, and 782.33 +/- 269.4 mm, respectively). There was negligible to moderate levels of association between insertion torque and lateral displacement in the Blossomt, fluted and non-fluted design groups, respectively. Conclusion: Insertion torque was reduced in implant macrodesigns that incorporated cutting edges, and lesser insertion torque was generally associated with decreased micromovement. However, insertion torque and micromotion were unrelated within implant designs, particularly for those designs showing the least insertion torque.

Evaluation of a nanometer roughness scale resorbable media-processed surface: a study in dogs

Objectives: This study compared the biomechanical fixation and bone-to-implant contact (BIC) of implants with different surfaces treatment (experimental resorbable blasting media-processed nanometer roughness scale surface, and control dual acid-etched) in a dog model. Material and methods: Surface characterization was made in six implants by means of scanning electron microscopic imaging, atomic force microscopy to evaluate roughness parameters, and X-ray photoelectron spectroscopy (XPS) for chemical assessment. The animal model comprised the bilateral placement of control (n = 24) and experimental surface (n = 24) implants along the proximal tibiae of six mongrel dogs, which remained in place for 2 or 4 weeks. Half of the specimens were biomechanically tested (torque), and the other half was subjected to histomorphologic/ morphometric evaluation. BIC and resistance to failure measures were each evaluated as a function of time and surface treatment in a mixed model ANOVA. Results: Surface texturing was significantly higher for the experimental compared with the control surface. The survey XPS spectra detected O, C, Al, and Ti at the control group, and Ca (similar to 0.2-0.9%) and P (similar to 1.7-4.1%) besides O, C, Al, and Ti at experimental surfaces. While no statistical difference in BIC was found between experimental and control surfaces or between 2 and 4 weeks in vivo, both longer time and use of experimental surface significantly increased resistance to failure. Conclusions: The experimental surface resulted in enhanced biomechanical fixation but comparable BIC relative to control, suggesting higher bone mechanical properties around the experimental implants.

Diminished Mycophenolic Acid Exposure Caused by Omeprazole May Be Clinically Relevant in the First Week Posttransplantation

Background: Some studies have reported a decreased absorption of mycophenolic acid (MPA) from mycophenolate mofetil (MMF) in renal transplanted (RTx) patients under proton-pump inhibitors (PPIs). There is still a lack of information regarding (1) whether this effect occurs when MMF is administered with either tacrolimus or cyclosporine A [calcineurin inhibitors (CNIs)], (2) whether the effect has the same amplitude during the first year after RTx, and finally (3) whether this decrease in exposure is clinically relevant. Methods: We retrospectively analyzed the omeprazole effect in 348 12-hour pharmacokinetic samplings [area under the curve (AUC) 0-12h] performed on days 7, 14, 30, 60, 180, and 360 after RTx in 77 patients who participated in previous trials. Results: For all periods, the groups with and without PPI did not differ in all variables. By mixed-model analysis of variance, PPI reduced the MPA AUC(0-12h) (P < 0.0008) in the patients under both CNIs mainly due to decreased absorption (P = 0.049). In the tacrolimus group, a lower exposure seemed also due to a decreased MPA reabsorption at 10-12 hours. The PPI effect remains throughout the first year but was clinically more important on day 7. By Cox analysis, the use of PPI was associated with a 25% less chance of being adequately exposed to MPA (95% confidence interval 0.58-0.99, P = 0.04). Similarly, the number of patients underexposed to MPA (AUC < 30 ng.h/mL) was higher at most periods in the PPI group but again not statistically significant. Conclusions: These data indicate that PPI decreases the MPA exposure when associated with both CNIs but particularly in the first week after RTx. In this period, the MMF dose should be increased. This effect lasts throughout the first year but does not seem to be clinically relevant after the first week.

Multilevel Analysis in Household Survey: An Application to Health Condition Data

The aim of this thesis is to apply multilevel regression model in context of household surveys. Hierarchical structure in this type of data is characterized by many small groups. In last years comparative and multilevel analysis in the field of perceived health have grown in size. The purpose of this thesis is to develop a multilevel analysis with three level of hierarchy for Physical Component Summary outcome to: evaluate magnitude of within and between variance at each level (individual, household and municipality); explore which covariates affect on perceived physical health at each level; compare model-based and design-based approach in order to establish informativeness of sampling design; estimate a quantile regression for hierarchical data. The target population are the Italian residents aged 18 years and older. Our study shows a high degree of homogeneity within level 1 units belonging from the same group, with an intraclass correlation of 27% in a level-2 null model. Almost all variance is explained by level 1 covariates. In fact, in our model the explanatory variables having more impact on the outcome are disability, unable to work, age and chronic diseases (18 pathologies). An additional analysis are performed by using novel procedure of analysis :"Linear Quantile Mixed Model", named "Multilevel Linear Quantile Regression", estimate. This give us the possibility to describe more generally the conditional distribution of the response through the estimation of its quantiles, while accounting for the dependence among the observations. This has represented a great advantage of our models with respect to classic multilevel regression. The median regression with random effects reveals to be more efficient than the mean regression in representation of the outcome central tendency. A more detailed analysis of the conditional distribution of the response on other quantiles highlighted a differential effect of some covariate along the distribution.

On the car sequencing problem: analysis and solution methods

This work deals with the car sequencing (CS) problem, a combinatorial optimization problem for sequencing mixed-model assembly lines. The aim is to find a production sequence for different variants of a common base product, such that work overload of the respective line operators is avoided or minimized. The variants are distinguished by certain options (e.g., sun roof yes/no) and, therefore, require different processing times at the stations of the line. CS introduces a so-called sequencing rule H:N for each option, which restricts the occurrence of this option to at most H in any N consecutive variants. It seeks for a sequence that leads to no or a minimum number of sequencing rule violations. In this work, CS’ suitability for workload-oriented sequencing is analyzed. Therefore, its solution quality is compared in experiments to the related mixed-model sequencing problem. A new sequencing rule generation approach as well as a new lower bound for the problem are presented. Different exact and heuristic solution methods for CS are developed and their efficiency is shown in experiments. Furthermore, CS is adjusted and applied to a resequencing problem with pull-off tables.

Counterbalancing effects of competition for resources and facilitation against grazing in alpine snowbed communities

Alpine snowbeds are habitats where the major limiting factors for plant growth are herbivory and a small time window for growth due to late snowmelt. Despite these limitations, snowbed vegetation usually forms a dense carpet of palatable plants due to favourable abiotic conditions for plant growth within the short growing season. These environmental characteristics make snowbeds particularly interesting to study the interplay of facilitation and competition. We hypothesised an interplay between resource competition and facilitation against herbivory. Further, we investigated whether these predicted neighbour effects were species-specific and/or dependent on ontogeny, and whether the balance of positive and negative plant–plant interactions shifted along a snowmelt gradient. We determined the neighbour effects by means of neighbour removal experiments along the snowmelt gradient, and linear mixed model analyses. The results showed that the effects of neighbour removal were weak but generally consistent among species and snowmelt dates, and depended on whether biomass production or survival was considered. Higher total biomass and increased fruiting in removal plots indicated that plants competed for nutrients, water, and light, thereby supporting the hypothesis of prevailing competition for resources in snowbeds. However, the presence of neighbours reduced herbivory and thereby also facilitated survival. For plant growth the facilitative effects against herbivores in snowbeds counterbalanced competition for resources, leading to a weak negative net effect. Overall the neighbour effects were not species-specific and did not change with snowmelt date. Our finding of counterbalancing effects of competition and facilitation within a plant community is of special theoretical value for species distribution models and can explain the success of models that give primary importance to abiotic factors and tend to overlook interrelations between biotic and abiotic effects on plants.

Coping personality type and environmental enrichment affect aggression at weaning in pigs

This study investigated the effects of different environmental treatments and personality types on aggression at mixing of newly weaned domestic piglets. From birth to weaning, 16 litters were housed with their dams in either barren (B) or larger, substrate-enriched (E) environments. At 15 days old, piglets were classified as 'high' (HR) or low resistant' (LR) in a manual restraint test (backtest), which is thought to identify proactive (HR) and reactive (LR) stress coping strategies that may reflect different personality types. At 30 days old, 128 piglets were weaned, relocated and mixed into 32 pens comprising two HR and two LR unfamiliar pigs, balanced for sex and weaning weight. Eight B and eight E groups changed environmental condition whereas the others remained in the same type of environment. Number and duration of fights. fight outcomes and unilateral fighting were scored for 5 h post-mixing and skin lesions were counted before and 5 h, 1 day and 2 days after mixing. On the day following weaning, fighting and also exploratory and oral manipulative behaviours were measured for 6 h. Generalized Linear Mixed Model analyses suggested interactions between pre-weaning environment, post-weaning environment and personality type. Overall, pre-weaning E pigs had longer fights at weaning and mixing (P=0.01) and fought for longer on the next day (P=0.02) than pre-weaning B pigs, and inflicted more skin lesions (P=0.02). Post-weaning enrichment did not affect fighting at mixing but reduced the time spent fighting the next day (P=0.03). Personality had subtle and environment-dependent effects on fighting, and influenced the "structure" rather than the amount of aggressive behaviour. HR pigs, for instance, bullied (i.e. chased surrendering pigs) more often (P=0.009) and their fighting behaviour was less affected by their relative body weight than that of LR pigs. Post-weaning E pigs showed relatively higher levels of exploratory behaviour (P=0.02) and less oral manipulative behaviour (P=0.04) than post-weaning B pigs. In particular, switching from a good quality environment (E) to a worse quality one (B) at weaning decreased exploratory behaviour on the next day, especially for LR pigs, who also tended to fight with and orally manipulate their pen mates more in that condition, and seemed to be more affected by a deterioration of the environment. Overall, pre-weaning enrichment increased aggression after weaning whereas post-weaning enrichment reduced it, and personality type related to some aspects of fighting behaviour. (C) 2011 Elsevier B.V. All rights reserved.

Mapping genetic variants associated with beta-adrenergic responses in inbred mice

β-blockers and β-agonists are primarily used to treat cardiovascular diseases. Inter-individual variability in response to both drug classes is well recognized, yet the identity and relative contribution of the genetic players involved are poorly understood. This work is the first genome-wide association study (GWAS) addressing the values and susceptibility of cardiovascular-related traits to a selective β(1)-blocker, Atenolol (ate), and a β-agonist, Isoproterenol (iso). The phenotypic dataset consisted of 27 highly heritable traits, each measured across 22 inbred mouse strains and four pharmacological conditions. The genotypic panel comprised 79922 informative SNPs of the mouse HapMap resource. Associations were mapped by Efficient Mixed Model Association (EMMA), a method that corrects for the population structure and genetic relatedness of the various strains. A total of 205 separate genome-wide scans were analyzed. The most significant hits include three candidate loci related to cardiac and body weight, three loci for electrocardiographic (ECG) values, two loci for the susceptibility of atrial weight index to iso, four loci for the susceptibility of systolic blood pressure (SBP) to perturbations of the β-adrenergic system, and one locus for the responsiveness of QTc (p<10(-8)). An additional 60 loci were suggestive for one or the other of the 27 traits, while 46 others were suggestive for one or the other drug effects (p<10(-6)). Most hits tagged unexpected regions, yet at least two loci for the susceptibility of SBP to β-adrenergic drugs pointed at members of the hypothalamic-pituitary-thyroid axis. Loci for cardiac-related traits were preferentially enriched in genes expressed in the heart, while 23% of the testable loci were replicated with datasets of the Mouse Phenome Database (MPD). Altogether these data and validation tests indicate that the mapped loci are relevant to the traits and responses studied.

A prospective study of the impact of air pollution on respiratory symptoms and infections in infants

Rationale: There is increasing evidence that short-term exposure to air pollution has a detrimental effect on respiratory health, but data from healthy populations, particularly infants, are scarce. Objectives: To assess the association of air pollution with frequency and severity of respiratory symptoms and infections measured weekly in healthy infants. Methods: In a prospective birth cohort of 366 infants of unselected mothers, respiratory health was assessed weekly by telephone interviews during the first year of life (19,106 total observations). Daily mean levels of particulate matter (PM10), nitrogen dioxide (NO2), and ozone (O3) were obtained from local monitoring stations. We determined the association of the preceding week's pollutant levels with symptom scores and respiratory tract infections using a generalized additive mixed model with an autoregressive component. In addition, we assessed whether neonatal lung function influences this association and whether duration of infectious episodes differed between weeks with normal PM10 and weeks with elevated levels. Measurements and Main Results: We found a significant association between air pollution and respiratory symptoms, particularly in the week after respiratory tract infections (risk ratio, 1.13 [1.02-1.24] per 10 μg/m(3) PM10 levels) and in infants with premorbid lung function. During times of elevated PM10 (>33.3 μg/m(3)), duration of respiratory tract infections increased by 20% (95% confidence interval, 2-42%). Conclusions: Exposure to even moderate levels of air pollution was associated with increased respiratory symptoms in healthy infants. Particularly in infants with premorbid lung function and inflammation, air pollution contributed to longer duration of infectious episodes with a potentially large socioeconomic impact.

Clinical and molecular characterization of the potential CF disease modifier syntaxin 1A

Cystic fibrosis (CF) is caused by mutations in the CF transmembrane conductance regulator gene (CFTR). Disease severity in CF varies greatly, and sibling studies strongly indicate that genes other than CFTR modify disease outcome. Syntaxin 1A (STX1A) has been reported as a negative regulator of CFTR and other ion channels. We hypothesized that STX1A variants act as a CF modifier by influencing the remaining function of mutated CFTR. We identified STX1A variants by genomic resequencing patients from the Bernese CF Patient Data Registry and applied linear mixed model analysis to establish genotype-phenotype correlations, revealing STX1A rs4363087 (c.467-38A>G) to significantly influence lung function. The same STX1A risk allele was recognized in the European CF Twin and Sibling Study (P=0.0027), demonstrating that the genotype-phenotype association of STX1A to CF disease severity is robust enough to allow replication in two independent CF populations. rs4363087 is in linkage disequilibrium to the exonic variant rs2228607 (c.204C>T). Considering that neither rs4363087 nor rs2228607 changes the amino-acid sequence of STX1A, we investigated their effects on mRNA level. We show that rs2228607 reinforces aberrant splicing of STX1A mRNA, leading to nonsense-mediated mRNA decay. In conclusion, we demonstrate the clinical relevance of STX1A variants in CF, and evidence the functional relevance of STX1A variant rs2228607 at molecular level. Our findings show that genes interacting with CFTR can modify CF disease progression.European Journal of Human Genetics advance online publication, 10 April 2013; doi:10.1038/ejhg.2013.57.

Computational Techniques for Spatial Logistic Regression with Large Datasets

In epidemiological work, outcomes are frequently non-normal, sample sizes may be large, and effects are often small. To relate health outcomes to geographic risk factors, fast and powerful methods for fitting spatial models, particularly for non-normal data, are required. We focus on binary outcomes, with the risk surface a smooth function of space. We compare penalized likelihood models, including the penalized quasi-likelihood (PQL) approach, and Bayesian models based on fit, speed, and ease of implementation. A Bayesian model using a spectral basis representation of the spatial surface provides the best tradeoff of sensitivity and specificity in simulations, detecting real spatial features while limiting overfitting and being more efficient computationally than other Bayesian approaches. One of the contributions of this work is further development of this underused representation. The spectral basis model outperforms the penalized likelihood methods, which are prone to overfitting, but is slower to fit and not as easily implemented. Conclusions based on a real dataset of cancer cases in Taiwan are similar albeit less conclusive with respect to comparing the approaches. The success of the spectral basis with binary data and similar results with count data suggest that it may be generally useful in spatial models and more complicated hierarchical models.

FAST ADAPTIVE PENALIZED SPLINES

This paper proposes a numerically simple routine for locally adaptive smoothing. The locally heterogeneous regression function is modelled as a penalized spline with a smoothly varying smoothing parameter modelled as another penalized spline. This is being formulated as hierarchical mixed model, with spline coe±cients following a normal distribution, which by itself has a smooth structure over the variances. The modelling exercise is in line with Baladandayuthapani, Mallick & Carroll (2005) or Crainiceanu, Ruppert & Carroll (2006). But in contrast to these papers Laplace's method is used for estimation based on the marginal likelihood. This is numerically simple and fast and provides satisfactory results quickly. We also extend the idea to spatial smoothing and smoothing in the presence of non normal response.