RESUMO
This paper proposes a flexible modeling approach for so-called comet assay data regularly encountered in preclinical research. While such data consist of non-Gaussian outcomes in a multilevel hierarchical structure, traditional analyses typically completely or partly ignore this hierarchical nature by summarizing measurements within a cluster. Non-Gaussian outcomes are often modeled using exponential family models. This is true not only for binary and count data, but also for, example, time-to-event outcomes. Two important reasons for extending this family are for (1) the possible occurrence of overdispersion, meaning that the variability in the data may not be adequately described by the models, which often exhibit a prescribed mean-variance link, and (2) the accommodation of a hierarchical structure in the data, owing to clustering in the data. The first issue is dealt with through so-called overdispersion models. Clustering is often accommodated through the inclusion of random subject-specific effects. Though not always, one conventionally assumes such random effects to be normally distributed. In the case of time-to-event data, one encounters, for example, the gamma frailty model (Duchateau and Janssen, 2007 ). While both of these issues may occur simultaneously, models combining both are uncommon. Molenberghs et al. ( 2010 ) proposed a broad class of generalized linear models accommodating overdispersion and clustering through two separate sets of random effects. Here, we use this method to model data from a comet assay with a three-level hierarchical structure. Although a conjugate gamma random effect is used for the overdispersion random effect, both gamma and normal random effects are considered for the hierarchical random effect. Apart from model formulation, we place emphasis on Bayesian estimation. Our proposed method has an upper hand over the traditional analysis in that it (1) uses the appropriate distribution stipulated in the literature; (2) deals with the complete hierarchical nature; and (3) uses all information instead of summary measures. The fit of the model to the comet assay is compared against the background of more conventional model fits. Results indicate the toxicity of 1,2-dimethylhydrazine dihydrochloride at different dose levels (low, medium, and high).
Assuntos
Teorema de Bayes , Ensaio Cometa/estatística & dados numéricos , Algoritmos , Análise de Variância , Animais , Análise por Conglomerados , Técnicas Citológicas , Dano ao DNA , Interpretação Estatística de Dados , Dimetilidrazinas/toxicidade , Relação Dose-Resposta a Droga , Avaliação Pré-Clínica de Medicamentos , Fígado/citologia , Fígado/efeitos dos fármacos , Masculino , Modelos Estatísticos , Ratos , Resultado do TratamentoRESUMO
Several pharmacological studies involve experiments aimed at testing for a difference between experimental groups wherein the data are longitudinal in nature, frequently with long sequences per subject. Oftentimes, treatment effect, if present, is not constant over time. In such situations, imposing a parametric mean structure can be too complicated and/or restrictive. A more flexible approach is to model the mean using a semiparametric smooth function, estimated using, for example, penalized smoothing splines. We formulate a series of models exhibiting how the group-specific mean profiles could possibly differ. Once an appropriate model is chosen, interest lies in identifying specific time points where the groups differ. For this purpose, we propose the use of simultaneous confidence bands around the fitted models wherein the bands take into account within and between-subject variability, as well as variability arising from smoothing.