A landscape-scale test of the predictive ability of a spatially explicit model for population viability analysis

1. Although population viability analysis (PVA) is widely employed, forecasts from PVA models are rarely tested. This study in a fragmented forest in southern Australia contrasted field data on patch occupancy and abundance for the arboreal marsupial greater glider Petauroides volans with predictions from a generic spatially explicit PVA model. This work represents one of the first landscape-scale tests of its type. 2. Initially we contrasted field data from a set of eucalypt forest patches totalling 437 ha with a naive null model in which forecasts of patch occupancy were made, assuming no fragmentation effects and based simply on remnant area and measured densities derived from nearby unfragmented forest. The naive null model predicted an average total of approximately 170 greater gliders, considerably greater than the true count (n = 81). 3. Congruence was examined between field data and predictions from PVA under several metapopulation modelling scenarios. The metapopulation models performed better than the naive null model. Logistic regression showed highly significant positive relationships between predicted and actual patch occupancy for the four scenarios (P = 0.001-0.006). When the model-derived probability of patch occupancy was high (0.50-0.75, 0.75-1.00), there was greater congruence between actual patch occupancy and the predicted probability of occupancy. 4. For many patches, probability distribution functions indicated that model predictions for animal abundance in a given patch were not outside those expected by chance. However, for some patches the model either substantially over-predicted or under-predicted actual abundance. Some important processes, such as inter-patch dispersal, that influence the distribution and abundance of the greater glider may not have been adequately modelled. 5. Additional landscape-scale tests of PVA models, on a wider range of species, are required to assess further predictions made using these tools. This will help determine those taxa for which predictions are and are not accurate and give insights for improving models for applied conservation management.

Production under uncertainty and choice under uncertainty in the emergence of generalized expected utility theory

This paper presents a personal view of the interaction between the analysis of choice under uncertainty and the analysis of production under uncertainty. Interest in the foundations of the theory of choice under uncertainty was stimulated by applications of expected utility theory such as the Sandmo model of production under uncertainty. This interest led to the development of generalized models including rank-dependent expected utility theory. In turn, the development of generalized expected utility models raised the question of whether such models could be used in the analysis of applied problems such as those involving production under uncertainty. Finally, the revival of the state-contingent approach led to the recognition of a fundamental duality between choice problems and production problems.

Multicollinearity and financial constraint in investment decisions: a Bayesian generalized ridge regression

This paper addresses the investment decisions considering the presence of financial constraints of 373 large Brazilian firms from 1997 to 2004, using panel data. A Bayesian econometric model was used considering ridge regression for multicollinearity problems among the variables in the model. Prior distributions are assumed for the parameters, classifying the model into random or fixed effects. We used a Bayesian approach to estimate the parameters, considering normal and Student t distributions for the error and assumed that the initial values for the lagged dependent variable are not fixed, but generated by a random process. The recursive predictive density criterion was used for model comparisons. Twenty models were tested and the results indicated that multicollinearity does influence the value of the estimated parameters. Controlling for capital intensity, financial constraints are found to be more important for capital-intensive firms, probably due to their lower profitability indexes, higher fixed costs and higher degree of property diversification.

Characterization of pore wall heterogeneity in nanoporous carbons using adsorption: the slit pore model revisited

In this paper, we propose a new nonlocal density functional theory characterization procedure, the finite wall thickness model, for nanoporous carbons, whereby heterogeneity of pore size and pore walls in the carbon is probed simultaneously. We determine the pore size distributions and pore wall thickness distributions of several commercial activated carbons and coal chars, with good correspondence with X-ray diffraction. It is shown that the conventional infinite wall thickness approach overestimates the pore size slightly. Pore-pore correlation has been shown to have a negligible effect on prediction of pore size and pore wall thickness distributions for small molecules such as argon used in characterization. By utilizing the structural parameters (pore size and pore wall thickness distribution) in the generalized adsorption isotherm (GAI) we are able to predict adsorption uptake of supercritical gases in BPL and Norit RI Extra carbons, in excellent agreement with experimental adsorption uptake data up to 60 MPa. The method offers a useful technique for probing features of the solid skeleton, hitherto studied by crystallographic methods.

Use of pleopod morphology to determine sexual dimorphism and maturity in hermit crabs: Isocheles sawayai as a model

