Methods of Estimation in Multiple Linear Regression: Application to Clinical Data

Nesse artigo, tem-se o interesse em avaliar diferentes estratégias de estimação de parâmetros para um modelo de regressão linear múltipla. Para a estimação dos parâmetros do modelo foram utilizados dados de um ensaio clínico em que o interesse foi verificar se o ensaio mecânico da propriedade de força máxima (EM-FM) está associada com a massa femoral, com o diâmetro femoral e com o grupo experimental de ratas ovariectomizadas da raça Rattus norvegicus albinus, variedade Wistar. Para a estimação dos parâmetros do modelo serão comparadas três metodologias: a metodologia clássica, baseada no método dos mínimos quadrados; a metodologia Bayesiana, baseada no teorema de Bayes; e o método Bootstrap, baseado em processos de reamostragem.

Improved score tests in symmetric linear regression models

The class of symmetric linear regression models has the normal linear regression model as a special case and includes several models that assume that the errors follow a symmetric distribution with longer-than-normal tails. An important member of this class is the t linear regression model, which is commonly used as an alternative to the usual normal regression model when the data contain extreme or outlying observations. In this article, we develop second-order asymptotic theory for score tests in this class of models. We obtain Bartlett-corrected score statistics for testing hypotheses on the regression and the dispersion parameters. The corrected statistics have chi-squared distributions with errors of order O(n(-3/2)), n being the sample size. The corrections represent an improvement over the corresponding original Rao`s score statistics, which are chi-squared distributed up to errors of order O(n(-1)). Simulation results show that the corrected score tests perform much better than their uncorrected counterparts in samples of small or moderate size.

Regression analysis of mean quality-adjusted survival time based on pseudo-observations

Regression models for the mean quality-adjusted survival time are specified from hazard functions of transitions between two states and the mean quality-adjusted survival time may be a complex function of covariates. We discuss a regression model for the mean quality-adjusted survival (QAS) time based on pseudo-observations, which has the advantage of directly modeling the effect of covariates in the QAS time. Both Monte Carlo Simulations and a real data set are studied. Copyright (C) 2009 John Wiley & Sons, Ltd.

Nutritional intervention programme among a Japanese-Brazilian community: procedures and results according to gender

Objective: To describe the results of a nutritional intervention programme among Japanese-Brazilians according to gender. Design: A non-controlled experimental study. Setting: The research included three points of clinical, nutritional and physical activity evaluation: at baseline (in 2005), after the first year and at the end of the second year (in 2007). The paired Student t test and multiple linear regression analysis were used to evaluate changes in the subjects` profile (clinical, nutritional and physical activity variables). Subjects: Japanese-Brazilians (n 575) of both genders, aged over 30 years. Results: We verified statistically significant reductions in body weight (0.9 kg), waist circumference (2.9 cm), blood pressure, fasting blood glucose (>3 mg/dl) and total cholesterol (>20 mg/dl) and its fractions, in both genders. We also found reductions in intake of energy (among men), protein (among women) and fat (both genders) and increases in intake of total fibre (among women) and carbohydrate (among men). Conclusions: The intervention programme indicated meaningful benefits for the intervention subjects, with changes in their habits that led to a `healthier` lifestyle positively impacting their nutritional and metabolic profile.

Iron Concentrations in Breast Milk and Selected Maternal Factors of Human Milk Bank Donors

The aim of this study was to evaluate the relationship between iron concentration in mature breast milk and characteristics of 136 donors of a Brazilian milk bank. Iron, vitamin A, zinc, and copper concentrations were assessed in human milk and maternal blood. Data were collected on maternal anthropometrics, obstetric, socioeconomic, demographic, and lifestyle factors. Iron, zinc, and copper in milk and zinc and copper in blood were detected by spectrophotometry. Vitamin A in milk and blood was determined by high-performance liquid chromatography. Hemoglobin was measured by electronic counting and serum iron and ferritin by colorimetry and chemoluminescence, respectively. Transferrin and ceruloplasmin were determined by nephelometry. According to multivariate linear regression analysis, iron in milk was positively associated with vitamin A in milk and with smoking but negatively associated with timing of breast milk donation (P < .001). These results indicate that iron concentration in milk of Brazilian donors may be influenced by nutritional factors and smoking. J Hum Lact. 26(2):175-179

