如何按Session与Trial Index修改DataFrame中label=d的time值？

Got it, let's work through this problem step by step. You need to update the time value for the "d" label in each group (grouped by session and trial_index) by subtracting the total of time values for "a", "b", and "c" in the same group—while leaving "a", "b", "c" completely untouched. Here are two reliable, easy-to-follow approaches:

1. Using dplyr (Tidyverse Approach) #

This is the most readable and intuitive method for grouped operations in R. First, make sure you have the dplyr package installed (if not, run install.packages("dplyr") first).

library(dplyr)

# Your original dataset
session <- c(rep(34,8), rep(28,8))
trial_index <- c(rep(2,4),rep(5,4),rep(6,4),rep(8,4))
label <- c(rep(c("a","b","c","d"),4))
time <- c(10,2,7,40,4,3,6,20,5,3,5,15,4,2,3,17)
data <- data.frame(session, trial_index,label,time)

# Modify the time values for "d" labels
updated_data <- data %>%
  group_by(session, trial_index) %>%
  mutate(
    time = ifelse(label == "d", 
                  time - sum(time[label %in% c("a", "b", "c")]),
                  time)
  ) %>%
  ungroup() # Optional: clears grouping if you don't need it for further work

# Check the result for your example group (session 34, trial_index 2)
updated_data %>% filter(session == 34, trial_index == 2)

Breakdown of this method: #

• group_by(session, trial_index): Groups the data so all calculations stay within each unique session-trial pair.
• mutate(): Updates the time column conditionally. For rows where label is "d", it subtracts the sum of time values from "a", "b", "c" in the same group. All other rows keep their original time value.
• ungroup(): Removes the grouping structure (optional, but good practice if you're done with grouped operations).

2. Base R Approach (No External Packages) #

If you prefer sticking to base R without installing extra packages, you can use the ave() function to handle grouped calculations:

# Your original dataset (same as above)
session <- c(rep(34,8), rep(28,8))
trial_index <- c(rep(2,4),rep(5,4),rep(6,4),rep(8,4))
label <- c(rep(c("a","b","c","d"),4))
time <- c(10,2,7,40,4,3,6,20,5,3,5,15,4,2,3,17)
data <- data.frame(session, trial_index,label,time)

# Calculate the sum of a/b/c times for each group (repeats the sum for every row in the group)
group_sum <- ave(data$time, 
                 list(data$session, data$trial_index), 
                 FUN = function(x) sum(x[data$label[seq_along(x)] %in% c("a","b","c")]))

# Update only the "d" label time values
data$time[data$label == "d"] <- data$time[data$label == "d"] - group_sum[data$label == "d"]

Breakdown of this method: #

• ave() computes the sum of "a"/"b"/"c" time values for each group, then repeats that sum for every row in the group so we can match it to the "d" rows.
• We then subset the time column to only "d" rows and subtract the corresponding group sum values.

Quick Verification #

For your example group (session 34, trial_index 2):

• Original "d" time: 40
• Sum of "a"+"b"+"c": 10+2+7=19
• Updated "d" time: 40-19=21, which you'll see in the output of either method.

内容的提问来源于stack exchange，提问作者unomas83