Hi,
I'd really appreciate if someone could give me some help or advice about
this - I've tried everything I know and am clueless about how to proceed!
I've written a script to import ASCII data of raster maps, bin them into
categories, calculate the mean values within these bins and plot the two in
a simple graph. I'm running into problems with my x axis, as R cannot add
the bin categories as labels and simply outputs the factor values of the
number of bins.
In much more detail:
One set of data is fractional values for the amount of cropland in a
designated area (raster cell). The second is annual average temperature
data, in the same format. I import these using:
xRaw<-scan(".../Temp.txt")
yRaw<-scan(".../Cropland.txt")
and then designate cells with values of -9999 as NAs:
x<-ifelse(xRaw==-9999,NA,xRaw)
y<-ifelse(yRaw==-9999,NA,yRaw)
Then, I compile them into a 2 column data frame, and exclude the NA values:
zRaw<-data.frame(x,y)
z<-na.exclude(zRaw)
I am then left with data frame 'z' which is a refined list of data. Using
the 'cut' function, I assign each row into one of 69 different bins:
z$bins<-cut(z$x, breaks=c(-36:33), include.lowest=TRUE)
Each bin is now 1 degree Celsius wide. Within each one, I calculate the mean
value of the corresponding fractional cropland data:
a<-tapply(z$y, z$bins, mean)
Object 'a' now contains each bin category and the associated mean cropland
values, so 69 points in total.
Now, a problem arises when I try and plot these values using:
plot(a, xlab="Temperature", ylab="Fraction of Cropland", pch=20)
What happens is that R intuitively labels my x-axis from 1-70, which
corresponds to the bin factor values that were returned when the 'cut'
function was used. I would like to simply label the x-axis with my
temperature values, from -36 degrees to 33 degrees, instead of having
numbers 1-69. How can I do this??
Thanks,
Jeremy Proville
M.Sc. Candidate
Bioresource Engineering, McGill University
--
View this message in context:
http://www.nabble.com/Bin-Category-Labels-on-Axis-tp24042325p24042325.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