In the Anomura, studies on growth patterns are infrequent, possibly because the heterogeneity of the group, especially in terms of morphology, makes it difficult to construct generalized growth models. Particularly hermit crabs are an interesting group to evaluate aspects of growth, because of their unique body. Isocheles sawayai, a hermit crab found only in the western Atlantic Ocean, poorly known with respect to its sexual dimorphism and maturity, was investigated here based on morphometry. Monthly collections (July 2001 through June 2003) were made from a shrimp fishing boat in the Caraguatatuba region on the northern coast of the state of SA o pound Paulo, Brazil. The specimens were measured and weighed, and had their sex checked. Throughout the sampling period, 374 specimens of I. sawayai were collected (11.23% nonovigerous females, 6.69% ovigerous females, 79.41% males and 2.67% intersexes). The size at which morphological sexual maturity was reached by both sexes ranged from 4.0 to 4.3 mm shield length, according to the relative growth and the size of the smallest ovigerous female. Sexual dimorphism was shown by males, which were significantly larger than females, and by differences in growth pattern between the sexes, especially for relationships that involved the pleopods, which is related to their different functions in males and females. The present study is one of the first to use pleopod morphometry to determine sexual maturity and dimorphism in hermit crabs, especially for species with intersexuality such as I. sawayai.

Derivation and external validation of a simple prediction model for the diagnosis of type 2 Diabetes Mellitus in the Brazilian urban population

A risk score model was developed based in a population of 1,224 individuals from the general population without known diabetes aging 35 years or more from an urban Brazilian population sample in order to select individuals who should be screened in subsequent testing and improve the efficacy of public health assurance. External validation was performed in a second, independent, population from a different city ascertained through a similar epidemiological protocol. The risk score was developed by multiple logistic regression and model performance and cutoff values were derived from a receiver operating characteristic curve. Model`s capacity of predicting fasting blood glucose levels was tested analyzing data from a 5-year follow-up protocol conducted in the general population. Items independently and significantly associated with diabetes were age, BMI and known hypertension. Sensitivity, specificity and proportion of further testing necessary for the best cutoff value were 75.9, 66.9 and 37.2%, respectively. External validation confirmed the model`s adequacy (AUC equal to 0.72). Finally, model score was also capable of predicting fasting blood glucose progression in non-diabetic individuals in a 5-year follow-up period. In conclusion, this simple diabetes risk score was able to identify individuals with an increased likelihood of having diabetes and it can be used to stratify subpopulations in which performing of subsequent tests is necessary and probably cost-effective.

Early identification of leptospirosis-associated pulmonary hemorrhage syndrome by use of a validated prediction model

Objective: To identify prediction factors for the development of leptospirosis-associated pulmonary hemorrhage syndrome (LPHS). Methods: We conducted a prospective cohort study. The study comprised of 203 patients, aged >= 14 years, admitted with complications of the severe form of leptospirosis at the Emilio Ribas Institute of Infectology (Sao Paulo, Brazil) between 1998 and 2004. Laboratory and demographic data were obtained and the severity of illness score and involvement of the lungs and others organs were determined. Logistic regression was performed to identify independent predictors of LPHS. A prospective validation cohort of 97 subjects with severe form of leptospirosis admitted at the same hospital between 2004 and 2006 was used to independently evaluate the predictive value of the model. Results: The overall mortality rate was 7.9%. Multivariate logistic regression revealed that five factors were independently associated with the development of LPHS: serum potassium (mmol/L) (OR = 2.6; 95% CI = 1.1-5.9); serum creatinine (mmol/L) (OR = 1.2; 95% CI = 1.1-1.4); respiratory rate (breaths/min) (OR = 1.1; 95% CI = 1.1-1.2); presenting shock (OR = 69.9; 95% CI = 20.1-236.4), and Glasgow Coma Scale Score (GCS) < 15 (OR = 7.7; 95% CI = 1.3-23.0). We used these findings to calculate the risk of LPHS by the use of a spreadsheet. In the validation cohort, the equation classified correctly 92% of patients (Kappa statistic = 0.80). Conclusions: We developed and validated a multivariate model for predicting LPHS. This tool should prove useful in identifying LPHS patients, allowing earlier management and thereby reducing mortality. (C) 2009 The British Infection Society. Published by Elsevier Ltd. All rights reserved.

Involvement of median raphe nucleus 5-HT1A receptors in the regulation of generalized anxiety-related defensive behaviours in rats

