Modelos para dados grupados e sensurados aplicados à área biológica: comparação usando fator de Bayes

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

Efeito da quantidade de habitat, do tamanho do fragmento e da posição dos recursos sobre a frugivoria em ambientes fragmentados

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

Efeito do gradiente interior-borda-matriz, e do tipo de matriz sobre a entomofaúna edáfica em fragmentos de cerrado

edge effect. Thus, under the influence of the adjacent matrix, fragments undergo microclimatic alterations that accentuate changes in species composition and community structure. In order to better understand edge and matrix effects on the richness and abundance of edaphic arthropods, this study assessed: (a) the difference between habitat (fragment) and non-habitat (matrix); (b) whether there is a continuous interior-edge-matrix gradient; and (c) the difference between matrices for arthropod orders richness and abundance. We selected 15 landscapes, 5 of which contained a cerrado fragment surrounded by sugarcane cultivation, 5 with a cerrado fragment within eucalyptus and 5 with a cerrado fragment within pasture. In each landscape the soil fauna was collected along with the soil and then extracted with the aid of the modified Berlese-Tullgren funnel. We chose the orders Coleoptera, Collembola, Mesostigmata and Oribatida for analysis, and after separation of the individuals we used model selection analysis via AIC. The model type fragment x matrix was the most likely to explain richness, total and relative abundances of the four orders (wAICc between 0,6623 and 1,0). The model of edge distance (edge effect) was plausible to total abundance and relative abundance of Mesostigmata order (wAICc=0,2717 and 0,186). Local environmental variables (soil texture, temperature and relative humidity), and fragment size were also measured to avoid confounding factors and were not presented as plausible models to explain the patterns. So edaphic arthropods, despite protecting themselves under the ground, are extremely sensitive to fragmentation, even with the replacement of natural habitat by agricultural use, such as sugarcane, pasture and eucalyptus. This group should be studied environmental impact assessments because provides important ecosystem se ravincde s inacnludd eisd ainn efficient bio-indicator

Efeito da quantidade de habitat, do tamanho do fragmento e da posição dos recursos sobre a frugivoria em ambientes fragmentados

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

Efeito do gradiente interior-borda-matriz, e do tipo de matriz sobre a entomofaúna edáfica em fragmentos de cerrado

edge effect. Thus, under the influence of the adjacent matrix, fragments undergo microclimatic alterations that accentuate changes in species composition and community structure. In order to better understand edge and matrix effects on the richness and abundance of edaphic arthropods, this study assessed: (a) the difference between habitat (fragment) and non-habitat (matrix); (b) whether there is a continuous interior-edge-matrix gradient; and (c) the difference between matrices for arthropod orders richness and abundance. We selected 15 landscapes, 5 of which contained a cerrado fragment surrounded by sugarcane cultivation, 5 with a cerrado fragment within eucalyptus and 5 with a cerrado fragment within pasture. In each landscape the soil fauna was collected along with the soil and then extracted with the aid of the modified Berlese-Tullgren funnel. We chose the orders Coleoptera, Collembola, Mesostigmata and Oribatida for analysis, and after separation of the individuals we used model selection analysis via AIC. The model type fragment x matrix was the most likely to explain richness, total and relative abundances of the four orders (wAICc between 0,6623 and 1,0). The model of edge distance (edge effect) was plausible to total abundance and relative abundance of Mesostigmata order (wAICc=0,2717 and 0,186). Local environmental variables (soil texture, temperature and relative humidity), and fragment size were also measured to avoid confounding factors and were not presented as plausible models to explain the patterns. So edaphic arthropods, despite protecting themselves under the ground, are extremely sensitive to fragmentation, even with the replacement of natural habitat by agricultural use, such as sugarcane, pasture and eucalyptus. This group should be studied environmental impact assessments because provides important ecosystem se ravincde s inacnludd eisd ainn efficient bio-indicator

A model of road effect using line integrals and a test of the performance of two new road indices using the distribution of small mammals in an Atlantic Forest landscape

