Aspectes metodològics i aplicacions de la modelització del temps de supervivència multivariant mitjançant models mixtes

Els estudis de supervivència s'interessen pel temps que passa des de l'inici de l'estudi (diagnòstic de la malaltia, inici del tractament,...) fins que es produeix l'esdeveniment d'interès (mort, curació, millora,...). No obstant això, moltes vegades aquest esdeveniment s'observa més d'una vegada en un mateix individu durant el període de seguiment (dades de supervivència multivariant). En aquest cas, és necessari utilitzar una metodologia diferent a la utilitzada en l'anàlisi de supervivència estàndard. El principal problema que l'estudi d'aquest tipus de dades comporta és que les observacions poden no ser independents. Fins ara, aquest problema s'ha solucionat de dues maneres diferents en funció de la variable dependent. Si aquesta variable segueix una distribució de la família exponencial s'utilitzen els models lineals generalitzats mixtes (GLMM); i si aquesta variable és el temps, variable amb una distribució de probabilitat no pertanyent a aquesta família, s'utilitza l'anàlisi de supervivència multivariant. El que es pretén en aquesta tesis és unificar aquests dos enfocs, és a dir, utilitzar una variable dependent que sigui el temps amb agrupacions d'individus o d'observacions, a partir d'un GLMM, amb la finalitat d'introduir nous mètodes pel tractament d'aquest tipus de dades.

Markov models of aging: theory and practice

We review and structure some of the mathematical and statistical models that have been developed over the past half century to grapple with theoretical and experimental questions about the stochastic development of aging over the life course. We suggest that the mathematical models are in large part addressing the problem of partitioning the randomness in aging: How does aging vary between individuals, and within an individual over the lifecourse? How much of the variation is inherently related to some qualities of the individual, and how much is entirely random? How much of the randomness is cumulative, and how much is merely short-term flutter? We propose that recent lines of statistical inquiry in survival analysis could usefully grapple with these questions, all the more so if they were more explicitly linked to the relevant mathematical and biological models of aging. To this end, we describe points of contact among the various lines of mathematical and statistical research. We suggest some directions for future work, including the exploration of information-theoretic measures for evaluating components of stochastic models as the basis for analyzing experiments and anchoring theoretical discussions of aging.

ESTIMATING SENSORY SHELF LIFE OF CHOCOLATE AND CARROT CUPCAKES USING ACCEPTANCE TESTS

This study determined the sensory shelf life of a commercial brand of chocolate and carrot cupcakes, aiming at increasing the current 120 days of shelf life to 180. Appearance, texture, flavor and overall quality of cakes stored at six different storage times were evaluated by 102 consumers. The data were analyzed by analysis of variance and linear regression. For both flavors, the texture presented a greater loss in acceptance during the storage period, showing an acceptance mean close to indifference on the hedonic scale at 120 days. Nevertheless, appearance, flavor and overall quality stayed acceptable up to 150 days. The end of shelf life was estimated at about 161 days for chocolate cakes and 150 days for carrot cakes. This study showed that the current 120 days of shelf life can be extended to 150 days for carrot cake and to 160 days for chocolate cake. However, the 180 days of shelf life desired by the company were not achieved. PRACTICAL APPLICATIONS This research shows the adequacy of using sensory acceptance tests to determine the shelf life of two food products (chocolate and carrot cupcakes). This practical application is useful because the precise determination of the shelf life of a food product is of vital importance for its commercial success. The maximum storage time should always be evaluated in the development or reformulation of new products, changes in packing or storage conditions. Once the physical-chemical and microbiological stability of a product is guaranteed, sensorial changes that could affect consumer acceptance will determine the end of the shelf life of a food product. Thus, the use of sensitive and reliable methods to estimate the sensory shelf life of a product is very important. Findings show the importance of determining the shelf life of each product separately and to avoid using the shelf time estimated for a specific product on other, similar products.

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.

The complementary exponential geometric distribution: Model, properties, and a comparison with its counterpart

