Ludmilla, bom dia!Não cheguei a rodar seu script, mas percebo que as formulas que você está gerando não estão com a referência correta. Os colnames ("resp_1", "pred_1",) estão diferentes da notação da formula (resp1~pred1).# formulas<-paste("resp",1:37,"~pred",1:6,sep='')formulas<-paste("resp_",1:37,"~pred_",1:6,sep='') ### alterar para essa formaOutra coisa é verificar se a função está entendendo as strings como formulas. Pra testar, se a primeira forma não rodar, tente a segunda.glmmadmb(formulas[1],data=data,random=1|~landscape,family="nbinom",zeroInflation=FALSE)glmmadmb(as.formula(formulas[1]),data=data,random=1|~landscape,family="nbinom",zeroInflation=FALSE)Espero que ajude,Em 28 de março de 2015 18:44, Ludmila Rattis <ludmilarattis@gmail.com> escreveu:_______________________________________________Ola,Tenho uma serie de modelos para rodar e quando tento colocar a funcao glmmadmb, pertencente ao pacote glmmADMB, num looping, retorna o seguinte erro:Error in if (rchar == "") rchar <- "1" : missing value where TRUE/FALSE neededProcurei por esse erro no google e nao encontrei nada semelhante. Quando tento fora do looping, a funcao roda normalmente. Eis um codigo de exemplo:install.packages("glmmADMB", repos="http://r-forge.r-project.org", type="source")data<-matrix(0,180,43);for(i in 1:43){data[,i]<-rnorm(180,mean=0,sd=1) }landscape<-factor(sort(rep(1:30,6)))resp<-paste("resp_",1:37,sep='')pred<-paste("pred_",1:6,sep='')data<-cbind(data,landscape)colnames(data)<-c(resp,pred,"landscape")formulas<-paste("resp",1:37,"~pred",1:6,sep='')fitted.glmm<-list()for (i in 1:length(formulas)){fitted.glmm[[]]<-glmmadmb(formulas[i],data=data,random=1|~landscape,family="nbinom",zeroInflation=FALSE)}Alguem poderia me dizer o que esta acontecendo?Muito obrigada, Lud--Ludmila Rattis
Programa de Pós-graduação em Ecologia/UNICAMP
Geomatics and Landscape Ecology Research Laboratory (Carleton University - Canada)
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.
_______________________________________________
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.