Changes in 5-HT1A receptor-mediated neurotransmission at the level of the median raphe nucleus (MRN) are reported to affect the expression of defensive responses that are associated with generalized anxiety disorder (e.g. inhibitory avoidance) but not with panic (e.g. escape). The objective of this study was to further explore the involvement of MRN 5-HT1A receptors in the regulation of generalized anxiety-related behaviours. Results of experiment 1 showed that intra-MRN injection of the 5-HT1A/7 receptor agonist 8-OH-DPAT (0.6 nmol) in male Wistar rats impaired the acquisition of inhibitory avoidance, without interfering with the performance of escape in the elevated T-maze test of anxiety. Pre-treatment with the 5-HT1A receptor antagonist WAY-100635 (0.18 nmol) fully blocked this anxiolytic-like effect. As revealed by experiment 2, intra-MRN injection of 8-OH-DPAT (0.6, 3 or 15 nmol) also caused anxiolytic effect in rats submitted to the light-dark transition test, another animal model that has been associated with generalized anxiety. In the same test, intra-MRN injection of WAY-100635 (0.18, 0.37 or 0.74 nmol) caused the opposite effect. Overall, the current findings support the view that MRN 5-HT neurons, through the regulation of 5-HT1A somatodendritic autoreceptors, are implicated in the regulation of generalized anxiety-associated behaviours. (C) 2008 Elsevier Ireland Ltd. All rights reserved.

Involvement of the lateral habenula in the regulation of generalized anxiety- and panic-related defensive responses in rats

Recently obtained evidence points to the involvement of the lateral habenular nuclei (LHb) in the mediation of coping defensive responses to threatening/stressful stimuli. Nevertheless, the role of this brain area in the regulation of defensive responses that have been associated with specific subtypes of anxiety disorders recognized in clinical settings is presently unknown. To address this question, we investigated the effects of either electrolytic lesions or chemical stimulation of the LHb on the defensive behaviors generated in rats by the elevated T-maze. This experimental model allows the measurement, in a same rat, of two defensive behaviors, inhibitory avoidance and escape, that have been related in terms of psychopathology to generalized anxiety and panic disorders, respectively. Bilateral electrolytic lesions of the LHb (1 mA, 10 s) impaired inhibitory avoidance acquisition and facilitated escape performance. On the other hand, chemical stimulation of the LHb by bilateral microinjection of kainic acid (30-60 pmol/0.2 mu L) had the opposite effect, i.e., facilitated inhibitory avoidance and impaired escape. The present results indicate that the LHb exerts an opposed regulatory control on generalized anxiety- and panic-related defensive responses in rats. (c) 2008 Elsevier Inc. All rights reserved.

Estimating Income Inequality in China Using Grouped Data and the Generalized Beta Distribution

There are two main types of data sources of income distributions in China: household survey data and grouped data. Household survey data are typically available for isolated years and individual provinces. In comparison, aggregate or grouped data are typically available more frequently and usually have national coverage. In principle, grouped data allow investigation of the change of inequality over longer, continuous periods of time, and the identification of patterns of inequality across broader regions. Nevertheless, a major limitation of grouped data is that only mean (average) income and income shares of quintile or decile groups of the population are reported. Directly using grouped data reported in this format is equivalent to assuming that all individuals in a quintile or decile group have the same income. This potentially distorts the estimate of inequality within each region. The aim of this paper is to apply an improved econometric method designed to use grouped data to study income inequality in China. A generalized beta distribution is employed to model income inequality in China at various levels and periods of time. The generalized beta distribution is more general and flexible than the lognormal distribution that has been used in past research, and also relaxes the assumption of a uniform distribution of income within quintile and decile groups of populations. The paper studies the nature and extent of inequality in rural and urban China over the period 1978 to 2002. Income inequality in the whole of China is then modeled using a mixture of province-specific distributions. The estimated results are used to study the trends in national inequality, and to discuss the empirical findings in the light of economic reforms, regional policies, and globalization of the Chinese economy.

Long-term survivor mixture model with random effects: application to a multi-centre clinical trial of carcinoma

