Modelos log-lineares em tabelas de contingência

Os modelos log-lineares permitem enriquecer bastante a análise e interpretação das tabelas de contingência. Embora a nível teórico a sua importância já tenha sido reconhecida há bastante tempo, a nível da sua aplicação prática só há relativamente pouco tempo é que foi reconhecida, devido, sobretudo, às dificuldades de cálculo que lhe são inerentes e que só se resolveram completamente com o desenvolvimento dos computadores e do software adequado. Neste trabalho apresentam-se os métodos básicos da análise log-linear de tabelas de contingência bidimensionais e tridimensionais.

A Log-Linear Regression Model for the Beta-Weibull Distribution

We introduce the log-beta Weibull regression model based on the beta Weibull distribution (Famoye et al., 2005; Lee et al., 2007). We derive expansions for the moment generating function which do not depend on complicated functions. The new regression model represents a parametric family of models that includes as sub-models several widely known regression models that can be applied to censored survival data. We employ a frequentist analysis, a jackknife estimator, and a parametric bootstrap for the parameters of the proposed model. We derive the appropriate matrices for assessing local influences on the parameter estimates under different perturbation schemes and present some ways to assess global influences. Further, for different parameter settings, sample sizes, and censoring percentages, several simulations are performed. In addition, the empirical distribution of some modified residuals are displayed and compared with the standard normal distribution. These studies suggest that the residual analysis usually performed in normal linear regression models can be extended to a modified deviance residual in the proposed regression model applied to censored data. We define martingale and deviance residuals to evaluate the model assumptions. The extended regression model is very useful for the analysis of real data and could give more realistic fits than other special regression models.

Análise multivariada: um exemplo usando modelo log-linear

Apresenta-se de forma resumida análise multivariada de dados categóricos, usando modelo log-linear para a situação de uma tabela de contingência 2 x 2 x 2.

A Log-linear model of grain size influence on the geochemistry of sediments

Sediment composition is mainly controlled by the nature of the source rock(s), and chemical (weathering) and physical processes (mechanical crushing, abrasion, hydrodynamic sorting) during alteration and transport. Although the factors controlling these processes are conceptually well understood, detailed quantification of compositional changes induced by a single process are rare, as are examples where the effects of several processes can be distinguished. The present study was designed to characterize the role of mechanical crushing and sorting in the absence of chemical weathering. Twenty sediment samples were taken from Alpine glaciers that erode almost pure granitoid lithologies. For each sample, 11 grain-size fractions from granules to clay (ø grades &-1 to &9) were separated, and each fraction was analysed for its chemical composition.The presence of clear steps in the box-plots of all parts (in adequate ilr and clr scales) against ø is assumed to be explained by typical crystal size ranges for the relevant mineral phases. These scatter plots and the biplot suggest a splitting of the full grain size range into three groups: coarser than ø=4 (comparatively rich in SiO2, Na2O, K2O, Al2O3, and dominated by “felsic” minerals like quartz and feldspar), finer than ø=8 (comparatively rich in TiO2, MnO, MgO, Fe2O3, mostly related to “mafic” sheet silicates like biotite and chlorite), and intermediate grains sizes (4≤ø &8; comparatively rich in P2O5 and CaO, related to apatite, some feldspar).To further test the absence of chemical weathering, the observed compositions were regressed against three explanatory variables: a trend on grain size in ø scale, a step function for ø≥4, and another for ø≥8. The original hypothesis was that the trend could be identified with weathering effects, whereas each step function would highlight those minerals with biggest characteristic size at its lower end. Results suggest that this assumption is reasonable for the step function, but that besides weathering some other factors (different mechanical behavior of minerals) have also an important contribution to the trend.Key words: sediment, geochemistry, grain size, regression, step function

A Log-linear model of grain size influence on the geochemistry of sediments

