Странное поведение фильтрации / сопоставления в data.frame в R

Question

У меня есть data.frame с 31 столбцом и 31 наблюдением.Все данные в df являются числовыми, с положительными и отрицательными значениями.

Из первых 13 столбцов моего фрейма данных я хочу извлечь индекс положения наблюдения с максимальным абсолютным значением.

Мой код ведет себя странно, и я не могу понять, почему.

Имя базы данных: loads.P_Y_wolf_df_2012_pca

> head(loads.P_Y_wolf_df_2012_pca)


               PC1          PC2         PC3          PC4         PC5         PC6
Mol_17  0.20956737 -0.148512375  0.24587841 -0.269043967  0.18686684 -0.17759436
Mol_18  0.16206391 -0.137572327  0.37520832 -0.243564027  0.18204151 -0.11839496
Mol_20  0.15833954  0.009470871 -0.04321647  0.094148199  0.36981320 -0.34281987
Mol_23 -0.07904978  0.034748195 -0.13504722 -0.271102820  0.03052371  0.15275525
Mol_30 -0.14805422 -0.333887605 -0.05379184  0.003507335  0.19356645  0.01928621
Mol_33 -0.22142249  0.064829612 -0.16049180 -0.278572239 -0.16350937 -0.16854468

Я копирую только первые шесть переменных здесь вместе с первыми 5 наблюдениями.

Ниже приводится цикл, который я пытаюсь, и относительная ошибка.

> indexvec <- vector()
> for (i in 1:13) {
    indexvec[i] <- which(loads.P_Y_wolf_df_2012_pca[, i] ==     max(abs(loads.P_Y_wolf_df_2012_pca[ , i])))
  }

Error in indexvec[i] <- which(loads.P_Y_wolf_df_2012_pca[, i] == max(abs(loads.P_Y_wolf_df_2012_pca[,  : 
  replacement has length zero

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

> which(loads.P_Y_wolf_df_2012_pca[, 1] == max(abs(loads.P_Y_wolf_df_2012_pca[ , 1])))
[1] 15

Итак, следующим шагом было изменить что-то в цикле for и напечатать вывод, и здесь это странно:

> for (i in 1:13) {
    print(which(loads.P_Y_wolf_df_2012_pca[, i] == max(abs(loads.P_Y_wolf_df_2012_pca[ , i]))))
}

[1] 15
[1] 28
[1] 22
[1] 21
[1] 3
integer(0)
integer(0)
integer(0)
[1] 20
integer(0)
[1] 13
[1] 18
integer(0)

Кажется, что код работает на некоторых столбцах, но не на других!И это странно, потому что вот значения из столбца 6.

> loads.P_Y_wolf_df_2012_pca[, 6]
 [1] -0.177594359 -0.118394964 -0.342819874  0.152755248  0.019286210 -0.168544683
 [7] -0.232132045 -0.037674359 -0.164510267 -0.123064519 -0.427363710  0.249120454
[13]  0.189818509 -0.125318890 -0.090817036  0.132368681  0.137849210 -0.179389033
[19] -0.251168502  0.284717837 -0.203070939 -0.025723435  0.161473077  0.072130066
[25]  0.001923286  0.028732683 -0.100890265 -0.091261136  0.216944725  0.085688749
[31]  0.151103335

И если я просто ищу значение, я получаю соответствующий ответ:

> max(abs(loads.P_Y_wolf_df_2012_pca[ , 6]))
[1] 0.4273637

Любая помощь, чтобы позволить мнепонять, в чем проблема, будет высоко ценится.

Answer 1

Мы можем использовать apply, чтобы избежать цикла следующим образом:

apply(df,2,function(x) which.max(abs(x)))

Если мы хотим использовать цикл (не рекомендуется в большинстве случаев по вычислительным соображениям):

 res<-vector()
   for(i in 1:ncol(df)){
  res[i]<-which.max(abs(df[,i]))
}
res

Вариант для цикла:

for(i in 1:ncol(df)){
  res[i]<-which(abs(df[,i])==max(abs(df[,i])))
}
res

С sapply:

sapply(df,function(x) which.max(abs(x)))

Как подсказывает @akrun, мы также можем использовать max.col

Результаты:apply (более информативно):

PC1 PC2 PC3 PC4 PC5 PC6 
  6   5   2   6   3   3 

Явный цикл:

[1] 6 5 2 6 3 3

С max.col:

max.col(t(abs(df)), 'first')
[1] 6 5 2 6 3 3

С sapply:

PC1 PC2 PC3 PC4 PC5 PC6 
  6   5   2   6   3   3 

С purrr:

purrr::map_dbl(df,function(x) which.max(abs(x)))
PC1 PC2 PC3 PC4 PC5 PC6 
  6   5   2   6   3   3