R: L oop для фильтрации кадра данных и вывода среднего значения столбца отфильтрованных данных в новый кадр данных - PullRequest
Новая
       30

R: L oop для фильтрации кадра данных и вывода среднего значения столбца отфильтрованных данных в новый кадр данных

0 голосов
/ 22 апреля 2020

У меня есть этот фрейм данных с именем PlantGroNoT.RCP2.first.0 (случайная выборка его с помощью dput ()). 

structure(list(Year = c(2040L, 2068L, 2096L, 2049L, 2072L, 2042L, 
2087L, 2047L, 2076L, 2075L, 2052L, 2054L, 2034L, 2060L, 2041L, 
2031L, 2067L, 2080L, 2055L, 2067L, 2095L, 2097L, 2097L, 2077L, 
2057L, 2038L, 2052L, 2095L, 2043L, 2075L), DOY = c(78L, 13L, 
20L, 364L, 23L, 14L, 352L, 72L, 54L, 2L, 1L, 5L, 53L, 43L, 63L, 
26L, 21L, 30L, 360L, 22L, 46L, 62L, 355L, 4L, 61L, 355L, 75L, 
51L, 46L, 39L), DAP = c(93L, 28L, 35L, 14L, 38L, 29L, 2L, 87L, 
69L, 17L, 16L, 20L, 68L, 58L, 79L, 41L, 36L, 45L, 10L, 37L, 61L, 
78L, 5L, 20L, 77L, 5L, 90L, 66L, 61L, 54L), NSTD = c(0.244, 0, 
0.023, 0, 0, 0, 0, 0.547, 0, 0, 0, 0, 0, 0, 0.148, 0, 0, 0.001, 
0, 0, 0.186, 0.443, 0, 0, 0.278, 0, 0.403, 0.16, 0.085, 0.069
), RCP = c("2", "4", "4", "4", "2", "4", "8", "4", "8", "4", 
"4", "4", "8", "4", "4", "2", "8", "8", "8", "4", "8", "2", "4", 
"8", "4", "8", "4", "4", "8", "2")), row.names = c(5634575L, 
1218676L, 788139L, 7051284L, 6262001L, 1122307L, 3560772L, 3925983L, 
2584130L, 3401700L, 1217186L, 1455727L, 2993448L, 597952L, 5893799L, 
6696664L, 2345802L, 3105731L, 6347255L, 5005396L, 3330449L, 780805L, 
2522272L, 3668256L, 2099405L, 818891L, 1123257L, 5237622L, 1259779L, 
1623133L), class = "data.frame")

Я хочу отфильтровать его по RCP, DAP и Year. Из результатов этого фильтра я хотел бы взять среднее значение NSTD и вывести его в столбец NSTDmean кадра данных HO1_2020_RCP2_stress

DAP в диапазоне от 0 до 114, поэтому в первая ячейка HO1_2020_RCP2_stress$WSPDmean Я хотел бы иметь значение NSTD для DAP == 0, во второй ячейке среднее значение для DAP == 1 и т. д. c.

Я создал этот l oop, чтобы сделать работу. 

i = -1
  repeat{
    i= i+1
    PlantGroNoT.RCP2.first.0 <- filter(PlantGroNoT_1, RCP==2 & DAP==i & Year<=2060)
    H01_2020_RCP2_stress$NSTDmean <- mean(PlantGroNoT.RCP2.first.0$NSTD, na.rm=TRUE)
    if (i>114)
      break
  }

Сам фильтр и вычисление среднего значения работают.

Однако l oop неправильно возвращает NSTDmean в столбце HO1_2020_RCP2_stress$WSPDmean. Он всегда печатает значение для последнего i (т.е. здесь DAS == 114) для всего столбца NSTDmean.

1 Ответ

1 голос
/ 22 апреля 2020

Надеюсь, я не пойму неправильно, что вы хотите сделать, но я думаю, что использование group_by является правильным подходом: 

df %>%
  filter(RCP == 2 & Year <= 2060) %>% # These are the only conditions and they are fixed
  group_by(DAP) %>%
  summarise(NTSDT = mean(NSTD,na.rm =TRUE)) %>%
  arrange(DAP)

Это отфильтрует фрейм данных только по строкам с RCP == 2 и год <= 2060. Затем он суммирует среднее значение NSTD для всех значений DAP. </p>

Если вам нужно добавить этот результат в исходный фрейм данных, вы можете просто присоединиться к нему или использовать mutate вместо обобщать.

...