Perguntas com a marcação «mixed-model»

Suponhamos que tem alguns alguma variável resposta yijyijy_{ij} que foi medido a partir de jjj th irmão em th família. Além disso, alguns dados comportamentais foram coletados ao mesmo tempo de cada sujeito. Estou tentando analisar a situação com o seguinte modelo linear de efeitos mistos:iiixijxijx_{ij} yij=α0+α1xij+δ1ixij+εijyij=α0+α1xij+δ1ixij+εijy_{ij} = \alpha_0 + …

14 mixed-model  lme4-nlme  covariance-matrix  non-independent 

Alguém pode me dizer como R estimar o ponto de interrupção em um modelo linear por partes (como um parâmetro fixo ou aleatório), quando eu também precisar estimar outros efeitos aleatórios? Incluí um exemplo de brinquedo abaixo que se encaixa em uma regressão de taco de hóquei / taco quebrado …

14 r  mixed-model  lme4-nlme  change-point  piecewise-linear 

Douglas Bates afirma que os seguintes modelos são equivalentes "se a matriz de variância-covariância para efeitos aleatórios com valor vetorial tiver uma forma especial, chamada simetria composta" ( slide 91 desta apresentação ): m1 <- lmer(y ~ factor + (0 + factor|group), data) m2 <- lmer(y ~ factor + (1|group) …

13 r  anova  mixed-model  repeated-measures  lme4-nlme 

Vejo mencionado em vários lugares que a ANOVA faz sua estimativa usando o método dos momentos. Estou confuso com essa afirmação porque, embora não esteja familiarizado com o método dos momentos, meu entendimento é que é algo diferente e não equivalente ao método da máxima probabilidade; por outro lado, a …

13 anova  mixed-model  maximum-likelihood  method-of-moments 

Aprendi nas estatísticas elementares que, com um modelo linear geral, para que as inferências sejam válidas, as observações devem ser independentes. Quando ocorre o agrupamento, a independência não pode mais ser mantida, levando a inferência inválida, a menos que isso seja explicado. Uma maneira de explicar esse cluster é usar …

13 r  mixed-model  inference  independence 

Vamos considerar este conjunto de dados hipotéticos: set.seed(12345) num.subjects <- 10 dose <- rep(c(1,10,50,100), num.subjects) subject <- rep(1:num.subjects, each=4) group <- rep(1:2, each=num.subjects/2*4) response <- dose*dose/10 * group + rnorm(length(dose), 50, 30) df <- data.frame(dose=dose, response=response, subject=subject, group=group) podemos usar lmepara modelar a resposta com um modelo de efeito aleatório: …

13 mixed-model  linear-model  prediction 

Suponha que eu tenha medido alguma variável em irmãos, que estão aninhados nas famílias. A estrutura de dados fica assim: valor do irmão da família ------ ------- ----- 1 1 y_11 1 2 y_12 2 1 y_21 2 2 y_22 2 3 y_23 ... ... ... Quero saber a correlação …

13 mixed-model  intraclass-correlation 

Alguém pode me dizer a diferença entre usar aov()e lme()analisar dados longitudinais e como interpretar os resultados desses dois métodos? Abaixo, eu analisar o mesmo conjunto de dados usando aov()e lme()e tem 2 resultados diferentes. Com aov(), obtive um resultado significativo no tempo por interação com o tratamento, mas, ao …

13 r  mixed-model  repeated-measures  panel-data 

Ouvi dizer que os modelos LME são mais sólidos na análise de dados de precisão (ou seja, em experimentos de psicologia), na medida em que podem trabalhar com distribuições binomiais e outras não-normais que as abordagens tradicionais (por exemplo, ANOVA) não podem. Qual é a base matemática dos modelos LME …

13 mixed-model 

Eu tenho lido sobre o cálculo de valores de R2R2R^2 em modelos mistos e depois de ler as perguntas frequentes do R-sig, outras postagens neste fórum (eu vincularia algumas, mas não tenho reputação suficiente) e várias outras referências que eu entendo que usando valores no contexto de modelos mistos é …

13 r  mixed-model  r-squared  lme4-nlme 

Estive analisando alguns dados usando a modelagem linear de efeitos mistos em R. Estou planejando fazer um pôster com os resultados e fiquei imaginando se alguém com modelos de efeitos mistos poderia sugerir quais gráficos usar para ilustrar os resultados do modelo. Eu estava pensando em gráficos residuais, gráficos de …

13 r  data-visualization  mixed-model 

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 …

12 r  mixed-model  confidence-interval  lme4-nlme 

Na modelagem multinível, a correlação intraclasse geralmente é calculada a partir de uma ANOVA de efeitos aleatórios yeu j= γ00+ uj+ eeu jyEuj=γ00+vocêj+eEuj y_{ij} = \gamma_{00} + u_j + e_{ij} onde são os resíduos de nível 2 e são os resíduos de nível 1. Em seguida, obtemos estimativas e para …

12 correlation  mixed-model  intraclass-correlation 

TL; DR: a lme4otimização parece linear no número de parâmetros do modelo por padrão e é muito mais lenta que um glmmodelo equivalente com variáveis ​​dummy para grupos. Existe algo que eu possa fazer para acelerar isso? Estou tentando ajustar um modelo de logit hierárquico bastante grande (~ 50k linhas, …

12 r  mixed-model  optimization  lme4-nlme  numerics 

Visão geral da pergunta Aviso: Esta pergunta requer muita configuração. Por favor, tenha paciência comigo. Um colega meu e eu estamos trabalhando em um experimento. O design deve contornar um grande número de restrições, que listarei abaixo. Eu desenvolvi um design que satisfaz as restrições e que nos fornece estimativas …

12 mixed-model  experiment-design  simulation  confounding