Hello,
I remember this dataset from not a long time ago. Please use dput() to
post your data and do not use formatting to post code.
As for the question, though your example is not reproducible (it doesn't
run) see
http://gallery.r-enthusiasts.com/graph/contour_plot_%3Cbr%3EMaunga_Whau_Volcano,21
Hope this helps,
Rui Barradas
Em 19-09-2012 16:55, Tinus Sonnekus escreveu:
Hi All,
I am trying to reproduce this using R instead.
[image: Full-size image (38 K)]
I tried using the following code
*SChla <- read.csv("SM_Chla_data.csv")*
*Atlantis <- SChla[16:66,]*
*head(Atlantis)*
*
*
Seamount Station Depth Pico Nano Micro Total_Ch dbar Latitude
Longitud
16 Atlantis 1217 Surface 0.0639 0.1560 0.0398 0.2597 2.082 -32.71450
57.29733
17 Atlantis 1217 Shallow 0.0305 0.1250 0.0740 0.2295 24.524 -32.71450
57.29733
18 Atlantis 1217 Deep 0.1660 0.3560 0.0734 0.5954 49.573 -32.71450
57.29733
19 Atlantis 1217 Fmax 0.1740 0.5200 0.1830 0.8770 79.595 -32.71450
57.29733
20 Atlantis 1217 Below 0.0235 0.0486 0.0242 0.0963 199.067 -32.71450
57.29733
21 Atlantis 1225 Surface 0.0764 0.1490 0.0333 0.2587 0.409 -32.71133
57.27150
*# Obtaining longitude and depth data*
*Long <- SMfrac[,6]*
*dbar <- -SMfrac[,5] #reverse depth*
*
*
*#converting data to percentage of total chl a*
*
*
*picoFrac <- Atlantis$Pico/Atlantis$Total_Ch*100*
*nanoFrac <- Atlantis$Nano/Atlantis$Total_Ch*100*
*microFrac <- Atlantis$Micro/Atlantis$Total_Ch*100*
*
*
*
*
*library(akima)*
*
*
*par(mfrow = c(3,1)) *
*
*
*contour(interp(Long, dbar, picoFrac),nlevels = 6,*
* xlab ="Longitude", main= "Pico", *
* ylab ="Depth", ylim = c(-200,0), method = "simple",*
* col = "light green")*
*
*
*contour(interp(Long, dbar, nanoFrac), nlevels = 6,*
* xlab ="Longitude", main= "Nano",*
* ylab ="Depth",ylim = c(-200,0), method = "simple",*
* col = "green")*
*
*
*contour(interp(Long, dbar, microFrac),nlevels = 6,*
* xlab ="Longitude", main= "Micro",*
* ylab ="Depth",ylim = c(-200,0), *
* lylab ="Depth", method = "simple",*
* col = "dark green")*
I can get the contours plotted. Not sure if this is 100 % correct but it
looks as expected. Can you guys please help me to add/display my data
points. I would also like to fill the contours in grayscale/colour (e.g.
levels 0, 20, 40, 60, 80. 100 %) and add a legend to the graph. As of yet I
haven't tried to figure out how to add the stations as a secondary axis to
the plot.
Thanks for the help and advice.
Tinus
______________________________________________
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.
