An R implementation for generalized Birnbaum-Saunders distributions

The Birnbaum-Saunders (BS) model is a positively skewed statistical distribution that has received great attention in recent decades. A generalized version of this model was derived based on symmetrical distributions in the real line named the generalized BS (GBS) distribution. The R package named gbs was developed to analyze data from GBS models. This package contains probabilistic and reliability indicators and random number generators from GBS distributions. Parameter estimates for censored and uncensored data can also be obtained by means of likelihood methods from the gbs package. Goodness-of-fit and diagnostic methods were also implemented in this package in order to check the suitability of the GBS models. in this article, the capabilities and features of the gbs package are illustrated by using simulated and real data sets. Shape and reliability analyses for GBS models are presented. A simulation study for evaluating the quality and sensitivity of the estimation method developed in the package is provided and discussed. (C) 2008 Elsevier B.V. All rights reserved.

A new class of survival regression models with heavy-tailed errors: robustness and diagnostics

Birnbaum-Saunders models have largely been applied in material fatigue studies and reliability analyses to relate the total time until failure with some type of cumulative damage. In many problems related to the medical field, such as chronic cardiac diseases and different types of cancer, a cumulative damage caused by several risk factors might cause some degradation that leads to a fatigue process. In these cases, BS models can be suitable for describing the propagation lifetime. However, since the cumulative damage is assumed to be normally distributed in the BS distribution, the parameter estimates from this model can be sensitive to outlying observations. In order to attenuate this influence, we present in this paper BS models, in which a Student-t distribution is assumed to explain the cumulative damage. In particular, we show that the maximum likelihood estimates of the Student-t log-BS models attribute smaller weights to outlying observations, which produce robust parameter estimates. Also, some inferential results are presented. In addition, based on local influence and deviance component and martingale-type residuals, a diagnostics analysis is derived. Finally, a motivating example from the medical field is analyzed using log-BS regression models. Since the parameter estimates appear to be very sensitive to outlying and influential observations, the Student-t log-BS regression model should attenuate such influences. The model checking methodologies developed in this paper are used to compare the fitted models.

Beta-binomial/Poisson regression models for repeated bivariate counts

