R-如何在循環中保存循環結果？

Question

好的，所以目標是通過一個循環讀取包含特定於氣象站的氣象信息的csv文件，該循環提供開始日期，結束日期以及干球溫度和濕度的百分數以及站號。 到目前為止，我已經能夠制作一個循環，為我提供一個csv文件的信息以及6-20個電台的數據，因此效果很好。 現在，我試圖在第一個循環中設置第一個循環，該循環從一個文件夾運行多頭csv文件，但是我只能得到工作站ID，而不能得到其他四個結果。 我想知道如何以列表，單獨的向量或數據框的形式保存所有這些數據。

這是代碼。

files <- list.files("R_Test_Nicholas", pattern="*.csv", full.names=TRUE)
l <- length(files)
r <- NULL
first <- NULL
for(k in 1:l) {
  first <- read.csv(files[k])
  library(openair)
  station <-unique(first$STN_ID)
  n = length(station)
  first$date <- paste(first$LOCAL_YEAR, first$LOCAL_MONTH, first$LOCAL_DAY, first$LOCAL_HOUR, sep=" ")
  first$date <- as.POSIXct(first$date, format="%Y %m %d %H", "UTC")

  result7 <- NULL
  for (i in 1:n) {
    curr_station <- subset(first, STN_ID==station[i], select=c(date, DRY_BULB_TEMP, HUMIDEX, LOCAL_YEAR))
    test <- cutData(curr_station, type = "season")
    timers <- selectByDate(test, start = "YYYY-mm-dd", year = 1971:2000)
    trial <- subset(timers, season == "summer (JJA)")
    Ptile <- quantile(trial$DRY_BULB_TEMP, c(0.95), na.rm=T)
    Htile <- quantile(trial$HUMIDEX, c(0.95), na.rm=T)
    date_start <- min(trial$LOCAL_YEAR, na.rm=T)
    date_last <- max(trial$LOCAL_YEAR, na.rm=T)
    result7$stn[i] <- station[i]
    result7$Ptile[i] <- Ptile
    result7$Htile[i] <- Htile
    result7$date_start[i] <- date_start
    result7$date_last[i] <- date_last
  }
  r$STN_ID[k] <- result7
  r$Ptile[k] <- result7$Ptile[i]
  r$Htile[k] <- result7$Htile[i]
  r$date_start[k] <- result7$date_start[i]
  r$date_last[k] <-result7$date_last[i]
}
r
$STN_ID
$STN_ID[[1]]
[1] 6130 6137 6140 6157 6205 6207 6250 6256 6915 6916 6918 6919 7026 7558 9821

$STN_ID[[2]]
[1] 10808 10981 26823 26968 27746 29493 30309 30598 42103 43323 43383 45090 48568 48768 50309 50310 51537

$STN_ID[[3]]
[1] 6312 6330 6339 6345 6354 6356 6358 6369 6399 6442 6454 6465 6468 6486 6491 6501 6516 6923 7103 7162 7169 7173 8990 9033 9833

$STN_ID[[4]]
[1] 10078 10661 10792 10848 10859 10914 10936 10945 10969 10970 26824 26864 27141 27223 27592 27600 27868 30326 30668 31829 41575 42083 42243 43123 43124 43183 43403 43404 43405 43406 44363
[32] 44503 46007 47187 48668 49748 50133 50408 50620

$STN_ID[[5]]
[1] 6526 6547 7177

$STN_ID[[6]]
[1] 10800 10814 27846 30308 31029 41903 50621


$Ptile
[1] 23.9   NA 19.5   NA 21.2   NA

$Htile
[1] 29.35    NA 24.34    NA 26.96    NA

$date_start
[1] 1986  Inf 1986  Inf 1978  Inf

$date_last
[1] 2000 -Inf 1996 -Inf 2000 -Inf

Answer 1

您可以使用[[運算符來創建和訪問list元素。 由於缺少可復制的數據，因此將mtcars數據集用於此演示。

data(mtcars)
NCol = ncol(mtcars)
ObjList=list()

for(i in 1:NCol){

ObjList[[i]]=list()
ObjList[[i]]['max']=max(mtcars[,i])
ObjList[[i]]['min']=min(mtcars[,i])
ObjList[[i]]['mean']=mean(mtcars[,i])
ObjList[[i]]['sd']=sd(mtcars[,i])
ObjList[[i]][["head"]]=head(mtcars)
ObjList[[i]][["tail"]]=tail(mtcars)


}

max_vec=as.vector(do.call(cbind,lapply(ObjList,function(x) x[["max"]])))
#[1]  33.900   8.000 472.000 335.000   4.930   5.424  22.900   1.000   1.000   5.000   8.000


max_orig=as.vector(apply(mtcars,2,max))
#[1]  33.900   8.000 472.000 335.000   4.930   5.424  22.900   1.000   1.000   5.000   8.000

setdiff(max_vec,max_orig)
#numeric(0)

Answer 2

沒有簡短的數據集很難重現代碼並理解您想要做什么，但是我看到您應該使用矩陣來收集數據

result7 <-matrix（ncol =？，nrow =？））＃而不是NULL

for（i in 1：n）{...（其余的代碼，然后得到結果時，將它們收集：） result7 [i，1] <-結果的計算結果將在第i行，第一列result7 [i，2] <-計算結果B將位於i行和第二列}

您可以將其轉換為數據框，然后為各列命名：

result7 <-data.frame（result7）
colnames（georef.data）<-c（“此處的列名）

當您完成所有數據幀時，只需使用rbind按行連接數據幀即可，這些行可以由ID定義。

Answer 3

r$STN_ID[k] <- result7
  r$Ptile[k] <- result7[2]
  r$Htile[k] <- result7[3]
  r$date_start[k] <- result7[4]
  r$date_last[k] <-result7[5]

這就是我最終所做的，這使我在列表中得到了一堆列表，然后將其更改為數據幀。 感謝那些花時間回答的人。

Answer 4

這是一個關於在循環中將新數據連接到現有數據幀的常見問題的非常具體的示例。 為了完整起見，我想分享一種通用的方法。

這個想法是從一個空的數據幀開始，在一個循環中進行計算，然后使用rbind()添加數據。 所以...

df <- NULL
for (n in seq(1,10)){

df <- rbind(df,
            data.frame(n = n,
                       n.squared = n)
}

在此示例中，您將獲得一個10行數據幀，其中包含列n和n.squared 。