A mixture model incorporating long-term survivors has been adopted in the field of biostatistics where some individuals may never experience the failure event under study. The surviving fractions may be considered as cured. In most applications, the survival times are assumed to be independent. However, when the survival data are obtained from a multi-centre clinical trial, it is conceived that the environ mental conditions and facilities shared within clinic affects the proportion cured as well as the failure risk for the uncured individuals. It necessitates a long-term survivor mixture model with random effects. In this paper, the long-term survivor mixture model is extended for the analysis of multivariate failure time data using the generalized linear mixed model (GLMM) approach. The proposed model is applied to analyse a numerical data set from a multi-centre clinical trial of carcinoma as an illustration. Some simulation experiments are performed to assess the applicability of the model based on the average biases of the estimates formed. Copyright (C) 2001 John Wiley & Sons, Ltd.

Three-dimensional orthotropic viscoelastic finite element model of a human ligament

Ligaments undergo finite strain displaying hyperelastic behaviour as the initially tangled fibrils present straighten out, combined with viscoelastic behaviour (strain rate sensitivity). In the present study the anterior cruciate ligament of the human knee joint is modelled in three dimensions to gain an understanding of the stress distribution over the ligament due to motion imposed on the ends, determined from experimental studies. A three dimensional, finite strain material model of ligaments has recently been proposed by Pioletti in Ref. [2]. It is attractive as it separates out elastic stress from that due to the present strain rate and that due to the past history of deformation. However, it treats the ligament as isotropic and incompressible. While the second assumption is reasonable, the first is clearly untrue. In the present study an alternative model of the elastic behaviour due to Bonet and Burton (Ref. [4]) is generalized. Bonet and Burton consider finite strain with constant modulii for the fibres and for the matrix of a transversely isotropic composite. In the present work, the fibre modulus is first made to increase exponentially from zero with an invariant that provides a measure of the stretch in the fibre direction. At 12% strain in the fibre direction, a new reference state is then adopted, after which the material modulus is made constant, as in Bonet and Burton's model. The strain rate dependence can be added, either using Pioletti's isotropic approximation, or by making the effect depend on the strain rate in the fibre direction only. A solid model of a ligament is constructed, based on experimentally measured sections, and the deformation predicted using explicit integration in time. This approach simplifies the coding of the material model, but has a limitation due to the detrimental effect on stability of integration of the substantial damping implied by the nonlinear dependence of stress on strain rate. At present, an artificially high density is being used to provide stability, while the dynamics are being removed from the solution using artificial viscosity. The result is a quasi-static solution incorporating the effect of strain rate. Alternate approaches to material modelling and integration are discussed, that may result in a better model.

Integrable coupling in a model for Josephson tunneling between non-identical BCS systems

We extend a recent construction for an integrable model describing Josephson tunneling between identical BCS systems to the case where the BCS systems have different single particle energy levels. The exact solution of this generalized model is obtained through the Bethe ansatz.

A note on quasi-stationary distributions of birth-death processes and the SIS logistic epidemic

For Markov processes on the positive integers with the origin as an absorbing state, Ferrari, Kesten, Martinez and Picco studied the existence of quasi-stationary and limiting conditional distributions by characterizing quasi-stationary distributions as fixed points of a transformation Phi on the space of probability distributions on {1, 2,.. }. In the case of a birth-death process, the components of Phi(nu) can be written down explicitly for any given distribution nu. Using this explicit representation, we will show that Phi preserves likelihood ratio ordering between distributions. A conjecture of Kryscio and Lefevre concerning the quasi-stationary distribution of the SIS logistic epidemic follows as a corollary.

Bridging the gap between different statistical approaches: An integrated framework for modelling

This paper proposes a template for modelling complex datasets that integrates traditional statistical modelling approaches with more recent advances in statistics and modelling through an exploratory framework. Our approach builds on the well-known and long standing traditional idea of 'good practice in statistics' by establishing a comprehensive framework for modelling that focuses on exploration, prediction, interpretation and reliability assessment, a relatively new idea that allows individual assessment of predictions. The integrated framework we present comprises two stages. The first involves the use of exploratory methods to help visually understand the data and identify a parsimonious set of explanatory variables. The second encompasses a two step modelling process, where the use of non-parametric methods such as decision trees and generalized additive models are promoted to identify important variables and their modelling relationship with the response before a final predictive model is considered. We focus on fitting the predictive model using parametric, non-parametric and Bayesian approaches. This paper is motivated by a medical problem where interest focuses on developing a risk stratification system for morbidity of 1,710 cardiac patients given a suite of demographic, clinical and preoperative variables. Although the methods we use are applied specifically to this case study, these methods can be applied across any field, irrespective of the type of response.