Steve,
I use filled.contour but have a semicolon, ";", between the two
axis calls. Do you need one after
"axis(1, seq(450000 , 580000, by = 10000))"
?
Clint
--
Clint Bowman INTERNET: cl...@ecy.wa.gov
Air Quality Modeler INTERNET: cl...@math.utah.edu
Department of Ecology VOICE: (360) 407-6815
PO Box 47600 FAX: (360) 407-7534
Olympia, WA 98504-7600
USPS: PO Box 47600, Olympia, WA 98504-7600
Parcels: 300 Desmond Drive, Lacey, WA 98503-1274
On Fri, 1 Apr 2011, steve_fried...@nps.gov wrote:
Hi Peter,
Thanks for taking the time to consider the problem. I realize the procedure
is not reproducible. I use a 4 dimensional netCDF file (4.2 GB) in size to
pull data into this process. Nobody in their right mind should work with
such things. It's a spatial temporal database with 10 years of daily data
in a irregular area spanning 405 x 287 cells.
Anyway, I tried your suggestion adding a { in front of the for loop and a
closing } following the last line.
It did not work.
Steve Friedman Ph. D.
Ecologist / Spatial Statistical Analyst
Everglades and Dry Tortugas National Park
950 N Krome Ave (3rd Floor)
Homestead, Florida 33034
steve_fried...@nps.gov
Office (305) 224 - 4282
Fax (305) 224 - 4147
Peter Ehlers
<ehlers@ucalgary.
ca> To
"steve_fried...@nps.gov"
04/01/2011 05:08 <steve_fried...@nps.gov>
PM cc
"r-help@r-project.org"
<r-help@r-project.org>
Subject
Re: [R] filled contour plot with
contour lines
Aren't you missing a set of parentheses?
I can't run your code since it's not reproducible, but to
my aging eyes it seems that you need a set of '{}' around
the contents of your loop:
for(j in 1:(varsize[4]-1)) { loop stuff }
Peter Ehlers
On 2011-04-01 11:48, steve_fried...@nps.gov wrote:
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_fried...@nps.gov
Office (305) 224 - 4282
Fax (305) 224 - 4147
______________________________________________
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.
______________________________________________
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.
______________________________________________
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.
