Sampling, WLS, and Mixed Models

Mixed models may be defined with or without reference to sampling, and can be used to predict realized random effects, as when estimating the latent values of study subjects measured with response error. When the model is specified without reference to sampling, a simple mixed model includes two random variables, one stemming from an exchangeable distribution of latent values of study subjects and the other, from the study subjects` response error distributions. Positive probabilities are assigned to both potentially realizable responses and artificial responses that are not potentially realizable, resulting in artificial latent values. In contrast, finite population mixed models represent the two-stage process of sampling subjects and measuring their responses, where positive probabilities are only assigned to potentially realizable responses. A comparison of the estimators over the same potentially realizable responses indicates that the optimal linear mixed model estimator (the usual best linear unbiased predictor, BLUP) is often (but not always) more accurate than the comparable finite population mixed model estimator (the FPMM BLUP). We examine a simple example and provide the basis for a broader discussion of the role of conditioning, sampling, and model assumptions in developing inference.

Predicting random effects with an expanded finite population mixed model

Prediction of random effects is an important problem with expanding applications. In the simplest context, the problem corresponds to prediction of the latent value (the mean) of a realized cluster selected via two-stage sampling. Recently, Stanek and Singer [Predicting random effects from finite population clustered samples with response error. J. Amer. Statist. Assoc. 99, 119-130] developed best linear unbiased predictors (BLUP) under a finite population mixed model that outperform BLUPs from mixed models and superpopulation models. Their setup, however, does not allow for unequally sized clusters. To overcome this drawback, we consider an expanded finite population mixed model based on a larger set of random variables that span a higher dimensional space than those typically applied to such problems. We show that BLUPs for linear combinations of the realized cluster means derived under such a model have considerably smaller mean squared error (MSE) than those obtained from mixed models, superpopulation models, and finite population mixed models. We motivate our general approach by an example developed for two-stage cluster sampling and show that it faithfully captures the stochastic aspects of sampling in the problem. We also consider simulation studies to illustrate the increased accuracy of the BLUP obtained under the expanded finite population mixed model. (C) 2007 Elsevier B.V. All rights reserved.

The Relationship between Intraocular Pressure Reduction and Rates of Progressive Visual Field Loss in Eyes with Optic Disc Hemorrhage

Purpose: To evaluate rates of visual field progression in eyes with optic disc hemorrhages and the effect of intraocular pressure (IOP) reduction on these rates. Design: Observational cohort study. Participants: The study included 510 eyes of 348 patients with glaucoma who were recruited from the Diagnostic Innovations in Glaucoma Study (DIGS) and followed for an average of 8.2 years. Methods: Eyes were followed annually with clinical examination, standard automated perimetry visual fields, and optic disc stereophotographs. The presence of optic disc hemorrhages was determined on the basis of masked evaluation of optic disc stereophotographs. Evaluation of rates of visual field change during follow-up was performed using the visual field index (VFI). Main Outcome Measures: The evaluation of the effect of optic disc hemorrhages on rates of visual field progression was performed using random coefficient models. Estimates of rates of change for individual eyes were obtained by best linear unbiased prediction (BLUP). Results: During follow-up, 97 (19%) of the eyes had at least 1 episode of disc hemorrhage. The overall rate of VFI change in eyes with hemorrhages was significantly faster than in eyes without hemorrhages (-0.88%/year vs. -0.38%/year, respectively, P < 0.001). The difference in rates of visual field loss pre- and post-hemorrhage was significantly related to the reduction of IOP in the post-hemorrhage period compared with the pre-hemorrhage period (r = -0.61; P < 0.001). Each 1 mmHg of IOP reduction was associated with a difference of 0.31%/year in the rate of VFI change. Conclusions: There was a beneficial effect of treatment in slowing rates of progressive visual field loss in eyes with optic disc hemorrhage. Further research should elucidate the reasons why some patients with hemorrhages respond well to IOP reduction and others seem to continue to progress despite a significant reduction in IOP levels. Financial Disclosure(s): Proprietary or commercial disclosure may be found after the references. Ophthalmology 2010; 117: 2061-2066 (C) 2010 by the American Academy of Ophthalmology.

Peak torque and rate of torque development in elderly with and without fall history

Background Falls are one of the greatest concerns among the elderly A number of studies have described peak torque as one of the best fall-related predictor. No studies have comprehensively focused on the rate of torque development of the lower limb muscles among elderly fallers. Then, the aim of this study was to determine the relationship between muscle peak torque and rate of torque development of the lower limb joints in elderly with and without fall history It was also aimed to determine whether these parameters of muscle performance (i e, peak torque and rate of torque development) are related to the number of falls. Methods: Thirty-one women volunteered to participate in the study and were assigned in one of the groups according to the number of falls over the 12 months that preceded the present Then, participants with no fall history (Cl; n = 13; 67.6[7.5] years-old), one fall (GII; n = 8, 66 0[4 91 years-old) and two or more falls (GIII, n = 10; 67.8[8.8] years-old) performed a number of lower limb maximal isometric voluntary contractions from which peak torque and rate of torque development were quantified Findings. Primary outcomes indicated no peak torque differences between experimental groups in any lower limb joint. The rate of torque development of the knee flexor muscles observed in the non-fallers (Cl) was greater than that observed in the fallers (P < 0.05) and had a significant relationship with the number of falls (P < 0 05) Interpretation. The greater knee flexor muscles` rate of torque development found in the non-fallers in comparison to the fallers indicated that the ability of the elderly to rapidly reorganise the arrangement of the lower limb may play a significant role in allowing the elderly to recover balance after a trip. Thus, training stimulus aimed to improve the rate of torque development may be more beneficial to prevent falls among the elderly than other training stimulus, which are not specifically designed to improve the ability to rapidly produce large amounts of torque (C) 2010 Published by Elsevier Ltd

Anti-mullerian hormone is the best predictor of poor response in ICSI cycles of patients with endometriosis

Purpose: To correlate ovarian reserve (OR) markers with response in assisted reproduction techniques (ART) and determine their ability to predict poor response among patients with endometriosis (EDT). Methods: We evaluated ART cycles of 27 women with EDT and 50 with exclusive male factor. Basal follicle stimulating hormone (FSH) and anti-mullerian hormone (AMH) levels were determined. Ovarian response to gonadotropin stimulation was assessed and correlation coefficients calculated between the variables and reserve markers. Areas under the curve (AUC) determined ability of tests to predict poor response. Results: AMH was significantly correlated with response in both groups and it was the only marker with significant discriminative capacity to predict poor response among EDT (AUC = 0.842; 95% CI: 0.651-0.952) and control group (AUC = 0.869; 95% CI: 0.743-0.947). Conclusion: Infertile patients with endometriosis can benefit from the pre-therapeutic assessment of OR markers. However, regardless of disease presence, only AMH predicts poor response to stimulus.

Linear dimensional changes in plaster die models using different elastomeric materials

Dental impression is an important step in the preparation of prostheses since it provides the reproduction of anatomic and surface details of teeth and adjacent structures. The objective of this study was to evaluate the linear dimensional alterations in gypsum dies obtained with different elastomeric materials, using a resin coping impression technique with individual shells. A master cast made of stainless steel with fixed prosthesis characteristics with two prepared abutment teeth was used to obtain the impressions. References points (A, B, C, D, E and F) were recorded on the occlusal and buccal surfaces of abutments to register the distances. The impressions were obtained using the following materials: polyether, mercaptan-polysulfide, addition silicone, and condensation silicone. The transfer impressions were made with custom trays and an irreversible hydrocolloid material and were poured with type IV gypsum. The distances between identified points in gypsum dies were measured using an optical microscope and the results were statistically analyzed by ANOVA (p < 0.05) and Tukey's test. The mean of the distances were registered as follows: addition silicone (AB = 13.6 µm, CD=15.0 µm, EF = 14.6 µm, GH=15.2 µm), mercaptan-polysulfide (AB = 36.0 µm, CD = 36.0 µm, EF = 39.6 µm, GH = 40.6 µm), polyether (AB = 35.2 µm, CD = 35.6 µm, EF = 39.4 µm, GH = 41.4 µm) and condensation silicone (AB = 69.2 µm, CD = 71.0 µm, EF = 80.6 µm, GH = 81.2 µm). All of the measurements found in gypsum dies were compared to those of a master cast. The results demonstrated that the addition silicone provides the best stability of the compounds tested, followed by polyether, polysulfide and condensation silicone. No statistical differences were obtained between polyether and mercaptan-polysulfide materials.

Modelo de programação linear para otimização econômica do projeto de irrigação Baixo Acaraú - CE

