Generalized beta-generated distributions

This article introduces generalized beta-generated (GBG) distributions. Sub-models include all classical beta-generated, Kumaraswamy-generated and exponentiated distributions. They are maximum entropy distributions under three intuitive conditions, which show that the classical beta generator skewness parameters only control tail entropy and an additional shape parameter is needed to add entropy to the centre of the parent distribution. This parameter controls skewness without necessarily differentiating tail weights. The GBG class also has tractable properties: we present various expansions for moments, generating function and quantiles. The model parameters are estimated by maximum likelihood and the usefulness of the new class is illustrated by means of some real data sets. (c) 2011 Elsevier B.V. All rights reserved.

The log-exponentiated generalized gamma regression model for censored data

For the first time, we introduce a generalized form of the exponentiated generalized gamma distribution [Cordeiro et al. The exponentiated generalized gamma distribution with application to lifetime data, J. Statist. Comput. Simul. 81 (2011), pp. 827-842.] that is the baseline for the log-exponentiated generalized gamma regression model. The new distribution can accommodate increasing, decreasing, bathtub- and unimodal-shaped hazard functions. A second advantage is that it includes classical distributions reported in the lifetime literature as special cases. We obtain explicit expressions for the moments of the baseline distribution of the new regression model. The proposed model can be applied to censored data since it includes as sub-models several widely known regression models. It therefore can be used more effectively in the analysis of survival data. We obtain maximum likelihood estimates for the model parameters by considering censored data. We show that our extended regression model is very useful by means of two applications to real data.

A generalized finite element method with global-local enrichment functions for confined plasticity problems

The main feature of partition of unity methods such as the generalized or extended finite element method is their ability of utilizing a priori knowledge about the solution of a problem in the form of enrichment functions. However, analytical derivation of enrichment functions with good approximation properties is mostly limited to two-dimensional linear problems. This paper presents a procedure to numerically generate proper enrichment functions for three-dimensional problems with confined plasticity where plastic evolution is gradual. This procedure involves the solution of boundary value problems around local regions exhibiting nonlinear behavior and the enrichment of the global solution space with the local solutions through the partition of unity method framework. This approach can produce accurate nonlinear solutions with a reduced computational cost compared to standard finite element methods since computationally intensive nonlinear iterations can be performed on coarse global meshes after the creation of enrichment functions properly describing localized nonlinear behavior. Several three-dimensional nonlinear problems based on the rate-independent J (2) plasticity theory with isotropic hardening are solved using the proposed procedure to demonstrate its robustness, accuracy and computational efficiency.

Generalized Metropolis dynamics with a generalized master equation: An approach for time-independent and time-dependent Monte Carlo simulations of generalized spin systems

The extension of Boltzmann-Gibbs thermostatistics, proposed by Tsallis, introduces an additional parameter q to the inverse temperature beta. Here, we show that a previously introduced generalized Metropolis dynamics to evolve spin models is not local and does not obey the detailed energy balance. In this dynamics, locality is only retrieved for q = 1, which corresponds to the standard Metropolis algorithm. Nonlocality implies very time-consuming computer calculations, since the energy of the whole system must be reevaluated when a single spin is flipped. To circumvent this costly calculation, we propose a generalized master equation, which gives rise to a local generalized Metropolis dynamics that obeys the detailed energy balance. To compare the different critical values obtained with other generalized dynamics, we perform Monte Carlo simulations in equilibrium for the Ising model. By using short-time nonequilibrium numerical simulations, we also calculate for this model the critical temperature and the static and dynamical critical exponents as functions of q. Even for q not equal 1, we show that suitable time-evolving power laws can be found for each initial condition. Our numerical experiments corroborate the literature results when we use nonlocal dynamics, showing that short-time parameter determination works also in this case. However, the dynamics governed by the new master equation leads to different results for critical temperatures and also the critical exponents affecting universality classes. We further propose a simple algorithm to optimize modeling the time evolution with a power law, considering in a log-log plot two successive refinements.

Generalized non-commutative inflation

