Preciso plotar um gráfico de variáveis climáticas sendo o eixo Y esquerdo referente à precipitação e o eixo Y direito referente às temperaturas min, med e max. A precipitação deve ser apresentada em barras e as temperaturas em linhas. Como posso fazer isso? x=c(seq(1:12)) rainfall=c(17.54,7.68,9.87,5.45,0.15,1.61,0.01,0.33,1.53,4.85,3.77,9.35)##barras tmax=c(30.21,30.71,27.23,28.23,25.19,24.39,26.59,28.39,29.77,28.49,28.61,29.79)##linhas tmin=c(16.82,16.36,15.91,15.8,11.57,9.07,13.07,12.65,12.98,15.47,16.21,17.26)##linhas tmed=c(23.51,23.53,21.57,22.02,18.38,16.73,19.83,20.52,21.38,21.98,22.41,23.53)##linhas Muito obrigado pela ajuda. Abraço. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil
Bom, de inicio acredito que vai ser algo assim: barras<-barplot(rainfall,ylim=c(0,35)) points(tmed~barras,type="l") points(tmin~barras,type="l",lty=2) points(tmax~barras,type="l",lty=2) axis(4) #vc vai salvar a posição das barras no vertor barras #ai depois note que vc esta fazendo as linhas # e por ultimo use o comando axes pra fazer o segundo ou mais eixos, se quiser algo no eixo X tb #leia a documentação pra deixa do jeito que gostar ?axis #agora não sei se entra nas suas necessidades mas #de uma olhada nesse pacote e nos graficos que da pra fazer com ele #é exatamente o tipo de dado que vc ilustrou #por exemplo: install.packages("climatol") require(climatol) data(datcli) #exemplo dele datcli #seus dados t(data.frame(rainfall,tmax,tmin,tmed)) #exemplo do comando diagwl(datcli,est="Example station",alt=100,per="1961-90",mlab="en") #seus dados de exemplo diagwl(t(data.frame(rainfall,tmax,tmin,tmed)),est="Example station",alt=100,per="1961-90",mlab="en") Em 30 de julho de 2012 15:18, Marcelo Claro de Souza <marcelo_claro@yahoo.com.br> escreveu:
Preciso plotar um gráfico de variáveis climáticas sendo o eixo Y esquerdo referente à precipitação e o eixo Y direito referente às temperaturas min, med e max. A precipitação deve ser apresentada em barras e as temperaturas em linhas. Como posso fazer isso?
x=c(seq(1:12)) rainfall=c(17.54,7.68,9.87,5.45,0.15,1.61,0.01,0.33,1.53,4.85,3.77,9.35)##barras tmax=c(30.21,30.71,27.23,28.23,25.19,24.39,26.59,28.39,29.77,28.49,28.61,29.79)##linhas tmin=c(16.82,16.36,15.91,15.8,11.57,9.07,13.07,12.65,12.98,15.47,16.21,17.26)##linhas tmed=c(23.51,23.53,21.57,22.02,18.38,16.73,19.83,20.52,21.38,21.98,22.41,23.53)##linhas
Muito obrigado pela ajuda. Abraço.
Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil
_______________________________________________ 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.
-- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056
participantes (2)
-
Augusto Ribas -
Marcelo Claro de Souza