[R] filled contour plot with contour lines

From: <Steve_Friedman_at_nps.gov>
Date: Fri, 01 Apr 2011 14:48:31 -0400

I'm stumped, can anyone find my error in this sequence.

 for(j in 1:(varsize[4]-1))

      temp <- get.var.ncdf(nc=input,
varid="p_foraging",c(1,1,j),c(varsize[1],varsize[2],1))

         filled.contour(x, y, temp, color = terrain.colors,
             plot.title = title(main = paste("Everglades Wood Stork
Foraging Potential \nYear", (2000+j)),
             xlab = "UTM East", ylab = "UTM North") ,
             plot.axes = { contour(temp, add=T)
                axis(1, seq(450000 , 580000, by = 10000))
                axis(2, seq(2800000,4000000, by = 10000)) },
             key.title = title(main="Probability") ,
             key.axes = axis(4, seq(0 , 1 , by = 0.1))


The routine will work in a modified form without adding the coordinates (the axis lines) but when I include these the routine produces various errors, such as "dimension mismatch", or "unexpected end encountered."

I tried to follow the example on filled.contour help page.

Thanks in advance

   Steve

Steve Friedman Ph. D.
Ecologist / Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor)
Homestead, Florida 33034

Steve_Friedman_at_nps.gov
Office (305) 224 - 4282
Fax (305) 224 - 4147



R-help_at_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. Received on Fri 01 Apr 2011 - 18:51:45 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Fri 01 Apr 2011 - 21:50:26 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive