Relatórios de Pesquisa

11/2014 Robust Mixture Regression Modeling Based on Scale Mixtures of Skew-Normal Distributions
Camila Borelli Zeller, Celso R. B. Cabral, Víctor H. Lachos

The traditional estimation of mixture regression models is based on the assumption of normality (symmetry) of component errors and thus is sensitive to outliers, heavy-tailed errors and/or asymmetric errors. In this work we present a proposal to deal with these issues simultaneously in the context of the mixture regression by extending the classic normal model byassuming that the random errors follow a scale mixtures of skew-normal distributions. This approach allows us to model data with great flexibility, accommodating skewness and heavy tails. The main virtue of considering the mixture regression models under the class of scale mixtures of skew-normal distributions is that they have a nice hierarchical representation whichallows easy implementation of inference. We develop a simple EM-type algorithm to perform maximum likelihood inference of the parameters of the proposed model. In order to examine the robust aspect of this flexible model against outlying observations, some simulation studies are also been presented. Finally, a real data set is analyzed, illustrating the usefulness of the proposed method.

PDF icon rp-2014-11.pdf
10/2014 Bifurcations of mutually coupled equations in random graphs
Eduardo Garibaldi, Tiago Pereira

We study the behavior of solutions of mutually coupled equations in heterogeneous random graphs. Heterogeneity means that some equations receive many inputs whereas most of the equations are given only with a few connections. Starting from a situation where the isolated equations are unstable, we prove that the a heterogeneous interaction structure leads to the appearance of stable subspaces of solutions. Moreover, we show that, for certain classes of heterogeneous networks, increasing the strength of interaction leads to a cascade of bifurcations in which the dimension of the stable subspace of solu-tions increases. We explicitly determine the bifurcation scenario in terms of the graph structure.

PDF icon rp-2014-10.pdf
9/2014 A Note on Certain Maximal Curves
Saeed Tafazolian, Fernando Torres

We characterize certain maximal curves over finite fields whose plane models are of Hurwitz type, namely x m y a + y n + x b = 0. We also consider maximal hyperelliptic curves of maximal genus. Finally, we discuss maximal curves of type y q + y = x m via class field theory.

PDF icon rp-2014-9.pdf
8/2014 Lyapunov Graphs for Circle Valued Morse Functions
Ketty A. de Rezende, Guido G. E. Ledesma, Oziride Manzoli Neto, Gioia M. Vago

Within the context of Novikov theory, the Conley index is to used to obtain results for circular Morse flows on compact n-manifold. Examples are provided for orientable and non-orientable surfaces via a complete characterization of circular Morse digraphs. A generalization of these results is presented for smooth flows associated to circle valued Lyapunov functions.

PDF icon rp-2014-8.pdf
7/2014 Augmented mixed models for clustered proportion data
Dipankar Bandyopadhyay, Diana M. Galvis, Víctor H. Lachos

Often in biomedical research, we deal with continuous (clustered) proportion responses ranging between zero and one quantifying the disease status of the cluster units. Interestingly, the study population might also consist of relatively disease-free as well as highly diseased subjects, contributing to proportion values in the interval [0, 1]. Regression on a variety of parametric densities with support lying in (0, 1), such as beta regression, can assess important covariate effects. However, they are deemed inappropriate due the presence of zeros and/or ones. To evade this, we introduce a class of general proportion density (GPD), and further augment the probabilities of zero and one to this GPD, controlling for the clustering. Our approach is Bayesian, and presents a computationally convenient framework amenable to available freeware.Bayesian case-deletion influence diagnostics based on q-divergence measures are automatic from the MCMC output. The methodology is illustrated using both simulation studies and application to a real dataset from a clinical periodontology study.

PDF icon rp-2014-7.pdf
6/2014 Censored Linear Regression Models for Irregularly Observed Longitudinal Data using the Multivariate-t Distribution
Aldo M. Garay, Luis M. Castro, Jacek Leskow, Víctor H. Lachos