In this paper, we proposed a new two-parameter lifetime distribution with increasing failure rate, the complementary exponential geometric distribution, which is complementary to the exponential geometric model proposed by Adamidis and Loukas (1998). The new distribution arises on a latent complementary risks scenario, in which the lifetime associated with a particular risk is not observable; rather, we observe only the maximum lifetime value among all risks. The properties of the proposed distribution are discussed, including a formal proof of its probability density function and explicit algebraic formulas for its reliability and failure rate functions, moments, including the mean and variance, variation coefficient, and modal value. The parameter estimation is based on the usual maximum likelihood approach. We report the results of a misspecification simulation study performed in order to assess the extent of misspecification errors when testing the exponential geometric distribution against our complementary one in the presence of different sample size and censoring percentage. The methodology is illustrated on four real datasets; we also make a comparison between both modeling approaches. (C) 2011 Elsevier B.V. All rights reserved.

Log-Burr XII regression models with censored data

In survival analysis applications, the failure rate function may frequently present a unimodal shape. In such case, the log-normal or log-logistic distributions are used. In this paper, we shall be concerned only with parametric forms, so a location-scale regression model based on the Burr XII distribution is proposed for modeling data with a unimodal failure rate function as an alternative to the log-logistic regression model. Assuming censored data, we consider a classic analysis, a Bayesian analysis and a jackknife estimator for the parameters of the proposed model. For different parameter settings, sample sizes and censoring percentages, various simulation studies are performed and compared to the performance of the log-logistic and log-Burr XII regression models. Besides, we use sensitivity analysis to detect influential or outlying observations, and residual analysis is used to check the assumptions in the model. Finally, we analyze a real data set under log-Buff XII regression models. (C) 2008 Published by Elsevier B.V.

Cutpoint selection for discretizing a continuous covariate for generalized estimating equations

We consider consider the problem of dichotomizing a continuous covariate when performing a regression analysis based on a generalized estimation approach. The problem involves estimation of the cutpoint for the covariate and testing the hypothesis that the binary covariate constructed from the continuous covariate has a significant impact on the outcome. Due to the multiple testing used to find the optimal cutpoint, we need to make an adjustment to the usual significance test to preserve the type-I error rates. We illustrate the techniques on one data set of patients given unrelated hematopoietic stem cell transplantation. Here the question is whether the CD34 cell dose given to patient affects the outcome of the transplant and what is the smallest cell dose which is needed for good outcomes. (C) 2010 Elsevier BM. All rights reserved.

Effectiveness of artemether-lumefantrine provided by community health workers in under-five children with uncomplicated malaria in rural Tanzania : an open label prospective study

Background: Home-management of malaria (HMM) strategy improves early access of anti-malarial medicines to high-risk groups in remote areas of sub-Saharan Africa. However, limited data are available on the effectiveness of using artemisinin-based combination therapy (ACT) within the HMM strategy. The aim of this study was to assess the effectiveness of artemether-lumefantrine (AL), presently the most favoured ACT in Africa, in under-five children with uncomplicated Plasmodium falciparum malaria in Tanzania, when provided by community health workers (CHWs) and administered unsupervised by parents or guardians at home. Methods: An open label, single arm prospective study was conducted in two rural villages with high malaria transmission in Kibaha District, Tanzania. Children presenting to CHWs with uncomplicated fever and a positive rapid malaria diagnostic test (RDT) were provisionally enrolled and provided AL for unsupervised treatment at home. Patients with microscopy confirmed P. falciparum parasitaemia were definitely enrolled and reviewed weekly by the CHWs during 42 days. Primary outcome measure was PCR corrected parasitological cure rate by day 42, as estimated by Kaplan-Meier survival analysis. This trial is registered with ClinicalTrials.gov, number NCT00454961. Results: A total of 244 febrile children were enrolled between March-August 2007. Two patients were lost to follow up on day 14, and one patient withdrew consent on day 21. Some 141/241 (58.5%) patients had recurrent infection during follow-up, of whom 14 had recrudescence. The PCR corrected cure rate by day 42 was 93.0% (95% CI 88.3%-95.9%). The median lumefantrine concentration was statistically significantly lower in patients with recrudescence (97 ng/mL [IQR 0-234]; n = 10) compared with reinfections (205 ng/mL [114-390]; n = 92), or no parasite reappearance (217 [121-374] ng/mL; n = 70; p <= 0.046). Conclusions: Provision of AL by CHWs for unsupervised malaria treatment at home was highly effective, which provides evidence base for scaling-up implementation of HMM with AL in Tanzania.