Effects of roads on wildlife and its habitat have been measured using metrics, such as the nearest road distance, road density, and effective mesh size. In this work we introduce two new indices: (1) Integral Road Effect (IRE), which measured the sum effects of points in a road at a fixed point in the forest; and (2) Average Value of the Infinitesimal Road Effect (AVIRE), which measured the average of the effects of roads at this point. IRE is formally defined as the line integral of a special function (the infinitesimal road effect) along the curves that model the roads, whereas AVIRE is the quotient of IRE by the length of the roads. Combining tools of ArcGIS software with a numerical algorithm, we calculated these and other road and habitat cover indices in a sample of points in a human-modified landscape in the Brazilian Atlantic Forest, where data on the abundance of two groups of small mammals (forest specialists and habitat generalists) were collected in the field. We then compared through the Akaike Information Criterion (AIC) a set of candidate regression models to explain the variation in small mammal abundance, including models with our two new road indices (AVIRE and IRE) or models with other road effect indices (nearest road distance, mesh size, and road density), and reference models (containing only habitat indices, or only the intercept without the effect of any variable). Compared to other road effect indices, AVIRE showed the best performance to explain abundance of forest specialist species, whereas the nearest road distance obtained the best performance to generalist species. AVIRE and habitat together were included in the best model for both small mammal groups, that is, higher abundance of specialist and generalist small mammals occurred where there is lower average road effect (less AVIRE) and more habitat. Moreover, AVIRE was not significantly correlated with habitat cover of specialists and generalists differing from the other road effect indices, except mesh size, which allows for separating the effect of roads from the effect of habitat on small mammal communities. We suggest that the proposed indices and GIS procedures could also be useful to describe other spatial ecological phenomena, such as edge effect in habitat fragments. (C) 2012 Elsevier B.V. All rights reserved.

Height-diameter relationships of tropical Atlantic moist forest trees in southeastern Brazil

Site-specific height-diameter models may be used to improve biomass estimates for forest inventories where only diameter at breast height (DBH) measurements are available. In this study, we fit height-diameter models for vegetation types of a tropical Atlantic forest using field measurements of height across plots along an altitudinal gradient. To fit height-diameter models, we sampled trees by DBH class and measured tree height within 13 one-hectare permanent plots established at four altitude classes. To select the best model we tested the performance of 11 height-diameter models using the Akaike Information Criterion (AIC). The Weibull and Chapman-Richards height-diameter models performed better than other models, and regional site-specific models performed better than the general model. In addition, there is a slight variation of height-diameter relationships across the altitudinal gradient and an extensive difference in the stature between the Atlantic and Amazon forests. The results showed the effect of altitude on tree height estimates and emphasize the need for altitude-specific models that produce more accurate results than a general model that encompasses all altitudes. To improve biomass estimation, the development of regional height-diameter models that estimate tree height using a subset of randomly sampled trees presents an approach to supplement surveys where only diameter has been measured.

The importance of plant life form on spatial associations along a subtropical coastal dune gradient

Questions Does the spatial association between isolated adult trees and understorey plants change along a gradient of sand dunes? Does this association depend on the life form of the understorey plant? Location Coastal sand dunes, southeast Brazil. Methods We recorded the occurrence of understorey plant species in 100 paired 0.25 m2 plots under adult trees and in adjacent treeless sites along an environmental gradient from beach to inland. Occurrence probabilities were modelled as a function of the fixed variables of the presence of a neighbour, distance from the seashore and life form, and a random variable, the block (i.e. the pair of plots). Generalized linear mixed models (GLMM) were fitted in a backward step-wise procedure using Akaike's information criterion (AIC) for model selection. Results The occurrence of understorey plants was affected by the presence of an adult tree neighbour, but the effect varied with the life form of the understorey species. Positive spatial association was found between isolated adult neighbour and young trees, whereas a negative association was found for shrubs. Moreover, a neutral association was found for lianas, whereas for herbs the effect of the presence of an adult neighbour ranged from neutral to negative, depended on the subgroup considered. The strength of the negative association with forbs increased with distance from the seashore. However, for the other life forms, the associational pattern with adult trees did not change along the gradient. Conclusions For most of the understorey life forms there is no evidence that the spatial association between isolated adult trees and understorey plants changes with the distance from the seashore, as predicted by the stress gradient hypothesis, a common hypothesis in the literature about facilitation in plant communities. Furthermore, the positive spatial association between isolated adult trees and young trees identified along the entire gradient studied indicates a positive feedback that explains the transition from open vegetation to forest in subtropical coastal dune environments.

Drag reduction phenomenon in viscous oil-water dispersed pipe flow: Experimental investigation and phenomenological modeling

An experimental study on drag-reduction phenomenon in dispersed oil-water flow has been performed in a 26-mm-i.d. Twelve meter long horizontal glass pipe. The flow was characterized using a novel wire-mesh sensor based on capacitance measurements and high-speed video recording. New two-phase pressure gradient, volume fraction, and phase distribution data have been used in the analysis. Drag reduction and slip ratio were detected at oil volume fractions between 10 and 45% and high mixture Reynolds numbers, and with water as the dominant phase. Phase-fraction distribution diagrams and cross-sectional imaging of the flow suggested the presence of a higher amount of water near to the pipe wall. Based on that, a phenomenology for explaining drag reduction in dispersed flow in a flow situation where slip ratio is significant is proposed. A simple phenomenological model is developed and the agreement between model predictions and data, including data from the literature, is encouraging. (c) 2011 American Institute of Chemical Engineers AIChE J, 2012

