
Olhe novamente para o seu exemplo. O fator 'Animal' está com seis observações!! Deveria estar como mesmo número de observações das demais fontes de variação e variáveis, ou seja, cinco observações. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman@yahoo.com.br Tel: (35)3826-6608/9900-2924

Os dados da variável animal está maior do que as demais. O erro que está mostrando é sobre isso. Abraços, Victor Eduardo Em 20 de maio de 2011 08:59, Ivan Bezerra Allaman <ivanalaman@yahoo.com.br>escreveu:
Olhe novamente para o seu exemplo. O fator 'Animal' está com seis observações!! Deveria estar como mesmo número de observações das demais fontes de variação e variáveis, ou seja, cinco observações.
Allaman (S,f,P)
*M.Sc Ivan Bezerra Allaman* Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA **email e msn - ivanalaman@yahoo.com.br Tel: (35)3826-6608/9900-2924
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

Certo Ivan. Corrigido o problema de digitação. Mas o erro persiste.??????? Rode os dados abaixo: Herd<-as.factor(c(1,2,1,2,1)) Animal<-as.factor(c(1,2,3,4,5)) Sire<-c(0,0,0,1,4) Dam<-c(0,0,0,2,3) Phenotyp<-c(100,130,120,110,140) require(pedigreemm) p<- pedigree(sire = Sire, dam = Dam, label= 1:5) model<-pedigreemm(Phenotyp~ (1|Animal)+ Herd, pedigree=list(Animal=p)) --- Em sex, 20/5/11, Ivan Bezerra Allaman <ivanalaman@yahoo.com.br> escreveu: De: Ivan Bezerra Allaman <ivanalaman@yahoo.com.br> Assunto: Re: [R-br] Modelo Animal - Pedigreemm Para: "R Brasil" <r-br@listas.c3sl.ufpr.br> Data: Sexta-feira, 20 de Maio de 2011, 11:59 Olhe novamente para o seu exemplo. O fator 'Animal' está com seis observações!! Deveria estar como mesmo número de observações das demais fontes de variação e variáveis, ou seja, cinco observações. Allaman(S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman@yahoo.com.br Tel: (35)3826-6608/9900-2924 -----Anexo incorporado----- _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

Eu só consegui reproduzir o exemplo depois de ler o artigo. Recomendo o mesmo. São dois artigos citados na library. Veja na documentação! Valeu!!! Fábio Mathias Corrêa Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e da Terra - DCET Campus Soane Nazaré de Andrade, km 16 Rodovia Ilhéus-Itabuna CEP 45662-900. Ilhéus-Bahia Tel.: 73-3680-5076 ________________________________ De: Marcelo Cardoso mello <marcelo_cm32@yahoo.com.br> Para: r-br@listas.c3sl.ufpr.br; Ivan Bezerra Allaman <ivanalaman@yahoo.com.br> Enviadas: Sexta-feira, 20 de Maio de 2011 13:04 Assunto: Re: [R-br] Modelo Animal - Pedigreemm Certo Ivan. Corrigido o problema de digitação. Mas o erro persiste. ??????? Rode os dados abaixo: Herd<-as.factor(c(1,2,1,2,1)) Animal<-as.factor(c(1,2,3,4,5)) Sire<-c(0,0,0,1,4) Dam<-c(0,0,0,2,3) Phenotyp<-c(100,130,120,110,140) require(pedigreemm) p<- pedigree(sire = Sire, dam = Dam, label= 1:5) model<-pedigreemm(Phenotyp~ (1|Animal)+ Herd, pedigree=list(Animal=p)) --- Em sex, 20/5/11, Ivan Bezerra Allaman <ivanalaman@yahoo.com.br> escreveu:
De: Ivan Bezerra Allaman <ivanalaman@yahoo.com.br> Assunto: Re: [R-br] Modelo Animal - Pedigreemm Para: "R Brasil" <r-br@listas.c3sl.ufpr.br> Data: Sexta-feira, 20 de Maio de 2011, 11:59
Olhe novamente para o seu exemplo. O fator 'Animal' está com seis observações!! Deveria estar como mesmo número de observações das demais fontes de variação e variáveis, ou seja, cinco observações.
Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman@yahoo.com.br Tel: (35)3826-6608/9900-2924 -----Anexo incorporado-----
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

Ok Fábio, verifiquei no Technical note: An R package for fitting generalized linear mixed models in animal. Os exemplos tem dados do próprio pacote e consegui reproduzir. A questão é que estes exemplos tem mais de uma observação por indivíduo e no exemplo que estou pedindo um help tenho indivíduos sem repetição. Este é o problema que gostaria de contornar. Se tiver como resolver e um material sobre esta solução específica me indique. Se não tiver solução via R me avise também. Grato. Marcelo. --- Em sex, 20/5/11, Fabio Mathias Corrêa <fabio.ufla@yahoo.com.br> escreveu: De: Fabio Mathias Corrêa <fabio.ufla@yahoo.com.br> Assunto: Re: [R-br] Modelo Animal - Pedigreemm Para: "r-br@listas.c3sl.ufpr.br" <r-br@listas.c3sl.ufpr.br> Data: Sexta-feira, 20 de Maio de 2011, 16:08 Eu só consegui reproduzir o exemplo depois de ler o artigo. Recomendo o mesmo. São dois artigos citados na library. Veja na documentação! Valeu!!! Fábio Mathias Corrêa Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e da Terra - DCET Campus Soane Nazaré de Andrade, km 16 Rodovia Ilhéus-Itabuna CEP 45662-900. Ilhéus-Bahia Tel.: 73-3680-5076 De: Marcelo Cardoso mello <marcelo_cm32@yahoo.com.br> Para: r-br@listas.c3sl.ufpr.br; Ivan Bezerra Allaman <ivanalaman@yahoo.com.br> Enviadas: Sexta-feira, 20 de Maio de 2011 13:04 Assunto: Re: [R-br] Modelo Animal - Pedigreemm Certo Ivan. Corrigido o problema de digitação. Mas o erro persiste.??????? Rode os dados abaixo: Herd<-as.factor(c(1,2,1,2,1)) Animal<-as.factor(c(1,2,3,4,5)) Sire<-c(0,0,0,1,4) Dam<-c(0,0,0,2,3) Phenotyp<-c(100,130,120,110,140) require(pedigreemm) p<- pedigree(sire = Sire, dam = Dam, label= 1:5) model<-pedigreemm(Phenotyp~ (1|Animal)+ Herd, pedigree=list(Animal=p)) --- Em sex, 20/5/11, Ivan Bezerra Allaman <ivanalaman@yahoo.com.br> escreveu: De: Ivan Bezerra Allaman <ivanalaman@yahoo.com.br> Assunto: Re: [R-br] Modelo Animal - Pedigreemm Para: "R Brasil" <r-br@listas.c3sl.ufpr.br> Data: Sexta-feira, 20 de Maio de 2011, 11:59 Olhe novamente para o seu exemplo. O fator 'Animal' está com seis observações!! Deveria estar como mesmo número de observações das demais fontes de variação e variáveis, ou seja, cinco observações. Allaman(S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman@yahoo.com.br Tel: (35)3826-6608/9900-2924 -----Anexo incorporado----- _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br -----Anexo incorporado----- _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

Qual o tipo de variável resposta que vc tem? O fato de vc ter apenas uma repetição por indivíduo pode ser o problema! A solução que posso oferecer é com o uso de Inferência Bayesiana. Tem uma library muito boa para glmm utilizando bayesiana, ela se chama MCMCglmm. Vc pode colocar a estrutura de covariância que desejar. Ela é muito rápida, pelo menos nos exemplos que testei!!! Qualquer coisa é só falar! Valeu!!! Fábio Mathias Corrêa Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e da Terra - DCET Campus Soane Nazaré de Andrade, km 16 Rodovia Ilhéus-Itabuna CEP 45662-900. Ilhéus-Bahia Tel.: 73-3680-5076 ________________________________ De: Marcelo Cardoso mello <marcelo_cm32@yahoo.com.br> Para: r-br@listas.c3sl.ufpr.br; Fabio Mathias Corrêa <fabio.ufla@yahoo.com.br> Enviadas: Domingo, 22 de Maio de 2011 15:39 Assunto: Re: [R-br] Modelo Animal - Pedigreemm Ok Fábio, verifiquei no Technical note: An R package for fitting generalized linear mixed models in animal. Os exemplos tem dados do próprio pacote e consegui reproduzir. A questão é que estes exemplos tem mais de uma observação por indivíduo e no exemplo que estou pedindo um help tenho indivíduos sem repetição. Este é o problema que gostaria de contornar. Se tiver como resolver e um material sobre esta solução específica me indique. Se não tiver solução via R me avise também. Grato. Marcelo. --- Em sex, 20/5/11, Fabio Mathias Corrêa <fabio.ufla@yahoo.com.br> escreveu:
De: Fabio Mathias Corrêa <fabio.ufla@yahoo.com.br> Assunto: Re: [R-br] Modelo Animal - Pedigreemm Para: "r-br@listas.c3sl.ufpr.br" <r-br@listas.c3sl.ufpr.br> Data: Sexta-feira, 20 de Maio de 2011, 16:08
Eu só consegui reproduzir o exemplo depois de ler o artigo. Recomendo o mesmo.
São dois artigos citados na library. Veja na documentação!
Valeu!!! Fábio Mathias Corrêa
Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e da Terra - DCET
Campus Soane Nazaré de Andrade, km 16 Rodovia Ilhéus-Itabuna CEP 45662-900. Ilhéus-Bahia
Tel.: 73-3680-5076
________________________________ De: Marcelo Cardoso mello <marcelo_cm32@yahoo.com.br> Para: r-br@listas.c3sl.ufpr.br; Ivan Bezerra Allaman <ivanalaman@yahoo.com.br> Enviadas: Sexta-feira, 20 de Maio de 2011 13:04 Assunto: Re: [R-br] Modelo Animal - Pedigreemm
Certo Ivan. Corrigido o problema de digitação. Mas o erro persiste. ???????
Rode os dados abaixo:
Herd<-as.factor(c(1,2,1,2,1)) Animal<-as.factor(c(1,2,3,4,5)) Sire<-c(0,0,0,1,4) Dam<-c(0,0,0,2,3) Phenotyp<-c(100,130,120,110,140) require(pedigreemm) p<- pedigree(sire = Sire, dam = Dam, label= 1:5) model<-pedigreemm(Phenotyp~ (1|Animal)+ Herd, pedigree=list(Animal=p)) --- Em sex, 20/5/11, Ivan Bezerra Allaman <ivanalaman@yahoo.com.br> escreveu:
De: Ivan Bezerra Allaman <ivanalaman@yahoo.com.br> Assunto: Re: [R-br] Modelo Animal - Pedigreemm Para: "R Brasil" <r-br@listas.c3sl.ufpr.br> Data: Sexta-feira, 20 de Maio de 2011, 11:59
Olhe novamente para o seu exemplo. O fator 'Animal' está com seis observações!! Deveria estar como mesmo número de observações das demais fontes de variação e variáveis, ou seja, cinco observações.
Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman@yahoo.com.br Tel: (35)3826-6608/9900-2924 -----Anexo incorporado-----
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
-----Anexo incorporado-----
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

Ok Fábio, vou conferir a análise bayesiana quando tiver um tempo. Quanto a solução do problema, encontrei um texto na web que pode ajudar, mas não estou conseguindo aplicar. Parece que tem haver com alteração nas funções. Mas não entendo muito desta parte de programação. Segue abaixo. Se vc enxergar uma luz e puder aplicar neste exemplo avise. Obrigado. pedigreemm with one observation per individual Here is a somewhat crude fix; it adds an additional parameter (checklevels) to the lmer_finalize and glmer_finalize functions in lme4 that turns off the unwanted sanity check. The default is checklevels=TRUE, which maintains the previous behavior. The pedigreemm function is then changed to set checklevels=FALSE before calling the appropriate finalizing function. I'm running these changes locally and they seem to work fine. Comments or suggestions are welcome, as is consideration of this change or some variant on it for inclusion in the pedigreemm and lme4 packages. ----- in lmer.R in the lme4 package ----- lmer_finalize <- function(fr, FL, start, REML, verbose, checklevels=TRUE) { ... ### This checks that the number of levels in a grouping factor < n ### Only need to check the first factor because it is the one with ### the most levels. if (checklevels & !(length(levels(dm$flist[[1]])) < length(Y))) ... glmer_finalize <- function(fr, FL, glmFit, start, nAGQ, verbose, checklevels=TRUE) { ... ### This checks that the number of levels in a grouping factor < n ### Only need to check the first factor because it is the one with ### the most levels. if (checklevels & !(length(levels(dm$flist[[1]])) < ncol(dm$Zt))) ... ----- in pedigree.R in the pedigreemm package ----- pedigreemm <- function(formula, data, family = NULL, REML = TRUE, pedigree = list(), ... lmf$checklevels <- FALSE; ans <- do.call(if (!is.null(lmf$glmFit)) lme4:::glmer_finalize else lme4:::lmer_finalize, lmf) ... --- Em dom, 22/5/11, Fabio Mathias Corrêa <fabio.ufla@yahoo.com.br> escreveu: De: Fabio Mathias Corrêa <fabio.ufla@yahoo.com.br> Assunto: Re: [R-br] Modelo Animal - Pedigreemm Para: "r-br@listas.c3sl.ufpr.br" <r-br@listas.c3sl.ufpr.br> Data: Domingo, 22 de Maio de 2011, 22:22 Qual o tipo de variável resposta que vc tem? O fato de vc ter apenas uma repetição por indivíduo pode ser o problema! A solução que posso oferecer é com o uso de Inferência Bayesiana. Tem uma library muito boa para glmm utilizando bayesiana, ela se chama MCMCglmm. Vc pode colocar a estrutura de covariância que desejar. Ela é muito rápida, pelo menos nos exemplos que testei!!! Qualquer coisa é só falar! Valeu!!! Fábio Mathias Corrêa Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e da Terra - DCET Campus Soane Nazaré de Andrade, km 16 Rodovia Ilhéus-Itabuna CEP 45662-900. Ilhéus-Bahia Tel.: 73-3680-5076 De: Marcelo Cardoso mello <marcelo_cm32@yahoo.com.br> Para: r-br@listas.c3sl.ufpr.br; Fabio Mathias Corrêa <fabio.ufla@yahoo.com.br> Enviadas: Domingo, 22 de Maio de 2011 15:39 Assunto: Re: [R-br] Modelo Animal - Pedigreemm Ok Fábio, verifiquei no Technical note: An R package for fitting generalized linear mixed models in animal. Os exemplos tem dados do próprio pacote e consegui reproduzir. A questão é que estes exemplos tem mais de uma observação por indivíduo e no exemplo que estou pedindo um help tenho indivíduos sem repetição. Este é o problema que gostaria de contornar. Se tiver como resolver e um material sobre esta solução específica me indique. Se não tiver solução via R me avise também. Grato. Marcelo. --- Em sex, 20/5/11, Fabio Mathias Corrêa <fabio.ufla@yahoo.com.br> escreveu: De: Fabio Mathias Corrêa <fabio.ufla@yahoo.com.br> Assunto: Re: [R-br] Modelo Animal - Pedigreemm Para: "r-br@listas.c3sl.ufpr.br" <r-br@listas.c3sl.ufpr.br> Data: Sexta-feira, 20 de Maio de 2011, 16:08 Eu só consegui reproduzir o exemplo depois de ler o artigo. Recomendo o mesmo. São dois artigos citados na library. Veja na documentação! Valeu!!! Fábio Mathias Corrêa Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e da Terra - DCET Campus Soane Nazaré de Andrade, km 16 Rodovia Ilhéus-Itabuna CEP 45662-900. Ilhéus-Bahia Tel.: 73-3680-5076 De: Marcelo Cardoso mello <marcelo_cm32@yahoo.com.br> Para: r-br@listas.c3sl.ufpr.br; Ivan Bezerra Allaman <ivanalaman@yahoo.com.br> Enviadas: Sexta-feira, 20 de Maio de 2011 13:04 Assunto: Re: [R-br] Modelo Animal - Pedigreemm Certo Ivan. Corrigido o problema de digitação. Mas o erro persiste.??????? Rode os dados abaixo: Herd<-as.factor(c(1,2,1,2,1)) Animal<-as.factor(c(1,2,3,4,5)) Sire<-c(0,0,0,1,4) Dam<-c(0,0,0,2,3) Phenotyp<-c(100,130,120,110,140) require(pedigreemm) p<- pedigree(sire = Sire, dam = Dam, label= 1:5) model<-pedigreemm(Phenotyp~ (1|Animal)+ Herd, pedigree=list(Animal=p)) --- Em sex, 20/5/11, Ivan Bezerra Allaman <ivanalaman@yahoo.com.br> escreveu: De: Ivan Bezerra Allaman <ivanalaman@yahoo.com.br> Assunto: Re: [R-br] Modelo Animal - Pedigreemm Para: "R Brasil" <r-br@listas.c3sl.ufpr.br> Data: Sexta-feira, 20 de Maio de 2011, 11:59 Olhe novamente para o seu exemplo. O fator 'Animal' está com seis observações!! Deveria estar como mesmo número de observações das demais fontes de variação e variáveis, ou seja, cinco observações. Allaman(S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman@yahoo.com.br Tel: (35)3826-6608/9900-2924 -----Anexo incorporado----- _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br -----Anexo incorporado----- _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br -----Anexo incorporado----- _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
participantes (4)
-
Fabio Mathias Corrêa
-
Ivan Bezerra Allaman
-
Marcelo Cardoso mello
-
Victor Eduardo