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")