Non-commutative geometry indicates a deformation of the energy-momentum dispersion relation f (E) = E/pc (not equal 1) for massless particles. This distorted energy-momentum relation can affect the radiation-dominated phase of the universe at sufficiently high temperature. This prompted the idea of non-commutative inflation by Alexander et al (2003 Phys. Rev. D 67 081301) and Koh and Brandenberger (2007 JCAP06(2007) 021 and JCAP11(2007) 013). These authors studied a one-parameter family of a non-relativistic dispersion relation that leads to inflation: the a family of curves f (E) = 1 + (lambda E)(alpha). We show here how the conceptually different structure of symmetries of non-commutative spaces can lead, in a mathematically consistent way, to the fundamental equations of non-commutative inflation driven by radiation. We describe how this structure can be considered independently of (but including) the idea of non-commutative spaces as a starting point of the general inflationary deformation of SL(2, C). We analyze the conditions on the dispersion relation that leads to inflation as a set of inequalities which plays the same role as the slow-roll conditions on the potential of a scalar field. We study conditions for a possible numerical approach to obtain a general one-parameter family of dispersion relations that lead to successful inflation.

ON THE CONSISTENCY OF TWISTED GENERALIZED WEYL ALGEBRAS

A twisted generalized Weyl algebra A of degree n depends on a. base algebra R, n commuting automorphisms sigma(i) of R, n central elements t(i) of R and on some additional scalar parameters. In a paper by Mazorchuk and Turowska, it is claimed that certain consistency conditions for sigma(i) and t(i) are sufficient for the algebra to be nontrivial. However, in this paper we give all example which shows that this is false. We also correct the statement by finding a new set of consistency conditions and prove that the old and new conditions together are necessary and sufficient for the base algebra R to map injectively into A. In particular they are sufficient for the algebra A to be nontrivial. We speculate that these consistency relations may play a role in other areas of mathematics, analogous to the role played by the Yang-Baxter equation in the theory of integrable systems.

Variance Effective Population Size for Dioecious Species

The concept of effective population size (N(e)) is an important measure of representativeness in many areas. In this research, we consider the statistical properties of the number of contributed gametes under practical situations by adapting Crow and Denninston's (1988) N(e) formulas for dioecious species. Three sampling procedures were considered. In all circumstances, results show that as the offspring sex ratio (r) deviates from 0.5, N(e) values become smaller, and the efficiency of gametic control for increasing N(e) is reduced. For finite populations, where all individuals are potentially functional parents, the reduction in N(e) due to an unequal sex ratio can be compensated for through female gametic control when 0.28 <= r <= 0.72. This outcome is important when r is unknown. When only a fraction of the individuals in a population is taken for reproduction, N(e) is meaningful only if the size of the reference population is clearly defined. Gametic control is a compensating factor in accession regeneration when the viability of the accession is around 70 or 75%. For germ-plasm collection, when parents are a very small fraction of the population, maximum N(e) will be approximately 47 and 57% of the total number of offspring sampled, with female gametic control, r varying between 0.3 and 0.5, and being constant over generations.

GENERALIZED FINITE ELEMENT METHOD ON NONCONVENTIONAL HYBRID-MIXED FORMULATION

The generalized finite element method (GFEM) is applied to a nonconventional hybrid-mixed stress formulation (HMSF) for plane analysis. In the HMSF, three approximation fields are involved: stresses and displacements in the domain and displacement fields on the static boundary. The GFEM-HMSF shape functions are then generated by the product of a partition of unity associated to each field and the polynomials enrichment functions. In principle, the enrichment can be conducted independently over each of the HMSF approximation fields. However, stability and convergence features of the resulting numerical method can be affected mainly by spurious modes generated when enrichment is arbitrarily applied to the displacement fields. With the aim to efficiently explore the enrichment possibilities, an extension to GFEM-HMSF of the conventional Zienkiewicz-Patch-Test is proposed as a necessary condition to ensure numerical stability. Finally, once the extended Patch-Test is satisfied, some numerical analyses focusing on the selective enrichment over distorted meshes formed by bilinear quadrilateral finite elements are presented, thus showing the performance of the GFEM-HMSF combination.

A Single Statistic for Monitoring the Covariance Matrix of Bivariate Processes

In this article, we present a new control chart for monitoring the covariance matrix in a bivariate process. In this method, n observations of the two variables were considered as if they came from a single variable (as a sample of 2n observations), and a sample variance was calculated. This statistic was used to build a new control chart specifically as a VMIX chart. The performance of the new control chart was compared with its main competitors: the generalized sampled variance chart, the likelihood ratio test, Nagao's test, probability integral transformation (v(t)), and the recently proposed VMAX chart. Among these statistics, only the VMAX chart was competitive with the VMIX chart. For shifts in both variances, the VMIX chart outperformed VMAX; however, VMAX showed better performance for large shifts (higher than 10%) in one variance.

Optimal mean-variance control for discrete-time linear systems with Markovian jumps and multiplicative noises

