变形虫-细菌互作建模：R语言模型构建困境及求解

一、当前模型的核心逻辑问题 #

1. Holling函数误用：你用的(θ * M_P * P_M )/(1 + θ * H_t*P_M)*K4项物理意义模糊。Holling函数通常描述捕食/感染的速率，而这里是胞内细菌杀死宿主（感染变形虫P_M）的过程，应该和单个变形虫内的细菌负荷（M_P/P_M）强相关，而非两者的乘积。
2. 变量转化逻辑混乱：
   • dM中保留了β * P * M的消耗项，但dP_M的感染项和dM_P的增殖项重复关联β * P * M，没有明确摄入的游离菌如何转化为胞内细菌。
   • dM_P的增殖限制项1 - M_P/K4 * P_M错误，K4作为胞内细菌承载量，应该是每个变形虫的承载上限，合理的限制应为1 - M_P/(K4 * P_M)。
3. 释放项不合理：R * M_P * P_M是二次项，不符合“胞内细菌释放为游离菌”的生物逻辑——释放量应与胞内细菌总量M_P直接相关，而非两者乘积。

二、过程拆解与方程修正 #

先明确完整互作链：

未感染变形虫（P）摄入游离细菌（M）→ 转化为感染变形虫（P_M），同时游离菌转化为胞内细菌（M_P）→ 胞内细菌在P_M内增殖→ 胞内负荷过高杀死P_M→ 胞内细菌释放为游离菌（M）

1. 感染与胞内细菌初始形成 #

• 游离菌M因被摄入减少：dM需减去β * P * M
• 未感染变形虫P转化为感染变形虫P_M：dP_M增加β * P * M
• 摄入的游离菌转化为胞内细菌：dM_P增加β * P * M * η（η为转化效率，100%转化则取1）

2. 胞内细菌杀死宿主的过程 #

替换原Holling项，用负荷驱动的死亡率更符合生物逻辑：

• 变形虫死亡率：δ * (M_P / P_M) / (K_δ + M_P / P_M) * P_M
  • δ：变形虫最大死亡率（胞内负荷极高时）
  • K_δ：半饱和常数（死亡率达最大值一半时的胞内负荷）
• 对应变化：
  • dP_M减去该死亡率项（感染变形虫死亡）
  • dM加上δ * (M_P / P_M) / (K_δ + M_P / P_M) * P_M * ε（ε为死亡后释放为游离菌的比例）
  • dM_P减去δ * (M_P / P_M) / (K_δ + M_P / P_M) * P_M * (1 - ε)（未释放的胞内细菌降解/死亡）

如果坚持用Holling型，应调整为针对胞内负荷的攻击率：

# 基于Holling II型的宿主死亡率
kill_rate <- (θ_k * M_P) / (1 + θ_k * H_k * (M_P / P_M))

θ_k为单个胞内细菌的攻击率，H_k为攻击导致死亡的处理时间，避免原方程中无意义的乘积项。

3. 修正后完整方程示例 #

# 游离细菌（M）的变化
dM <- M * g_M * (1 - M / K3) - πM * M - (M / V_trough) * D * NC - β * P * M + 
      # 感染变形虫死亡释放的细菌
      δ * (M_P / P_M) / (K_δ + M_P / P_M) * P_M * ε + 
      # 主动释放的细菌（若存在该过程）
      R * M_P  

# 感染变形虫（P_M）的变化
dP_M <- β * P * M + g_P * P_M * (1 - (P + P_M) / K_2) - πP * P_M - 
        (P_M / V_trough) * D * NC - 
        # 胞内细菌导致的死亡
        δ * (M_P / P_M) / (K_δ + M_P / P_M) * P_M  

# 胞内细菌（M_P）的变化
dM_P <- g_MP * M_P * (1 - M_P / (K4 * P_M)) +  # 基于单变形虫承载量的增殖
        β * P * M * η -  # 摄入游离菌转化为胞内菌
        πMP * M_P + e * M_PC - n * M_P - 
        δ * (M_P / P_M) / (K_δ + M_P / P_M) * P_M * (1 - ε) -  # 未释放的胞内菌损耗
        R * M_P  # 主动释放的胞内菌

三、后续优化建议 #

• 变量重命名：将M/M_P、P/P_M改为更清晰的名称（如M_free/M_intra、A_uninf/A_inf），避免混淆。
• 参数校准：θ、H_t需调整为对应胞内攻击的参数，建议参考变形虫-分枝杆菌互作的已发表文献获取合理范围。
• 敏感性分析：测试δ、K_δ、η等关键参数对模型输出的影响，验证是否符合生物预期（如胞内负荷越高，变形虫死亡率越高）。
• 模块验证：单独模拟胞内细菌杀死宿主的过程，确认当M_P增加时，P_M的下降速率和M的上升速率匹配生物逻辑。

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