The impact of maternal cortisol concentrations on child arterial elasticity

Epidemiological studies suggest that glucocorticoid excess in the fetus may contribute to the pathophysiology of cardiovascular diseases in adulthood. However, the impact of maternal glucocorticoid on the cardiovascular system of the offspring has not been much explored in studies involving humans, especially in childhood. The objective of this study was to assess the influence of maternal cortisol concentrations on child arterial elasticity. One hundred and thirty pregnant women followed from 1997 to 2000, and respective children 5-7 years of age followed from 2004 to 2006 were included in the study. Maternal cortisol was determined in saliva by an enzyme immunoassay utilizing the mean concentration of nine samples of saliva. Arterial elasticity was assessed by the large artery elasticity index (LAEI; the capacitive elasticity of large arteries) by recording radial artery pulse wave, utilizing the equipment HDI/PulseWave CR-2000 Cardiovascular Profiling System (R). The nutritional status of the children was determined by the body mass index (BMI). Insulin concentration was assessed by chemiluminescence, and insulin resistance by the homeostasis model assessment. Blood glucose, total cholesterol and fractions (LDL-c and HDL-c) and triglyceride concentrations were determined by automated enzymatic methods. The association between maternal cortisol and child arterial elasticity was assessed by multivariate linear regression analysis. There was a statistically significant association between maternal cortisol and LAEI (P=0.02), controlling for birth weight, age, BMI and HDL-c of the children. This study suggests that exposure to higher glucocorticoid concentrations in the prenatal period is associated to lower arterial elasticity in childhood, an earlier cardiovascular risk marker.

A Longitudinal Three-Center Study of Dental Arch Relationship in Patients With Bilateral Cleft Lip and Palate

Objective: To compare and evaluate longitudinally the dental arch relationships from 4.5 to 13.5 years of age with the Bauru-BCLP Yardstick in a large sample of patients with bilateral cleft lip and palate (BCLP). Design: Retrospective longitudinal intercenter outcome study. Patients: Dental casts of 204 consecutive patients with complete BCLP were evaluated at 6, 9, and 12 years of age. All models were identified only by random identification numbers. Setting: Three cleft palate centers with different treatment protocols. Main Outcome Measures: Dental arch relationships were categorized with the Bauru-BCLP yardstick. Increments for each interval (from 6 to 9 years, 6 to 12 years, and 9 to 12 years) were analyzed by logistic and linear regression models. Results: There were no significant differences in outcome measures between the centers at age 12 or at age 9. At age 6, center B showed significantly better results (p = .027), but this difference diminished as the yardstick score for this group increased over time (linear regression analysis), the difference with the reference category (center C, boys) for the intervals 6 to 12 and 9 to 12 years being 10.4% (p = .041) and 12.9% (p = .009), respectively. Conclusions: Despite different treatment protocols, dental arch relationships in the three centers were comparable in final scores at age 9 and 12 years. Delaying hard palate closure and employing infant orthopedics did not appear to be advantageous in the long run. Premaxillary osteotomy employed in center B appeared to be associated with less favorable development of the dental arch relationship between 9 and 12 years.

A reduction of CETP activity, not an increase, is associated with modestly impaired postprandial lipemia and increased HDL-Cholesterol in adult asymptomatic women