A prediction model for assessing residential radon concentration in Switzerland

Indoor radon is regularly measured in Switzerland. However, a nationwide model to predict residential radon levels has not been developed. The aim of this study was to develop a prediction model to assess indoor radon concentrations in Switzerland. The model was based on 44,631 measurements from the nationwide Swiss radon database collected between 1994 and 2004. Of these, 80% randomly selected measurements were used for model development and the remaining 20% for an independent model validation. A multivariable log-linear regression model was fitted and relevant predictors selected according to evidence from the literature, the adjusted R², the Akaike's information criterion (AIC), and the Bayesian information criterion (BIC). The prediction model was evaluated by calculating Spearman rank correlation between measured and predicted values. Additionally, the predicted values were categorised into three categories (50th, 50th-90th and 90th percentile) and compared with measured categories using a weighted Kappa statistic. The most relevant predictors for indoor radon levels were tectonic units and year of construction of the building, followed by soil texture, degree of urbanisation, floor of the building where the measurement was taken and housing type (P-values <0.001 for all). Mean predicted radon values (geometric mean) were 66 Bq/m³ (interquartile range 40-111 Bq/m³) in the lowest exposure category, 126 Bq/m³ (69-215 Bq/m³) in the medium category, and 219 Bq/m³ (108-427 Bq/m³) in the highest category. Spearman correlation between predictions and measurements was 0.45 (95%-CI: 0.44; 0.46) for the development dataset and 0.44 (95%-CI: 0.42; 0.46) for the validation dataset. Kappa coefficients were 0.31 for the development and 0.30 for the validation dataset, respectively. The model explained 20% overall variability (adjusted R²). In conclusion, this residential radon prediction model, based on a large number of measurements, was demonstrated to be robust through validation with an independent dataset. The model is appropriate for predicting radon level exposure of the Swiss population in epidemiological research. Nevertheless, some exposure misclassification and regression to the mean is unavoidable and should be taken into account in future applications of the model.

Fitting of Mixtures with Unspecified Number of Components Using Cross Validation Distance Estimate

Estimation of the number of mixture components (k) is an unsolved problem. Available methods for estimation of k include bootstrapping the likelihood ratio test statistics and optimizing a variety of validity functionals such as AIC, BIC/MDL, and ICOMP. We investigate the minimization of distance between fitted mixture model and the true density as a method for estimating k. The distances considered are Kullback-Leibler (KL) and “L sub 2”. We estimate these distances using cross validation. A reliable estimate of k is obtained by voting of B estimates of k corresponding to B cross validation estimates of distance. This estimation methods with KL distance is very similar to Monte Carlo cross validated likelihood methods discussed by Smyth (2000). With focus on univariate normal mixtures, we present simulation studies that compare the cross validated distance method with AIC, BIC/MDL, and ICOMP. We also apply the cross validation estimate of distance approach along with AIC, BIC/MDL and ICOMP approach, to data from an osteoporosis drug trial in order to find groups that differentially respond to treatment.

On the Behaviour of Marginal and Conditional Akaike Information Criteria in Linear Mixed Models

In linear mixed models, model selection frequently includes the selection of random effects. Two versions of the Akaike information criterion (AIC) have been used, based either on the marginal or on the conditional distribution. We show that the marginal AIC is no longer an asymptotically unbiased estimator of the Akaike information, and in fact favours smaller models without random effects. For the conditional AIC, we show that ignoring estimation uncertainty in the random effects covariance matrix, as is common practice, induces a bias that leads to the selection of any random effect not predicted to be exactly zero. We derive an analytic representation of a corrected version of the conditional AIC, which avoids the high computational cost and imprecision of available numerical approximations. An implementation in an R package is provided. All theoretical results are illustrated in simulation studies, and their impact in practice is investigated in an analysis of childhood malnutrition in Zambia.

IDH/MGMT-driven molecular classification of low-grade glioma is a strong predictor for long-term survival

