R перекрывает несколько участков в цикле

Question

Итак, я создал цикл, который составляет 10 отдельных графиков:

for (k in 1:nrow(sites)) {
  temp_title <- paste("site",k, "county", sites[k,2],"site",sites[k,3])
  l <- which(hourly_nj_table$County.Code==sites[k,2]&hourly_nj_table$Site.Num==sites[k,3])#grab data for each site individually
  temp_filename <- paste("/Users/bob111higgins/Documents/School/College/Rutgers/Atmospheric Research",temp_title,".pdf")
  PM_site <- hourly_nj_table[l,]
  PM_site$realTime <- as.numeric(PM_site$Time.Local)
  PM_mean_site <- aggregate(PM_site, by=list(PM_site$Time.Local),FUN="mean",na.rm=TRUE)  #Make it average by time of day so can make time series plots. 
  plot(PM_mean_site$realTime,PM_mean_site$Sample.Measurement, type="l",lwd=10,main=paste(temp_title),xlab="LocalTime",ylab="Ozone (ppm)")#,ylim=c(0,0.05))
}

Но хотелось бы посмотреть, как они сравниваются по одной оси. Обычно (если я просто жестко его кодирую), я бы добавил новый параметр, а затем создал следующий график, но я не уверен, как включить это в цикл.

Все данные поступают из одного CSV-файла, если это помогает..

Спасибо!

Answer 1

Ты действительно очень близко. Plot() получает мяч катится, lines() позволит вам нарисовать внутри сюжета:

for (k in 1:nrow(sites)) {
  temp_title <- paste("site",k, "county", sites[k,2],"site",sites[k,3])
  l <- which(hourly_nj_table$County.Code==sites[k,2]&hourly_nj_table$Site.Num==sites[k,3])#grab data for each site individually
  temp_filename <- paste("/Users/bob111higgins/Documents/School/College/Rutgers/Atmospheric Research",temp_title,".pdf")
  PM_site <- hourly_nj_table[l,]
  PM_site$realTime <- as.numeric(PM_site$Time.Local)
  PM_mean_site <- aggregate(PM_site, by=list(PM_site$Time.Local),FUN="mean",na.rm=TRUE)  #Make it average by time of day so can make time series plots. 
  ifesle(k ==1 , 
  plot(PM_mean_site$realTime,PM_mean_site$Sample.Measurement, type="l",lwd=10,main=paste(temp_title),xlab="LocalTime",ylab="Ozone (ppm)")#,ylim=c(0,0.05)),
  lines(PM_mean_site$realTime,PM_mean_site$Sample.Measurement, lwd=10))
}

Я уверен, что есть лучшие способы сделать это, но так я делал это в прошлом.