The present work had as objective uses a model of lineal programming algorithm to optimize the use of the water in the District of Irrigation Baixo Acarau-CE proposing the best combination of crop types and areas established of 8,0 ha. The model aim maximize the net benefit of small farmer, incorporating the constraints in water and land availability, and constraints on the market. Considering crop types and the constraints, the study lead to the following conclusions: 1. The water availability in the District was not a limiting resources, while all available land was assigned in six of the seven cultivation plans analyzed. Furthermore, water availability was a restrictive factor as compared with land only when its availability was made to reduce to 60% of its actual value; 2. The combination of soursop and melon plants was the one that presented the largest net benefit, corresponding to R$ 5,250.00/ha/yr. The planting area for each crop made up to 50% of the area of the plot; 3. The plan that suggests the substitution of the cultivation of the soursop, since a decrease in annual net revenue of 5.87%. However, the plan that contemplates the simultaneous substitution of both soursop and melon produced the lowest liquid revenue, with reduction of 33.8%.

Adaptive Carrier Tracking for Mars to Earth Communications During Entry, Descent, and Landing

We propose a robust and low complexity scheme to estimate and track carrier frequency from signals traveling under low signal-to-noise ratio (SNR) conditions in highly nonstationary channels. These scenarios arise in planetary exploration missions subject to high dynamics, such as the Mars exploration rover missions. The method comprises a bank of adaptive linear predictors (ALP) supervised by a convex combiner that dynamically aggregates the individual predictors. The adaptive combination is able to outperform the best individual estimator in the set, which leads to a universal scheme for frequency estimation and tracking. A simple technique for bias compensation considerably improves the ALP performance. It is also shown that retrieval of frequency content by a fast Fourier transform (FFT)-search method, instead of only inspecting the angle of a particular root of the error predictor filter, enhances performance, particularly at very low SNR levels. Simple techniques that enforce frequency continuity improve further the overall performance. In summary we illustrate by extensive simulations that adaptive linear prediction methods render a robust and competitive frequency tracking technique.

Development of an IgG4-Based Predictor of Endemic Pemphigus Foliaceus (Fogo Selvagem)

Fogo selvagem (FS) is mediated by pathogenic, predominantly IgG4, anti-desmoglein 1 (Dsg1) autoantibodies and is endemic in Limao Verde, Brazil. IgG and IgG subclass autoantibodies were tested in a sample of 214 FS patients and 261 healthy controls by Dsg1 ELISA. For model selection, the sample was randomly divided into training (50%), validation (25%), and test (25%) sets. Using the training and validation sets, IgG4 was chosen as the best predictor of FS, with index values above 6.43 classified as FS. Using the test set, IgG4 has sensitivity of 92% (95% confidence interval (95% CI): 82-95%), specificity of 97% (95% CI: 89-100%), and area under the curve of 0.97 ( 95% CI: 0.94-1.00). The IgG4 positive predictive value (PPV) in Limao Verde (3% FS prevalence) was 49%. The sensitivity, specificity, and PPV of IgG anti-Dsg1 were 87, 91, and 23%, respectively. The IgG4-based classifier was validated by testing 11 FS patients before and after clinical disease and 60 Japanese pemphigus foliaceus patients. It classified 21 of 96 normal individuals from a Limao Verde cohort as having FS serology. On the basis of its PPV, half of the 21 individuals may currently have preclinical FS and could develop clinical disease in the future. Identifying individuals during preclinical FS will enhance our ability to identify the etiological agent(s) triggering FS.

Bowel dilation as a predictor of adverse outcome in isolated fetal gastroschisis

Objectives This study was designed to evaluate bowel diameter as a predictor of adverse outcome in isolated fetal gastroschisis Methods Retrospective study involving 94 singleton pregnancies Ultrasound measurements of herniated bowel transverse diameter (BTD) were performed up to 3 weeks before delivery Adverse outcome was intrauterine/neonatal death and/or bowel complications Results Last BTD was recorded at 35 6 +/- 1 6 weeks and mean interval to delivery was 6 2 +/- 5 0 days Intrauterine/neonatal death occurred in 10 (10 6%) cases, bowel complications were observed in 8 (8 5%) BTD >= 15, >= 20, >= 25, and >= 30 mm were found in 87, 46, 13, and 4% of pregnancies with a favorable outcome. respectively BTD >= 25 mm sensitivity was 38%. and positive and negative predictive values were 38 and 87% For BTD >= 30 mm. the values were 19, 50, and 85% Observed/expected BTD ROC curve showed an area of 0 67, best cut-off value at 1 39, prediction values were similar to those for BTD >= 25 mm Bowel dilatation was also significantly associated with lower rate of primary surgical closure. longer period to full oral feeding, and prolonged hospital stay Conclusions Bowel dilatation demonstrated up to 3 weeks before delivery is a predictor of intestinal complications and is associated with lower late of primary surgical closure, longer period to achieve full oral feeding. and hospital stay Copyright (C) 2010 John Wiley & Sons, Ltd

