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

Modelos mistos (também conhecidos como multiníveis ou hierárquicos) são modelos lineares que incluem efeitos fixos e efeitos aleatórios. Eles são usados ​​para modelar dados longitudinais ou aninhados.

1
Extraindo inclinações para casos de um modelo de efeitos mistos (lme4)
Gostaria de extrair as inclinações para cada indivíduo em um modelo de efeito misto, conforme descrito no parágrafo a seguir Modelos de efeitos mistos foram usados ​​para caracterizar caminhos individuais de mudança nas medidas de resumo cognitivo, incluindo termos para idade, sexo e anos de educação como efeitos fixos (Laird …
12 r  mixed-model 



1
Calcular probabilidade logarítmica "à mão" para regressão não-linearizada dos mínimos quadrados generalizada (nlme)
Estou tentando calcular a probabilidade de log para uma regressão de mínimos quadrados não linear generalizada para a função f(x)=β1(1+xβ2)β3f(x)=β1(1+xβ2)β3f(x)=\frac{\beta_1}{(1+\frac x\beta_2)^{\beta_3}}otimizado pelagnlsfunção no pacote Rnlme, usando a matriz de covariância de variância gerada pelas distâncias em uma árvore filogenética assumindo movimento browniano (corBrownian(phy=tree)daapeembalagem). O seguinte código R reprodutível se ajusta …



2
Modelo misto com 1 observação por nível
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 de nível superior (norte, sul, leste, …

1
Reconciliando notações para modelos mistos
Estou familiarizado com notações como: ondeβ0j=β0+uj, eyij=β0+βixij+uj+eij=β0j+βixij+eijyij=β0+βixij+uj+eij=β0j+βixij+eij\begin{align} y_{ij} &= \beta_0 + \beta_i x_{ij} + u_j + e_{ij}\\ &= \beta_{0j} + \beta_i x_{ij} + e_{ij} \end{align}β0j=β0+ujβ0j=β0+uj\beta_{0j}=\beta_{0}+u_j ondeβ0j=β0+u0jeβ1j=β1+u1jyij=β0+β1xij+u0j+u1jxij+eij=β0j+β1jxij+eijyij=β0+β1xij+u0j+u1jxij+eij=β0j+β1jxij+eij\begin{align} y_{ij} &= \beta_0 + \beta_1 x_{ij} + u_{0j} + u_{1j} x_{ij} + e_{ij} \\ &= \beta_{0j} + \beta_{1j} x_{ij} + e_{ij} \end{align}β0j=β0+u0jβ0j=β0+u0j\beta_{0j}=\beta_{0}+u_{0j}β1j=β1+u1jβ1j=β1+u1j\beta_{1j}=\beta_1+u_{1j} para um …

1
Regressão linear com medidas repetidas em R
Não consegui descobrir como executar a regressão linear em R in para um design de medida repetida. Em uma pergunta anterior (ainda sem resposta), me foi sugerido não usar, lmmas usar modelos mistos. Eu usei lmda seguinte maneira: lm.velocity_vs_Velocity_response <- lm(Velocity_response~Velocity*Subject, data=mydata) (mais detalhes sobre o conjunto de dados podem …

1
Como otimizar a eficiência computacional ao ajustar um modelo complexo a um grande conjunto de dados repetidamente?
Estou tendo problemas de desempenho usando o MCMCglmmpacote no R para executar um modelo de efeitos mistos. O código fica assim: MC1<-MCMCglmm(bull~1,random=~school,data=dt,family="categorical" , prior=list(R=list(V=1,fix=1), G=list(G1=list(V=1, nu=0))) , slice=T, nitt=iter, ,burnin=burn, verbose=F) Existem cerca de 20.000 observações nos dados e elas estão agrupadas em cerca de 200 escolas. Soltei todas as …

1

2
Como se pode fazer um teste de hipótese do MCMC em um modelo de regressão de efeito misto com inclinações aleatórias?
A biblioteca languageR fornece um método (pvals.fnc) para realizar testes de significância do MCMC dos efeitos fixos em um modelo de regressão de efeito misto ajustado usando o lmer. No entanto, pvals.fnc comete um erro quando o modelo lmer inclui inclinações aleatórias. Existe uma maneira de fazer um teste de …




Ao utilizar nosso site, você reconhece que leu e compreendeu nossa Política de Cookies e nossa Política de Privacidade.
Licensed under cc by-sa 3.0 with attribution required.