BACKGROUND Low-grade gliomas (LGGs) are rare brain neoplasms, with survival spanning up to a few decades. Thus, accurate evaluations on how biomarkers impact survival among patients with LGG require long-term studies on samples prospectively collected over a long period. METHODS The 210 adult LGGs collected in our databank were screened for IDH1 and IDH2 mutations (IDHmut), MGMT gene promoter methylation (MGMTmet), 1p/19q loss of heterozygosity (1p19qloh), and nuclear TP53 immunopositivity (TP53pos). Multivariate survival analyses with multiple imputation of missing data were performed using either histopathology or molecular markers. Both models were compared using Akaike's information criterion (AIC). The molecular model was reduced by stepwise model selection to filter out the most critical predictors. A third model was generated to assess for various marker combinations. RESULTS Molecular parameters were better survival predictors than histology (ΔAIC = 12.5, P< .001). Forty-five percent of studied patients died. MGMTmet was positively associated with IDHmut (P< .001). In the molecular model with marker combinations, IDHmut/MGMTmet combined status had a favorable impact on overall survival, compared with IDHwt (hazard ratio [HR] = 0.33, P< .01), and even more so the triple combination, IDHmut/MGMTmet/1p19qloh (HR = 0.18, P< .001). Furthermore, IDHmut/MGMTmet/TP53pos triple combination was a significant risk factor for malignant transformation (HR = 2.75, P< .05). CONCLUSION By integrating networks of activated molecular glioma pathways, the model based on genotype better predicts prognosis than histology and, therefore, provides a more reliable tool for standardizing future treatment strategies.

The Application of Comorbidity Indices to Predict Early Postoperative Outcomes After Laparoscopic Roux-en-Y Gastric Bypass: A Nationwide Comparative Analysis of Over 70,000 Cases

BACKGROUND: Patients undergoing laparoscopic Roux-en-Y gastric bypass (LRYGB) often have substantial comorbidities, which must be taken into account to appropriately assess expected postoperative outcomes. The Charlson/Deyo and Elixhauser indices are widely used comorbidity measures, both of which also have revised algorithms based on enhanced ICD-9-CM coding. It is currently unclear which of the existing comorbidity measures best predicts early postoperative outcomes following LRYGB. METHODS: Using the Nationwide Inpatient Sample, patients 18 years or older undergoing LRYGB for obesity between 2001 and 2008 were identified. Comorbidities were assessed according to the original and enhanced Charlson/Deyo and Elixhauser indices. Using multivariate logistic regression, the following early postoperative outcomes were assessed: overall postoperative complications, length of hospital stay, and conversion to open surgery. Model performance for the four comorbidity indices was assessed and compared using C-statistics and the Akaike's information criterion (AIC). RESULTS: A total of 70,287 patients were included. Mean age was 43.1 years (SD, 10.8), 81.6 % were female and 60.3 % were White. Both the original and enhanced Elixhauser indices modestly outperformed the Charlson/Deyo in predicting the surgical outcomes. All four models had similar C-statistics, but the original Elixhauser index was associated with the smallest AIC for all of the surgical outcomes. CONCLUSIONS: The original Elixhauser index is the best predictor of early postoperative outcomes in our cohort of patients undergoing LRYGB. However, differences between the Charlson/Deyo and Elixhauser indices are modest, and each of these indices provides clinically relevant insight for predicting early postoperative outcomes in this high-risk patient population.

Combination of expression levels of miR-21 and miR-126 is associated with cancer-specific survival in clear-cell renal cell carcinoma

BACKGROUND Renal cell carcinoma (RCC) is marked by high mortality rate. To date, no robust risk stratification by clinical or molecular prognosticators of cancer-specific survival (CSS) has been established for early stages. Transcriptional profiling of small non-coding RNA gene products (miRNAs) seems promising for prognostic stratification. The expression of miR-21 and miR-126 was analysed in a large cohort of RCC patients; a combined risk score (CRS)-model was constructed based on expression levels of both miRNAs. METHODS Expression of miR-21 and miR-126 was evaluated by qRT-PCR in tumour and adjacent non-neoplastic tissue in n = 139 clear cell RCC patients. Relation of miR-21 and miR-126 expression with various clinical parameters was assessed. Parameters were analysed by uni- and multivariate COX regression. A factor derived from the z-score resulting from the COX model was determined for both miRs separately and a combined risk score (CRS) was calculated multiplying the relative expression of miR-21 and miR-126 by this factor. The best fitting COX model was selected by relative goodness-of-fit with the Akaike information criterion (AIC). RESULTS RCC with and without miR-21 up- and miR-126 downregulation differed significantly in synchronous metastatic status and CSS. Upregulation of miR-21 and downregulation of miR-126 were independently prognostic. A combined risk score (CRS) based on the expression of both miRs showed high sensitivity and specificity in predicting CSS and prediction was independent from any other clinico-pathological parameter. Association of CRS with CSS was successfully validated in a testing cohort containing patients with high and low risk for progressive disease. CONCLUSIONS A combined expression level of miR-21 and miR-126 accurately predicted CSS in two independent RCC cohorts and seems feasible for clinical application in assessing prognosis.