Improved likelihood inference for the shape parameter in Weibull regression

We obtain adjustments to the profile likelihood function in Weibull regression models with and without censoring. Specifically, we consider two different modified profile likelihoods: (i) the one proposed by Cox and Reid [Cox, D.R. and Reid, N., 1987, Parameter orthogonality and approximate conditional inference. Journal of the Royal Statistical Society B, 49, 1-39.], and (ii) an approximation to the one proposed by Barndorff-Nielsen [Barndorff-Nielsen, O.E., 1983, On a formula for the distribution of the maximum likelihood estimator. Biometrika, 70, 343-365.], the approximation having been obtained using the results by Fraser and Reid [Fraser, D.A.S. and Reid, N., 1995, Ancillaries and third-order significance. Utilitas Mathematica, 47, 33-53.] and by Fraser et al. [Fraser, D.A.S., Reid, N. and Wu, J., 1999, A simple formula for tail probabilities for frequentist and Bayesian inference. Biometrika, 86, 655-661.]. We focus on point estimation and likelihood ratio tests on the shape parameter in the class of Weibull regression models. We derive some distributional properties of the different maximum likelihood estimators and likelihood ratio tests. The numerical evidence presented in the paper favors the approximation to Barndorff-Nielsen`s adjustment.

Inference for a skew extension of the Grubbs model

In this paper, we discuss inferential aspects for the Grubbs model when the unknown quantity x (latent response) follows a skew-normal distribution, extending early results given in Arellano-Valle et al. (J Multivar Anal 96:265-281, 2005b). Maximum likelihood parameter estimates are computed via the EM-algorithm. Wald and likelihood ratio type statistics are used for hypothesis testing and we explain the apparent failure of the Wald statistics in detecting skewness via the profile likelihood function. The results and methods developed in this paper are illustrated with a numerical example.

Prediction in Multilevel Logistic Regression

The purpose of this article is to present a new method to predict the response variable of an observation in a new cluster for a multilevel logistic regression. The central idea is based on the empirical best estimator for the random effect. Two estimation methods for multilevel model are compared: penalized quasi-likelihood and Gauss-Hermite quadrature. The performance measures for the prediction of the probability for a new cluster observation of the multilevel logistic model in comparison with the usual logistic model are examined through simulations and an application.

Improved testing inference in mixed linear models

Mixed linear models are commonly used in repeated measures studies. They account for the dependence amongst observations obtained from the same experimental unit. Often, the number of observations is small, and it is thus important to use inference strategies that incorporate small sample corrections. In this paper, we develop modified versions of the likelihood ratio test for fixed effects inference in mixed linear models. In particular, we derive a Bartlett correction to such a test, and also to a test obtained from a modified profile likelihood function. Our results generalize those in [Zucker, D.M., Lieberman, O., Manor, O., 2000. Improved small sample inference in the mixed linear model: Bartlett correction and adjusted likelihood. Journal of the Royal Statistical Society B, 62,827-838] by allowing the parameter of interest to be vector-valued. Additionally, our Bartlett corrections allow for random effects nonlinear covariance matrix structure. We report simulation results which show that the proposed tests display superior finite sample behavior relative to the standard likelihood ratio test. An application is also presented and discussed. (C) 2008 Elsevier B.V. All rights reserved.

Reduction of Models in the Presence of Nuisance Parameters

In many statistical inference problems, there is interest in estimation of only some elements of the parameter vector that defines the adopted model. In general, such elements are associated to measures of location and the additional terms, known as nuisance parameters, to control the dispersion and asymmetry of the underlying distributions. To estimate all the parameters of the model and to draw inferences only on the parameters of interest. Depending on the adopted model, this procedure can be both algebraically is common and computationally very costly and thus it is convenient to reduce it, so that it depends only on the parameters of interest. This article reviews estimation methods in the presence of nuisance parameters and consider some applications in models recently discussed in the literature.

Quasi-elastic barrier distribution of the (17)O+(64)Zn system and the derivation of the (17)O nuclear matter diffuseness

The quasi-elastic excitation function for the (17)O+(64)Zn system was measured at energies near and below the Coulomb barrier, at the backward angle theta(lab) = 161 degrees. The corresponding quasi-elastic barrier distribution was derived. The excitation function for the neutron stripping reactions was also measured, at the same angle and energies, and the experimental values of the spectroscopic factors were deduced by fitting the data. A reasonably good agreement was obtained between the experimental quasi-elastic barrier distribution with the coupled-channel calculations including a very large number of channels. Of the channels investigated, three dominated the coupling matrix: two inelastic channels, (64)Zn(2(1)(+)) and (17)O(1/(+)(2)), and one-neutron transfer channel, particularly the first one. On the other hand, a very good agreement is obtained when we use a nuclear diffuseness for the (17)O nucleus larger than the one for (16)O. We verify that quasi-elastic barrier distribution is a sensitive tool for determining nuclear matter diffuseness.

Structured minimal-memory inexact quasi-Newton method and secant preconditioners for augmented Lagrangian optimization

Augmented Lagrangian methods for large-scale optimization usually require efficient algorithms for minimization with box constraints. On the other hand, active-set box-constraint methods employ unconstrained optimization algorithms for minimization inside the faces of the box. Several approaches may be employed for computing internal search directions in the large-scale case. In this paper a minimal-memory quasi-Newton approach with secant preconditioners is proposed, taking into account the structure of Augmented Lagrangians that come from the popular Powell-Hestenes-Rockafellar scheme. A combined algorithm, that uses the quasi-Newton formula or a truncated-Newton procedure, depending on the presence of active constraints in the penalty-Lagrangian function, is also suggested. Numerical experiments using the Cute collection are presented.

Self-perception and malocclusion and their relation to oral appearance and function

The aim of this study was to evaluate the relationship between malocclusion and self-perception of oral appearance/function, in 12/15-year-old Brazilian adolescents. The cluster sample consisted of 717 teenagers attending 24 urban public (n=611) and 5 rural public (n=107) schools in Maringá/PR. Malocclusion was measured using the Dental Aesthetic Index (DAI), in accordance with WHO recommendations. A parental questionnaire was applied to collect information on esthetic perception level and oral variables related to oral health. Univariate and multiple logistic regression analyses were performed. Multiple logistic regression confirmed that for 12-year-old, missing teeth (OR=2.865) and presence of openbite (open occlusal relationship) (OR=2.865) were risk indicators for speech capability. With regard to 15-year-old, presence of mandibular overjet (horizontal overlap) (OR=4.016) was a risk indicator for speech capability and molar relationship (OR=1.661) was a risk indicator for chewing capability. The impact of malocclusion on adolescents' life was confirmed in this study. Speech and chewing capability were associated with orthodontic deviations, which should be taken into consideration in oral health planning, to identify risk groups and improve community health services.

Richards growth model and viability indicators for populations subject to interventions

In this work we study the problem of modeling identification of a population employing a discrete dynamic model based on the Richards growth model. The population is subjected to interventions due to consumption, such as hunting or farming animals. The model identification allows us to estimate the probability or the average time for a population number to reach a certain level. The parameter inference for these models are obtained with the use of the likelihood profile technique as developed in this paper. The identification method here developed can be applied to evaluate the productivity of animal husbandry or to evaluate the risk of extinction of autochthon populations. It is applied to data of the Brazilian beef cattle herd population, and the the population number to reach a certain goal level is investigated.

Flow cytometry as a tool in the evaluation of blood leukocyte function in Chelonia mydas (Linnaeus, 1758) (Testudines, Cheloniidae)

Chelonia mydas is a sea turtle that feeds and nests on the Brazilian coast and a disease called fibropapillomatosis is a threat to this species. Because of this, it is extremely necessary to determine a methodology that would enable the analysis of blood leukocyte function in these sea turtles. In order to achieve this aim, blood samples were collected from C. mydas with or without fibropapillomas captured on the São Paulo north coast. Blood samples were placed in tubes containing sodium heparin and were transported under refrigeration to the laboratory in sterile RPMI 1640 cell culture medium. Leukocytes were separated by density gradient using Ficoll-PaqueTM Plus, Amershan Biociences®. The following stimuli were applied in the assessment of leukocyte function: Phorbol Miristate-Acetate (PMA) for oxidative burst activity evaluation and Zymosan A (Saccharomyces cerevisiae) Bio Particles®, Alexa Fluor® 594 conjugate for phagocytosis evaluation. Three cell populations were identified: heterophils, monocytes and lymphocytes. Monocytes were the cells responsible for phagocytosis and oxidative burst.

Eccentric Hip Muscle Function in Females With and Without Patellofemoral Pain Syndrome

Context: Patellofemoral pain syndrome (PFPS) is a common knee condition in athletes. Recently, researchers have indicated that factors proximal to the knee, including hip muscle weakness and motor control impairment, contribute to the development of PFPS. However, no investigators have evaluated eccentric hip muscle function in people with PFPS. Objective: To compare the eccentric hip muscle function between females with PFPS and a female control group. Design: Cross-sectional study. Setting: Musculoskeletal laboratory. Patients or Other Participants: Two groups of females were studied: a group with PFPS (n = 10) and a group with no history of lower extremity injury or surgery (n = 10). Intervention(s): Eccentric torque of the hip musculature was evaluated on an isokinetic dynamometer. Main Outcome Measure(s): Eccentric hip abduction, adduction, and external and internal rotation peak torque were measured and expressed as a percentage of body mass (Nm/kg x 100). We also evaluated eccentric hip adduction to abduction and internal to external rotation torque ratios. The peak torque value of 5 maximal eccentric contractions was used for calculation. Two-tailed, independent-samples t tests were used to compare torque results between groups. Results: Participants with PFPS exhibited much lower eccentric hip abduction (t(18) = -2.917, P = .008) and adduction (t(18) = -2.764, P =.009) peak torque values than did their healthy counterparts. No differences in eccentric hip external (t(18) = 0.45, P = .96) or internal (t(18) = -0.742, P =.47) rotation peak torque values were detected between the groups. The eccentric hip adduction to abduction torque ratio was much higher in the PFPS group than in the control group (t(18) = 2.113, P = .04), but we found no difference in the eccentric hip internal to external rotation torque ratios between the 2 groups (t(18) = -0.932, P = .36). Conclusions: Participants with PFPS demonstrated lower eccentric hip abduction and adduction peak torque and higher eccentric adduction to abduction torque ratios when compared with control participants. Thus, clinicians should consider eccentric hip abduction strengthening exercises when developing rehabilitation programs for females with PFPS.

Genomic islands in the pathogenic filamentous fungus Aspergillus fumigatus

We present the genome sequences of a new clinical isolate of the important human pathogen, Aspergillus fumigatus, A1163, and two closely related but rarely pathogenic species, Neosartorya fischeri NRRL181 and Aspergillus clavatus NRRL1. Comparative genomic analysis of A1163 with the recently sequenced A. fumigatus isolate Af293 has identified core, variable and up to 2% unique genes in each genome. While the core genes are 99.8% identical at the nucleotide level, identity for variable genes can be as low 40%. The most divergent loci appear to contain heterokaryon incompatibility ( het) genes associated with fungal programmed cell death such as developmental regulator rosA. Cross-species comparison has revealed that 8.5%, 13.5% and 12.6%, respectively, of A. fumigatus, N. fischeri and A. clavatus genes are species-specific. These genes are significantly smaller in size than core genes, contain fewer exons and exhibit a subtelomeric bias. Most of them cluster together in 13 chromosomal islands, which are enriched for pseudogenes, transposons and other repetitive elements. At least 20% of A. fumigatus-specific genes appear to be functional and involved in carbohydrate and chitin catabolism, transport, detoxification, secondary metabolism and other functions that may facilitate the adaptation to heterogeneous environments such as soil or a mammalian host. Contrary to what was suggested previously, their origin cannot be attributed to horizontal gene transfer ( HGT), but instead is likely to involve duplication, diversification and differential gene loss (DDL). The role of duplication in the origin of lineage-specific genes is further underlined by the discovery of genomic islands that seem to function as designated ""gene dumps'' and, perhaps, simultaneously, as "" gene factories''.

The four hexamerin genes in the honey bee: structure, molecular evolution and function deduced from expression patterns in queens, workers and drones

Background: Hexamerins are hemocyanin-derived proteins that have lost the ability to bind copper ions and transport oxygen; instead, they became storage proteins. The current study aimed to broaden our knowledge on the hexamerin genes found in the honey bee genome by exploring their structural characteristics, expression profiles, evolution, and functions in the life cycle of workers, drones and queens. Results: The hexamerin genes of the honey bee (hex 70a, hex 70b, hex 70c and hex 110) diverge considerably in structure, so that the overall amino acid identity shared among their deduced protein subunits varies from 30 to 42%. Bioinformatics search for motifs in the respective upstream control regions (UCRs) revealed six overrepresented motifs including a potential binding site for Ultraspiracle (Usp), a target of juvenile hormone (JH). The expression of these genes was induced by topical application of JH on worker larvae. The four genes are highly transcribed by the larval fat body, although with significant differences in transcript levels, but only hex 110 and hex 70a are re-induced in the adult fat body in a caste-and sex-specific fashion, workers showing the highest expression. Transcripts for hex 110, hex 70a and hex70b were detected in developing ovaries and testes, and hex 110 was highly transcribed in the ovaries of egg-laying queens. A phylogenetic analysis revealed that HEX 110 is located at the most basal position among the holometabola hexamerins, and like HEX 70a and HEX 70c, it shares potential orthology relationship with hexamerins from other hymenopteran species. Conclusions: Striking differences were found in the structure and developmental expression of the four hexamerin genes in the honey bee. The presence of a potential binding site for Usp in the respective 5' UCRs, and the results of experiments on JH level manipulation in vivo support the hypothesis of regulation by JH. Transcript levels and patterns in the fat body and gonads suggest that, in addition to their primary role in supplying amino acids for metamorphosis, hexamerins serve as storage proteins for gonad development, egg production, and to support foraging activity. A phylogenetic analysis including the four deduced hexamerins and related proteins revealed a complex pattern of evolution, with independent radiation in insect orders.

Comparison of the effects of carbon dioxide and helium pneumoperitoneum on renal function

Background and Purpose: Carbon dioxide pneumoperitoneum is associated with significant hypercarbia and acidosis. The aim of this study is to evaluate the effects of carbon dioxide and helium pneumoperitoneum on renal function. Materials and Methods: Thirty adult dogs were put randomly into one of three groups ( n = 10 animals each): group A - pneumoperitoneum not performed; group B - CO2 pneumoperitoneum; and group C - helium pneumoperitoneum. The groups were analyzed with consideration given to body weight, hematologic values, hemodynamic parameters ( heart rate, mean arterial pressure, central venous pressure, cardiac output, stroke volume, systemic vascular resistance, pulmonary vascular resistance, left cardiac work index, cardiac index, mean pulmonary artery pressure, and pulmonary capillary wedge pressure), and renal function ( plasma renin activity, urinary output, creatinine clearance, and sodium excretory fraction). Results: An accentuated decrease in urinary output was observed during pneumoperitoneum in groups B and C compared to the control group. In groups B and C, creatinine clearance declined significantly during pneumoperitoneum in comparison to group A, but after deflation a faster recovery of glomerular filtration was noticed for group C, and a significant increase in sodium excretory fraction was seen for group B. On the other hand, in comparison to the control group, group B had a significant increase in plasma renin activity, with late recovery of glomerular function. Conclusion: Helium ameliorates renal alterations when used for pneumoperitoneum, and it might be used for patients with compromised renal function who have to undergo laparoscopic surgery.

Aerobic capacity and skeletal muscle function in children with asthma

Background Peripheral muscle strength and endurance are decreased in patients with chronic pulmonary diseases and seem to contribute to patients' exercise intolerance. However, the authors are not aware of any studies evaluating peripheral muscle function in children with asthma. It seems to be implied that children with asthma have lower aerobic fitness, but there are limited studies comparing the aerobic capacity of children with and without asthma. The present study aimed to evaluate muscle strength and endurance in children with persistent asthma and their association with aerobic capacity and inhaled corticosteroid consumption. Methods Forty children with mild persistent asthma (MPA) or severe persistent asthma (SPA) (N=20 each) and 20 children without asthma (control group) were evaluated. Upper (pectoralis and latissimus dorsi) and lower (quadriceps) muscle strength and endurance were assessed, and cardiopulmonary exercise testing was performed. Inhaled corticosteroid consumption during the last 6 and 24 months was also quantified. Results Children with SPA presented a reduction in peak oxygen consumption (VO(2)) (28.2 +/- 8.1 vs 34.7 +/- 6.9 ml/kg/min; p<0.01) and quadriceps endurance (43.1 +/- 6.7 vs 80.9 +/- 11.9 repetitions; p<0.05) compared with the control group, but not the MPA group (31.5 +/- 6.1 ml/kg/min and 56.7 +/- 47.7 repetitions respectively; p>0.05). Maximal upper and lower muscle strength was preserved in children with both mild and severe asthma (p>0.05). Finally, the authors observed that lower muscle endurance weakness was not associated with reductions in either peak VO(2) (r=0.22, p>0.05) or corticosteroid consumption (r=-0.31, p>0.05) in children with asthma. Conclusion The findings suggest that cardiopulmonary exercise and lower limb muscle endurance should be a priority during physical training programs for children with severe asthma.