Estou ajustando um modelo de efeitos aleatórios glmera alguns dados comerciais. O objetivo é analisar o desempenho de vendas por distribuidor, levando em consideração a variação regional. Eu tenho as seguintes variáveis: distcode: ID do distribuidor, com cerca de 800 níveis region: ID geográfico...
Fiquei com a impressão de que a função lmer()no lme4pacote não produzia valores de p (veja lmervalores de p e tudo isso ). Eu tenho usado valores MCMC gerados p vez como por esta pergunta: efeito significativo no lme4modelo misto e esta pergunta: Não é possível encontrar p-valores na saída do...
Eu tenho um conjunto de dados muito pequeno sobre a abundância de abelhas solitárias que estou tendo problemas para analisar. São dados de contagem, e quase todas as contagens estão em um tratamento com a maioria dos zeros no outro tratamento. Também existem alguns valores muito altos (um em cada...
Meus dados são descritos aqui O que pode causar um "erro () é um erro singular" no aov ao ajustar uma ANOVA de medidas repetidas? Estou tentando ver o efeito de uma interação usando, lmerentão meu caso base é: my_null.model <- lmer(value ~ Condition+Scenario+ (1|Player)+(1|Trial), data = my,...
Eu gostaria de obter intervalos de confiança de 95% nas previsões de um nlmemodelo misto não linear . Como nada é fornecido para fazer isso dentro de nlmemim, eu me perguntava se é correto usar o método de "intervalos de previsão da população", conforme descrito no capítulo do livro de Ben Bolker...
Eu sou um usuário mais familiarizado com R e tenho tentado estimar inclinações aleatórias (coeficientes de seleção) há cerca de 35 indivíduos com mais de 5 anos para quatro variáveis de habitat. A variável de resposta é se um local foi habitat "usado" (1) ou "disponível" (0) ("uso"...
O mgcvpacote para Rpossui duas funções para ajustar as interações do produto tensorial: te()e ti(). Entendo a divisão básica do trabalho entre os dois (ajustando uma interação não linear versus decompondo essa interação em efeitos principais e uma interação). O que não entendo é o porquê te(x1,...
Recentemente me disseram que não era possível incorporar covariáveis variáveis no tempo em modelos mistos longitudinais sem introduzir um intervalo de tempo para essas covariáveis. Você pode confirmar / negar isso? Você tem alguma referência sobre esta situação? Proponho uma situação simples...
Como se pode obter pesos de regressão padronizados (efeito fixo) a partir de uma regressão multinível? E, como um "complemento": Qual é a maneira mais fácil de obter esses pesos padronizados de um merobjeto-(da lmerfunção do lme4pacote em
Estou trabalhando em um conjunto de dados para avaliar o impacto da secagem nas atividades microbianas dos sedimentos. O objetivo é determinar se o impacto da secagem varia entre os tipos e / ou profundidade dos sedimentos. O projeto experimental é o seguinte: O primeiro fator Sedimento...
Eu ajustei um modelo mais recente com o seguinte (embora composto): Random effects: Groups Name Std.Dev. day:sample (Intercept) 0.09 sample (Intercept) 0.42 Residual 0.023 Eu realmente gostaria de criar um intervalo de confiança para cada efeito usando a seguinte
Estou executando um glmm com uma variável de resposta binomial e um preditor categórico. O efeito aleatório é dado pelo design aninhado usado para a coleta de dados. Os dados são assim: m.gen1$treatment [1] sucrose control protein control no_injection ..... Levels: no_injection control sucrose...
Minha pergunta é baseada nesta resposta que mostrou qual lme4::lmermodelo corresponde a uma ANOVA de medidas repetidas nos dois sentidos: require(lme4) set.seed(1234) d <- data.frame( y = rnorm(96), subject = factor(rep(1:12, 4)), a = factor(rep(1:2, each=24)), b = factor(rep(rep(1:2,...
Estou terminando algumas análises em um grande conjunto de dados. Eu gostaria de pegar o modelo linear usado na primeira parte do trabalho e reajustá-lo usando um modelo linear misto (LME). O LME seria muito semelhante, exceto que uma das variáveis usadas no modelo seria usada como efeito...
Dadas três variáveis, ye x, que são positivas contínuas e z, que são categóricas, tenho dois modelos candidatos dados por: fit.me <- lmer( y ~ 1 + x + ( 1 + x | factor(z) ) ) e fit.fe <- lm( y ~ 1 + x ) Espero comparar esses modelos para determinar qual modelo é mais apropriado. Parece-me...
Quero modelar duas variáveis de tempo diferentes, algumas das quais são fortemente colineares nos meus dados (idade + coorte = período). Ao fazer isso, tive alguns problemas com ee lmerinterações poly(), mas provavelmente não está limitado a lmer, obtive os mesmos resultados com o...
Estou realizando testes post-hoc em um modelo linear de efeitos mistos em R( lme4pacote). Estou usando o multcomppacote ( glht()função) para executar os testes post-hoc. Meu projeto experimental é de medidas repetidas, com um efeito de bloqueio aleatório. Os modelos são especificados como: mymod...
Para um exemplo simples, assuma que existem dois modelos de regressão linear Modelo 1 tem três preditores, x1a, x2b, ex2c O modelo 2 possui três preditores do modelo 1 e dois preditores adicionais x2aex2b Existe uma equação de regressão populacional em que a variação populacional explicada é...
Eu tenho o seguinte modelo m_plotequipado com lme4::lmerefeitos aleatórios cruzados para participantes ( lfdn) e itens ( content): Random effects: Groups Name Variance Std.Dev. Corr lfdn (Intercept) 172.173 13.121 role1 62.351 7.896 0.03 inference1 24.640 4.964 0.08 -0.30 inference2...
Como posso obter efeitos aleatórios agrupados para lmer após imputação múltipla? Estou usando mouses para atribuir múltiplos quadros de dados. E lme4 para um modelo misto com interceptação aleatória e inclinação aleatória. O pool de pool funciona bem, exceto que ele não combina os efeitos...