Sediment composition is mainly controlled by the nature of the source rock(s), and chemical (weathering) and physical processes (mechanical crushing, abrasion, hydrodynamic sorting) during alteration and transport. Although the factors controlling these processes are conceptually well understood, detailed quantification of compositional changes induced by a single process are rare, as are examples where the effects of several processes can be distinguished. The present study was designed to characterize the role of mechanical crushing and sorting in the absence of chemical weathering. Twenty sediment samples were taken from Alpine glaciers that erode almost pure granitoid lithologies. For each sample, 11 grain-size fractions from granules to clay (ø grades <-1 to >9) were separated, and each fraction was analysed for its chemical composition. The presence of clear steps in the box-plots of all parts (in adequate ilr and clr scales) against ø is assumed to be explained by typical crystal size ranges for the relevant mineral phases. These scatter plots and the biplot suggest a splitting of the full grain size range into three groups: coarser than ø=4 (comparatively rich in SiO2, Na2O, K2O, Al2O3, and dominated by “felsic” minerals like quartz and feldspar), finer than ø=8 (comparatively rich in TiO2, MnO, MgO, Fe2O3, mostly related to “mafic” sheet silicates like biotite and chlorite), and intermediate grains sizes (4≤ø <8; comparatively rich in P2O5 and CaO, related to apatite, some feldspar). To further test the absence of chemical weathering, the observed compositions were regressed against three explanatory variables: a trend on grain size in ø scale, a step function for ø≥4, and another for ø≥8. The original hypothesis was that the trend could be identified with weathering effects, whereas each step function would highlight those minerals with biggest characteristic size at its lower end. Results suggest that this assumption is reasonable for the step function, but that besides weathering some other factors (different mechanical behavior of minerals) have also an important contribution to the trend. Key words: sediment, geochemistry, grain size, regression, step function

Análisis y aplicación del modelo log-linear a la investigación educativa.

Analizar el modelo log-linear y sus posibilidades de aplicación en la investigación educativa. Al mismo tiempo se hace un estudio de un medio didáctico: el libro de texto, analizando su relación con la función docente del profesor. 308 profesores de EGB de las dos provincias canarias. El muestreo fue ocasional. Se trata de un diseño ex post facto. Se procedió a la aplicación piloto de un cuestionario en una muestra reducida y a su análisis por parte de un grupo de expertos. Después de las modificaciones oportunas, se procedió a su aplicación final con la colaboración de diversos encuestadores para su distribución y recogida. Las variables principales fueron las siguientes. Rreferidas al profesor: años de experiencia, grado de dependencia del libro de texto y ciclo. Referidas a medios didácticos: frecuencia de uso, razones de uso, finalidad didáctica y dimensiones más valoradas para la enseñanza. Cuestionario 'Uso de medios en la enseñanza'. No existe entre el profesorado, considerado globalmente, una tendencia mayoritaria por la dependencia e independencia hacia el libro de texto. Los profesores veteranos tienden a ser dependientes del libro de texto. El resto de los profesores no se inclinan por la dependencia o independencia. La relación del texto con los programas oficiales sólo es valorada por los profesores del ciclo medio. La dimensión curricular centrada en la metodología que más importancia merece son las actividades que propone el texto, seguida del planteamiento metodológico que se desprende de la guía didáctica. La dimensión más valorada es la adecuación del texto, seguida del planteamiento metodológico que se desprende de la guía didáctica. La dimensión más valorada es la adecuación del texto al nivel de conocimientos de los alumnos. Le sigue en importancia el lenguaje utilizado y, finalmente, los aspectos formales del texto (colorido, tamaño, ilustraciones, etc.). Se pone de manifiesto el papel del libro de texto como un medio destinado básicamente a uso del alumnado. Su uso para el profesor se limita a servir de apoyo en sus explicaciones, en tanto que motivar y evaluar el aprendizaje son funciones con las cuales parece incompatible el uso del libro de texto. El análisis log-linear constituye un poderoso instrumento de análisis de variables nominales, con un grado de sofisticación estadística solo disponible hasta ahora para variables continuas. La abundancia de variables nominales en la investigación educativa, le hace especialmente apropiado para nuestro campo. Las ventajas del análisis log-linear dependen de: la naturaleza de las variables, mínimo número de categorías si se incluyen datos continuos, puntos de corte, estrategias de muestreo, etc. Sigue sin disponerse de criterios claros con respecto al tamaño de la muestra y la interpretación de la intensidad de los parámetros. Tampoco se ha desarrollado un sistema de representación gráfica con esta técnica de análisis.