We analyze data obtained from a study designed to evaluate training effects on the performance of certain motor activities of Parkinson`s disease patients. Maximum likelihood methods were used to fit beta-binomial/Poisson regression models tailored to evaluate the effects of training on the numbers of attempted and successful specified manual movements in 1 min periods, controlling for disease stage and use of the preferred hand. We extend models previously considered by other authors in univariate settings to account for the repeated measures nature of the data. The results suggest that the expected number of attempts and successes increase with training, except for patients with advanced stages of the disease using the non-preferred hand. Copyright (c) 2008 John Wiley & Sons, Ltd.

MODELING THE EVOLUTION OF COMPLEX NETWORKS THROUGH THE PATH-STAR TRANSFORMATION AND OPTIMAL MULTIVARIATE METHODS

The topology of real-world complex networks, such as in transportation and communication, is always changing with time. Such changes can arise not only as a natural consequence of their growth, but also due to major modi. cations in their intrinsic organization. For instance, the network of transportation routes between cities and towns ( hence locations) of a given country undergo a major change with the progressive implementation of commercial air transportation. While the locations could be originally interconnected through highways ( paths, giving rise to geographical networks), transportation between those sites progressively shifted or was complemented by air transportation, with scale free characteristics. In the present work we introduce the path-star transformation ( in its uniform and preferential versions) as a means to model such network transformations where paths give rise to stars of connectivity. It is also shown, through optimal multivariate statistical methods (i.e. canonical projections and maximum likelihood classification) that while the US highways network adheres closely to a geographical network model, its path-star transformation yields a network whose topological properties closely resembles those of the respective airport transportation network.

A fast and robust statistical test based on likelihood ratio with Bartlett correction to identify Granger causality between gene sets

We propose a likelihood ratio test ( LRT) with Bartlett correction in order to identify Granger causality between sets of time series gene expression data. The performance of the proposed test is compared to a previously published bootstrapbased approach. LRT is shown to be significantly faster and statistically powerful even within non- Normal distributions. An R package named gGranger containing an implementation for both Granger causality identification tests is also provided.

Comparasion of two minimally invasive methods on the longevity of glass ionomer cement restorations: short-term results of a pilot study

The purpose of this study was to evaluate the clinical performance of glass ionomer cement (GIC) restorations comparing two minimally invasive methods in permanent teeth after 12 months. Fifty pregnant women (second trimester of pregnancy), mean age 22 ± 5.30 years, were treated by two previously trained operators. The treatment approaches tested were: chemomechanical method (CarisolvTM; MediTeam) and atraumatic restorative treatment (ART). A split-mouth study design was used in which the two treatments were randomly placed in 50 matched pairs of permanent teeth. The chemomechanical method (CM) was the test group and the ART was the control group. The treatments were performed in Public Health Centers. The tested restorative material was a high-strength GIC (Ketac Molar; 3M/ESPE). The restorations were placed according to the ART guidelines. Two calibrated independent examiners evaluated the restorations in accordance with ART criteria. The inter-examiner kappa was 0.97. Data were analyzed using 95% confidence interval on the binomial distribution and Fisher's exact test at 5% significance level. In a 12-month follow-up, 86% of the restorations were evaluated. In the test group (CM), 100% (CI=93.3-100%) of the restorations were considered successful. In the control group (ART) 97.6% (CI=87.4-99.9%) of the restorations were considered successful and 2.4% unsuccessful (marginal defect >0.5 mm). There was no statistically significant difference between the 12-mounth success rate for both groups (Fisher's exact test: P=0.49) and between the two operators (Fisher's exact test: P=1.00). Both minimally invasive methods, chemomechanical method and ART, showed a similar clinical performance after 12 months of follow up.

Histopathological evaluation of different methods of experimental induction of periapical periodontitis

This study evaluated histopathologically different methods of experimental induction of periapical periodontitis. The radiographic and microbiological evaluations have been performed in a previous investigation. Fifty-seven root canals from dogs' teeth were assigned to 4 groups. In GI (n=14) and GII (n=14), the root canals were exposed to oral environment for 180 days; in GIII (n=14) and GIV (n=15) the root canals were exposed for 7 days and then the access cavities were restored and remained sealed for 53 days. The root apices of GI and GIII were perforated, whilst those of GII and GIV remained intact. After induction of periapical periodontitis, the dogs were euthanized. Serial sections were obtained and stained with hematoxylin and eosin. Data of the histopathological evaluation were submitted to Kruskal-Wallis and Dunn's tests at 5% significance level. The inflammatory periapical reaction and resorption of mineralized tissues were less intense in GII than in the other groups (p<0.05). There was no histopathological difference among the experimentally induced periapical lesions in the teeth with coronal sealing. On the other hand, when coronal sealing was not performed, greater intensity of induced periapical periodontitis was observed in the teeth with apical perforation.

Microbial distribution in the root canal system after periapical lesion induction using different methods

The aim of this study was to evaluate the microbial distribution in the root canal system after periapical lesion induction in dogs' teeth using different methods. Fifty-two root canals were assigned to 4 groups (n=13). Groups I and II: root canals were exposed to the oral cavity for 180 days; groups III and IV: root canals were exposed for 7 days and then the coronal openings were sealed for 53 days. The root apices of groups I and III were perforated, while those of groups II and IV remained intact. After the experimental periods, the animals were euthanized and the anatomic pieces containing the roots were processed and stained with the Brown & Brenn method to assess the presence and distribution of microorganisms. The incidence of microorganisms at different sites of the roots and periapical lesions was analyzed statistically by the chi-square test at 5% significance level. All groups presented microorganisms in the entire root canal system. A larger number of microorganisms was observed on the root canal walls, apical delta and dentinal tubules (p<0.05), followed by cementum and cemental resorption areas. In spite of the different periods of exposure to the oral environment, the methods used for induction of periapical periodontitis yielded similar distribution of microorganisms in the root canal system.

Comparative study of two commercially pure titanium casting methods

The interest in using titanium to fabricate removable partial denture (RPD) frameworks has increased, but there are few studies evaluating the effects of casting methods on clasp behavior. OBJECTIVE: This study compared the occurrence of porosities and the retentive force of commercially pure titanium (CP Ti) and cobalt-chromium (Co-Cr) removable partial denture circumferential clasps cast by induction/centrifugation and plasma/vacuum-pressure. MATERIAL AND METHODS: 72 frameworks were cast from CP Ti (n=36) and Co-Cr alloy (n=36; control group). For each material, 18 frameworks were casted by electromagnetic induction and injected by centrifugation, whereas the other 18 were casted by plasma and injected by vacuum-pressure. For each casting method, three subgroups (n=6) were formed: 0.25 mm, 0.50 mm, and 0.75 mm undercuts. The specimens were radiographed and subjected to an insertion/removal test simulating 5 years of framework use. Data were analyzed by ANOVA and Tukey's to compare materials and cast methods (α=0.05). RESULTS: Three of 18 specimens of the induction/centrifugation group and 9 of 18 specimens of plasma/vacuum-pressure cast presented porosities, but only 1 and 7 specimens, respectively, were rejected for simulation test. For Co-Cr alloy, no defects were found. Comparing the casting methods, statistically significant differences (p<0.05) were observed only for the Co-Cr alloy with 0.25 mm and 0.50 mm undercuts. Significant differences were found for the 0.25 mm and 0.75 mm undercuts dependent on the material used. For the 0.50 mm undercut, significant differences were found when the materials were induction casted. CONCLUSION: Although both casting methods produced satisfactory CP Ti RPD frameworks, the occurrence of porosities was greater in the plasma/vacuum-pressure than in the induction/centrifugation method, the latter resulting in higher clasp rigidity, generating higher retention force values.

Comparative study of intermaxillary relationships of manual and swallowing methods

This study compared the mandibular displacement from three methods of centric relation record using an anterior jig associated with (A) chin point guidance, (B) swallowing (control group) and (C) bimanual manipulation. Ten patients aged 25-39 years were selected if they met the following inclusion criteria: complete dentition (up to the second molars), Angle class I and absence of signs and symptoms of temporomandibular disorders and diagnostic casts showing stability in the maximum intercuspation (MI) position. Impressions of maxillary and mandibular arches were made with an irreversible hydrocolloid impression material. Master casts of each patient were obtained, mounted on a microscope table in MI as a reference position and 5 records of each method were made per patient. The mandibular casts were then repositioned with records interposed and new measurements were obtained. The difference between the two readings allowed measuring the displacement of the mandible in the anteroposterior and lateral axes. Data were analyzed statistically by ANOVA and Tukey's test at 5% significance level. There was no statistically significant differences (p>0.05) among the three methods for measuring lateral displacement (A=0.38 ± 0.26, B=0.32 ± 0.25 and C=0.32 ± 0.23). For the anteroposterior displacement (A=2.76 ± 1.43, B=2.46 ± 1.48 and C=2.97 ± 1.51), the swallowing method (B) differed significantly from the others (p<0.05), but no significant difference (p>0.05) was found between chin point guidance (A) and bimanual manipulation (C). In conclusion, the swallowing method produced smaller mandibular posterior displacement than the other methods.

Validity and reliability of methods for the detection of secondary caries around amalgam restorations in primary teeth

Secondary caries has been reported as the main reason for restoration replacement. The aim of this in vitro study was to evaluate the performance of different methods - visual inspection, laser fluorescence (DIAGNOdent), radiography and tactile examination - for secondary caries detection in primary molars restored with amalgam. Fifty-four primary molars were photographed and 73 suspect sites adjacent to amalgam restorations were selected. Two examiners evaluated independently these sites using all methods. Agreement between examiners was assessed by the Kappa test. To validate the methods, a caries-detector dye was used after restoration removal. The best cut-off points for the sample were found by a Receiver Operator Characteristic (ROC) analysis, and the area under the ROC curve (Az), and the sensitivity, specificity and accuracy of the methods were calculated for enamel (D2) and dentine (D3) thresholds. These parameters were found for each method and then compared by the McNemar test. The tactile examination and visual inspection presented the highest inter-examiner agreement for the D2 and D3 thresholds, respectively. The visual inspection also showed better performance than the other methods for both thresholds (Az = 0.861 and Az = 0.841, respectively). In conclusion, the visual inspection presented the best performance for detecting enamel and dentin secondary caries in primary teeth restored with amalgam.

Reliability and discriminatory power of methods for dental plaque quantification

OBJECTIVE: This in situ study evaluated the discriminatory power and reliability of methods of dental plaque quantification and the relationship between visual indices (VI) and fluorescence camera (FC) to detect plaque. MATERIAL AND METHODS: Six volunteers used palatal appliances with six bovine enamel blocks presenting different stages of plaque accumulation. The presence of plaque with and without disclosing was assessed using VI. Images were obtained with FC and digital camera in both conditions. The area covered by plaque was assessed. Examinations were done by two independent examiners. Data were analyzed by Kruskal-Wallis and Kappa tests to compare different conditions of samples and to assess the inter-examiner reproducibility. RESULTS: Some methods presented adequate reproducibility. The Turesky index and the assessment of area covered by disclosed plaque in the FC images presented the highest discriminatory powers. CONCLUSION: The Turesky index and images with FC with disclosing present good reliability and discriminatory power in quantifying dental plaque.

New methods to quantify NH3 volatilization from fertilized surface soil with urea

Gaseous N losses from soil are considerable, resulting mostly from ammonia volatilization linked to agricultural activities such as pasture fertilization. The use of simple and accessible measurement methods of such losses is fundamental in the evaluation of the N cycle in agricultural systems. The purpose of this study was to evaluate quantification methods of NH3 volatilization from fertilized surface soil with urea, with minimal influence on the volatilization processes. The greenhouse experiment was arranged in a completely randomized design with 13 treatments and five replications, with the following treatments: (1) Polyurethane foam (density 20 kg m-3) with phosphoric acid solution absorber (foam absorber), installed 1, 5, 10 and 20 cm above the soil surface; (2) Paper filter with sulfuric acid solution absorber (paper absorber, 1, 5, 10 and 20 cm above the soil surface); (3) Sulfuric acid solution absorber (1, 5 and 10 cm above the soil surface); (4) Semi-open static collector; (5) 15N balance (control). The foam absorber placed 1 cm above the soil surface estimated the real daily rate of loss and accumulated loss of NH3N and proved efficient in capturing NH3 volatized from urea-treated soil. The estimates based on acid absorbers 1, 5 and 10 cm above the soil surface and paper absorbers 1 and 5 cm above the soil surface were only realistic for accumulated N-NH3 losses. Foam absorbers can be indicated to quantify accumulated and daily rates of NH3 volatilization losses similarly to an open static chamber, making calibration equations or correction factors unnecessary.

Critical applications of SW 846 US EPA methods to evaluation of marine samples quality

Technical evaluation of analytical data is of extreme relevance considering it can be used for comparisons with environmental quality standards and decision-making as related to the management of disposal of dredged sediments and the evaluation of salt and brackish water quality in accordance with CONAMA 357/05 Resolution. It is, therefore, essential that the project manager discusses the environmental agency's technical requirements with the laboratory contracted for the follow-up of the analysis underway and even with a view to possible re-analysis when anomalous data are identified. The main technical requirements are: (1) method quantitation limits (QLs) should fall below environmental standards; (2) analyses should be carried out in laboratories whose analytical scope is accredited by the National Institute of Metrology (INMETRO) or qualified or accepted by a licensing agency; (3) chain of custody should be provided in order to ensure sample traceability; (4) control charts should be provided to prove method performance; (5) certified reference material analysis or, if that is not available, matrix spike analysis, should be undertaken and (6) chromatograms should be included in the analytical report. Within this context and with a view to helping environmental managers in analytical report evaluation, this work has as objectives the discussion of the limitations of the application of SW 846 US EPA methods to marine samples, the consequences of having data based on method detection limits (MDL) and not sample quantitation limits (SQL), and present possible modifications of the principal method applied by laboratories in order to comply with environmental quality standards.