比較嵌套循環

Question

我有一個面板，看起來像

df <-read.table(text=" 
preis10_6 akt_datum10_6 preis11_6 akt_datum11_6 
1.55  10.06.2014 1.45  10.06.2014 
1.56  09.06.2014 1.49  11.06.2014 
",header=TRUE,sep="") 

在奇數列（preise）包含收集在那一天，而偶數列（akt_datum）包含信息的價格與列標題進入細胞關於數據的話題性。因此，當df[2, 2]是09.06.2014這意味着df[1:2, 2]中的信息是從前一天開始的。我想設置這些案例NA。

這是我迄今所做的：

# Instrall stringr for the function str_sub 
require(stringr) || install.packages("stringr") 

# Get indices for the columns with topicality information 
spalten <- seq(2, length(df), 2) 

# Loop over these columns 
for (spalte in spalten) { 

    # Construct the benchmark date from the column name 
    splitter <- str_sub(names(df)[spalte], 10, -1) 
    splitter <- strsplit(splitter, "_") 
    # Account for the case where the column name is in short time format (no trailing 0) 
    splitter[[1]][1] <- ifelse(nchar(splitter[[1]][1])==1, 
          paste0("0", splitter[[1]][1]), 
          splitter[[1]][1] 
) 
    splitter[[1]][2] <- ifelse(nchar(splitter[[1]][2])==1, 
          paste0("0", splitter[[1]][2]), 
          splitter[[1]][2] 
) 
    date <- paste(splitter[[1]][1], splitter[[1]][2], "2014", sep=".") 

    # Loop over all rows in the actual column 
    for (zeile in 1:nrow(df)) { 
    # and set the cell and the one before equal to NA 
    ifelse(df[zeile, spalte]!=date, df[zeile, spalte] <- df[zeile, spalte-1] <- NA, "") 
    } 

} 

這工作，但需要年齡，因爲我有事實933x324面板。也就是說，我每天有933個價格和時事信息，總共162天，這使得162個價格熱門對= 324列。

我該如何使此過程更快？

Answer 1

這是做你想做的嗎？至少它可以產生與你的代碼相同的輸出（以換算的日期爲模）。

# Read data example 
df <-read.table(text=" 
preis10_6 akt_datum10_6 preis11_6 akt_datum11_6 
1.55  10.06.2014 1.45  10.06.2014 
1.56  09.06.2014 1.49  11.06.2014 
",header=TRUE,sep="") 

# Convert to date (to allow for comparison between dates) 
df$akt_datum10_6 <- as.Date(df$akt_datum10_6, "%d.%m.%Y") 
df$akt_datum11_6 <- as.Date(df$akt_datum11_6, "%d.%m.%Y") 

# Check which date is first, and substitute `NA`s 
first <- df$akt_datum10_6 < df$akt_datum11_6 
df[first, 1:2] <- NA 
df[!first, 3:4] <- NA 
print(df) 
#preis10_6 akt_datum10_6 preis11_6 akt_datum11_6 
#1  1.55 2014-06-10  NA   <NA> 
#2  NA   <NA>  1.49 2014-06-11