## Scripts para a aula de gráficos # ################################## # Tipos de funções gráficas no R # ################################## # plot(galhadores$Altura.folha, galhadores$N.galhas) # nível alto, desenha o gráfico abline(2, 0.05) # nível baixo, adiciona linha de tendênciia # # Feche a janela gráfica e experimente o comando: abline(2, 0.05) # # Clique com o mouse sobre o gráfico. Clique com o botão direito para terminar. locator() # comando iterativo, "lê" o dispositivo grafico locator(3) # lê tres cliques # # ################### # Função plot ( ) # ################### # plot(aves.c$urubu, aves.c$carcara) # ou plot(carcara ~ urubu, data=aves.c) # ou plot(aves.c[,2:3]) # matriz de duas colunas # x<-ts(rpois(30, 5)) # cria 30 valores de uma serie temporal com dados de Poisson plot(x) # gráfico de série temporal # plot(aves.c[,1:2], notch=T) # box-plot # plot(aves.c) # matriz de dispersão: todas as variáveis plot(~ urubu + carcara + seriema, data=aves.c) # variáveis da fórmula x11(); plot(seriema ~ urubu + carcara, data=aves.c) # 'seriema' vs. as demais # # Voltar ao PPT # # ######################### # Dispositivos gráficos # ######################### # pdf("figura1.pdf") plot(aves.c[,1:2]) dev.off() # tiff("figura2.tif") plot(aves.c[,1:2]) dev.off() # jpeg("figura3_boa.jpg", quality = 80) plot(aves.c[,1:2]) dev.off() # jpeg("figura3_ruim.jpg", quality = 10) plot(aves.c[,1:2]) dev.off() # x11() jpeg("figura%02d_boa.jpg", quality = 80) plot(aves.c[,1:2]) plot(aves.c[,2:3]) dev.list() #lista dispositivos abertos dev.cur() # mostra dispositivo ativo dev.set(2) # indica o dispositivo ativo plot(aves.c[,2:3]) dev.cur() dev.off() # fecha dispositivo # # ############################################## # Criando e formatando gráficos de qualidade # ############################################## # # Nesse exercício desmonstrarei os recursos básicos para criação e formatação de gráficos # para apresentação final (publicação, palestra etc). # Vamos usar o conjunto de dados "vespas.csv", o qual consiste no número de vespas, sementes e flores # produzidos em figos de duas árvores de Ficus. # vespas<-read.csv2("vespas.csv") str(vespas) # # Vamos criar um gráfico de dispersão do número de sementes vs. vespas das árvores A e B (simbolos diferentes) # ########################### # Dispersão # plot(Vespas ~ Sementes, data=vespas) plot(Vespas ~ Sementes, data=vespas, xlab="Nº de Sementes", ylab="Nº de Vespas") plot(Vespas ~ Sementes, data=vespas, xlab="Nº de Sementes", ylab="Nº de Vespas", col="red") plot(Vespas ~ Sementes, data=vespas, subset=Arvore=="A", xlim=range(Sementes), ylim=range(Vespas), xlab="Nº de Sementes", ylab="Nº de Vespas", col="red") plot(Vespas ~ Sementes, data=vespas, subset=Arvore=="A", xlim=range(Sementes), ylim=range(Vespas), xlab="Nº de Sementes", ylab="Nº de Vespas", col="red", bty="l") plot(Vespas ~ Sementes, data=vespas, subset=Arvore=="A", xlim=range(Sementes), ylim=range(Vespas), xlab="Nº de Sementes", ylab="Nº de Vespas", col="red", bty="l", las=1) # acrescenta a nova sequencia points(Vespas ~ Sementes, data=vespas, subset=Arvore=="B", pch=15, col="blue") # acrescenta legenda legend(locator(1),legend=c("Arvore A", "Arvore B"),pch=c(1,15),col=c("red","blue"), bty="n") # acrescenta um polígono polygon(locator()) # # Mudar parâmetros globais par(family="serif") # muda a fonte boxplot(Vespas ~ Arvore, data=vespas, xlab="Árvore", ylab="Nº de vespas/figo", main="Diagrama de caixas") # par(mar=c(3,2,2,1)) # modifica as margens, o padrão é 5,4,4,2 plot(Vespas ~ Sementes, data=vespas, xlab="Nº de Sementes", ylab="Nº de Vespas", main="Diagrama de caixas") # par(cex=1.5) plot(Vespas ~ Sementes, data=vespas, xlab="Nº de Sementes", ylab="Nº de Vespas", main="Diagrama de caixas") # par(cex.lab=1.5) plot(Vespas ~ Sementes, data=vespas, xlab="Nº de Sementes", ylab="Nº de Vespas", main="Diagrama de caixas") # par(cex.main=1.5) plot(Vespas ~ Sementes, data=vespas, xlab="Nº de Sementes", ylab="Nº de Vespas", main="Diagrama de caixas") # par(mar=c(4,4,2,1), family="serif", col="red", cex=1.5) plot(Vespas ~ Sementes, data=vespas, xlab="Nº de Sementes", ylab="Nº de Vespas", main="Diagrama de caixas") # # # Dividir a tela em 4... par(mfrow=c(2,2)) # Box-plot boxplot(Vespas ~ Arvore, data=vespas, xlab="Árvore", ylab="Nº de vespas/figo", main="Diagrama de caixas") # # qqnorm de vespas, árvore A qqnorm(vespas$Vespas[vespas$Arvore=="A"], xlab="Quantis teóricos", ylab="Quantis de vespas", main="Árvore A") qqline(vespas$Vespas[vespas$Arvore=="A"]) # # qqnorm de vespas, árvore B qqnorm(vespas$Vespas[vespas$Arvore=="B"], xlab="Quantis teóricos", ylab="Quantis de vespas", main="Árvore B") qqline(vespas$Vespas[vespas$Arvore=="B"]) # # Dispersão plot(Vespas ~ Sementes, data=vespas, subset=Arvore=="A", xlim=range(Sementes), ylim=range(Vespas), xlab="Nº de Sementes", ylab="Nº de Vespas", col="red", main="Disperão") points(Vespas ~ Sementes, data=vespas, subset=Arvore=="B", pch=15, col="blue") par(mfrow=c(1,1))