Background: The relationship between CETP and postprandial hyperlipemia is still unclear. We verified the effects of varying activities of plasma CETP on postprandial lipemia and precocious atherosclerosis in asymptomatic adult women. Methods: Twenty-eight women, selected from a healthy population sample (n = 148) were classified according to three CETP levels, all statistically different: CETP deficiency (CETPd <= 4.5%, n = 8), high activity (CETPi >= 23.8, n = 6) and controls (CTL, CETP >= 4.6% and <= 23.7%, n = 14). After a 12 h fast they underwent an oral fat tolerance test (40 g of fat/m(2) of body surface area) for 8 hours. TG, TG-rich-lipoproteins (TRL), cholesterol and TRL-TG measurements (AUC, AUIC, AR, RR and late peaks) and comparisons were performed on all time points. Lipases and phospholipids transfer protein (PLTP) were determined. Correlation between carotid atherosclerosis (c-IMT) and postprandial parameters was determined. CETP TaqIB and I405V and ApoE-epsilon 3/epsilon 2/epsilon 4 polymorphisms were examined. To elucidate the regulation of increased lipemia in CETPd a multiple linear regression analysis was performed. Results: In the CETPi and CTL groups, CETP activity was respectively 9 and 5.3 higher compared to the CETPd group. Concentrations of all HDL fractions and ApoA-I were higher in the CETPd group and clearance was delayed, as demonstrated by modified lipemia parameters (AUC, AUIC, RR, AR and late peaks and meal response patterns). LPL or HL deficiencies were not observed. No genetic determinants of CETP deficiency or of postprandial lipemia were found. Correlations with c-IMT in the CETPd group indicated postprandial pro-atherogenic associations. In CETPd the regression multivariate analysis (model A) showed that CETP was largely and negatively predicted by VLDL-C lipemia (R(2) = 92%) and much less by TG, LDL-C, ApoAI, phospholipids and non-HDL-C. CETP (model B) influenced mainly the increment in ApoB-100 containing lipoproteins (R(2) = 85% negatively) and phospholipids (R(2) = 13%), at the 6(th)h point. Conclusion: The moderate CETP deficiency phenotype included a paradoxically high HDL-C and its sub fractions (as earlier described), positive associations with c-IMT, a postprandial VLDL-C increment predicting negatively CETP activity and CETP activity regulating inversely the increment in ApoB100-containing lipoproteins. We hypothesize that the enrichment of TG content in triglyceride-rich ApoB-containing lipoproteins and in TG rich remnants increases lipoproteins` competition to active lipolysis sites, reducing their catabolism and resulting on postprandial lipemia with atherogenic consequences.

The log-exponentiated Weibull regression model for interval-censored data

In interval-censored survival data, the event of interest is not observed exactly but is only known to occur within some time interval. Such data appear very frequently. In this paper, we are concerned only with parametric forms, and so a location-scale regression model based on the exponentiated Weibull distribution is proposed for modeling interval-censored data. We show that the proposed log-exponentiated Weibull regression model for interval-censored data represents a parametric family of models that include other regression models that are broadly used in lifetime data analysis. Assuming the use of interval-censored data, we employ a frequentist analysis, a jackknife estimator, a parametric bootstrap and a Bayesian analysis for the parameters of the proposed model. We derive the appropriate matrices for assessing local influences on the parameter estimates under different perturbation schemes and present some ways to assess global influences. Furthermore, for different parameter settings, sample sizes and censoring percentages, various simulations are performed; in addition, the empirical distribution of some modified residuals are displayed and compared with the standard normal distribution. These studies suggest that the residual analysis usually performed in normal linear regression models can be straightforwardly extended to a modified deviance residual in log-exponentiated Weibull regression models for interval-censored data. (C) 2009 Elsevier B.V. All rights reserved.

Deviance residuals in generalised log-gamma regression models with censored observations

In this article, we compare three residuals based on the deviance component in generalised log-gamma regression models with censored observations. For different parameter settings, sample sizes and censoring percentages, various simulation studies are performed and the empirical distribution of each residual is displayed and compared with the standard normal distribution. For all cases studied, the empirical distributions of the proposed residuals are in general symmetric around zero, but only a martingale-type residual presented negligible kurtosis for the majority of the cases studied. These studies suggest that the residual analysis usually performed in normal linear regression models can be straightforwardly extended for the martingale-type residual in generalised log-gamma regression models with censored data. A lifetime data set is analysed under log-gamma regression models and a model checking based on the martingale-type residual is performed.

Influence diagnostics in Birnbaum-Saunders nonlinear regression models

We consider the issue of assessing influence of observations in the class of Birnbaum-Saunders nonlinear regression models, which is useful in lifetime data analysis. Our results generalize those in Galea et al. [8] which are confined to Birnbaum-Saunders linear regression models. Some influence methods, such as the local influence, total local influence of an individual and generalized leverage are discussed. Additionally, the normal curvatures for studying local influence are derived under some perturbation schemes. We also give an application to a real fatigue data set.

Effect of inbreeding on growth and reproductive traits of Nellore cattle in Brazil

The objective of this study was to identify and quantify the influence of F (inbreeding coefficient) on weaning weight (WW), weight gain from weaning to 18 months of age (WG345), finishing visual score (precocity) at 18 months of age, muscling visual score at 18 months of age (MUS), hip height (HH), scrotal circumference at 18 months of age (SC), heifer probability of pregnancy at 14 months of age (PP14), and stayability (STAY) in Brazilian Nellore cattle. The complete pedigree included 417,552 animals born between 1984 and 2007 on 12 farms located in the states of Mato Grosso do Sul, Sao Paulo and Bahia. Following the observation of a statistically significant effect (P<0.05) of the covariates individual inbreeding coefficient (F) and maternal inbreeding coefficient, regression analysis of each trait, adjusted for all other effects, was performed as a function of the linear and quadratic effect of F and maternal F (when significant). Inbreeding negatively affected all traits studied (P<0.05), except for muscling. A quadratic effect of individual F on WW, WG345, HH and PP14, and a quadratic effect of maternal F on WG345 and HH were observed. Levels of inbreeding higher than 7-11% affected negatively growth and reproductive performance of Nellore cattle. Therefore, inbreeding should be avoided, except for purposes of genetic breeding whose main objective is the fixation of certain alleles in the population. (C) 2010 Elsevier B.V. All rights reserved.

ESTIMATING SENSORY SHELF LIFE OF CHOCOLATE AND CARROT CUPCAKES USING ACCEPTANCE TESTS

This study determined the sensory shelf life of a commercial brand of chocolate and carrot cupcakes, aiming at increasing the current 120 days of shelf life to 180. Appearance, texture, flavor and overall quality of cakes stored at six different storage times were evaluated by 102 consumers. The data were analyzed by analysis of variance and linear regression. For both flavors, the texture presented a greater loss in acceptance during the storage period, showing an acceptance mean close to indifference on the hedonic scale at 120 days. Nevertheless, appearance, flavor and overall quality stayed acceptable up to 150 days. The end of shelf life was estimated at about 161 days for chocolate cakes and 150 days for carrot cakes. This study showed that the current 120 days of shelf life can be extended to 150 days for carrot cake and to 160 days for chocolate cake. However, the 180 days of shelf life desired by the company were not achieved. PRACTICAL APPLICATIONS This research shows the adequacy of using sensory acceptance tests to determine the shelf life of two food products (chocolate and carrot cupcakes). This practical application is useful because the precise determination of the shelf life of a food product is of vital importance for its commercial success. The maximum storage time should always be evaluated in the development or reformulation of new products, changes in packing or storage conditions. Once the physical-chemical and microbiological stability of a product is guaranteed, sensorial changes that could affect consumer acceptance will determine the end of the shelf life of a food product. Thus, the use of sensitive and reliable methods to estimate the sensory shelf life of a product is very important. Findings show the importance of determining the shelf life of each product separately and to avoid using the shelf time estimated for a specific product on other, similar products.

Ratios of UV, PAR and NIR components to global solar radiation measured at Botucatu site in Brazil

The relationships between the four radiant fluxes are analyzed based on a 4 year data archive of hourly and daily global ultraviolet (I(UV)), photosynthetically active-PAR (I(PAR)), near infrared (I(NIR)) and broadband global solar radiation (I(G)) collected at Botucatu, Brazil. These data are used to establish both the fractions of spectral components to global solar radiation and the proposed linear regression models. Verification results indicated that the proposed regression models predict accurately the spectral radiant fluxes at least for the Brazilian environment. Finally, results obtained in this analysis agreed well with most published results in the literature. (c) 2010 Elsevier Ltd. All rights reserved.

Modeling hourly and daily fractions of UV, PAR and NIR to global solar radiation under various sky conditions at Botucatu, Brazil

In this analysis, using available hourly and daily radiometric data performed at Botucatu, Brazil, several empirical models relating ultraviolet (UV), photosynthetically active (PAR) and near infrared (NIR) solar global components with solar global radiation (G) are established. These models are developed and discussed through clearness index K(T) (ratio of the global-to-extraterrestrial solar radiation). Results obtained reveal that the proposed empirical models predict hourly and daily values accurately. Finally. the overall analysis carried Out demonstrates that the sky conditions are more important in developing correlation models between the UV component and the global solar radiation. The linear regression models derived to estimate PAR and NIR components may be obtained without sky condition considerations within a maximum variation of 8%. In the case of UV, not taking into consideration the sky condition may cause a discrepancy of up to 18% for hourly values and 15% for daily values. (C) 2008 Elsevier Ltd. All rights reserved.