Default mercado de crédito parcelado para bens duráveis : veículos automotores

The present work aims to study the macroeconomic factors influence in credit risk for installment autoloans operations. The study is based on 4.887 credit operations surveyed in the Credit Risk Information System (SCR) hold by the Brazilian Central Bank. Using Survival Analysis applied to interval censured data, we achieved a model to estimate the hazard function and we propose a method for calculating the probability of default in a twelve month period. Our results indicate a strong time dependence for the hazard function by a polynomial approximation in all estimated models. The model with the best Akaike Information Criteria estimate a positive effect of 0,07% for males over de basic hazard function, and 0,011% for the increasing of ten base points on the operation annual interest rate, toward, for each R$ 1.000,00 on the installment, the hazard function suffer a negative effect of 0,28% , and an estimated elevation of 0,0069% for the same amount added to operation contracted value. For de macroeconomics factors, we find statistically significant effects for the unemployment rate (-0,12%) , for the one lag of the unemployment rate (0,12%), for the first difference of the industrial product index(-0,008%), for one lag of inflation rate (-0,13%) and for the exchange rate (-0,23%). We do not find statistic significant results for all other tested variables.

Falência bancária e capital regulatório: evidência para o Brasil

Este trabalho tem como proposta estudar se o nível de capital sobre os ativos ponderados pelo risco, o Índice de Basileia, mantido pelos intermediários financeiros, pode servir como preditor de falência dos intermediários financeiros. Um dos desafios apresentados foi o fato do Índice de Basileia reportado para cada instituição estar disponível ao público apenas a partir de 2009, ao passo que grande parte das falências bancárias no Brasil ocorreram no período entre 1995 e 2005. Dessa forma, construindo um Índice de Basileia Sintético (IBS) para o período de dezembro de 1995 a dezembro de 2014, testou-se a hipótese de que um nível mais alto de capital em relação aos ativos de risco diminui a probabilidade de falência da instituição. São utilizados modelos logit com variável binária discreta e análise survival, possibilitando estimar o quanto que um aumento no nível de capital proporciona em tempo de vida para a instituição. A amostra estudada é composta por 313 intermediários financeiros atuando no Brasil, tanto de controle público quanto privado, com dados semestrais. Em linha com estudos anteriores, foi encontrada evidência empírica apontando para uma relação inversa entre nível de capital e probabilidade de falência, tanto com o emprego de logit como de survival.

Atividade antiplasmódica e toxicológica de plantas medicinais usadas popularmente no Brasil : uma abordagem etnobotânica

