Generalized linear mixed models for the genetic evaluation of binary reproductive traits: a simulation study

The objective of this study was to evaluate the use of probit and logit link functions for the genetic evaluation of early pregnancy using simulated data. The following simulation/analysis structures were constructed: logit/logit, logit/probit, probit/logit, and probit/probit. The percentages of precocious females were 5, 10, 15, 20, 25 and 30% and were adjusted based on a change in the mean of the latent variable. The parametric heritability (h²) was 0.40. Simulation and genetic evaluation were implemented in the R software. Heritability estimates (ĥ²) were compared with h² using the mean squared error. Pearson correlations between predicted and true breeding values and the percentage of coincidence between true and predicted ranking, considering the 10% of bulls with the highest breeding values (TOP10) were calculated. The mean ĥ² values were under- and overestimated for all percentages of precocious females when logit/probit and probit/logit models used. In addition, the mean squared errors of these models were high when compared with those obtained with the probit/probit and logit/logit models. Considering ĥ², probit/probit and logit/logit were also superior to logit/probit and probit/logit, providing values close to the parametric heritability. Logit/probit and probit/logit presented low Pearson correlations, whereas the correlations obtained with probit/probit and logit/logit ranged from moderate to high. With respect to the TOP10 bulls, logit/probit and probit/logit presented much lower percentages than probit/probit and logit/logit. The genetic parameter estimates and predictions of breeding values of the animals obtained with the logit/logit and probit/probit models were similar. In contrast, the results obtained with probit/logit and logit/probit were not satisfactory. There is need to compare the estimation and prediction ability of logit and probit link functions.

Logistic regression analysis of potential prognostic factors for pulmonary thromboembolism

Objective: To identify potential prognostic factors for pulmonary thromboembolism (PTE), establishing a mathematical model to predict the risk for fatal PTE and nonfatal PTE.Method: the reports on 4,813 consecutive autopsies performed from 1979 to 1998 in a Brazilian tertiary referral medical school were reviewed for a retrospective study. From the medical records and autopsy reports of the 512 patients found with macroscopically and/or microscopically,documented PTE, data on demographics, underlying diseases, and probable PTE site of origin were gathered and studied by multiple logistic regression. Thereafter, the jackknife method, a statistical cross-validation technique that uses the original study patients to validate a clinical prediction rule, was performed.Results: the autopsy rate was 50.2%, and PTE prevalence was 10.6%. In 212 cases, PTE was the main cause of death (fatal PTE). The independent variables selected by the regression significance criteria that were more likely to be associated with fatal PTE were age (odds ratio [OR], 1.02; 95% confidence interval [CI], 1.00 to 1.03), trauma (OR, 8.5; 95% CI, 2.20 to 32.81), right-sided cardiac thrombi (OR, 1.96; 95% CI, 1.02 to 3.77), pelvic vein thrombi (OR, 3.46; 95% CI, 1.19 to 10.05); those most likely to be associated with nonfatal PTE were systemic arterial hypertension (OR, 0.51; 95% CI, 0.33 to 0.80), pneumonia (OR, 0.46; 95% CI, 0.30 to 0.71), and sepsis (OR, 0.16; 95% CI, 0.06 to 0.40). The results obtained from the application of the equation in the 512 cases studied using logistic regression analysis suggest the range in which logit p > 0.336 favors the occurrence of fatal PTE, logit p < - 1.142 favors nonfatal PTE, and logit P with intermediate values is not conclusive. The cross-validation prediction misclassification rate was 25.6%, meaning that the prediction equation correctly classified the majority of the cases (74.4%).Conclusions: Although the usefulness of this method in everyday medical practice needs to be confirmed by a prospective study, for the time being our results suggest that concerning prevention, diagnosis, and treatment of PTE, strict attention should be given to those patients presenting the variables that are significant in the logistic regression model.

Choosing between Cox proportional hazards and logistic models for interval-censored data via bootstrap

This work develops a new methodology in order to discriminate models for interval-censored data based on bootstrap residual simulation by observing the deviance difference from one model in relation to another, according to Hinde (1992). Generally, this sort of data can generate a large number of tied observations and, in this case, survival time can be regarded as discrete. Therefore, the Cox proportional hazards model for grouped data (Prentice & Gloeckler, 1978) and the logistic model (Lawless, 1982) can befitted by means of generalized linear models. Whitehead (1989) considered censoring to be an indicative variable with a binomial distribution and fitted the Cox proportional hazards model using complementary log-log as a link function. In addition, a logistic model can be fitted using logit as a link function. The proposed methodology arises as an alternative to the score tests developed by Colosimo et al. (2000), where such models can be obtained for discrete binary data as particular cases from the Aranda-Ordaz distribution asymmetric family. These tests are thus developed with a basis on link functions to generate such a fit. The example that motivates this study was the dataset from an experiment carried out on a flax cultivar planted on four substrata susceptible to the pathogen Fusarium oxysoprum. The response variable, which is the time until blighting, was observed in intervals during 52 days. The results were compared with the model fit and the AIC values.

Late administration of a specific COX-2 inhibitor does not treat and/or prevent progression of gastric tumors in rats submitted to duodenogastric reflux procedure

PURPOSE:To assess whether late introduction of a specific COX-2 inhibitor (Meloxicam) can treat and/or prevent the progression of tumors in the stomach of rats submitted to duodenogastric reflux. METHODS: Seventy five male Wistar rats, weighing 150 grams, were submitted to the induction of duodenogastric reflux through the pylorus. At 36 weeks of follow-up were established three experimental groups: DGR36 sacrificed immediately, DGR54 and DGR54MLX both sacrificed at 54th week of follow-up . The animals of the latter group were fed with a rat chow premixed with Meloxicam (2.0 mg/ kg feed; 0.3 mg / kg bw / day) and the other two with standard rat chow. The lesions found in the pyloric mucosa and gastrojejunal anastomosis were analyzed macroscopically and histologically. For statistical analysis was adjusted a generalized linear model assuming a binomial distribution with LOGIT link function. RESULTS: No significant differences were found when comparing the incidences of benign tumor lesions (Adenomatous Hyperplasia), p=0.4915, or malignant (Mucinous Adenocarcinoma), p=0.2731, among groups. CONCLUSION: Late introduction of specific COX-2 inhibitor (Meloxicam) did not treat and was not able to prevent the progression of tumoral lesions induced by duodenogastric reflux in the rat stomachs.

Modelos lineares generalizados mistos na avaliação genética da prenhez precoce na raça Nelore

This study aimed to use the generalized linear models with probit and logit link function to evaluate early pregnancy, and to observe the effects on genetic variability and on sire selection when different ages are adopted in the definition of this trait. Early pregnancy was studied at 15 (EP15), and 21 (EP21) months. The analysis was done in R software. Pearson correlations (PC), between genetic predicted values and percentage of bulls in common considering only 10% of bulls with higher genetic values (TOP 10), between classification by logit and probit models and in each model among EP15 and EP21, were calculated. The heritability for EP15 and EP21 were close between models, except for EP15 using probit link function. PC and TOP10 among models were high. The Akaike and Bayesian criteria reported was similar between models. TOP10, considering the same model, among EP15-EP21 were moderated between EP15-EP21.

Determining the season of death from the family composition of insects infesting carrion

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Um perfil de concluintes do curso superior com base no ENADE (2005)

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)