On the use of correlated beta random variables with animal population modelling

We give reasons why demographic parameters such as survival and reproduction rates are often modelled well in stochastic population simulation using beta distributions. In practice, it is frequently expected that these parameters will be correlated, for example with survival rates for all age classes tending to be high or low in the same year. We therefore discuss a method for producing correlated beta random variables by transforming correlated normal random variables, and show how it can be applied in practice by means of a simple example. We also note how the same approach can be used to produce correlated uniform triangular, and exponential random variables. (C) 2008 Elsevier B.V. All rights reserved.

General results for the beta-modified Weibull distribution

We study in detail the so-called beta-modified Weibull distribution, motivated by the wide use of the Weibull distribution in practice, and also for the fact that the generalization provides a continuous crossover towards cases with different shapes. The new distribution is important since it contains as special sub-models some widely-known distributions, such as the generalized modified Weibull, beta Weibull, exponentiated Weibull, beta exponential, modified Weibull and Weibull distributions, among several others. It also provides more flexibility to analyse complex real data. Various mathematical properties of this distribution are derived, including its moments and moment generating function. We examine the asymptotic distributions of the extreme values. Explicit expressions are also derived for the chf, mean deviations, Bonferroni and Lorenz curves, reliability and entropies. The estimation of parameters is approached by two methods: moments and maximum likelihood. We compare by simulation the performances of the estimates from these methods. We obtain the expected information matrix. Two applications are presented to illustrate the proposed distribution.

Exponential Rates of Convergence in the Ergodic Theorem: A Constructive Approach

We prove that, once an algorithm of perfect simulation for a stationary and ergodic random field F taking values in S(Zd), S a bounded subset of R(n), is provided, the speed of convergence in the mean ergodic theorem occurs exponentially fast for F. Applications from (non-equilibrium) statistical mechanics and interacting particle systems are presented.

Expectation Propagation with Factorizing Distributions: A Gaussian Approximation and Performance Results for Simple Models

We discuss the expectation propagation (EP) algorithm for approximate Bayesian inference using a factorizing posterior approximation. For neural network models, we use a central limit theorem argument to make EP tractable when the number of parameters is large. For two types of models, we show that EP can achieve optimal generalization performance when data are drawn from a simple distribution.

Microscopic origin of non-Gaussian distributions of financial returns

In this paper we study the possible microscopic origin of heavy-tailed probability density distributions for the price variation of financial instruments. We extend the standard log-normal process to include another random component in the so-called stochastic volatility models. We study these models under an assumption, akin to the Born-Oppenheimer approximation, in which the volatility has already relaxed to its equilibrium distribution and acts as a background to the evolution of the price process. In this approximation, we show that all models of stochastic volatility should exhibit a scaling relation in the time lag of zero-drift modified log-returns. We verify that the Dow-Jones Industrial Average index indeed follows this scaling. We then focus on two popular stochastic volatility models, the Heston and Hull-White models. In particular, we show that in the Hull-White model the resulting probability distribution of log-returns in this approximation corresponds to the Tsallis (t-Student) distribution. The Tsallis parameters are given in terms of the microscopic stochastic volatility model. Finally, we show that the log-returns for 30 years Dow Jones index data is well fitted by a Tsallis distribution, obtaining the relevant parameters. (c) 2007 Elsevier B.V. All rights reserved.

