Assessment of variance components in nonlinear mixed-effects elliptical models

The issue of assessing variance components is essential in deciding on the inclusion of random effects in the context of mixed models. In this work we discuss this problem by supposing nonlinear elliptical models for correlated data by using the score-type test proposed in Silvapulle and Silvapulle (1995). Being asymptotically equivalent to the likelihood ratio test and only requiring the estimation under the null hypothesis, this test provides a fairly easy computable alternative for assessing one-sided hypotheses in the context of the marginal model. Taking into account the possible non-normal distribution, we assume that the joint distribution of the response variable and the random effects lies in the elliptical class, which includes light-tailed and heavy-tailed distributions such as Student-t, power exponential, logistic, generalized Student-t, generalized logistic, contaminated normal, and the normal itself, among others. We compare the sensitivity of the score-type test under normal, Student-t and power exponential models for the kinetics data set discussed in Vonesh and Carter (1992) and fitted using the model presented in Russo et al. (2009). Also, a simulation study is performed to analyze the consequences of the kurtosis misspecification.

On Wald Residuals in Generalized Linear Models

A rigorous asymptotic theory for Wald residuals in generalized linear models is not yet available. The authors provide matrix formulae of order O(n(-1)), where n is the sample size, for the first two moments of these residuals. The formulae can be applied to many regression models widely used in practice. The authors suggest adjusted Wald residuals to these models with approximately zero mean and unit variance. The expressions were used to analyze a real dataset. Some simulation results indicate that the adjusted Wald residuals are better approximated by the standard normal distribution than the Wald residuals.

Growth curve in Santa Ines sheep

The objective of this research was to use non-linear models to describe the growth pattern in Santa Ines sheep and to study the influence of environmental effects on curve parameters with the best-fit model. The models included the Brody, Richards, Von Bertalanffy, Gompertz, and Logistic models. We used 773 field reports on 162 animals ranging in age from 120 to 774 days, including 46 males and 116 females. The statistics used to evaluate the quality of fit included RMS (residual mean square), C% (percentage of convergence), R-2 (adjusted determination coefficient) and MAD (mean absolute deviation). Of the fixed effects studied, the only significant relationship was the effect of sex on parameter A. The Richards model was problematic during the process of convergence. Considering all studied criteria, the Logistic model presented the best fit in describing the growth pattern in Santa Ines sheep. (C) 2011 Elsevier B.V. All rights reserved.

Dynamics of adventitious rooting in mini-cuttings of Eucalyptus benthamii x Eucalyptus dunnii

It is possible to determine the optimum time for permanence of vegetative propagules (mini-cuttings) inside a greenhouse for rooting, and this value can be used to optimize the structure of the nursery. The aim of this study was to determine the dynamics of adventitious rooting in mini-cuttings of three clones of Eucalyptus benthamii x Eucalyptus dunnii. Sprouts of H12, H19 and H20 clones were collected from mini-stumps that were planted in gutters containing sand and grown in a semi-hydroponic system. The basal region of the mini-cuttings was immersed in 2,000 mg L-1 indole-3-butyric acid (IBA) solution for 10 seconds. The rooting percentage of the mini-cuttings, the total length of the root system and the rooting rate per mini-cutting were also evaluated at 0 (time of planting), 7, 14, 21, 28, 35, 42, 49 and 56 days. We used logistic and exponential regression to mathematically model the speed of rhizogenesis. The rooting percentage was best represented as a logistic model, and the total length of the root system was best represented as an exponential model. The clones had different speeds of adventitious rooting. The optimum time for permanence of the mini-cuttings inside the greenhouse for rooting was between 35 and 42 days, and varied depending on the genetic material.

Risk factors for hepatitis C virus infection in Inland Brazil: An analysis of pooled epidemiological sectional studies

