如何使用for循环计算蛙类个体每9分钟间隔的水分流失速率？

你现有代码的计算逻辑是取单个个体0分钟和45分钟的重量差除以总时长，得到全实验周期的平均失水率，之后将该值赋值给该个体的所有行，所以才会出现同一个个体所有区间速率相同的情况，和你需要的单9分钟间隔速率计算逻辑不符。

推荐解法（dplyr包，代码简洁不易出错） #

按个体分组后取上一行重量计算单间隔速率即可，自带的lag()函数可以直接取同分组内前一行的数值：

library(dplyr)
EWL <- EWL %>%
  arrange(Animal, Time) %>% # 先按个体+时间排序，避免原始数据顺序混乱导致计算错误
  group_by(Animal) %>%
  mutate(interval_WR = (Weight - lag(Weight))/9) %>% # 每个间隔时长为9分钟
  ungroup()

计算结果中每个个体Time=0的行没有前序数据，速率会返回NA，其余行对应当前9分钟区间的失水速率，符合你的需求。

基础包for循环解法 #

如果不想依赖第三方包，可以用嵌套for循环实现，注意去掉你原代码中完全没用的外层i循环：

# 先初始化间隔速率列
EWL$interval_WR <- NA
# 遍历所有个体
for (j in unique(EWL$Animal)) {
  # 提取当前个体的所有行索引
  idx <- which(EWL$Animal == j)
  # 按时间升序排列索引，避免原始数据顺序混乱
  idx_sort <- idx[order(EWL$Time[idx])]
  # 从第二个时间点开始计算每个间隔的速率
  for (k in 2:length(idx_sort)) {
    EWL$interval_WR[idx_sort[k]] <- (EWL$Weight[idx_sort[k]] - EWL$Weight[idx_sort[k-1]])/9
  }
}

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