Malaria is a disease of global distribution, recognized by governments around the world as a serious public health problem, affecting more than 109 countries and territories and endangering more than 3.3 billion people. The economic costs of this disease are also relevant: the African continent itself has malaria-related costs of about $ 12 billion annually. Nowadays, in addition to chloroquine, Plasmodium falciparum is resistant to many drugs used in the treatment of malaria, such as amodiaquine, mefloquine, quinine and sulfadoxine-pyrimethamine; resistance of Plasmodium vivax to treatments, although less studied, is also reported. Nature, in general, is responsible for the production of most known organic substances, and the plant kingdom is responsible for the most of the chemical diversity known and reported in the literature. Most medicinal plants commercialized in Brazil, however, are of exotic origin, which makes the search for endemic medicinal plants, besides a patent necessity, a fascinating subject of academic research and development. This study aimed to: (i) verify the antimalarial activity of ethanolic and hydroalcoholic extracts of Boerhavia paniculata Rich. And acetonic extract of Clethra scabra Pers. in Swiss albino mice infected by Plasmodium berghei NK65, (ii) observe possible combined effects between the course of infection by P. berghei NK65 and administration of these extracts in Swiss albino mice, and (iii) conduct a preliminary study of the acute toxicity of these extracts in Swiss albino mice. All extracts notable pharmacological activities - with parasite infections inhibitions ranging from 22% to 54%.These characteristics suggest that the activities are relevant, although comparatively lower than the activity displayed by the positive control group (always above 90%). The general framework of survival analysis demonstrates an overall reduction in survival times for all groups. Necroscopy has not pointed no change in color, shape, size and/or consistency in the evaluated organs - the only exception was the livers of rats submitted to treatment to hydroalcoholic extracts: these organs have been presented in a slightly congestive aspect with mass increasing roughly 28% higher than the other two groups and a p-value of 0.0365. The 250 mg/Kg ethanolic group has been pointed out by the Dunn s post test, as the only class with simultaneous inequalities (p<0.05) between positive and negative control groups. The extracts, notably ethanol extract, have, in fact, a vestigial antimalarial activity, although well below from the ones perceived to chloroquine-treated groups; nevertheless, the survival times of the animals fed with the extracts do not rise by presence of such therapy. Both the toxicopharmacological studies of the synergism between the clinical course of malaria and administration of extracts and the isolated evaluation of toxicity allow us to affirm the absence of toxicity of the extracts at the level of CNS and ANS, as well as their non-influence on food and water consumption patterns, until dosages of 500 mg/Kg. Necroscopic analysis leads us to deduct a possible hepatotoxic effect of hydroalcoholic extract at dosages of 500 mg/Kg, and an innocuous tissue activity of the ethanol extract, in the same dosage. We propose a continuation of the studies of these extracts, with protocol modifications capable of addressing more clearly and objectively their pharmacological and toxicological aspects

Influência de fatores biopsicossociais sobre a mortalidade de idosos no município de Santa Cruz/RN: um estudo prospectivo

The aging process if characterizes for a complex events network, from multidimensional nature, that encloses biological, social, psychic and functional aspects. The alteration of one or more aspects can speed up the aging process, anticipating limitations and until the death in the aged. For an adjusted confrontation of this question is necessary an interdisciplinary vision, in which the some areas of the knowledge can interact and with this to intervenes of the best possible form. Then, information derived from studies of aspects related to incidence, morbidity-mortality and transition patterns, involved in the health-illness process can more accurately identify risk groups thereby establishing links between social factors, illness, incapacity and death. Thus, this study aimed to identify, by a multidimensional vision, the risk factors of mortality in a coorth of elderly in a city in the interior of the state of Rio Grande do Norte (RN), Brazil. A prospective study carried out in Santa Cruz RN, where 310 elderly were randomly selected to form a baseline. The follow-up was 53 months. The predictive variables were divided into sociodemographic, physical health, neuropsychiatric and functional capacity. The statistical analysis carried out by bivariate analysis, survival analysis, followed by binary logistic regression and Cox regression, in the multivariate analysis, considering significant levels p < 0.05 and confidence interval (CI) of 95%. A total of 60 (19.3%) elderly died during the follow-up, where cardiovascular disease was the main cause. The survival was approximately 24.8 months. The study of general survival showed, at 12, 24, 36, and 48 months of observation, a survival rate of 97%, 54%, 31%, and 5% respectively, with a statistical difference in survival only observed for the variables of cognitive function and Basic Activities of Daily Living. In the logistic regression analysis, the risk factors identified were cognitive deficits (OR = 8.74), poor perception of health (OR = 3.89) and dependence for Basic Activities of Daily Living (OR = 3.96). In the Cox analysis, as well as dependence for Basic Activities of Daily Living (HR = 3.17), cognitive deficit (HR = 4.30) and stroke (CVA) (HR = 3.49) continued as independent risk factors for death. The risk factors found in the study can be interpreted as the primary predictors for death among elderly members of the community. Therefore, improvements in health conditions, with actions towards sustaining an autonomous life with special attention for elderly with cognitive impairment, could mean additional healthy quality of life, resulting in the reduction of premature mortality in this population

Metodologia estatística na solução do problema do caixeiro viajante e na avaliação de algoritmos : um estudo aplicado à transgenética computacional

