[R-sig-Geo] How to plot multiple variogram with different title in loop?

Uzzal uzzal at gist.ac.kr
Thu Sep 3 10:26:57 CEST 2015


I have a csv file contains hourly PM10 concentration of 1 march to 7 march. Please, download from here. I have already plotted all the hourly variogram (total 161) in a loop by automap package. Here is the coad, this code work properly     library(sp)
    library(gstat)
    library(rgdal)
    library(automap)
    library(latticeExtra)
    
    seoul1to7
    seoul1to7[seoul1to7==0] 
    seoul1to7 
    seoul1to7_split
    seq(seoul1to7_split)
    
        vars
    {
      dat
      coordinates(dat)
      proj4string(dat) 
      dat 
      variogram
      plot
      
      return(plot)
    })
    vars[[1]]
    vars[[2]] 
Here, I can get the individual plot by `vars[[1]], vars[[2]]...`etc command but all the variogram has a same title. Now, I want to get all the variogram  with different title in a loop. I want my variogram title be like, "Variogram for 2012-03-01 1.00",  "Variogram for 2012-03-01 2.00" .....etc. To do this I wrote this code but It doesn't work!     a
    a
    hours
    hours vars_hourlyvars_hourly[[1]]  I keep my desired different title in `hours` variable like "2012-03-01 01:00", "2012-03-01 02:00", "2012-03-01 03:00" ...etc. could you please help me to plot all the variograms (total 161)  with different title? I am badly struck with this problem. Thanks. Orpheus 




-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20150903/35015035/attachment.html>


More information about the R-sig-Geo mailing list