Recentemente, comecei a aprender sobre modelos mistos lineares generalizados e estava usando o R para explorar que diferença faz para tratar a associação ao grupo como efeito fixo ou aleatório. Em particular, estou analisando o exemplo de conjunto de dados discutido
Devo incluir efeitos aleatórios em um modelo, mesmo que não sejam estatisticamente significativos? Eu tenho um desenho experimental de medidas repetidas, no qual cada indivíduo experimenta três tratamentos diferentes em ordem aleatória. Eu gostaria de controlar quaisquer efeitos individuais e de...
Espero que alguém possa ajudar com o que considero uma pergunta relativamente simples, e acho que sei a resposta, mas sem confirmação ela se tornou algo que eu simplesmente não posso ter certeza. Eu tenho alguns dados de contagem como uma variável de resposta e quero medir como essa variável muda...
Introdução Tenho participantes que tocam repetidamente superfícies contaminadas com E. coli em duas condições ( A = usando luvas, B = sem luvas). Quero saber se há uma diferença entre a quantidade de bactérias na ponta dos dedos com e sem luvas, mas também entre o número de contatos. Ambos os...
Eu tenho algumas perguntas sobre especificação e interpretação de GLMMs. Definitivamente, três perguntas são estatísticas e duas são mais específicas sobre R. Estou publicando aqui porque, em última análise, acho que a questão é a interpretação dos resultados do GLMM. Atualmente, estou tentando...
Primeiro vamos simular alguns dados para uma regressão logística com partes fixas e aleatórias: set.seed(1) n <- 100 x <- runif(n) z <- sample(c(0,1), n, replace=TRUE) b <- rnorm(2) beta <- c(0.4, 0.8) X <- model.matrix(~x) Z <- cbind(z, 1-z) eta <- X%*%beta + Z%*%b pr...
Estou tendo dificuldades para acessar a mesma página que meu supervisor quando se trata de validar meu modelo. Analisei os resíduos (observados em relação aos valores ajustados) e usei isso como argumento para discutir os resultados obtidos pelo meu modelo, no entanto, meu supervisor insiste que a...
Estou usando o lme4pacote no R para fazer alguma modelagem logística de efeitos mistos. Meu entendimento era que a soma de cada efeito aleatório deveria ser zero. Quando eu uso modelos mistos lineares de brinquedo lmer, os efeitos aleatórios geralmente são < confirmando minha crença de que...
Considere o seguinte exemplo simples: library( rms ) library( lme4 ) params <- structure(list(Ns = c(181L, 191L, 147L, 190L, 243L, 164L, 83L, 383L, 134L, 238L, 528L, 288L, 214L, 502L, 307L, 302L, 199L, 156L, 183L), means = c(0.09, 0.05, 0.03, 0.06, 0.07, 0.07, 0.1, 0.1, 0.06, 0.11, 0.1,...
Eu tenho dados de coleta de longo prazo e gostaria de testar se o número de animais coletados é influenciado pelos efeitos climáticos. Meu modelo se parece abaixo: glmer(SumOfCatch ~ I(pc.act.1^2) +I(pc.act.2^2) + I(pc.may.1^2) + I(pc.may.2^2) + SampSize + as.factor(samp.prog) + (1|year/month),...
Eu li os ótimos comentários sobre como lidar com valores ausentes antes de aplicar o SVD, mas gostaria de saber como ele funciona com um exemplo simples: Movie1 Movie2 Movie3 User1 5 4 User2 2 5 5 User3 3 4 User4 1 5 User5 5 1 5 Dada a matriz acima, se eu remover os valores de NA, acabarei...
Meu design é o seguinte. yyy é a resposta de Bernoulli x1x1x_1 é uma variável contínua x2x2x_2 é uma variável categórica (fator) com dois níveis O experimento é completamente dentro dos sujeitos. Ou seja, cada sujeito recebe cada combinação de e .x 2x1x1x_1x2x2x_2 Esta é uma configuração de...
O que eu quero fazer é construir GLMMs para avaliar a seleção de recursos, e eu tenho um conjunto de variáveis (algumas representando distâncias e outras representando% da cobertura do solo). Posso testar a correlação entre variáveis antes de padronizá-las? Não tenho muita certeza do que devo...