[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