Linkage disequilibrium assessment via log-linear modeling of SNP haplotype frequencies

Analyses of high-density single-nucleotide polymorphism (SNP) data, such as genetic mapping and linkage disequilibrium (LD) studies, require phase-known haplotypes to allow for the correlation between tightly linked loci. However, current SNP genotyping technology cannot determine phase, which must be inferred statistically. In this paper, we present a new Bayesian Markov chain Monte Carlo (MCMC) algorithm for population haplotype frequency estimation, particulary in the context of LD assessment. The novel feature of the method is the incorporation of a log-linear prior model for population haplotype frequencies. We present simulations to suggest that 1) the log-linear prior model is more appropriate than the standard coalescent process in the presence of recombination (>0.02cM between adjacent loci), and 2) there is substantial inflation in measures of LD obtained by a "two-stage" approach to the analysis by treating the "best" haplotype configuration as correct, without regard to uncertainty in the recombination process. Genet Epidemiol 25:106-114, 2003. (C) 2003 Wiley-Liss, Inc.

Linear vs. log-linear unit-root specification: an application of mis-specification encompassing

The objective of this paper is to apply the mis-specification (M-S) encompassing perspective to the problem of choosing between linear and log-linear unit-root models. A simple M-S encompassing test, based on an auxiliary regression stemming from the conditional second moment, is proposed and its empirical size and power are investigated using Monte Carlo simulations. It is shown that by focusing on the conditional process the sampling distributions of the relevant statistics are well behaved under both the null and alternative hypotheses. The proposed M-S encompassing test is illustrated using US total disposable income quarterly data.

Estimativas de correlações genéticas entre escores visuais e características de carcaça medidas por ultrassonografia em bovinos Nelore utilizando modelos bayesianos linear-limiar

O objetivo neste estudo foi estimar as correlações genéticas entre escores visuais e características de carcaça medidas por ultrassonografia em bovinos da raça Nelore utilizando a estatística bayesiana por meio da Amostragem de Gibbs, sob modelo animal linear-limiar. Foram estudadas as características categóricas morfológicas de musculosidade, estrutura física, conformação e sacro, avaliadas aos 15 e 22 meses de idade. Para as características de carcaça, foram avaliadas as características área de olho-de-lombo, espessura de gordura subcutânea, espessura de gordura subcutânea na garupa e altura na garupa. Os escores visuais devem ser empregados como critérios de seleção para aumentar o progresso genético para a característica área de olhode-lombo e, consequentemente, melhorar o rendimento de carcaça. As estimativas de correlação genética obtidas para musculosidade com espessura de gordura subcutânea e espessura de gordura subcutânea na garupa indicaram que a seleção para musculosidade pode levar a animais com melhor acabamento de carcaça. A seleção para a estrutura física e conformação aos 15 e 22 meses de idade pode promover resposta correlacionada para o aumento da altura na garupa.

A log-linear regression model for the beta-Birnbaum-Saunders distribution with censored data

The beta-Birnbaum-Saunders (Cordeiro and Lemonte, 2011) and Birnbaum-Saunders (Birnbaum and Saunders, 1969a) distributions have been used quite effectively to model failure times for materials subject to fatigue and lifetime data. We define the log-beta-Birnbaum-Saunders distribution by the logarithm of the beta-Birnbaum-Saunders distribution. Explicit expressions for its generating function and moments are derived. We propose a new log-beta-Birnbaum-Saunders regression model that can be applied to censored data and be used more effectively in survival analysis. We obtain the maximum likelihood estimates of the model parameters for censored data and investigate influence diagnostics. The new location-scale regression model is modified for the possibility that long-term survivors may be presented in the data. Its usefulness is illustrated by means of two real data sets. (C) 2011 Elsevier B.V. All rights reserved.

Modeling multilevel sleep transitional data via Poisson log-linear multilevel models