In order to assess the contribution of different parenteral routes as risk exposure to the hepatitis C virus (HCV), samples from nine surveys or cross-sectional studies conducted in two Brazilian inland regions were pooled, including a total of 3,910 subjects. Heterogeneity among the study results for different risk factors was tested and the results were shown to be homogeneous. Anti-HCV antibodies were observed in 241 individuals, of which 146 (3.7%, 95% CI?=?3.24.4) had HCV exposure confirmed by immunoblot analysis or PCR test. After adjustment for relevant variables, a correlation between confirmed HCV exposure and injection drug use, tattooing, and advance age was observed. In a second logistic model that included exposures not searched in all nine studies, a smaller sample was analyzed, revealing an independent HCV association with past history of surgery and males who have sex with other males, in addition to repeated injection drug use. Overall, these analyses corroborate the finding that injection drug use is the main risk factor for HCV exposure and spread, in addition to other parenteral routes. J. Med. Virol. 84:756762, 2012. (C) 2012 Wiley Periodicals, Inc.

Spatial clusters of suicide in the municipality of Sao Paulo 1996-2005: an ecological study

Background: In a classical study, Durkheim mapped suicide rates, wealth, and low family density and realized that they clustered in northern France. Assessing others variables, such as religious society, he constructed a framework for the analysis of the suicide, which still allows international comparisons using the same basic methodology. The present study aims to identify possible significantly clusters of suicide in the city of Sao Paulo, and then, verify their statistical associations with socio-economic and cultural characteristics. Methods: A spatial scan statistical test was performed to analyze the geographical pattern of suicide deaths of residents in the city of Sao Paulo by Administrative District, from 1996 to 2005. Relative risks and high and/or low clusters were calculated accounting for gender and age as co-variates, were analyzed using spatial scan statistics to identify geographical patterns. Logistic regression was used to estimate associations with socioeconomic variables, considering, the spatial cluster of high suicide rates as the response variable. Drawing from Durkheim's original work, current World Health Organization (WHO) reports and recent reviews, the following independent variables were considered: marital status, income, education, religion, and migration. Results: The mean suicide rate was 4.1/100,000 inhabitant-years. Against this baseline, two clusters were identified: the first, of increased risk (RR = 1.66), comprising 18 districts in the central region; the second, of decreased risk (RR = 0.78), including 14 districts in the southern region. The downtown area toward the southwestern region of the city displayed the highest risk for suicide, and though the overall risk may be considered low, the rate climbs up to an intermediate level in this region. One logistic regression analysis contrasted the risk cluster (18 districts) against the other remaining 78 districts, testing the effects of socioeconomic-cultural variables. The following categories of proportion of persons within the clusters were identified as risk factors: singles (OR = 2.36), migrants (OR = 1.50), Catholics (OR = 1.37) and higher income (OR = 1.06). In a second logistic model, likewise conceived, the following categories of proportion of persons were identified as protective factors: married (OR = 0.49) and Evangelical (OR = 0.60). Conclusions: This risk/ protection profile is in accordance with the interpretation that, as a social phenomenon, suicide is related to social isolation. Thus, the classical framework put forward by Durkheim seems to still hold, even though its categorical expression requires re-interpretation.

Awareness and knowledge of HPV, cervical cancer, and vaccines in young women after first delivery in São Paulo, Brazil - a cross-sectional study

Abstract Background The success of HPV vaccination programs will require awareness regarding HPV associated diseases and the benefits of HPV vaccination for the general population. The aim of this study was to assess the level of awareness and knowledge of human papillomavirus (HPV) infection, cervical cancer prevention, vaccines, and factors associated with HPV awareness among young women after birth of the first child. Methods This analysis is part of a cross-sectional study carried out at Hospital Maternidade Leonor Mendes de Barros, a large public maternity hospital in Sao Paulo. Primiparous women (15-24 years) who gave birth in that maternity hospital were included. A questionnaire that included questions concerning knowledge of HPV, cervical cancer, and vaccines was applied. To estimate the association of HPV awareness with selected factors, prevalence ratios (PR) were estimated using a generalized linear model (GLM). Results Three hundred and one primiparous women were included; 37% of them reported that they "had ever heard about HPV", but only 19% and 7%, respectively, knew that HPV is a sexually transmitted infection (STI) and that it can cause cervical cancer. Seventy-four percent of interviewees mentioned the preventive character of vaccines and all participants affirmed that they would accept HPV vaccination after delivery. In the multivariate analysis, only increasing age (P for trend = 0.021) and previous STI (P < 0.001) were factors independently associated with HPV awareness ("had ever heard about HPV"). Conclusions This survey indicated that knowledge about the association between HPV and cervical cancer among primiparous young women is low. Therefore, these young low-income primiparous women could benefit greatly from educational interventions to encourage primary and secondary cervical cancer prevention programs.

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

