R：根据日期计算平均行。 tidyverse

Question

你们中的一些人通过基数R解决了这个问题，但是我现在正尝试使用整洁的数据来解决这个问题。这是我的数据集：

data <- tibble( "DATE_FIRE"= c("1989-07-31", "1989-07-31", "1989-07-31", "1989-07-31","1989-07-31","1989-08-31", "1989-08-31", "1989-08-31", "1989-08-31","1989-08-31"), 
       "FID_FIRE" = c(1,1,1,1,1,2,2,2,2,2),
       "date" = c("1988-01-01", "1989-01-01","1990-01-01","1991-06-07","1992-01-01", "1988-01-01", "1989-01-01","1990-01-01","1991-06-07","1992-01-01"),
       "NDVI" = c( 0.9, 0.8, 0.1, 0.2, 0.3, 0.8, 0.85, 0.15, 0.30, 0.50))
> data
# A tibble: 10 x 4
   DATE_FIRE  FID_FIRE date        NDVI
   <chr>         <dbl> <chr>      <dbl>
 1 1989-07-31        1 1988-01-01  0.9 
 2 1989-07-31        1 1989-01-01  0.8 
 3 1989-07-31        1 1990-01-01  0.1 
 4 1989-07-31        1 1991-06-07  0.2 
 5 1989-07-31        1 1992-01-01  0.3 
 6 1989-08-31        2 1988-01-01  0.8 
 7 198`enter code here`9-08-31        2 1989-01-01  0.85
 8 1989-08-31        2 1990-01-01  0.15
 9 1989-08-31        2 1991-06-07  0.3 
10 1989-08-31        2 1992-01-01  0.5 

如您所见，我有2个林地的日期，火灾1和火灾2（FID_FIRE列）。 DATE_FIRE告诉我们每场火灾发生的时间，在火灾发生之前和之后，我也有关于这些地块的NDVI（绿色参数）数据（请参阅date和NDVI列）。

我想为每幅图添加一列计算NDVI的均值 之前发生火灾。在FID_FIRE=1的情况下，它将是NDVI和row 1的2值的平均值。

输出应类似于：

> desired_output
# A tibble: 10 x 5
   DATE_FIRE  FID_FIRE date        NDVI meanPrefire
   <chr>         <dbl> <chr>      <dbl>       <dbl>
 1 1989-07-31        1 1988-01-01  0.9        0.85 
 2 1989-07-31        1 1989-01-01  0.8        0.85 
 3 1989-07-31        1 1990-01-01  0.1        0.85 
 4 1989-07-31        1 1991-06-07  0.2        0.85 
 5 1989-07-31        1 1992-01-01  0.3        0.85 
 6 1989-08-31        2 1988-01-01  0.8        0.825
 7 1989-08-31        2 1989-01-01  0.85       0.825
 8 1989-08-31        2 1990-01-01  0.15       0.825
 9 1989-08-31        2 1991-06-07  0.3        0.825
10 1989-08-31        2 1992-01-01  0.5        0.825

Answer 1

尝试：

library(dplyr)

data %>%
  group_by(FID_FIRE) %>%
  mutate(meanPrefire = mean(NDVI[date < DATE_FIRE], na.rm = TRUE))

输出：

# A tibble: 10 x 5
# Groups:   FID_FIRE [2]
   DATE_FIRE  FID_FIRE date        NDVI meanPrefire
   <chr>         <dbl> <chr>      <dbl>       <dbl>
 1 1989-07-31        1 1988-01-01  0.9        0.85 
 2 1989-07-31        1 1989-01-01  0.8        0.85 
 3 1989-07-31        1 1990-01-01  0.1        0.85 
 4 1989-07-31        1 1991-06-07  0.2        0.85 
 5 1989-07-31        1 1992-01-01  0.3        0.85 
 6 1989-08-31        2 1988-01-01  0.8        0.825
 7 1989-08-31        2 1989-01-01  0.85       0.825
 8 1989-08-31        2 1990-01-01  0.15       0.825
 9 1989-08-31        2 1991-06-07  0.3        0.825
10 1989-08-31        2 1992-01-01  0.5        0.825