Welcome as an Arctic Sea Ice observer. I am finding it incredibly interesting and a real data visualization challenge.
Here's a snippet of R code to help you display the month names rather than day of year in your trend chart.
#### To plot day of year axis ############################# ### Day of year axis setup ## Set up basic day of year vectors (mon_names, 1st day of mon, mo name plot pos) mon_names <- c("Jan", "Feb", "Mar", "April", "May", "June", "July", "Aug", "Sept", "Oct","Nov","Dec") mon_doy <- c(1,32,60,91,121,151,182,213,244,274,305,335) mon_pos <- c(16, 46, 75, 106,135, 165, 200, 228, 255, 289, 320, 355)
## Example DOY plot using month labels plot(DOY, sie, type="l", col = "red", axes=F) axis(side = 1, at=mon_doy, labels=F) axis(side=1, at= mon_pos, labels=mon_names, tick=F, line=F) axis(side=2, at=NULL, labels=T)
Here's how my SIE chart looks using the month names:
Thanks, Kelly I should have checked your site before embarking on this. I've added it to the blogroll to make sure I do so in the future. I'll experiment with your script.
Nick
ReplyDeleteWelcome as an Arctic Sea Ice observer. I am finding it incredibly interesting and a real data visualization challenge.
Here's a snippet of R code to help you display the month names rather than day of year in your trend chart.
#### To plot day of year axis #############################
### Day of year axis setup
## Set up basic day of year vectors (mon_names, 1st day of mon, mo name plot pos)
mon_names <- c("Jan", "Feb", "Mar", "April", "May", "June", "July", "Aug", "Sept", "Oct","Nov","Dec")
mon_doy <- c(1,32,60,91,121,151,182,213,244,274,305,335)
mon_pos <- c(16, 46, 75, 106,135, 165, 200, 228, 255, 289, 320, 355)
## Example DOY plot using month labels
plot(DOY, sie, type="l", col = "red", axes=F)
axis(side = 1, at=mon_doy, labels=F)
axis(side=1, at= mon_pos, labels=mon_names, tick=F, line=F)
axis(side=2, at=NULL, labels=T)
Here's how my SIE chart looks using the month names:
http://chartsgraphs.files.wordpress.com/2010/06/jaxa_sie_07_10_diff_7_06.png
Kelly
Thanks, Kelly
ReplyDeleteI should have checked your site before embarking on this. I've added it to the blogroll to make sure I do so in the future. I'll experiment with your script.