Ajuda para plotar um gráfico

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

Daniele, Forneça sempre código mínimo reproduzível. Segue uma possível solução http://www.leg.ufpr.br/doku.php/ridiculas?grafico_de_valores_observados_e_... Consulte a página de documentação das funções. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

Obrigada Walmes, o grafico deu certo. Abraços Daniele Portela de Oliveira Mestranda Produção Animal - UEM Cel: 44 9949 3525 67 9617 0965 ________________________________ De: Walmes Zeviani <walmeszeviani@gmail.com> Para: r-br@listas.c3sl.ufpr.br Enviadas: Segunda-feira, 4 de Julho de 2011 16:54:18 Assunto: Re: [R-br] Ajuda para plotar um gráfico Daniele, Forneça sempre código mínimo reproduzível. Segue uma possível solução http://www.leg.ufpr.br/doku.php/ridiculas?grafico_de_valores_observados_e_... Consulte a página de documentação das funções. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
participantes (2)
-
Daniele Portela de oliveira
-
Walmes Zeviani