Spontaneous Interblink Time Distributions in Patients with Graves` Orbitopathy and Normal Subjects

PURPOSE. To determine the shape of spontaneous interblink time interval distributions obtained in a long observation period in normal subjects and patients with Graves` orbitopathy. METHODS. The magnetic search coil technique was used to register the spontaneous blinking activity during 1 hour of video observation of two groups of 10 subjects each (normal controls aged 27-61 years, mean +/- SD = 46.0 +/- 13.6; patients with Graves` orbitopathy aged 33-61 years, mean +/- SD +/- 46.7 +/- 8.9). The spontaneous blink rate of each subject was calculated for the entire period of observation and for 56 five-minute bins. Histograms of the interblink time interval were plotted for each measurement of blink rate. RESULTS. Neither the overall mean blink rate (controls, 19.8 +/- 4.9; Graves`, 17.6 +/- 5.4) nor the interblink time (controls, 5.2 +/- 3.1, Graves`, 7.9 +/- 3.5) differed between the two groups. There was a large variation of both measurements when the 5-minute bins were considered. The interblink time distribution of all subjects was highly positively skewed when the 1-hour period was measured. A significant number of the 5-minute bin distributions deviated from the overall pattern and became symmetric. CONCLUSIONS. The normal blinking process is characterized by highly positively skewed interblink time distributions. This result means that most blinks have a short time interval, and occasionally a small number of blinks have long time intervals. The different patterns of distribution described in the early literature probably represent artifacts because of the small samples analyzed. (Invest Ophthalmol Vis Sci. 2011;52:3419-3424) DOI:10.1167/iovs.10-7060

Drosophila melanogaster lipins are tissue-regulated and developmentally regulated and present specific subcellular distributions

Lipins constitute a novel family of Mg2+-dependent phosphatidate phosphatases that catalyze the dephosphorylation of phosphatidic acid to yield diacylglycerol, an important intermediate in lipid metabolism and cell signaling. Whereas a single lipin is detected in less complex organisms, in mammals there are distinct lipin isoforms and paralogs that are differentially expressed among tissues. Compatible with organism tissue complexity, we show that the single Drosophila Lpin1 ortholog (CG8709, here named DmLpin) expresses at least three isoforms (DmLpinA, DmLpinK and DmLpinJ) in a temporal and spatially regulated manner. The highest levels of lipin in the fat body, where DmLpinA and DmLpinK are expressed, correlate with the highest levels of triacylglycerol (TAG) measured in this tissue. DmLpinK is the most abundant isoform in the central nervous system, where TAG levels are significantly lower than in the fat body. In the testis, where TAG levels are even lower, DmLpinJ is the predominant isoform. Together, these data suggest that DmLpinA might be the isoform that is mainly involved in TAG production, and that DmLpinK and DmLpinJ could perform other cellular functions. In addition, we demonstrate by immunofluorescence that lipins are most strongly labeled in the perinuclear region of the fat body and ventral ganglion cells. In visceral muscles of the larval midgut and adult testis, lipins present a sarcomeric distribution. In the ovary chamber, the lipin signal is concentrated in the internal rim of the ring canal. These specific subcellular localizations of the Drosophila lipins provide the basis for future investigations on putative novel cellular functions of this protein family.

The impact of smoke from forest fires on the spectral dispersion of cloud droplet size distributions in the Amazonian region

In this paper, the main microphysical characteristics of clouds developing in polluted and clean conditions in the biomass-burning season of the Amazon region are examined, with special attention to the spectral dispersion of the cloud droplet size distribution and its potential impact on climate modeling applications. The dispersion effect has been shown to alter the climate cooling predicted by the so-called Twomey effect. In biomass-burning polluted conditions, high concentrations of low dispersed cloud droplets are found. Clean conditions revealed an opposite situation. The liquid water content (0.43 +/- 0.19 g m(-3)) is shown to be uncorrelated with the cloud drop number concentration, while the effective radius is found to be very much correlated with the relative dispersion of the size distribution (R(2) = 0.81). The results suggest that an increase in cloud condensation nuclei concentration from biomass-burning aerosols may lead to an additional effect caused by a decrease in relative dispersion. Since the dry season in the Amazonian region is vapor limiting, the dispersion effect of cloud droplet size distributions could be substantially larger than in other polluted regions.

Dynamic picture of the inner asteroid belt: implications for the density, size and taxonomic distributions of real objects

In this paper, we construct a dynamic portrait of the inner asteroidal belt. We use information about the distribution of test particles, which were initially placed on a perfectly rectangular grid of initial conditions, after 4.2 Myr of gravitational interactions with the Sun and five planets, from Mars to Neptune. Using the spectral analysis method introduced by Michtchenko et al., the asteroidal behaviour is illustrated in detail on the dynamical, averaged and frequency maps. On the averaged and frequency maps, we superpose information on the proper elements and proper frequencies of real objects, extracted from the data base, AstDyS, constructed by Milani and Knezevic. A comparison of the maps with the distribution of real objects allows us to detect possible dynamical mechanisms acting in the domain under study; these mechanisms are related to mean-motion and secular resonances. We note that the two- and three-body mean-motion resonances and the secular resonances (strong linear and weaker non-linear) have an important role in the diffusive transportation of the objects. Their long-lasting action, overlaid with the Yarkovsky effect, may explain many observed features of the density, size and taxonomic distributions of the asteroids.

Caesium, potassium and ammonium distributions in different organs of tropical plants

In the present work the distribution of ions in aboveground plant parts was studied in order to establish the suitability of using radiocaesium as a tracer for the plant absorption of nutrients, such as potassium (K(+)) and ammonium (NH(4)(+)). We present the results for the distributions of (137)Cs, (40)K and NH(4)(+) from four tropical plant species: lemon (Citrus aurantifolia), orange (Citrus sinensis), guava (Psidium guajava) and chili pepper (Capsicum frutescens). Activity concentrations of (137)Cs and (40)K were measured by gamma spectrometry and concentrations of free NH(4)(+) ions by a colorimetric method. Similarly to potassium and ammonium, caesium showed a high mobility within the plants, exhibiting the highest values of concentration in the growing parts of the tree (fruits, new leaves, twigs, and barks). A significant correlation between activity concentrations of (137)Cs and (40)K was observed in these tropical plants. The K/Cs discrimination ratios were approximately equal to unity in different compartments of each individual plant, suggesting that caesium could be a good tracer for (40)K in tropical woody fruit species. Despite the similarity observed for the behaviour of caesium and ammonium in the newly grown plant compartments, (137)Cs was not well correlated with NH(4)(+). Significant temporal changes in the NH(4)(+) concentrations were observed during the development of fruits, while the (137)Cs activity concentration alterations were not of great importance, indicating, therefore, that Cs(+) and free NH(4)(+) ions could have distinct concentration ratios for each particular plant organ. (C) 2008 Elsevier B.V. All rights reserved.

INDEFINITE QUASILINEAR ELLIPTIC EQUATIONS IN EXTERIOR DOMAINS WITH EXPONENTIAL CRITICAL GROWTH

This paper is concerned with the existence of solutions for the quasilinear problem {-div(vertical bar del u vertical bar(N-2) del u) + vertical bar u vertical bar(N-2) u = a(x)g(u) in Omega u = 0 on partial derivative Omega, where Omega subset of R(N) (N >= 2) is an exterior domain; that is, Omega = R(N)\omega, where omega subset of R(N) is a bounded domain, the nonlinearity g(u) has an exponential critical growth at infinity and a(x) is a continuous function and changes sign in Omega. A variational method is applied to establish the existence of a nontrivial solution for the above problem.

A unified view on lifetime distributions arising from selection mechanisms

In this paper, we formulate a flexible density function from the selection mechanism viewpoint (see, for example, Bayarri and DeGroot (1992) and Arellano-Valle et al. (2006)) which possesses nice biological and physical interpretations. The new density function contains as special cases many models that have been proposed recently in the literature. In constructing this model, we assume that the number of competing causes of the event of interest has a general discrete distribution characterized by its probability generating function. This function has an important role in the selection procedure as well as in computing the conditional personal cure rate. Finally, we illustrate how various models can be deduced as special cases of the proposed model. (C) 2011 Elsevier B.V. All rights reserved.

A new family of generalized distributions

Kumaraswamy [Generalized probability density-function for double-bounded random-processes, J. Hydrol. 462 (1980), pp. 79-88] introduced a distribution for double-bounded random processes with hydrological applications. For the first time, based on this distribution, we describe a new family of generalized distributions (denoted with the prefix `Kw`) to extend the normal, Weibull, gamma, Gumbel, inverse Gaussian distributions, among several well-known distributions. Some special distributions in the new family such as the Kw-normal, Kw-Weibull, Kw-gamma, Kw-Gumbel and Kw-inverse Gaussian distribution are discussed. We express the ordinary moments of any Kw generalized distribution as linear functions of probability weighted moments (PWMs) of the parent distribution. We also obtain the ordinary moments of order statistics as functions of PWMs of the baseline distribution. We use the method of maximum likelihood to fit the distributions in the new class and illustrate the potentiality of the new model with an application to real data.

Bayesian nonlinear regression models with scale mixtures of skew-normal distributions: Estimation and case influence diagnostics

The purpose of this paper is to develop a Bayesian analysis for nonlinear regression models under scale mixtures of skew-normal distributions. This novel class of models provides a useful generalization of the symmetrical nonlinear regression models since the error distributions cover both skewness and heavy-tailed distributions such as the skew-t, skew-slash and the skew-contaminated normal distributions. The main advantage of these class of distributions is that they have a nice hierarchical representation that allows the implementation of Markov chain Monte Carlo (MCMC) methods to simulate samples from the joint posterior distribution. In order to examine the robust aspects of this flexible class, against outlying and influential observations, we present a Bayesian case deletion influence diagnostics based on the Kullback-Leibler divergence. Further, some discussions on the model selection criteria are given. The newly developed procedures are illustrated considering two simulations study, and a real data previously analyzed under normal and skew-normal nonlinear regression models. (C) 2010 Elsevier B.V. All rights reserved.

Robust linear mixed models with skew-normal independent distributions from a Bayesian perspective

Linear mixed models were developed to handle clustered data and have been a topic of increasing interest in statistics for the past 50 years. Generally. the normality (or symmetry) of the random effects is a common assumption in linear mixed models but it may, sometimes, be unrealistic, obscuring important features of among-subjects variation. In this article, we utilize skew-normal/independent distributions as a tool for robust modeling of linear mixed models under a Bayesian paradigm. The skew-normal/independent distributions is an attractive class of asymmetric heavy-tailed distributions that includes the skew-normal distribution, skew-t, skew-slash and the skew-contaminated normal distributions as special cases, providing an appealing robust alternative to the routine use of symmetric distributions in this type of models. The methods developed are illustrated using a real data set from Framingham cholesterol study. (C) 2009 Elsevier B.V. All rights reserved.