On beta regression residuals

We propose two new residuals for the class of beta regression models, and numerically evaluate their behaviour relative to the residuals proposed by Ferrari and Cribari-Neto. Monte Carlo simulation results and empirical applications using real and simulated data are provided. The results favour one of the residuals we propose.

Comparing diagnostic tests with missing data

When missing data occur in studies designed to compare the accuracy of diagnostic tests, a common, though naive, practice is to base the comparison of sensitivity, specificity, as well as of positive and negative predictive values on some subset of the data that fits into methods implemented in standard statistical packages. Such methods are usually valid only under the strong missing completely at random (MCAR) assumption and may generate biased and less precise estimates. We review some models that use the dependence structure of the completely observed cases to incorporate the information of the partially categorized observations into the analysis and show how they may be fitted via a two-stage hybrid process involving maximum likelihood in the first stage and weighted least squares in the second. We indicate how computational subroutines written in R may be used to fit the proposed models and illustrate the different analysis strategies with observational data collected to compare the accuracy of three distinct non-invasive diagnostic methods for endometriosis. The results indicate that even when the MCAR assumption is plausible, the naive partial analyses should be avoided.

Leverage analysis for linear mixed models

We consider a generalized leverage matrix useful for the identification of influential units and observations in linear mixed models and show how a decomposition of this matrix may be employed to identify high leverage points for both the marginal fitted values and the random effect component of the conditional fitted values. We illustrate the different uses of the two components of the decomposition with a simulated example as well as with a real data set.

On the skew-normal calibration model

In this article, we present the EM-algorithm for performing maximum likelihood estimation of an asymmetric linear calibration model with the assumption of skew-normally distributed error. A simulation study is conducted for evaluating the performance of the calibration estimator with interpolation and extrapolation situations. As one application in a real data set, we fitted the model studied in a dimensional measurement method used for calculating the testicular volume through a caliper and its calibration by using ultrasonography as the standard method. By applying this methodology, we do not need to transform the variables to have symmetrical errors. Another interesting aspect of the approach is that the developed transformation to make the information matrix nonsingular, when the skewness parameter is near zero, leaves the parameter of interest unchanged. Model fitting is implemented and the best choice between the usual calibration model and the model proposed in this article was evaluated by developing the Akaike information criterion, Schwarz`s Bayesian information criterion and Hannan-Quinn criterion.

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.

Bias-corrected Pearson estimating functions for Taylor`s power law applied to benthic macrofauna data

Estimation of Taylor`s power law for species abundance data may be performed by linear regression of the log empirical variances on the log means, but this method suffers from a problem of bias for sparse data. We show that the bias may be reduced by using a bias-corrected Pearson estimating function. Furthermore, we investigate a more general regression model allowing for site-specific covariates. This method may be efficiently implemented using a Newton scoring algorithm, with standard errors calculated from the inverse Godambe information matrix. The method is applied to a set of biomass data for benthic macrofauna from two Danish estuaries. (C) 2011 Elsevier B.V. All rights reserved.

A critical view on biaxial and short-beam uniaxial flexural strength tests applied to resin composites using Weibull, fractographic and finite element analyses

Objective. To evaluate the biaxial and short-beam uniaxial strength tests applied to resin composites based upon their Weibull parameters, fractographic features and stress distribution. Methods. Disk- (15 mm x 1 mm) and beam-shaped specimens (10 mm x 2 mm x 1 mm) of three commercial composites (Concept/Vigodent, CA; Heliomolar/Ivoclar-Vivadent, HE; Z250/3M ESPE, FZ) were prepared. After 48h dry storage at 37 degrees C, disks and beams were submitted to piston-on-three-balls (BI) and three-point bending (UNI) tests, respectively. Data were analyzed by Weibull statistics. Fractured surfaces were observed under stereomicroscope and scanning electron microscope. Maximum principal stress (sigma(1)) distribution was determined by finite element analysis (FEA). Maximum sigma(1-BI) and sigma(1-UNI) were compared to FZ strengths calculated by applying the average failure loads to the analytical equations (sigma(a-BI) and sigma(a-UNI)). Results. For BI, characteristic strengths were: 169.9a (FZ), 122.4b (CA) and 104.8c (HE), and for UNI were: 160.3a (FZ), 98.2b (CA) and 91.6b (HE). Weibull moduli ( m) were similar within the same test. CA and HE presented statistically higher m for BI. Surface pores ( BI) and edge flaws ( UNI) were the most frequent fracture origins. sigma(1-BI) was 14% lower than sigma(a-BI.) sigma(1-UNI) was 43% higher than sigma(a-UNI). Significance. Compared to the short-beam uniaxial test, the biaxial test detected more differences among composites and displayed less data scattering for two of the tested materials. Also, biaxial strength was closer to the material`s strength estimated by FEA. (C) 2009 Academy of Dental Materials. Published by Elsevier Ltd. All rights reserved.

