Estou usando o AIC (critério de informações de Akaike) para comparar modelos não lineares em R. É válido comparar os AICs de diferentes tipos de modelo? Especificamente, estou comparando um modelo ajustado pelo glm versus um modelo com um termo de efeito aleatório ajustado pelo glmer (lme4). Caso...
Como o consenso geral parece ser o uso de modelos mistos via lmer()R em vez da ANOVA clássica (pelas razões frequentemente citadas, como desenhos desequilibrados, efeitos aleatórios cruzados etc.), gostaria de tentar com meus dados. No entanto, estou preocupado com a possibilidade de "vender" essa...
O problema: Eu li em outras postagens que predictnão estão disponíveis para lmermodelos de efeitos mistos {lme4} no [R]. Tentei explorar esse assunto com um conjunto de dados de brinquedos ... Fundo: O conjunto de dados é adaptado desta fonte e está disponível como ... require(gsheet) data...
Vasculhei muitos sites de ajuda e ainda estou confuso sobre como especificar termos aninhados mais complicados em um modelo misto também. Eu também estou confuso como o uso de :e /e |em especificação de interações e de nidificação com fatores aleatórios utilizando lmer()no lme4pacote no R. Para o...
Eu tenho a seguinte saída: Generalized linear mixed model fit by the Laplace approximation Formula: aph.remain ~ sMFS2 +sAG2 +sSHDI2 +sbare +season +crop +(1|landscape) AIC BIC logLik deviance 4062 4093 -2022 4044 Random effects: Groups Name Variance Std.Dev. landscape (Intercept) 0.82453...
O lmerTestpacote fornece uma anova()função para modelos mistos lineares com a aproximação opcional de Satterthwaite (padrão) ou Kenward-Roger dos graus de liberdade (df). Qual é a diferença entre essas duas abordagens? Quando escolher
A cotação em bloco abaixo, dos líderes no campo da modelagem de efeitos mistos, afirma que as mudanças de coordenadas nos modelos com correlação zero entre efeitos aleatórios (modelos 'ZCP') alteram as previsões do modelo. Mas, alguém pode elaborar ou justificar ainda mais suas reivindicações? As...
Executei um projeto repetido, no qual testei 30 homens e 30 mulheres em três tarefas diferentes. Quero entender como o comportamento de homens e mulheres é diferente e como isso depende da tarefa. Eu usei o pacote lmer e o lme4 para investigar isso, no entanto, estou tentando verificar suposições...
Atualmente, estou executando alguns modelos lineares de efeito misto. Estou usando o pacote "lme4" em R. Meus modelos assumem a forma: model <- lmer(response ~ predictor1 + predictor2 + (1 | random effect)) Antes de executar meus modelos, verifiquei a possibilidade de multicolinearidade...
Fechadas. Esta questão está fora de tópico . No momento, não está aceitando respostas. Deseja melhorar esta pergunta? Atualize a pergunta para que ela esteja no tópico de Validação cruzada. Fechado há 6 meses . Eu estava trabalhando nos pacotes R nlme e...
Estou tentando deixar de usar o ezpacote lmepara medidas repetidas ANOVA (como espero poder usar contrastes personalizados lme). Seguindo o conselho desta postagem no blog, consegui configurar o mesmo modelo usando os dois aov(como faz ez, quando solicitado) e lme. No entanto, enquanto no exemplo...
Eu estive examinando essa visão geral das fórmulas lm / lmer R por @conjugateprior e fiquei confuso com a seguinte entrada: Agora assuma que A é aleatório, mas B é fixo e B está aninhado em A. aov(Y ~ B + Error(A/B), data=d) Abaixo, a fórmula do modelo misto análogo lmer(Y ~ B + (1 | A:B),...
Gostaria de saber se existem métodos para calcular o tamanho da amostra em modelos mistos? Estou usando lmerno R para ajustar os modelos (tenho inclinações e interceptações
Suponha que eu tenha alguma medida para cada assunto em cada site. Duas variáveis, assunto e site, são de interesse em termos de computação de valores de correlação intraclasse (ICC). Normalmente eu usaria a função lmerdo pacote R lme4e executava lmer(measurement ~ 1 + (1 | subject) + (1 | site),...
Executamos uma regressão logística de efeitos mistos usando a seguinte sintaxe; # fit model fm0 <- glmer(GoalEncoding ~ 1 + Group + (1|Subject) + (1|Item), exp0, family = binomial(link="logit")) # model output summary(fm0) Assunto e Item são os efeitos aleatórios. Estamos obtendo um resultado...
Tenho um experimento de medidas repetidas em que a variável dependente é uma porcentagem e tenho vários fatores como variáveis independentes. Eu gostaria de usar glmero pacote R lme4para tratá-lo como um problema de regressão logística (especificando family=binomial), pois parece acomodar essa...
Com o seguinte conjunto de dados, eu queria ver se a resposta (efeito) muda em relação a sites, temporada, duração e suas interações. Alguns fóruns on-line sobre estatísticas sugeriram que eu continuasse com os Modelos de efeitos mistos lineares, mas o problema é que, como as réplicas são...
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...
Eu tenho trabalhado com alguns dados que têm alguns problemas com medições repetidas. Ao fazer isso, notei um comportamento muito diferente entre lme()e lmer()usando meus dados de teste e quero saber o porquê. O conjunto de dados falsos que criei possui medidas de altura e peso para 10 indivíduos,...
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...