Question

最好通过首先描述我正在使用的系统和数据来解释我的编码问题。然后我会提出我的问题/问我的问题。

我正在尝试模拟种群中单个昆虫的生长和发育。发展很大程度上受温度驱动。因此，显影是以“累积热量”来衡量的，即，更多的热量等于更多的显影。但是，存在一个温度，在该温度以下不会发生显影（“基准”）。

每种昆虫在发育过程中都会经历多个阶段，每个阶段都有一个独特的基础温度。

最后，从一个阶段前进到下一个阶段所需的热量因人而异。

好的，一些示例数据：

# df1: Hourly temperatures
df1 = data.frame(DateTime = seq(
  from = as.POSIXct("1986-1-1 0:00"),
  to = as.POSIXct("1986-1-8 23:00"), 
  by = "hour"))

temp = c(5,5,5,6,7,8,9,10,13,17,22,25,26,28,26,25,25,22,19,14,10,8,5,5)

df1$temp <- temp 

# df2: Each row is an individual insect. 
#      s1_thresh is the number of degrees that need to accumulate for each 
#      individual to advance from stage 1 to stage 2. 

df2 <- data.frame(id = c("A", "B", "C", "D", "E"), 
              s1_thresh = c(21.5, 25.1, 19.9, 20.4, 21.4))

# Stage-specific base temperatures below which no degrees accumulate

base_s1 <- 10.5  # base temp for stage 1
base_s2 <- 8.6   # base temp for stage 2

# Temperature accumulation above base_s1 
df1$dd_s1 <- ifelse(temp >= base_s1, (temp - base_s1)/24, 0)
df1$cumdd_s1 <- cumsum(df1$dd_s1)

这是我的问题：由于热量需求的不均匀性，每个人都将在独立的时间进行过渡/提前阶段，因此当发生这种转变时，我该如何改变每个人的基本温度？这是df1中单个“ A”的理想结果（某种）。

# Example for single individual:
# Individual "A" has s1_thresh of 21.5, so a shift to base_2 occurs on 1986-01-04 16:00:00, row 89
df1$dd_s2 <-  ifelse(df1$cumdd_s1 > df2$s1_thresh[df2$id == "A"] & temp >= base_s2, (temp - base_s2)/24, 0) 
df1$cumdd_s2 <- cumsum(df1$dd_s2)

我试图避免为每个人设置多个温度累积列，但是了解在有限的时间范围内每个人的累积热量很重要。

非常感谢您的光临！

Answer 1

假设您只想为每个人获得一个累积2度日的向量：

# transformed temperatures series relative to stages
df1$temps_b1 <- pmax(df1$temp - base_s1, 0)
df1$temps_b2 <- pmax(df1$temp - base_s2, 0)

ddsum <- function(theid) {
  in_stage1 <- cumsum(df1$temps_b1/24) < df2$s1_thresh[df2$id==theid]
  s2_temp <- ifelse(in_stage1, 0, df1$temps_b2)
  return(cumsum(s2_temp/24))
}

res<-as.data.frame(sapply(df2$id, FUN=ddsum))
names(res) <- df2$id
res

R中基于个人的仿真

1 个答案: