Sabemos que um teste t emparelhado é apenas um caso especial de ANOVA unidirecional de medidas repetidas (ou dentro do sujeito), bem como modelo linear de efeitos mistos, que pode ser demonstrado com a função lme (), o pacote nlme em R como mostrado abaixo. #response data from 10 subjects under...
Eu estive analisando a modelagem de efeitos mistos usando o pacote lme4 em R. Estou usando principalmente o lmercomando, por isso vou fazer minha pergunta por meio de código que usa essa sintaxe. Suponho que possa ser uma pergunta fácil geral, é correto comparar dois modelos construídos no lmeruso...
A lógica da imputação múltipla (MI) é imputar os valores ausentes não uma vez, mas várias (normalmente M = 5) vezes, resultando em M conjuntos de dados concluídos. Os conjuntos de dados completos M são então analisados com métodos de dados completos nos quais as estimativas M e seus erros padrão...
Há uma distinção que está me atrapalhando com modelos mistos, e eu estou me perguntando se eu poderia ter alguma clareza sobre isso. Vamos supor que você tenha um modelo misto de dados de contagem. Há uma variável que você sabe que deseja como efeito fixo (A) e outra variável para o tempo (T),...
Digamos que estamos trabalhando em um modelo de efeitos aleatórios de alguns dados de contagem ao longo do tempo e queremos controlar algumas tendências. Normalmente, você faria algo como: lmer(counts ~ dependent_variable + (1+t+I(t^2)|ID), family="poisson") para incluir uma forma quadrática para...
Por favor, considere estes dados: dt.m <- structure(list(id = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12), occasion = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("g1", "g2"), class =...
Estou tentando produzir um modelo para o qual tenho uma variável de resposta que é uma proporção entre 0 e 1, isso inclui alguns 0s e 1s, mas também muitos valores no meio. Estou pensando em tentar uma regressão beta. O pacote que encontrei para R (betareg) permite apenas valores entre 0 e 1, mas...
Basicamente, o que eu quero saber é como as diferentes estruturas de covariância são aplicadas e como os valores dentro dessas matrizes são calculados. Funções como lme () permitem escolher qual estrutura gostaríamos, mas eu adoraria saber como elas são estimadas. Considere o modelo linear de...
Considere os seguintes dados de uma maneira bidirecional no design de assuntos: df <- "http://personality-project.org/r/datasets/R.appendix4.data" df <- read.table(df,header=T) head(df) Observation Subject Task Valence Recall 1 1 Jim Free Neg 8 2 2 Jim Free Neu 9 3 3 Jim Free Pos 5 4 4 Jim...
Em resposta a essa pergunta, sobre se meu projeto, onde eu aleatoriamente apresentei aos participantes fotos de diferentes categorias, era um exemplo em que eu deveria usar uma ANOVA de medidas repetidas, obtive a resposta de que deveria usar um modelo misto, com uma das opções razões são que eu...
Digamos que precisamos dos GLMMs mod1 <- glmer(y ~ x + A + (1|g), data = dat) mod2 <- glmer(y ~ x + B + (1|g), data = dat) Esses modelos não são aninhados no sentido usual de: a <- glmer(y ~ x + A + (1|g), data = dat) b <- glmer(y ~ x + A + B + (1|g), data = dat) então não podemos...
Estou tentando usar lmedo nlmepacote para replicar resultados de aovANOVAs de medidas repetidas. Eu fiz isso para um experimento de medidas repetidas de fator único e para um experimento de dois fatores com um fator entre sujeitos e um fator dentro dos sujeitos, mas estou tendo problemas para fazer...
Olá, tenho dois problemas que parecem candidatos naturais a modelos multiníveis / mistos, que nunca usei. O mais simples, e que espero tentar como introdução, é o seguinte: Os dados se parecem com muitas linhas do formulário x y innergroup outergroup onde x é uma covariável numérica na qual...
Quais são algumas das principais armadilhas do uso de modelos lineares de efeitos mistos? Quais são as coisas mais importantes a serem testadas / observadas na avaliação da adequação do seu modelo? Ao comparar modelos do mesmo conjunto de dados, quais são as coisas mais importantes a serem...
Eu sei que uma das vantagens dos modelos mistos é que eles permitem especificar a matriz de variância-covariância para os dados (simetria composta, auto-regressiva, não estruturada etc.). No entanto, a lmerfunção em R não permite uma especificação fácil dessa matriz. Alguém sabe qual estrutura...
Histórico: Nota: Meu conjunto de dados e código-r estão incluídos abaixo do texto Desejo usar o AIC para comparar dois modelos de efeitos mistos gerados usando o pacote lme4 em R. Cada modelo tem um efeito fixo e um efeito aleatório. O efeito fixo difere entre os modelos, mas o efeito aleatório...
Entendo que usamos modelos de efeitos aleatórios (ou efeitos mistos) quando acreditamos que alguns parâmetros do modelo variam aleatoriamente em algum fator de agrupamento. Desejo ajustar um modelo em que a resposta tenha sido normalizada e centralizada (não perfeitamente, mas bastante próxima) em...
Posso obter ajuda para concluir esta tentativa (em andamento) de tentar me orientar sobre os equivalentes da ANOVA e REGRESSION? Eu tenho tentado conciliar os conceitos, nomenclatura e sintaxe dessas duas metodologias. Há muitos posts sobre este local sobre a sua semelhança, por exemplo, este ou...
Percorri vários livros (Raudenbush & Bryk, Snijders & Bosker, Gelman & Hill, etc.) e vários artigos (Gelman, Jusko, Primo & Jacobsmeier, etc.), e ainda não entendi direito as principais diferenças entre o uso de erros padrão em cluster versam a modelagem multinível. Eu entendo as...
Em particular, como devem ser calculados os erros padrão dos efeitos fixos em um modelo linear de efeitos mistos (no sentido freqüentista)? Fui levado a acreditar que as estimativas típicas ( ), como as apresentadas em Laird e Ware [1982], darão aos SE que são subestimados em tamanho porque os...