The problems of combinatory optimization have involved a large number of researchers in search of approximative solutions for them, since it is generally accepted that they are unsolvable in polynomial time. Initially, these solutions were focused on heuristics. Currently, metaheuristics are used more for this task, especially those based on evolutionary algorithms. The two main contributions of this work are: the creation of what is called an -Operon- heuristic, for the construction of the information chains necessary for the implementation of transgenetic (evolutionary) algorithms, mainly using statistical methodology - the Cluster Analysis and the Principal Component Analysis; and the utilization of statistical analyses that are adequate for the evaluation of the performance of the algorithms that are developed to solve these problems. The aim of the Operon is to construct good quality dynamic information chains to promote an -intelligent- search in the space of solutions. The Traveling Salesman Problem (TSP) is intended for applications based on a transgenetic algorithmic known as ProtoG. A strategy is also proposed for the renovation of part of the chromosome population indicated by adopting a minimum limit in the coefficient of variation of the adequation function of the individuals, with calculations based on the population. Statistical methodology is used for the evaluation of the performance of four algorithms, as follows: the proposed ProtoG, two memetic algorithms and a Simulated Annealing algorithm. Three performance analyses of these algorithms are proposed. The first is accomplished through the Logistic Regression, based on the probability of finding an optimal solution for a TSP instance by the algorithm being tested. The second is accomplished through Survival Analysis, based on a probability of the time observed for its execution until an optimal solution is achieved. The third is accomplished by means of a non-parametric Analysis of Variance, considering the Percent Error of the Solution (PES) obtained by the percentage in which the solution found exceeds the best solution available in the literature. Six experiments have been conducted applied to sixty-one instances of Euclidean TSP with sizes of up to 1,655 cities. The first two experiments deal with the adjustments of four parameters used in the ProtoG algorithm in an attempt to improve its performance. The last four have been undertaken to evaluate the performance of the ProtoG in comparison to the three algorithms adopted. For these sixty-one instances, it has been concluded on the grounds of statistical tests that there is evidence that the ProtoG performs better than these three algorithms in fifty instances. In addition, for the thirty-six instances considered in the last three trials in which the performance of the algorithms was evaluated through PES, it was observed that the PES average obtained with the ProtoG was less than 1% in almost half of these instances, having reached the greatest average for one instance of 1,173 cities, with an PES average equal to 3.52%. Therefore, the ProtoG can be considered a competitive algorithm for solving the TSP, since it is not rare in the literature find PESs averages greater than 10% to be reported for instances of this size.

Modelos de sobrevivência com fração de cura e erro de medida nas covariáveis

In this work, we study the survival cure rate model proposed by Yakovlev et al. (1993), based on a competing risks structure concurring to cause the event of interest, and the approach proposed by Chen et al. (1999), where covariates are introduced to model the risk amount. We focus the measurement error covariates topics, considering the use of corrected score method in order to obtain consistent estimators. A simulation study is done to evaluate the behavior of the estimators obtained by this method for finite samples. The simulation aims to identify not only the impact on the regression coefficients of the covariates measured with error (Mizoi et al. 2007) but also on the coefficients of covariates measured without error. We also verify the adequacy of the piecewise exponential distribution to the cure rate model with measurement error. At the end, model applications involving real data are made

Modelos de sobrevivência com fração de cura e omissão nas covariáveis

In this work we study the survival cure rate model proposed by Yakovlev (1993) that are considered in a competing risk setting. Covariates are introduced for modeling the cure rate and we allow some covariates to have missing values. We consider only the cases by which the missing covariates are categorical and implement the EM algorithm via the method of weights for maximum likelihood estimation. We present a Monte Carlo simulation experiment to compare the properties of the estimators based on this method with those estimators under the complete case scenario. We also evaluate, in this experiment, the impact in the parameter estimates when we increase the proportion of immune and censored individuals among the not immune one. We demonstrate the proposed methodology with a real data set involving the time until the graduation for the undergraduate course of Statistics of the Universidade Federal do Rio Grande do Norte