2

私は、3 月 1 日から 3 月 7 日までの 1 時間ごとの PM10 濃度を含む csv ファイルを持っています。こちらからダウンロードしてください。automap パッケージにより、ループ内のすべてのバリオグラム (合計 161) をプロットしました。

library(sp)
library(gstat)
library(rgdal)
library(automap)
library(latticeExtra)

seoul1to7<-read.csv("seoul1to7.csv")
seoul1to7[seoul1to7==0] <-NA
seoul1to7 <- na.omit(seoul1to7)
seoul1to7_split<-split(seoul1to7,seoul1to7$time) 
seq(seoul1to7_split)



vars<-lapply(seq(seoul1to7_split), function(i)
{
  dat<-seoul1to7_split[[i]]  
  coordinates(dat)<-~LON+LAT
  proj4string(dat) <- "+proj=longlat +datum=WGS84" 
  dat <- spTransform(dat, CRS("+proj=utm +north +zone=52 +datum=WGS84"))
  variogram<-autofitVariogram(log(PM10)~1,dat, model="Sph")
  plot<- plot(variogram,plotit=FALSE, asp=1)

  return(plot)
})
vars[[1]]
vars[[2]]

ここでは、etc で個々のプロットを取得できますvars[[1]],vars[[2]]...が、すべてのバリオグラムには同じタイトルが付いています。ここで、ループ内の異なるタイトルのすべてのバリオグラム画像をプロットしたいと思います。バリオグラムのタイトルを「2012-03-01 1.00 のバリオグラム」、「2012-03-01 2.00 のバリオグラム」などにします。

a<-as.POSIXct(names(seoul1to7_split), format="%Y%m%d%H")
a
hours<-substr(a,1,16)
hours

hours「2012-03-01 01:00」、「2012-03-01 02:00」、「2012-03-01 03:00」など、目的の別のタイトルを変数に保持します。

ループを使用して、タイトルが異なるすべてのバリオグラム (合計 161) をプロットするにはどうすればよいですか?

4

1 に答える 1