Validation of body condition score as a predictor of subcutaneous fat in Nelore (Bos indicus) cows

The aim of the present study was to determine the relationship among body weight (BW), body condition score (BCS) and rump fat thickness (RFAT) measured by ultrasonography, and validate the relationship between BCS and RFAT over the time. Two hundred sixty and six Nelore cows had their BW, BCS and RFAT evaluated at five different moments during the production cycle: M1) weaning: M2) parturition, M3) 42 days post-partum; M4) 82 days postpartum and M5) 112 days post-partum. A BCS value was attributed for each cow following a I to 5 points scale. Ultrasonographic images for RFAT measurement were obtained using a 3.5 MHz linear transducer. Images were immediately analyzed as soon as they were formed and frozen. Body condition scores and ultrasound measurements were collected on the same day by a single trained technician. The relationship between BCS and RFAT values was investigated by regression models. The analysis of similarity among the five obtained models was performed using the proc MIXED from SAS and the correlations among variables were analyzed with proc CORR from SAS. The BCS was able to predict RFAT in Nelore cows in all different moments evaluated. Also, it was shown that BCS presented high correlation (r=0.82 to 0.93) and relationship (R(2) = 0.73 to 0.92) with RFAT. However, both BCS and RFAT showed low correlation (r=0.37 to 0.50) and relationship (R(2) = 0.13 to 0.25) with BW. The BCS classification by visual method using a 1 to 5 point scale, was able to predict the RFAT in Nelore cows over the time. (C) 2008 Elsevier B.V. All rights reserved.

Performance of balanced two-stage empirical predictors of realized cluster latent values from finite populations: A simulation study

