Existe uma razão intuitiva para que os efeitos aleatórios sejam reduzidos em relação ao seu valor esperado no modelo misto linear
Existe uma razão intuitiva para que os efeitos aleatórios sejam reduzidos em relação ao seu valor esperado no modelo misto linear
Estrutura de dados > str(data) 'data.frame': 6138 obs. of 10 variables: $ RT : int 484 391 422 516 563 531 406 500 516 578 ... $ ASCORE : num 5.1 4 3.8 2.6 2.7 6.5 4.9 2.9 2.6 7.2 ... $ HSCORE : num 6 2.1 7.9 1 6.9 8.9 8.2 3.6 1.7 8.6 ... $ MVMNT : Factor w/ 2 levels "_Withd","Appr": 2 2 1...
Eu tenho um GLMM do formulário: lmer(present? ~ factor1 + factor2 + continuous + factor1*continuous + (1 | factor3), family=binomial) Quando uso drop1(model, test="Chi"), obtenho resultados diferentes dos que utilizo Anova(model, type="III")na embalagem do carro ou summary(model). Estes dois...
Estou tentando ajustar um modelo misto multivariado (isto é, resposta múltipla) R. Além dos pacotes ASReml-re SabreR(que requerem software externo), parece que isso só é possível no MCMCglmm. No artigo que acompanha o MCMCglmmpacote (pp.6), Jarrod Hadfield descreve o processo de ajustar um modelo...
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...
Alguém sabe como calcular (ou extrair) a alavancagem e as distâncias de Cook para um merobjeto de classe (obtido através do lme4pacote)? Eu gostaria de plotá-los para uma análise de
Entendo o efeito de um efeito aleatório categórico em um modelo de efeitos mistos, na medida em que ele executa um agrupamento parcial das observações por nível no efeito aleatório, assumindo efetivamente que as observações não são independentes, mas apenas seus conjuntos parciais. Também para meu...
Estou trabalhando em um conjunto de dados. Depois de usar algumas técnicas de identificação de modelos, criei um modelo ARIMA (0,2,1). Usei a detectIOfunção no pacote TSAem R para detectar um outlier inovador (IO) na 48ª observação do meu conjunto de dados original. Como faço para incorporar esse...
Sempre fui ensinado que efeitos aleatórios influenciam apenas a variação (erro) e que efeitos fixos influenciam apenas a média. Mas encontrei um exemplo em que efeitos aleatórios também influenciam a média - a estimativa do coeficiente: require(nlme) set.seed(128) n <- 100 k <- 5 cat <-...
Há vários tópicos no R-sig-ME sobre a obtenção de intervalos de confiança para previsões usando lme4e nlmeem R. Por exemplo, aqui e aqui em 2010, incluindo alguns comentários de Dougals Bates, um dos autores de ambos os pacotes. Hesito em citá-lo literalmente, por medo de serem retirados do...
Esta pergunta / tópico surgiu em uma discussão com um colega e eu estava procurando algumas opiniões sobre isso: Estou modelando alguns dados usando uma regressão logística de efeitos aleatórios, mais precisamente uma regressão logística de interceptação aleatória. Para os efeitos fixos, tenho 9...
No modelo misto, assumimos que os efeitos aleatórios (parâmetros) são variáveis aleatórias que seguem distribuições normais. Parece muito semelhante ao método bayesiano, no qual todos os parâmetros são assumidos como aleatórios. Então, o modelo de efeito aleatório é um caso especial do método...
Estou analisando os resultados de um experimento de tempo de reação em R. Executei uma ANOVA de medidas repetidas (1 fator dentro do sujeito com 2 níveis e 1 fator entre sujeitos com 2 níveis). Eu executei um modelo misto linear semelhante e queria resumir os resultados anteriores na forma de...
Qual é a estrutura de variância-covariância padrão para efeitos aleatórios no glmerou lmerno lme4pacote? Como se especifica outra estrutura de variância-covariância para efeitos aleatórios no código? Não encontrei nenhuma informação sobre isso na
Eu estou analisando alguns dados comportamentais usando lme4em R, principalmente após excelentes tutoriais do Bodo de Inverno , mas eu não entendo se eu estou lidando com interações corretamente. Pior, ninguém mais envolvido nesta pesquisa usa modelos mistos, então estou um pouco à deriva quando se...
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...
Antes de fazer essa pergunta, pesquisei em nosso site e encontrei muitas perguntas semelhantes (como aqui , aqui e aqui ). Mas acho que essas perguntas relacionadas não foram bem respondidas ou discutidas, portanto, gostaria de levantar essa questão novamente. Eu acho que deve haver uma grande...
Espero sinceramente que tenha formulado esta pergunta de forma que ela possa ser respondida definitivamente - caso contrário, informe-me e tentarei novamente! Eu também acho que devo usar R para essas análises. Eu tenho várias medidas plant performance (Ys)que eu suspeito que foram influenciadas...
Estou confuso com declarações em uma página da UCLA sobre regressão logística de efeitos mistos. Eles mostram uma tabela de coeficientes de efeitos fixos ao se ajustar a esse modelo e o primeiro parágrafo abaixo parece interpretar os coeficientes exatamente como uma regressão logística normal. Mas...
Considere um experimento com vários participantes humanos, cada um medido várias vezes em duas condições. Um modelo de efeitos mistos pode ser formulado (usando a sintaxe lme4 ) como: fit = lmer( formula = measure ~ (1|participant) + condition ) Agora, digamos que eu queira gerar intervalos de...