In this paper, we consider the stochastic optimal control problem of discrete-time linear systems subject to Markov jumps and multiplicative noises under two criteria. The first one is an unconstrained mean-variance trade-off performance criterion along the time, and the second one is a minimum variance criterion along the time with constraints on the expected output. We present explicit conditions for the existence of an optimal control strategy for the problems, generalizing previous results in the literature. We conclude the paper by presenting a numerical example of a multi-period portfolio selection problem with regime switching in which it is desired to minimize the sum of the variances of the portfolio along the time under the restriction of keeping the expected value of the portfolio greater than some minimum values specified by the investor. (C) 2011 Elsevier Ltd. All rights reserved.

Bilateral subthalamic nucleus stimulation for generalized dystonia after bilateral pallidotomy

Background: Thalamotomies and pallidotomies were commonly performed before the deep brain stimulation (DBS) era. Although ablative procedures can lead to significant dystonia improvement, longer periods of analysis reveal disease progression and functional deterioration. Today, the same patients seek additional treatment possibilities. Methods: Four patients with generalized dystonia who previously had undergone bilateral pallidotomy came to our service seeking additional treatment because of dystonic symptom progression. Bilateral subthalamic nucleus DBS (B-STN-DBS) was the treatment of choice. The patients were evaluated with the BurkeFahnMarsden Dystonia Rating Scale (BFMDRS) and the Unified Dystonia Rating Scale (UDRS) before and 2 years after surgery. Results: All patients showed significant functional improvement, averaging 65.3% in BFMDRS (P = .014) and 69.2% in UDRS (P = .025). Conclusions: These results suggest that B-STN-DBS may be an interesting treatment option for generalized dystonia, even for patients who have already undergone bilateral pallidotomy. (c) 2012 Movement Disorder Society

Longitudinal familial analysis of blood pressure involving parametric (co)variance functions

Abstract Background For analyzing longitudinal familial data we adopted a log-linear form to incorporate heterogeneity in genetic variance components over the time, and additionally a serial correlation term in the genetic effects at different levels of ages. Due to the availability of multiple measures on the same individual, we permitted environmental correlations that may change across time. Results Systolic blood pressure from family members from the first and second cohort was used in the current analysis. Measures of subjects receiving hypertension treatment were set as censored values and they were corrected. An initial check of the variance and covariance functions proposed for analyzing longitudinal familial data, using empirical semi-variogram plots, indicated that the observed trait dispersion pattern follows the assumptions adopted. Conclusion The corrections for censored phenotypes based on ordinary linear models may be an appropriate simple model to correct the data, ensuring that the original variability in the data was retained. In addition, empirical semi-variogram plots are useful for diagnosis of the (co)variance model adopted.

Catch rates of artisanal fisheries targeting Dasyatis guttata (Bloch & Schneider, 1801) in Todos os Santos Bay, northeastern Brazil

Dasyatis guttata has been target of artisanal fisheries in the coast of Bahia (Northeast Brazil) mainly by “arraieira” (gillnet) and “grozeira” (bottom long-line), but until now there is no stock assessment study. One of the important data for this knowledge is reliable indices of abundance. The aims of the present work are to: (1) estimate the best predictor for relative abundance (catch-per-unit-of-effort, CPUE), examining whether catch (production – kg) was related to: soak time of the gear, size of the gillnet or number of hooks, applying generalized linear model (GLM); (2) estimate the annual CPUE (kg/hooks and kg/m) averaged by gear; and (3) assess the temporal CPUE variance. Based on monthly sampling between January 2012 and January 2013, 222 landings by grozeira and 76 by arraiaiera were recorded in the two landing sites in Todos os Santos Bay, Bahia. A total of 14,550 kg (average = 44 kg/month) of D. guttata was captured. Models for both gears were highly significant (P < 0.0001). The analysis indicated that the most appropriate variable for CPUE analysis was the size of the gillnet (P < 0.001) and the number of hooks (P < 0.0001). Soak time of the gear was not significant for both gears (P = 0.4). High residual deviance expresses the complexity of the relations between ecosystem factors and other fisheries factors affecting relative abundance, which were not considered in this study. The average CPUE by grozeira was 6.39 kg/100 hooks ± 8.89 and by arraieira, 1.47 kg/100 m ± 1.66 over the year. Kruskal-Wallis test showed effect of the month on the mean grozeira CPUE (P = <0.001), but no effect (P = 0.096) on the mean arraieira CPUE. Grozeira CPUE values were highest in December and March, and lowest between May to August