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

Uzzal uzzal at gist.ac.kr
Mon Aug 31 16:04:55 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.      a
    a
    hours
    hours 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. How can I plot all the variograms (total 161)  with different title?  Orpheus    
 




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


More information about the R-sig-Geo mailing list