Heteroscedastic controlled calibration model applied to analytical chemistry

In chemical analyses performed by laboratories, one faces the problem of determining the concentration of a chemical element in a sample. In practice, one deals with the problem using the so-called linear calibration model, which considers that the errors associated with the independent variables are negligible compared with the former variable. In this work, a new linear calibration model is proposed assuming that the independent variables are subject to heteroscedastic measurement errors. A simulation study is carried out in order to verify some properties of the estimators derived for the new model and it is also considered the usual calibration model to compare it with the new approach. Three applications are considered to verify the performance of the new approach. Copyright (C) 2010 John Wiley & Sons, Ltd.

Generalized Birnbaum-Saunders distributions applied to air pollutant concentration

The generalized Birnbaum-Saunders (GBS) distribution is a new class of positively skewed models with lighter and heavier tails than the traditional Birnbaum-Saunders (BS) distribution, which is largely applied to study lifetimes. However, the theoretical argument and the interesting properties of the GBS model have made its application possible beyond the lifetime analysis. The aim of this paper is to present the GBS distribution as a useful model for describing pollution data and deriving its positive and negative moments. Based on these moments, we develop estimation and goodness-of-fit methods. Also, some properties of the proposed estimators useful for developing asymptotic inference are presented. Finally, an application with real data from Environmental Sciences is given to illustrate the methodology developed. This example shows that the empirical fit of the GBS distribution to the data is very good. Thus, the GBS model is appropriate for describing air pollutant concentration data, which produces better results than the lognormal model when the administrative target is determined for abating air pollution. Copyright (c) 2007 John Wiley & Sons, Ltd.

Bayesian Analysis for the Generalized Lognormal Distribution Applied to Failure Time Analysis

There are several versions of the lognormal distribution in the statistical literature, one is based in the exponential transformation of generalized normal distribution (GN). This paper presents the Bayesian analysis for the generalized lognormal distribution (logGN) considering independent non-informative Jeffreys distributions for the parameters as well as the procedure for implementing the Gibbs sampler to obtain the posterior distributions of parameters. The results are used to analyze failure time models with right-censored and uncensored data. The proposed method is illustrated using actual failure time data of computers.

Oral health conditions among the elderly in Southeastern São Paulo State

OBJECTIVE: This study aimed to asses oral health conditions in a population aged 60 years and over living in Botucatu, Southeastern Brazil. A cross-sectional population-based study was carried out using a random sample (N=372) of the urban population aged 60 years and over from the city of Botucatu, in 2005. World Health Organization criteria and codes for oral health epidemiological surveys were used. Re-examination was carried out in 10% of individuals aiming to evaluate intra-examiner agreement. Statistical analysis was performed by one-way ANOVA or Kruskal-Wallis ANOVA, as applicable. Also, the t-test was used in the absence of homoscedasticity. Fisher's exact test was used for situations where the categories with less than five units were observed. Adjusted residuals and multiple-comparison analysis were conducted to identify associations between variable categories and subgroups. The intra-examiner agreement was 98% and Kappa statistics result was 0.95. Loss component represented 90.68% of DMF-T index, which was 29.85. The prevalence of edentulism was 63.17%. Upper and lower dentures were found in 80% and 58% respectively, with complete denture as the most commonly used. In those studied, 15% required upper and 38% lower dentures. There was more need for complete denture in both jaws. Approximately 20% had soft tissue alterations. For periodontal conditions, most sextants were excluded (81.81%). Periodontal pockets (4 - 5 mm) were seen in 11.29% of the examined individuals. The oral health status of the elderly population in Botucatu is poor, as well as in other Brazilian cities. The results of this study may help planning collective health actions, giving an accurate description of the oral problems among the elderly.

Evaluation of simulated bone lesion in the head of the mandible by using multislice computed tomography

Conventional radiography has shown limitation in acquiring image of the ATM region, thus, computed tomography (CT) scanning has been the best option to the present date for diagnosis, surgical planning and treatment of bone lesions, owing to its specific properties. OBJECTIVE: The aim of the study was to evaluate images of simulated bone lesions at the head of the mandible by multislice CT. MATERIAL AND METHODS: Spherical lesions were made with dental spherical drills (sizes 1, 3, and 6) and were evaluated by using multislice CT (64 rows), by two observers in two different occasions, deploying two protocols: axial, coronal, and sagittal images, and parasagittal images for pole visualization (anterior, lateral, posterior, medial and superior). Acquired images were then compared with those lesions in the dry mandible (gold standard) to evaluate the specificity and sensibility of both protocols. Statistical methods included: Kappa statistics, validity test and chi-square test. Results demonstrated the advantage of associating axial, coronal, and sagittal slices with parasagittal slices for lesion detection at the head of the mandible. RESULTS: There was no statistically significant difference between the types of protocols regarding a particular localization of lesions at the poles. CONCLUSIONS: Protocols for the assessment of the head of the mandible were established to improve the visualization of alterations of each of the poles of the mandible's head. The anterior and posterior poles were better visualized in lateral-medial planes while lateral, medial and superior poles were better visualized in the anterior-posterior plane.

