R help - Apr 2011 - filled contour plot with contour lines

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

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_Friedman at 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_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.