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)) 
} 

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