In AIDS studies it is quite common to observe viral load measurements collected irregularlyover time. Moreover, these measurements can be subjected to some upper and/or lower detection limitsdepending on the quantification assays. A complication arises when these continuous repeated measureshave a heavy-tailed behavior. For such data structures, we propose a robust structure for a censoredlinear model based on the multivariate Student-t distribution. To compensate for the autocorrelationexisting among irregularly observed measures, a damped exponential correlation structure is employed.An efficient EM-type algorithm is developed for computing the maximum likelihood estimates, obtaining as a by-product the standard errors of the fixed effects and the log-likelihood function. The proposed algorithm uses closed-form expressions at the E-step, that rely on formulas for the mean and variance of a truncated multivariate Student-t distribution. The methodology is illustrated through an application to an HIV-AIDS study and several simulation studies.

PDF icon rp-2014-6.pdf
5/2014 Further examples of maximal curves which cannot be covered by the Hermitian curve
Saeed Tafazolian, Arnoldo Teheran-Herrera, Fernando Torres

We construct examples of curves defined over the finite field Fq6 which are covered by the GK-curve. Thus such curves are maximal over Fq6 although they cannot be covered by the Hermitian curve for q > 2. We also give examples of maximal curvesthat cannot be Galois covered by the Hermitian curve over the finite field Fq2n with n > 3 odd and q > 2. We point out some applications to codes related to an array coming from telescopic semigroups.

PDF icon rp-2014-5.pdf
4/2014 Improving Shewhart-type Generalized Variance Control Charts for Multivariate Process Variability Monitoring using Cornish-Fisher Quantile Correction, Meijer-G Function and Other Tools
Emanuel P. Barbosa, Mario A. Gneri, Ariane Meneguetti

This paper presents an improved version of the Shewhart-type generalized variance |S| control chart for multivariate Gaussian process dispersion monitoring, based on the Cornish-Fisher quantile formula for non-normality correction of the traditional nor-mal based 3-sigma chart limits.Also, the exact sample distribution of |S| and its quantiles (chart exact limits) are obtained through the Meijer-G function (inverse Mellin-Barnes integral transform), and an auxiliary control chart based on the trace of the standardized S matrix is introduced in order to avoid non detection of certain changes in the process variance-covariance Σ matrix.The performance of the proposed CF-corrected control chart is compared, considering false alarm risk (using analytical and simulation tools), with the traditional normal based chart and with the exact distributed based chart (for dimensions d = 2 and d = 3). This study shows that the proposed control limit corrections do remove the drawback of excess of false alarm associated with the traditional normal based |S| control chart.The proposed new chart (with its corresponding auxiliary chart) is illustrated with two numerical examples.

PDF icon rp-2014-4.pdf
3/2014 On the curve Y n = X m + X over finite fields
Saeed Tafazolian, Fernando Torres

We show that a maximal curve over Fq2 defined by the affine equation y n = f (x), where f (x) ∈ Fq2 [x] has degree coprime to n, is such that n is a divisor of q + 1 if and only if f (x) has a root in Fq2 . In this case, all the roots of f (x) belong to Fq2 ;cf. Thm. 1.2, Thm. 4.3 in [J. Pure Appl. Algebra 212 (2008), 2513–2521]. In particular, we characterize certain maximal curves defined by equations of type y n = xm + x over finite fields.

PDF icon rp-2014-3.pdf
2/2014 A Mixed-Effect Model for Positive Responses Augmented by Zeros
Mariana R. Motta, Diana M. Galvis, Víctor H. Lachos, Filidor E. Vilca-Labra, Valéria Troncoso Baltar, Eliseu Verly Junior, Regina Mara Fisberg, Dirce Maria Lobo Marchioni

In this work we propose a model for positive and zero responses by means of a zero augmented mixed regression model. Under this class, we are particularly interested in studying positive responses whose distribution accommodates skewness. At the same time, responses can be zero and therefore we justify the use of a zero-augmented mixture model. We model the mean of the positive response in a logarithm scale and the mixture probability in a logit scale, both as afunction of fixed and random effects. Here, the random effects link the two random components through their joint distribution and incorporate within subject correlation due to repeated measurements and between-subject heterogeneity.An MCMC algorithm is tailored to obtain Bayesian posterior distributions of the unknown quantities of interest and Bayesian case-deletion influence diagnostics based on the q-divergence measure is performed. We motivate and illustrate the proposed methodology by means of a data set from a 24 hours dietary recall study obtained in the city of S ̃o Paulo, Brazil, and present a simulation study a to evaluate the performance of the proposed methods.Bayesian inference, gamma distribution, log-normal distribution, mixed models, random effects, usual intake, zero-augmented data.

PDF icon rp-2014-2.pdf