This paper proposes Poisson log-linear multilevel models to investigate population variability in sleep state transition rates. We specifically propose a Bayesian Poisson regression model that is more flexible, scalable to larger studies, and easily fit than other attempts in the literature. We further use hierarchical random effects to account for pairings of individuals and repeated measures within those individuals, as comparing diseased to non-diseased subjects while minimizing bias is of epidemiologic importance. We estimate essentially non-parametric piecewise constant hazards and smooth them, and allow for time varying covariates and segment of the night comparisons. The Bayesian Poisson regression is justified through a re-derivation of a classical algebraic likelihood equivalence of Poisson regression with a log(time) offset and survival regression assuming piecewise constant hazards. This relationship allows us to synthesize two methods currently used to analyze sleep transition phenomena: stratified multi-state proportional hazards models and log-linear models with GEE for transition counts. An example data set from the Sleep Heart Health Study is analyzed.

Perfil da demanda atendida em ambulatório de doenças profissionais e a presença de lesões por esforços repetitivos

OBJETIVO: Caracterizar o perfil do paciente atendido em ambulatório de doenças profissionais de um hospital geral e, em especial, o portador da condição patológica mais freqüente: as lesões por esforços repetitivos (LER). MÉTODOS: Foram analisadas as 565 fichas-resumo de atendimento ambulatorial em saúde do trabalhador, em 1996, em função de variáveis sociodemográficas, de características do trabalho e de diagnóstico. Foram adotados os seguintes procedimentos: descrição univariável de todas as variáveis; análise bivariada, utilizando-se o teste do qui-quadrado (chi²) para variáveis categóricas; comparação de médias ou amplitudes para variável categórica e outra quantitativa; análise multivariada por meio de modelos log-linear, com transformação logit, para detectar os grupos com maior proporção de LER. RESULTADOS: Houve predomínio no atendimento de pacientes do sexo feminino com menos de 40 anos, sendo que, nos casos de LER, esse predomínio mostrou-se mais acentuado. As ocupações agrupadas apareceram em ordem decrescente, na seqüência: administrativos, operários e prestadores de serviços, sendo 59,6% pertencentes ao ramo de atividade econômica de serviços. As LER corresponderam a 56% dos diagnósticos feitos nas primeiras consultas. CONCLUSÕES: Confirmou-se a tendência de mudança no perfil epidemiológico das demandas atendidas no ambulatório de doenças profissionais, no final da década de 80, quando emergiram as LER. Deslocaram-se para uma minoria as antigas doenças profissionais decorrentes da exposição a riscos químicos e físicos, em função daquelas relacionadas às condições ergonômicas inadequadas.

Tendência temporal e desigualdades na mortalidade por diarreias em menores de 5 anos

OBJETIVO: Analisar a tendência da mortalidade por diarreia entre menores de 5 anos, no município de Osasco (SP), entre 1980 e 2000. MÉTODOS: Trata-se de estudo observacional com dois delineamentos. Um descritivo, que toma o indivíduo como unidade do estudo, e outro ecológico, analisando agregado populacional que incluiu análise de séries temporais. A fonte de dados foi o sistema de informação de mortalidade do Estado de São Paulo e censos de 1980, 1991 e 2000. Descreveu-se a variação sazonal e para a análise de tendência aplicaram-se modelos log lineares de regressão polinomiais, utilizando-se variáveis sociodemográficas da criança e da mãe. Foram analisadas a evolução de indicadores sociodemográficos do município de 1980 a 2000, as taxas médias de mortalidade por diarreia nos menores de 5 anos e seus diferenciais por distrito nos anos 90. RESULTADOS: Dos 1.360 óbitos, 94,3 e 75,3% atingiram, respectivamente, menores de 1 ano e de 6 meses. O declínio da mortalidade foi de 98,3%, com deslocamento da sazonalidade do verão para o outono. A mediana da idade elevou-se de 2 meses nos primeiros períodos para 3 meses no último. O resíduo de óbitos manteve-se entre filhos de mães de 20 a 29 anos e escolaridade < 8 anos. O risco relativo entre o distrito mais atingido e a taxa média do município diminuiu de 3,4 para 1,3 do primeiro para o segundo quinquênio dos anos 90. CONCLUSÃO: Nossos resultados apontam uma elevação da idade mais vulnerável e a provável mudança do agente mais frequentemente associado ao óbito por diarreia.