Boa tarde colegas, estou recorrendo ao grupo pois estou com problemas para plotar um gráfico com dois tratamentos (Sol e Sombrite), e o tempo em função do ócio. Gostaria de colocar legenda também. Abaixo tem o modelo usando e também a minha rotina do grafico: mod.10<-glm(formula=ocio ~ tempo%in%trat + I(tempo^2)%in%trat+ I(tempo^3)%in%trat+ I(tempo^4)%in%trat, family=poisson(link=log), data = dados) summary(mod.10) pred.o4<-predict(mod.10) # gráfico pred<-cbind(pred.o4) trat1<-pred[1:280,] trat2<-pred[281:560,] grafico<-predict(mod.10) ordenax<- tempo[order(tempo)] ordenat1<-trat1[order(tempo)] ordenat2<-trat2[order(tempo)] plot(ordenax,ordenat1,col="red", xlim=c(8,18),type="l", xlab="Horas do Dia", ylab="Ócio") line(ordenax,ordenat2, col="blue",xlim=c(8,18),type="l") Daniele Portela de Oliveira Mestranda Produção Animal - UEM Cel: 44 9949 3525 67 9617 0965