Erro em função lmer do pacote nlme!
Bom dia amigos da lista. Estou tentando ajustar um modelo linear misto a um conjunto de dados (didático), em que Arvore foi considerado como aleatório e Local fixo. Veja os dados e CMR:
dados Altura Local Arvore 1 87 1 1 2 84 2 2 3 75 2 3 4 90 1 4 5 79 2 5
CMR: require(lme4) misto = lmer(Altura ~ Local +(1| Arvore), dados) Mas, a seguinte mensagem de erro aparece: Error in checkNlevels(reTrms$flist, n = n, control) : number of levels of each grouping factor must be < number of observations Alguém saberia dizer como corrigir este problema? Grato. Elton.
Acho que seu erro está no fato de você esta utilizando arvore como fator de agrupamento para o modelo de efeito misto. Você deve utilizar uma variável que você blocou. Pelo seu exemplo acho que árvore está dentro de local (você possui mais de uma árvore no mesmo local) , por isso acho que local deve ser o seu fator de agrupamento.Arvore seria se você avaliasse algo na mesma árvore (p.ex peso de insetos capturados na mesma arvore). Da forma como você fez você tem 5 grupos de árvores e 1 observação em cada grupo. É isso que esta dizendo. Se o fator de agrupamento for local você terá 2 grupos e n repetições. Altura Local Arvore 1 87 1 1 2 84 2 2 3 75 2 3 4 90 1 4 5 79 2 5 Em Seg 05 Mai 2014 14:43:09 BRT, Elton Gean Araújo Araújo escreveu:
Bom dia amigos da lista. Estou tentando ajustar um modelo linear misto a um conjunto de dados (didático), em que Arvore foi considerado como aleatório e Local fixo. Veja os dados e CMR:
dados Altura Local Arvore 1 87 1 1 2 84 2 2 3 75 2 3 4 90 1 4 5 79 2 5
CMR:
require(lme4) misto = lmer(Altura ~ Local +(1| Arvore), dados)
Mas, a seguinte mensagem de erro aparece:
Error in checkNlevels(reTrms$flist, n = n, control) : number of levels of each grouping factor must be < number of observations
Alguém saberia dizer como corrigir este problema? Grato.
Elton.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
Pelo visto só existe um registo de árvore dentro de cada local (n ávores = n registros = uma linha para cada árvore), sendo assim não tem como identificar o modelo que você tá especificando, ou seja, o componente de variância de árvore dentro de local declarado já é na realidade a variância residual. Esse modelo sé seria estimável ou identificável se você tivesse registros dentro das árvores que permitisse compor uma variância num nível de estratificação mais interno. À disposição. Walmes.
Ok Walmes! Obrigado! Elton. Em Segunda-feira, 5 de Maio de 2014 19:22, walmes . <walmeszeviani@gmail.com> escreveu: Pelo visto só existe um registo de árvore dentro de cada local (n ávores = n registros = uma linha para cada árvore), sendo assim não tem como identificar o modelo que você tá especificando, ou seja, o componente de variância de árvore dentro de local declarado já é na realidade a variância residual. Esse modelo sé seria estimável ou identificável se você tivesse registros dentro das árvores que permitisse compor uma variância num nível de estratificação mais interno. À disposição. Walmes. _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
participantes (3)
-
Elton Gean Araújo Araújo -
Fernando Souza -
walmes .