ODE模型随机参数多重复实验绘图异常问题求助

1. 先校准ODE求解逻辑 #

曲线反向回起点大概率是求解逻辑的问题，先确保你的ODE函数和时间序列设置正确。以下是示例ODE函数（替换成你实际的模型）：

# 替换为你的实际ODE方程
ode_model <- function(t, y, params) {
  At <- y[1]
  mu <- params$mu
  dAt <- -mu * At  # 根据你的模型调整符号或表达式
  list(dAt)
}

注意时间序列必须递增且无重复起始点，比如times <- seq(0, 100, by = 1)，不要出现c(0,100,0)这类错误序列。

2. 用replicate批量生成随机mu的模拟结果 #

replicate可以自动重复实验，每次生成0-0.04范围内的随机mu，代码如下：

library(deSolve)  # 假设你用deSolve包求解ODE

# 基础参数设置
n_reps <- 50  # 重复实验次数，按需调整
times <- seq(0, 100, by = 1)  # 你的时间范围
y0 <- c(At = 100)  # At的初始值，替换成你的实际初始条件

# 批量生成模拟结果
sim_results <- replicate(n_reps, {
  # 每次生成随机mu
  mu <- runif(1, min = 0, max = 0.04)
  # 求解ODE并转成数据框
  out <- ode(y = y0, times = times, func = ode_model, params = list(mu = mu))
  as.data.frame(out)
}, simplify = FALSE)  # 必须设为FALSE，保留每个实验的独立数据框

3. 绘制扇形分布曲线 #

提供两种绘图方案，都能实现扇形发散效果：

方案1：ggplot2绘图（推荐） #

library(ggplot2)

# 合并所有实验数据
combined_data <- do.call(rbind, sim_results)
combined_data$rep_id <- rep(1:n_reps, each = length(times))  # 给每组实验加标识

ggplot(combined_data, aes(x = time, y = At, group = rep_id)) +
  geom_line(color = "blue", alpha = 0.5) +  # alpha设置透明度，避免曲线重叠
  theme_minimal() +
  labs(x = "时间", y = "At值")

方案2：基础绘图 #

# 先绘制空图，提前设置y轴范围避免曲线被截断
plot(times, y0["At"], type = "n", xlab = "时间", ylab = "At值", 
     ylim = range(sapply(sim_results, function(x) x$At)))
# 循环绘制每条曲线
for (i in 1:n_reps) {
  lines(sim_results[[i]]$time, sim_results[[i]]$At, col = rgb(0, 0, 1, 0.5))
}

4. 异常排查 #

如果曲线仍反向回起点，检查这几点：

• ODE方程的符号是否正确：根据你的模型需求，调整dAt的表达式（比如是递增还是递减）
• 时间序列是否严格递增：确保times是从小到大排列的
• 初始值是否被错误重置：检查求解过程中有没有意外修改初始条件

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