Thermodynamics of Ising model with infinite-range interactions by generalized canonical ensemble

In this work we present the idea of how generalized ensembles can be used to simplify the operational study of non-additive physical systems. As alternative of the usual methods of direct integration or mean-field theory, we show how the solution of the Ising model with infinite-range interactions is obtained by using a generalized canonical ensemble. We describe how the thermodynamical properties of this model in the presence of an external magnetic field are founded by simple parametric equations. Without impairing the usual interpretation, we obtain an identical critical behaviour as observed in traditional approaches.

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 general threshold stress hybrid hazard model for lifetime data

In this paper we propose a hybrid hazard regression model with threshold stress which includes the proportional hazards and the accelerated failure time models as particular cases. To express the behavior of lifetimes the generalized-gamma distribution is assumed and an inverse power law model with a threshold stress is considered. For parameter estimation we develop a sampling-based posterior inference procedure based on Markov Chain Monte Carlo techniques. We assume proper but vague priors for the parameters of interest. A simulation study investigates the frequentist properties of the proposed estimators obtained under the assumption of vague priors. Further, some discussions on model selection criteria are given. The methodology is illustrated on simulated and real lifetime data set.

Connections between the Sznajd model with general confidence rules and graph theory

The Sznajd model is a sociophysics model that is used to model opinion propagation and consensus formation in societies. Its main feature is that its rules favor bigger groups of agreeing people. In a previous work, we generalized the bounded confidence rule in order to model biases and prejudices in discrete opinion models. In that work, we applied this modification to the Sznajd model and presented some preliminary results. The present work extends what we did in that paper. We present results linking many of the properties of the mean-field fixed points, with only a few qualitative aspects of the confidence rule (the biases and prejudices modeled), finding an interesting connection with graph theory problems. More precisely, we link the existence of fixed points with the notion of strongly connected graphs and the stability of fixed points with the problem of finding the maximal independent sets of a graph. We state these results and present comparisons between the mean field and simulations in Barabasi-Albert networks, followed by the main mathematical ideas and appendices with the rigorous proofs of our claims and some graph theory concepts, together with examples. We also show that there is no qualitative difference in the mean-field results if we require that a group of size q > 2, instead of a pair, of agreeing agents be formed before they attempt to convince other sites (for the mean field, this would coincide with the q-voter model).

Effective carrying capacity and analytical solution of a particular case of the Richards-like two-species population dynamics model

We consider a generalized two-species population dynamic model and analytically solve it for the amensalism and commensalism ecological interactions. These two-species models can be simplified to a one-species model with a time dependent extrinsic growth factor. With a one-species model with an effective carrying capacity one is able to retrieve the steady state solutions of the previous one-species model. The equivalence obtained between the effective carrying capacity and the extrinsic growth factor is complete only for a particular case, the Gompertz model. Here we unveil important aspects of sigmoid growth curves, which are relevant to growth processes and population dynamics. (C) 2011 Elsevier B.V. All rights reserved.

Generalized correlation functions for conductance fluctuations and the mesoscopic spin Hall effect

We study the spin Hall conductance fluctuations in ballistic mesoscopic systems. We obtain universal expressions for the spin and charge current fluctuations, cast in terms of current-current autocorrelation functions. We show that the latter are conveniently parametrized as deformed Lorentzian shape lines, functions of an external applied magnetic field and the Fermi energy. We find that the charge current fluctuations show quite unique statistical features at the symplectic-unitary crossover regime. Our findings are based on an evaluation of the generalized transmission coefficients correlation functions within the stub model and are amenable to experimental test. DOI: 10.1103/PhysRevB.86.235112

Boundedness of solutions of retarded functional differential equations with variable impulses via generalized ordinary differential equations

In this paper, we give sufficient conditions for the uniform boundedness and uniform ultimate boundedness of solutions of a class of retarded functional differential equations with impulse effects acting on variable times. We employ the theory of generalized ordinary differential equations to obtain our results. As an example, we investigate the boundedness of the solution of a circulating fuel nuclear reactor model.