Diagnostic agreement between panoramic radiographs and color doppler images of carotid atheroma

The aim of this study was to investigate the agreement between diagnoses of calcified atheroma seen on panoramic radiographs and color Doppler images. Our interest stems from the fact that panoramic images can show the presence of atheroma regardless of the level of obstruction detected by color Doppler images. Panoramic and color Doppler images of 16 patients obtained from the archives of the Health Department of the city of Valença, RJ, Brazil, were analyzed in this study. Both sides of each patient were observed on the images, with a total of 32 analyzed cervical regions. The level of agreement between diagnoses was analyzed using the Kappa statistics. There was a high level of agreement, with a Kappa value of 0.78. In conclusion, panoramic radiographs can help detecting calcifications in the cervical region of patients susceptible to vascular diseases predisposing to myocardial infarction and cerebrovascular accidents. If properly trained and informed, dentists can refer their patients to a physician for a cardiovascular evaluation in order to receive proper and timely medical treatment.

Evaluating resin-enamel bonds by microshear and microtensile bond strength tests: effects of composite resin

OBJECTIVES: The aims of this study were to evaluate the effect of resin composite (Filtek Z250 and Filtek Flow Z350) and adhesive system [(Solobond Plus, Futurabond NR (VOCO) and Adper Single Bond (3M ESPE)] on the microtensile (μTBS) and microshear bond strength (μSBS) tests on enamel, and to correlate the bond strength means between them. MATERIAL AND METHODS: Thirty-six extracted human molars were sectioned to obtain two tooth halves: one for μTBS and the other one for μSBS. Adhesive systems and resin composites were applied to the enamel ground surfaces and light-cured. After storage (37(0)C/24 h) specimens were stressed (0.5 mm/min). Fracture modes were analyzed under scanning electron microscopy. The data were analyzed using two-way ANOVA and Tukey's test (α=0.05). RESULTS: The correlation between tests was estimated with Pearson's product-moment correlation statistics (α =0.05). For both tests only the main factor resin composite was statistically significant (p<0.05). The correlation test detected a positive (r=0.91) and significant (p=0.01) correlation between the tests. CONCLUSIONS: The results were more influenced by the resin type than by the adhesives. Both microbond tests seem to be positive and linearly correlated and can therefore lead to similar conclusions.

Origem e dispersão dos Tupiguarani: o que diz a morfologia craniana?

A origem e a dispersão dos povos Tupiguarani têm sido intensamente debatidas entre arqueólogos e linguistas nas últimas cinco décadas. Em resumo, pode-se dizer que a ideia de que esses povos, que ocuparam grande parte do território brasileiro e parte da Bolívia, do Paraguai, do Uruguai e da Argentina, tiveram sua etnogênese na Amazônia e dali partiram para o leste e para o sul, por volta de 2.500 anos antes do presente, é bastante aceita entre os especialistas, embora uma dispersão no sentido oposto, isto é, do sul para o norte, com origem na bacia do Tietê-Paraná, não seja completamente descartada. Entre os arqueólogos que consideram a Amazônia como berço desses povos, alguns acreditam que esse surgimento se deu na Amazônia central. Outros acreditam que a etnogênese Tupiguarani ocorreu no sudoeste da Amazônia, onde hoje se concentra a maior diversidade linguística do tronco Tupi. Neste trabalho, a morfologia de 19 crânios associados à cerâmica Tupiguarani ou etnograficamente classificados como tais foram comparados a várias séries cranianas pré-históricas e etnográficas brasileiras por meio de estatísticas multivariadas. Duas técnicas multivariadas foram empregadas: Análise de Componentes Principais, aplicada sobre os centróides de cada série, e Distâncias de Mahalanobis, aplicadas aos dados individuais. Os resultados obtidos sugerem uma origem amazônica para os povos Tupiguarani, sobretudo pela forte associação encontrada entre crânios Tupi e Guarani do sudeste e do sul brasileiro e dos Tupi do norte do Brasil, com os espécimes provenientes da ilha de Marajó incluídos no estudo.