Predictors of random effects are usually based on the popular mixed effects (ME) model developed under the assumption that the sample is obtained from a conceptual infinite population; such predictors are employed even when the actual population is finite. Two alternatives that incorporate the finite nature of the population are obtained from the superpopulation model proposed by Scott and Smith (1969. Estimation in multi-stage surveys. J. Amer. Statist. Assoc. 64, 830-840) or from the finite population mixed model recently proposed by Stanek and Singer (2004. Predicting random effects from finite population clustered samples with response error. J. Amer. Statist. Assoc. 99, 1119-1130). Predictors derived under the latter model with the additional assumptions that all variance components are known and that within-cluster variances are equal have smaller mean squared error (MSE) than the competitors based on either the ME or Scott and Smith`s models. As population variances are rarely known, we propose method of moment estimators to obtain empirical predictors and conduct a simulation study to evaluate their performance. The results suggest that the finite population mixed model empirical predictor is more stable than its competitors since, in terms of MSE, it is either the best or the second best and when second best, its performance lies within acceptable limits. When both cluster and unit intra-class correlation coefficients are very high (e.g., 0.95 or more), the performance of the empirical predictors derived under the three models is similar. (c) 2007 Elsevier B.V. All rights reserved.

Abrasion wear resistance of different artificial teeth opposed to metal and composite antagonists

One of the most important properties of artificial teeth is the abrasion wear resistance, which is determinant in the maintenance of the rehabilitation's occlusal pattern. OBJECTIVES: This in vitro study aims to evaluate the abrasion wear resistance of 7 brands of artificial teeth opposed to two types of antagonists. MATERIAL AND METHODS: Seven groups were prepared with 12 specimens each (BIOLUX & BL, TRILUX & TR, BLUE DENT & BD, BIOCLER & BC, POSTARIS & PO, ORTHOSIT & OR, GNATHOSTAR & GN), opposed to metallic (M & nickel-chromium alloy), and to composite antagonists (C & Solidex indirect composite). A mechanical loading device was used (240 cycles/min, 4 Hz speed, 10 mm antagonist course). Initial and final contours of each specimen were registered with aid of a profile projector (20x magnification). The linear difference between the two profiles was measured and the registered values were subjected to ANOVA and Tukey's test. RESULTS: Regarding the antagonists, only OR (M = 10.45 ± 1.42 µm and C = 2.77 ± 0.69 µm) and BC (M = 6.70 ± 1.37 µm and C = 4.48 ± 0.80 µm) presented statistically significant differences (p < 0.05). Best results were obtained with PO (C = 2.33 ± 0.91 µm and M = 1.78 ± 0.42 µm), followed by BL (C = 3.70 ± 1.32 µm and M = 3.70 ± 0.61 µm), statistically similar for both antagonists (p>0.05). Greater result variance was obtained with OR, which presented the worse results opposed to Ni-Cr (10.45 ± 1.42 µm), and results similar to the best ones against composite (2.77 ± 0.69 µm). CONCLUSIONS: Within the limitations of this study, it may be concluded that the antagonist material is a factor of major importance to be considered in the choice of the artificial teeth to be used in the prosthesis.

Correlation between metal-ceramic bond strength and coefficient of linear thermal expansion difference

The purpose of this study was to evaluate the metal-ceramic bond strength (MCBS) of 6 metal-ceramic pairs (2 Ni-Cr alloys and 1 Pd-Ag alloy with 2 dental ceramics) and correlate the MCBS values with the differences between the coefficients of linear thermal expansion (CTEs) of the metals and ceramics. Verabond (VB) Ni-Cr-Be alloy, Verabond II (VB2), Ni-Cr alloy, Pors-on 4 (P), Pd-Ag alloy, and IPS (I) and Duceram (D) ceramics were used for the MCBS test and dilatometric test. Forty-eight ceramic rings were built around metallic rods (3.0 mm in diameter and 70.0 mm in length) made from the evaluated alloys. The rods were subsequently embedded in gypsum cast in order to perform a tensile load test, which enabled calculating the CMBS. Five specimens (2.0 mm in diameter and 12.0 mm in length) of each material were made for the dilatometric test. The chromel-alumel thermocouple required for the test was welded into the metal test specimens and inserted into the ceramics. ANOVA and Tukey's test revealed significant differences (p=0.01) for the MCBS test results (MPa), with PI showing higher MCBS (67.72) than the other pairs, which did not present any significant differences. The CTE (10-6 oC-1) differences were: VBI (0.54), VBD (1.33), VB2I (-0.14), VB2D (0.63), PI (1.84) and PD (2.62). Pearson's correlation test (r=0.17) was performed to evaluate of correlation between MCBS and CTE differences. Within the limitations of this study and based on the obtained results, there was no correlation between MCBS and CTE differences for the evaluated metal-ceramic pairs.

Previsões climáticas sazonais sobre o Brasil: avaliação do RegCM3 aninhado no modelo global CPTEC/COLA

Este trabalho avalia o desempenho de previsões sazonais do modelo climático regional RegCM3, aninhado ao modelo global CPTEC/COLA. As previsões com o RegCM3 utilizaram 60 km de resolução horizontal num domínio que inclui grande parte da América do Sul. As previsões do RegCM3 e CPTEC/COLA foram avaliadas utilizando as análises de chuva e temperatura do ar do Climate Prediction Center (CPC) e National Centers for Enviromental Prediction (NCEP), respectivamente. Entre maio de 2005 e julho de 2007, 27 previsões sazonais de chuva e temperatura do ar (exceto a temperatura do CPTEC/COLA, que possui 26 previsões) foram avaliadas em três regiões do Brasil: Nordeste (NDE), Sudeste (SDE) e Sul (SUL). As previsões do RegCM3 também foram comparadas com as climatologias das análises. De acordo com os índices estatísticos (bias, coeficiente de correlação, raiz quadrada do erro médio quadrático e coeficiente de eficiência), nas três regiões (NDE, SDE e SUL) a chuva sazonal prevista pelo RegCM3 é mais próxima da observada do que a prevista pelo CPTEC/COLA. Além disto, o RegCM3 também é melhor previsor da chuva sazonal do que da média das observações nas três regiões. Para temperatura, as previsões do RegCM3 são superiores às do CPTEC/COLA nas áreas NDE e SUL, enquanto o CPTEC/COLA é superior no SDE. Finalmente, as previsões de chuva e temperatura do RegCM3 são mais próximas das observações do que a climatologia observada. Estes resultados indicam o potencial de utilização do RegCM3 para previsão sazonal, que futuramente deverá ser explorado através de previsão por conjunto.