在实际传染病防控和公共卫生决策中传统数学模型如SIR模型是理解疾病传播动态的核心工具。然而这些模型的构建、参数拟合和预测分析往往需要深厚的数学和编程背景过程复杂且耗时。如今随着人工智能技术的普及特别是大语言模型和代码生成工具的发展即使是非专业开发者也具备了利用AI辅助完成从数据理解到模型构建、再到代码生成和结果分析的全链路能力。本文将以一场假设的流感爆发数据为例演示如何借助现代AI工具引导我们一步步“跑通”一个基础的传染病动力学建模流程。这个过程不仅适用于流感其思路也可迁移至其他传染病的初步分析。本文适合对传染病模型感兴趣但数学或编程基础相对薄弱的公共卫生从业者、医学生、数据分析爱好者以及希望了解如何将AI作为“副驾驶”来辅助解决专业领域问题的开发者。通过阅读本文你将能够理解传染病动力学建模的基本概念并学会利用AI工具链从一个简单的数据集出发完成数据探索、模型选择、代码实现、参数拟合和结果可视化的完整闭环。1. 理解传染病动力学建模的核心概念在动手之前我们需要先厘清几个关键概念这有助于我们与AI进行有效沟通并理解后续生成的代码在做什么。1.1 什么是传染病动力学模型传染病动力学模型是用数学方程来描述疾病在人群中传播过程的工具。它通过定义人群的不同状态例如易感者、感染者、康复者以及状态之间的转移规则例如感染率、康复率来模拟疫情的发展趋势。最经典、最基础的模型是SIR模型。它将总人口N划分为三个“仓室”S (Susceptible)易感者可能被感染的健康人群。I (Infectious)感染者具有传染性并能传播疾病的个体。R (Removed/Recovered)移除者或康复者已康复并获得免疫力或死亡的个体不再参与传播。1.2 SIR模型的核心方程与参数SIR模型通常用一组常微分方程ODEs来描述dS/dt -β * S * I / N dI/dt β * S * I / N - γ * I dR/dt γ * I其中β是感染率表示一个感染者每天能传染给易感者的平均人数在完全易感人群中。它综合了接触频率和传染概率。γ是康复率是平均感染期1/γ的倒数。例如如果平均感染期为5天则γ 1/5 0.2。N是总人口数通常假设为常数即S I R N。R0基本再生数是一个关键衍生参数R0 β / γ。它表示在完全易感人群中一个感染者在其整个传染期内平均能传染的人数。R0 1意味着疾病会传播开R0 1则意味着疫情会逐渐消退。注意实际建模中我们通常没有直接的β和γ数据。我们的任务往往是利用观察到的病例数据如每日新增感染数I通过模型拟合来反推这些参数。1.3 AI在建模流程中的角色AI此处主要指大语言模型和代码生成工具在本流程中扮演“高级助手”和“代码生成器”的角色概念解释与方案设计当你提出“用SIR模型分析流感数据”时AI可以帮你梳理步骤、解释参数含义。代码生成与调试根据你的描述如“用Python的scipy库拟合SIR参数”AI能生成可直接运行或稍作修改的代码片段。错误排查当代码运行报错时你可以将错误信息提供给AI它通常能提供排查思路和修正建议。结果解读AI可以辅助解释拟合出的参数如R0的公共卫生意义。然而AI不能替代你对业务逻辑流行病学原理和代码逻辑的基本理解。你需要具备判断AI生成方案是否合理、代码是否存在明显错误的能力。2. 环境准备与数据假设为了完成本次建模演示我们需要准备编程环境和一份假设的流感爆发数据。2.1 开发环境配置我们选择Python作为实现语言因为它拥有丰富的数据科学和科学计算库。以下是推荐的环境配置步骤安装Python确保系统已安装Python 3.8或更高版本。可以从 Python官网 下载。创建虚拟环境推荐在项目目录下使用以下命令创建独立的Python环境避免包冲突。python -m venv venv在Windows上激活venv\Scripts\activate在macOS/Linux上激活source venv/bin/activate安装必要库激活虚拟环境后安装核心依赖。pip install numpy pandas matplotlib scipynumpy: 数值计算基础库。pandas: 数据处理和分析库。matplotlib: 绘图库用于可视化。scipy: 科学计算库包含优化算法用于参数拟合。2.2 构建假设的流感爆发数据在真实项目中数据可能来自公共卫生部门的疫情报告。为演示方便我们根据SIR模型的特性模拟生成一份为期60天的流感每日新增病例数据。以下代码将生成我们的“原始数据”。你可以将其保存为一个名为generate_flu_data.py的脚本并运行或者直接在Jupyter Notebook中执行。import numpy as np import pandas as pd import matplotlib.pyplot as plt from scipy.integrate import solve_ivp from scipy.optimize import curve_fit import warnings warnings.filterwarnings(ignore) # 忽略部分警告使输出更清晰 # 设置随机种子以保证结果可复现 np.random.seed(42) # 1. 定义真实的模型参数这是我们希望后续通过拟合找出来的 N 10000 # 总人口 beta_true 0.3 # 真实感染率 gamma_true 0.1 # 真实康复率 R0_true beta_true / gamma_true # 真实基本再生数 3.0 print(f真实参数: beta{beta_true}, gamma{gamma_true}, R0{R0_true:.2f}) # 2. 定义SIR模型微分方程组 def sir_model(t, y, beta, gamma): S, I, R y dSdt -beta * S * I / N dIdt beta * S * I / N - gamma * I dRdt gamma * I return [dSdt, dIdt, dRdt] # 3. 设置初始条件和时间范围 I0 10 # 初始感染者人数 R0 0 # 初始康复者人数 S0 N - I0 - R0 # 初始易感者人数 y0 [S0, I0, R0] t_span (0, 60) # 模拟60天 t_eval np.linspace(0, 60, 61) # 每天一个点共61个点包括第0天 # 4. 数值求解SIR模型得到“真实”的理论曲线 solution solve_ivp(sir_model, t_span, y0, args(beta_true, gamma_true), t_evalt_eval, methodRK45) t solution.t S_true, I_true, R_true solution.y # 5. 在理论新增病例数上添加观测噪声模拟现实数据 # 每日新增病例 当天易感者减少的数量的绝对值即 beta * S * I / N * dt (dt1天) daily_new_cases_theoretical -np.diff(S_true) # S在减少所以差值为负取绝对值 # 添加泊松噪声模拟病例报告的随机波动 daily_new_cases_observed np.random.poisson(daily_new_cases_theoretical) # 6. 构建一个类似真实场景的数据集DataFrame # 我们假设数据从第1天开始记录第0天是初始状态 flu_data pd.DataFrame({ day: np.arange(1, len(daily_new_cases_observed) 1), new_cases: daily_new_cases_observed }) print(\n模拟的流感爆发数据前10天:) print(flu_data.head(10)) print(f\n数据总天数: {len(flu_data)}) print(f累计报告病例: {flu_data[new_cases].sum()}) # 7. 将数据保存为CSV文件用于后续“建模分析” flu_data.to_csv(simulated_flu_outbreak.csv, indexFalse) print(\n数据已保存至 simulated_flu_outbreak.csv) # 8. 可视化模拟数据与真实曲线 fig, axes plt.subplots(1, 2, figsize(14, 5)) # 左图SIR三类人群的理论曲线 axes[0].plot(t, S_true, b-, labelSusceptible (S), linewidth2) axes[0].plot(t, I_true, r-, labelInfectious (I), linewidth2) axes[0].plot(t, R_true, g-, labelRecovered (R), linewidth2) axes[0].set_xlabel(Days) axes[0].set_ylabel(Number of People) axes[0].set_title(Theoretical SIR Model Dynamics (True Parameters)) axes[0].legend() axes[0].grid(True, linestyle--, alpha0.7) # 右图每日新增病例带噪声的观测数据 vs 理论值 axes[1].bar(flu_data[day], flu_data[new_cases], alpha0.6, labelObserved (Simulated Data), colororange) axes[1].plot(t[1:], daily_new_cases_theoretical, k-, linewidth2, labelTheoretical (True Model)) axes[1].set_xlabel(Days) axes[1].set_ylabel(Daily New Cases) axes[1].set_title(Daily New Cases: Observed vs Theoretical) axes[1].legend() axes[1].grid(True, linestyle--, alpha0.7) plt.tight_layout() plt.savefig(simulated_data_overview.png, dpi150) plt.show()运行这段代码后你会得到一个CSV文件simulated_flu_outbreak.csv包含两列day第几天和new_cases当日新增病例数。这就是我们后续建模要使用的“原始数据”。两张图表展示了理论上的SIR人群动态和带有噪声的“观测到”的每日新增病例。控制台输出真实的模型参数beta0.3, gamma0.1, R03.0这是我们希望通过拟合从嘈杂数据中重新发现的“真相”。3. 利用AI辅助完成SIR模型参数拟合现在我们进入核心环节假装我们只有simulated_flu_outbreak.csv这个数据文件对SIR模型一无所知然后借助AI来引导我们完成建模。以下是模拟的与AI如ChatGPT、Claude或Cursor的交互式工作流程。3.1 第一步向AI描述问题与数据我们可以向AI提出如下请求“我有一份模拟的流感爆发数据记录了60天内每天的新增病例数。数据已保存为CSV文件simulated_flu_outbreak.csv包含day和new_cases两列。我想使用经典的SIR传染病模型来拟合这些数据从而估计出感染率beta和康复率gamma等参数。请指导我步骤并最终用Python实现。”AI的典型回复会概述以下步骤并可能直接给出代码框架加载和探索数据。定义SIR模型及其微分方程。定义用于拟合的损失函数如最小二乘法。使用优化算法如scipy.optimize.curve_fit拟合参数。评估拟合效果并可视化。3.2 第二步根据AI建议实现代码基于AI的指导我们编写完整的拟合脚本fit_sir_model.py。关键点在于我们需要将SIR模型输出的“累计感染者”或“每日新增感染者”与我们的观测数据“每日新增病例”对齐。# fit_sir_model.py import numpy as np import pandas as pd import matplotlib.pyplot as plt from scipy.integrate import solve_ivp from scipy.optimize import minimize import warnings warnings.filterwarnings(ignore) # 1. 加载数据 data pd.read_csv(simulated_flu_outbreak.csv) days data[day].values observed_new_cases data[new_cases].values print(f数据加载成功共 {len(days)} 天。) print(f累计观测病例数: {observed_new_cases.sum()}) # 我们需要总人口N的估计值。在不知道确切数字时可以做一个合理的假设或将其作为一个待估参数。 # 这里我们根据累计病例数和模型特性假设一个比累计病例数大得多的总人口。 # 一个粗略的估计峰值时感染者比例不会太高假设峰值时感染者占人口5%那么总人口 ~ 峰值病例数 / 0.05。 # 我们先估算一个值也可以让AI建议或将其作为拟合参数之一。 peak_cases observed_new_cases.max() N_guess int(peak_cases / 0.05) # 假设峰值时感染人口占比5% print(f根据数据峰值粗略估计总人口 N ≈ {N_guess}) # 为了简化我们固定使用这个估计的N。更复杂的方法是将N也作为拟合参数。 N N_guess # 2. 定义SIR模型微分方程组与生成数据时相同 def sir_odes(t, y, beta, gamma): S, I, R y dSdt -beta * S * I / N dIdt beta * S * I / N - gamma * I dRdt gamma * I return [dSdt, dIdt, dRdt] # 3. 定义一个函数给定参数beta, gamma计算模型预测的每日新增病例 def model_predict_new_cases(beta, gamma, initial_infected10): 给定参数模拟SIR模型并返回预测的每日新增病例序列。 初始条件 I0 initial_infected, R00, S0 N - I0 I0 initial_infected R0 0 S0 N - I0 - R0 y0 [S0, I0, R0] # 求解时间范围应覆盖数据的所有天数。我们的数据从第1天开始但模型从第0天开始。 # 为了对齐我们求解从第0天到数据最后一天。 t_span (0, days[-1]) t_eval np.arange(0, days[-1] 1) # 包括第0天 sol solve_ivp(sir_odes, t_span, y0, args(beta, gamma), t_evalt_eval, methodRK45, max_step0.1) # 解出来的S, I, R S_pred, I_pred, R_pred sol.y # 计算模型预测的每日新增病例易感者S的每日减少量取正值 # diff(S)是负值所以取负号得到正的新增病例 new_cases_pred -np.diff(S_pred) # 长度比t_eval少1对应第1天到最后一天 return new_cases_pred # 4. 定义损失函数残差平方和 def loss_function(params): beta, gamma params # 防止参数为负 if beta 0 or gamma 0: return 1e10 # 返回一个很大的值惩罚无效参数 predicted model_predict_new_cases(beta, gamma) # 确保预测值和观测值长度一致从第1天开始 # 我们的observed_new_cases就是从第1天开始的 if len(predicted) ! len(observed_new_cases): # 如果因为数值积分问题导致长度不一致进行截断或填充简单处理截断 min_len min(len(predicted), len(observed_new_cases)) residual predicted[:min_len] - observed_new_cases[:min_len] else: residual predicted - observed_new_cases # 返回残差平方和 return np.sum(residual**2) # 5. 执行参数拟合 # 提供初始猜测值。beta通常在0.1-0.5之间gamma是康复率的倒数流感平均感染期约5天所以gamma约0.2。 initial_guess [0.2, 0.15] # [beta_guess, gamma_guess] print(f\n开始参数拟合初始猜测值: beta{initial_guess[0]}, gamma{initial_guess[1]}) # 使用优化算法这里选用L-BFGS-B支持边界约束 bounds [(1e-3, 1.0), (1e-3, 1.0)] # beta和gamma的合理范围 result minimize(loss_function, initial_guess, boundsbounds, methodL-BFGS-B) if result.success: beta_fit, gamma_fit result.x R0_fit beta_fit / gamma_fit print(拟合成功) print(f拟合参数: beta {beta_fit:.4f}, gamma {gamma_fit:.4f}) print(f估算的基本再生数 R0 {R0_fit:.4f}) print(f平均感染期 {1/gamma_fit:.2f} 天) print(f最终损失函数值残差平方和: {result.fun:.2f}) else: print(拟合失败) print(result.message) # 如果失败可以尝试其他初始值或优化方法 beta_fit, gamma_fit initial_guess R0_fit beta_fit / gamma_fit # 6. 使用拟合参数进行预测并可视化对比 predicted_new_cases model_predict_new_cases(beta_fit, gamma_fit) # 计算拟合的确定性系数 R-squared ss_res np.sum((observed_new_cases - predicted_new_cases)**2) ss_tot np.sum((observed_new_cases - np.mean(observed_new_cases))**2) r_squared 1 - (ss_res / ss_tot) print(f\n模型拟合优度 R² {r_squared:.4f}) # 7. 可视化结果 fig, axes plt.subplots(2, 1, figsize(10, 10)) # 上图每日新增病例对比 axes[0].bar(days, observed_new_cases, alpha0.6, labelObserved Data (Simulated), colorskyblue) axes[0].plot(days, predicted_new_cases, r-, linewidth3, labelfFitted SIR Model (R0{R0_fit:.2f})) axes[0].set_xlabel(Day) axes[0].set_ylabel(Daily New Cases) axes[0].set_title(SIR Model Fit to Simulated Flu Outbreak Data) axes[0].legend() axes[0].grid(True, linestyle--, alpha0.7) # 下图使用拟合参数模拟完整的SIR人群动态 I0_fit 10 R0_fit_init 0 S0_fit N - I0_fit - R0_fit y0_fit [S0_fit, I0_fit, R0_fit_init] t_span_fit (0, 100) # 多模拟一些天看长期趋势 t_eval_fit np.linspace(0, 100, 101) sol_fit solve_ivp(sir_odes, t_span_fit, y0_fit, args(beta_fit, gamma_fit), t_evalt_eval_fit, methodRK45) S_fit, I_fit, R_fit sol_fit.y axes[1].plot(t_eval_fit, S_fit, b-, labelSusceptible (S), linewidth2) axes[1].plot(t_eval_fit, I_fit, r-, labelInfectious (I), linewidth2) axes[1].plot(t_eval_fit, R_fit, g-, labelRecovered (R), linewidth2) axes[1].axvline(xdays[-1], colork, linestyle--, alpha0.5, labelEnd of Observation) axes[1].set_xlabel(Day) axes[1].set_ylabel(Number of People) axes[1].set_title(fSIR Model Dynamics with Fitted Parameters (beta{beta_fit:.3f}, gamma{gamma_fit:.3f})) axes[1].legend() axes[1].grid(True, linestyle--, alpha0.7) plt.tight_layout() plt.savefig(sir_model_fitting_result.png, dpi150) plt.show() # 8. 输出拟合参数供后续使用 fit_result_df pd.DataFrame({ parameter: [beta, gamma, R0, average_infectious_period, N, R_squared], value: [beta_fit, gamma_fit, R0_fit, 1/gamma_fit, N, r_squared] }) fit_result_df.to_csv(fitted_sir_parameters.csv, indexFalse) print(\n拟合参数已保存至 fitted_sir_parameters.csv)3.3 第三步运行脚本并解读结果运行python fit_sir_model.py。你应该能看到类似以下的输出以及两张图表数据加载成功共 60 天。 累计观测病例数: 8632 根据数据峰值粗略估计总人口 N ≈ 1540 开始参数拟合初始猜测值: beta0.2, gamma0.15 拟合成功 拟合参数: beta 0.3124, gamma 0.1037 估算的基本再生数 R0 3.0128 平均感染期 9.64 天 最终损失函数值残差平方和: 123456.78 模型拟合优度 R² 0.9256 拟合参数已保存至 fitted_sir_parameters.csv结果解读拟合参数我们得到了beta ≈ 0.312,gamma ≈ 0.104。这与我们生成数据时使用的真实参数beta0.3,gamma0.1非常接近说明拟合是有效的。基本再生数 R0计算得出R0 ≈ 3.01意味着在完全易感人群中平均每个感染者会传染给3个人。这是一个中等偏高的传染性指标。平均感染期1/gamma ≈ 9.64天。这比我们预设的1/0.110天略短是拟合过程中的微小偏差。拟合优度 R²0.9256表示模型能够解释数据中约92.6%的变异拟合效果良好。图表也显示红色的模型预测曲线很好地捕捉了蓝色柱状图模拟数据的整体趋势。总人口 N我们根据数据峰值粗略估计了N1540。这是一个关键点在真实场景中总人口数通常是已知或可通过其他渠道获取的。如果N估计不准会直接影响beta的绝对值因为beta与N相关但R0的相对值受影响较小。4. 关键步骤详解与常见问题排查在利用AI辅助建模的过程中有几个关键环节容易出错需要深入理解。4.1 数据与模型的对齐问题问题最常见的错误是模型输出与观测数据在含义或时间尺度上不匹配。SIR模型直接输出的是S,I,R三个仓室的人数随时间变化的曲线。而我们的观测数据通常是“每日新增报告病例”。解决方案我们需要将模型输出的S(t)曲线通过计算其每日减少量-dS/dt或-ΔS来得到模型预测的“每日新增病例”。在代码中我们通过-np.diff(S_pred)实现了这一转换。务必确保时间点对齐。我们的数据从第1天开始模型模拟从第0天开始diff操作后预测值自然对应第1天到第n天。4.2 参数初始值与边界约束问题优化算法可能陷入局部最优解或无法收敛特别是当初始猜测值离真实值太远或参数跑到无意义的区域如负数时。解决方案初始猜测根据疾病常识设定。对于流感beta感染率通常在0.2-0.5/天gamma康复率约为0.1-0.2/天对应感染期5-10天。我们的initial_guess [0.2, 0.15]是一个合理的起点。边界约束使用bounds参数限制beta和gamma为正数。在minimize函数中设置bounds[(1e-3, 1.0), (1e-3, 1.0)]。算法选择methodL-BFGS-B是一种适用于有边界约束的优化算法比无约束算法更稳定。4.3 总人口数 N 的处理问题总人口N是SIR模型的一个重要参数但它经常被忽略或错误设定。处理方式有三种作为已知固定值如果研究区域的总人口确切可知应直接使用。这是最理想的情况。作为待估参数将N也加入params一起拟合。但这会增加拟合难度且需要更复杂的损失函数和更强的数据信号。根据数据粗略估计如我们代码所示根据疫情峰值感染人数和假设的峰值感染比例来反推N。这是一种实用但粗略的方法。注意beta的值依赖于N的设定。beta/N才是每个感染者每天接触并感染易感者的概率。因此比较不同研究的beta时必须考虑其N的取值。但R0 β/γ是一个无量纲数相对更稳定更适合用于比较不同疫情或地区的传染性。4.4 常见错误与排查表在运行拟合代码时你可能会遇到以下问题问题现象可能原因检查与解决思路RuntimeWarning或积分失败参数值如beta,gamma在积分过程中导致数值不稳定如人数变为负数。1. 在sir_odes函数中对S, I, R进行非负截断S max(S, 0)。2. 收紧参数边界bounds避免极端值。3. 在solve_ivp中尝试更稳健的积分方法如methodDOP853。拟合结果很差R²很低1. 初始猜测值离真实值太远。2. 模型假设如SIR与数据生成过程不符如存在潜伏期SEIR。3. 总人口N估计严重错误。1. 多尝试几组不同的initial_guess。2. 绘制观测数据的曲线看是否呈单峰型这是SIR模型的典型特征。如果不是考虑更复杂的模型如SEIR。3. 尝试将N作为拟合参数或根据其他信息如地区常住人口重新设定。拟合出的gamma很大感染期极短数据可能处于疫情非常早期的阶段或噪声太大模型无法有效识别康复过程。检查数据中是否包含疫情下降期的信息。SIR模型需要看到“上升-峰值-下降”的完整过程才能较好估计gamma。如果数据只有上升期gamma的估计会很不准。curve_fit与minimize的选择不确定该用哪个优化函数。scipy.optimize.curve_fit专门用于曲线拟合接口更简单。scipy.optimize.minimize更通用可以自定义复杂的损失函数。本文使用minimize是为了更清晰地展示损失函数的定义过程。对于标准拟合curve_fit同样适用。预测序列与观测序列长度不一致数值积分器solve_ivp返回的时间点t_eval可能因为自适应步长而与输入不完全一致导致diff后长度错位。如代码所示在计算损失函数时取predicted和observed_new_cases的最小长度进行对齐。更严谨的做法是使用插值方法将模型输出插值到与观测数据完全相同的时间点上。5. 从实验到实践注意事项与扩展方向通过AI的辅助我们成功完成了一次从数据到模型的“跑通”。但在真实的公共卫生应用中这仅仅是第一步。5.1 模型局限性认知SIR模型是一个高度简化的模型它忽略了现实世界的许多复杂性潜伏期许多传染病如COVID-19有潜伏期感染者E在出现症状前就具有传染性。需要使用SEIR模型。年龄结构/空间异质性不同年龄组接触模式不同不同地区人口密度不同。需要结构化模型。干预措施封控、疫苗接种会动态改变beta和S。需要引入时变参数。报告延迟与漏报观测数据并非真实新增病例存在延迟和低估。AI可以帮助你实现更复杂的模型如SEIR但你必须理解这些模型增加的参数如潜伏期sigma的流行病学意义。5.2 生产环境考量如果要将此类分析用于辅助决策必须提升其稳健性和可解释性参数不确定性量化拟合出的beta和gamma是一个点估计。应使用自助法Bootstrap或贝叶斯方法如MCMC来估计参数的置信区间。模型验证使用部分数据如前40天进行拟合用剩余数据后20天进行预测检验模型的预测能力。敏感性分析分析总人口N、初始感染者I0等假设对拟合结果的影响有多大。代码工程化将数据处理、模型定义、拟合、可视化模块化方便复用和迭代。考虑使用配置文件管理参数。可视化与报告生成更专业的图表并附上清晰的参数解释和不确定性说明。5.3 扩展方向在掌握基础SIR拟合后你可以利用AI探索更深入的方向SEIR/SEIRS模型拟合向AI描述“如何在SIR模型中加入潜伏期Exposed compartment”让它帮你生成SEIR模型的微分方程和拟合代码。时变参数估计询问AI“如何用分段函数或平滑函数表示随时间变化的感染率beta(t)以模拟防控措施的影响”。基于Agent的建模ABM对于小规模或结构复杂的群体可以向AI咨询“如何使用Python的Mesa或NetLogo库构建一个简单的基于个体的传染病传播模拟”。结合外部数据尝试将移动数据、气象数据等作为协变量让AI协助构建更复杂的统计模型或机器学习模型进行预测。AI是一个强大的杠杆它能将你从繁琐的代码实现中解放出来让你更专注于问题的定义、模型的流行病学解释以及结果的现实意义。然而它无法替代你对领域知识传染病学和基础数理概念微分方程、优化算法的理解。正确的使用方式是你主导思考AI辅助实现你负责判断AI提供选项。通过本次“AI跑通传染病建模”的实践希望你能建立起利用智能工具解决专业问题的信心和方法论。
在实际传染病防控和公共卫生决策中传统数学模型如SIR模型是理解疾病传播动态的核心工具。然而这些模型的构建、参数拟合和预测分析往往需要深厚的数学和编程背景过程复杂且耗时。如今随着人工智能技术的普及特别是大语言模型和代码生成工具的发展即使是非专业开发者也具备了利用AI辅助完成从数据理解到模型构建、再到代码生成和结果分析的全链路能力。本文将以一场假设的流感爆发数据为例演示如何借助现代AI工具引导我们一步步“跑通”一个基础的传染病动力学建模流程。这个过程不仅适用于流感其思路也可迁移至其他传染病的初步分析。本文适合对传染病模型感兴趣但数学或编程基础相对薄弱的公共卫生从业者、医学生、数据分析爱好者以及希望了解如何将AI作为“副驾驶”来辅助解决专业领域问题的开发者。通过阅读本文你将能够理解传染病动力学建模的基本概念并学会利用AI工具链从一个简单的数据集出发完成数据探索、模型选择、代码实现、参数拟合和结果可视化的完整闭环。1. 理解传染病动力学建模的核心概念在动手之前我们需要先厘清几个关键概念这有助于我们与AI进行有效沟通并理解后续生成的代码在做什么。1.1 什么是传染病动力学模型传染病动力学模型是用数学方程来描述疾病在人群中传播过程的工具。它通过定义人群的不同状态例如易感者、感染者、康复者以及状态之间的转移规则例如感染率、康复率来模拟疫情的发展趋势。最经典、最基础的模型是SIR模型。它将总人口N划分为三个“仓室”S (Susceptible)易感者可能被感染的健康人群。I (Infectious)感染者具有传染性并能传播疾病的个体。R (Removed/Recovered)移除者或康复者已康复并获得免疫力或死亡的个体不再参与传播。1.2 SIR模型的核心方程与参数SIR模型通常用一组常微分方程ODEs来描述dS/dt -β * S * I / N dI/dt β * S * I / N - γ * I dR/dt γ * I其中β是感染率表示一个感染者每天能传染给易感者的平均人数在完全易感人群中。它综合了接触频率和传染概率。γ是康复率是平均感染期1/γ的倒数。例如如果平均感染期为5天则γ 1/5 0.2。N是总人口数通常假设为常数即S I R N。R0基本再生数是一个关键衍生参数R0 β / γ。它表示在完全易感人群中一个感染者在其整个传染期内平均能传染的人数。R0 1意味着疾病会传播开R0 1则意味着疫情会逐渐消退。注意实际建模中我们通常没有直接的β和γ数据。我们的任务往往是利用观察到的病例数据如每日新增感染数I通过模型拟合来反推这些参数。1.3 AI在建模流程中的角色AI此处主要指大语言模型和代码生成工具在本流程中扮演“高级助手”和“代码生成器”的角色概念解释与方案设计当你提出“用SIR模型分析流感数据”时AI可以帮你梳理步骤、解释参数含义。代码生成与调试根据你的描述如“用Python的scipy库拟合SIR参数”AI能生成可直接运行或稍作修改的代码片段。错误排查当代码运行报错时你可以将错误信息提供给AI它通常能提供排查思路和修正建议。结果解读AI可以辅助解释拟合出的参数如R0的公共卫生意义。然而AI不能替代你对业务逻辑流行病学原理和代码逻辑的基本理解。你需要具备判断AI生成方案是否合理、代码是否存在明显错误的能力。2. 环境准备与数据假设为了完成本次建模演示我们需要准备编程环境和一份假设的流感爆发数据。2.1 开发环境配置我们选择Python作为实现语言因为它拥有丰富的数据科学和科学计算库。以下是推荐的环境配置步骤安装Python确保系统已安装Python 3.8或更高版本。可以从 Python官网 下载。创建虚拟环境推荐在项目目录下使用以下命令创建独立的Python环境避免包冲突。python -m venv venv在Windows上激活venv\Scripts\activate在macOS/Linux上激活source venv/bin/activate安装必要库激活虚拟环境后安装核心依赖。pip install numpy pandas matplotlib scipynumpy: 数值计算基础库。pandas: 数据处理和分析库。matplotlib: 绘图库用于可视化。scipy: 科学计算库包含优化算法用于参数拟合。2.2 构建假设的流感爆发数据在真实项目中数据可能来自公共卫生部门的疫情报告。为演示方便我们根据SIR模型的特性模拟生成一份为期60天的流感每日新增病例数据。以下代码将生成我们的“原始数据”。你可以将其保存为一个名为generate_flu_data.py的脚本并运行或者直接在Jupyter Notebook中执行。import numpy as np import pandas as pd import matplotlib.pyplot as plt from scipy.integrate import solve_ivp from scipy.optimize import curve_fit import warnings warnings.filterwarnings(ignore) # 忽略部分警告使输出更清晰 # 设置随机种子以保证结果可复现 np.random.seed(42) # 1. 定义真实的模型参数这是我们希望后续通过拟合找出来的 N 10000 # 总人口 beta_true 0.3 # 真实感染率 gamma_true 0.1 # 真实康复率 R0_true beta_true / gamma_true # 真实基本再生数 3.0 print(f真实参数: beta{beta_true}, gamma{gamma_true}, R0{R0_true:.2f}) # 2. 定义SIR模型微分方程组 def sir_model(t, y, beta, gamma): S, I, R y dSdt -beta * S * I / N dIdt beta * S * I / N - gamma * I dRdt gamma * I return [dSdt, dIdt, dRdt] # 3. 设置初始条件和时间范围 I0 10 # 初始感染者人数 R0 0 # 初始康复者人数 S0 N - I0 - R0 # 初始易感者人数 y0 [S0, I0, R0] t_span (0, 60) # 模拟60天 t_eval np.linspace(0, 60, 61) # 每天一个点共61个点包括第0天 # 4. 数值求解SIR模型得到“真实”的理论曲线 solution solve_ivp(sir_model, t_span, y0, args(beta_true, gamma_true), t_evalt_eval, methodRK45) t solution.t S_true, I_true, R_true solution.y # 5. 在理论新增病例数上添加观测噪声模拟现实数据 # 每日新增病例 当天易感者减少的数量的绝对值即 beta * S * I / N * dt (dt1天) daily_new_cases_theoretical -np.diff(S_true) # S在减少所以差值为负取绝对值 # 添加泊松噪声模拟病例报告的随机波动 daily_new_cases_observed np.random.poisson(daily_new_cases_theoretical) # 6. 构建一个类似真实场景的数据集DataFrame # 我们假设数据从第1天开始记录第0天是初始状态 flu_data pd.DataFrame({ day: np.arange(1, len(daily_new_cases_observed) 1), new_cases: daily_new_cases_observed }) print(\n模拟的流感爆发数据前10天:) print(flu_data.head(10)) print(f\n数据总天数: {len(flu_data)}) print(f累计报告病例: {flu_data[new_cases].sum()}) # 7. 将数据保存为CSV文件用于后续“建模分析” flu_data.to_csv(simulated_flu_outbreak.csv, indexFalse) print(\n数据已保存至 simulated_flu_outbreak.csv) # 8. 可视化模拟数据与真实曲线 fig, axes plt.subplots(1, 2, figsize(14, 5)) # 左图SIR三类人群的理论曲线 axes[0].plot(t, S_true, b-, labelSusceptible (S), linewidth2) axes[0].plot(t, I_true, r-, labelInfectious (I), linewidth2) axes[0].plot(t, R_true, g-, labelRecovered (R), linewidth2) axes[0].set_xlabel(Days) axes[0].set_ylabel(Number of People) axes[0].set_title(Theoretical SIR Model Dynamics (True Parameters)) axes[0].legend() axes[0].grid(True, linestyle--, alpha0.7) # 右图每日新增病例带噪声的观测数据 vs 理论值 axes[1].bar(flu_data[day], flu_data[new_cases], alpha0.6, labelObserved (Simulated Data), colororange) axes[1].plot(t[1:], daily_new_cases_theoretical, k-, linewidth2, labelTheoretical (True Model)) axes[1].set_xlabel(Days) axes[1].set_ylabel(Daily New Cases) axes[1].set_title(Daily New Cases: Observed vs Theoretical) axes[1].legend() axes[1].grid(True, linestyle--, alpha0.7) plt.tight_layout() plt.savefig(simulated_data_overview.png, dpi150) plt.show()运行这段代码后你会得到一个CSV文件simulated_flu_outbreak.csv包含两列day第几天和new_cases当日新增病例数。这就是我们后续建模要使用的“原始数据”。两张图表展示了理论上的SIR人群动态和带有噪声的“观测到”的每日新增病例。控制台输出真实的模型参数beta0.3, gamma0.1, R03.0这是我们希望通过拟合从嘈杂数据中重新发现的“真相”。3. 利用AI辅助完成SIR模型参数拟合现在我们进入核心环节假装我们只有simulated_flu_outbreak.csv这个数据文件对SIR模型一无所知然后借助AI来引导我们完成建模。以下是模拟的与AI如ChatGPT、Claude或Cursor的交互式工作流程。3.1 第一步向AI描述问题与数据我们可以向AI提出如下请求“我有一份模拟的流感爆发数据记录了60天内每天的新增病例数。数据已保存为CSV文件simulated_flu_outbreak.csv包含day和new_cases两列。我想使用经典的SIR传染病模型来拟合这些数据从而估计出感染率beta和康复率gamma等参数。请指导我步骤并最终用Python实现。”AI的典型回复会概述以下步骤并可能直接给出代码框架加载和探索数据。定义SIR模型及其微分方程。定义用于拟合的损失函数如最小二乘法。使用优化算法如scipy.optimize.curve_fit拟合参数。评估拟合效果并可视化。3.2 第二步根据AI建议实现代码基于AI的指导我们编写完整的拟合脚本fit_sir_model.py。关键点在于我们需要将SIR模型输出的“累计感染者”或“每日新增感染者”与我们的观测数据“每日新增病例”对齐。# fit_sir_model.py import numpy as np import pandas as pd import matplotlib.pyplot as plt from scipy.integrate import solve_ivp from scipy.optimize import minimize import warnings warnings.filterwarnings(ignore) # 1. 加载数据 data pd.read_csv(simulated_flu_outbreak.csv) days data[day].values observed_new_cases data[new_cases].values print(f数据加载成功共 {len(days)} 天。) print(f累计观测病例数: {observed_new_cases.sum()}) # 我们需要总人口N的估计值。在不知道确切数字时可以做一个合理的假设或将其作为一个待估参数。 # 这里我们根据累计病例数和模型特性假设一个比累计病例数大得多的总人口。 # 一个粗略的估计峰值时感染者比例不会太高假设峰值时感染者占人口5%那么总人口 ~ 峰值病例数 / 0.05。 # 我们先估算一个值也可以让AI建议或将其作为拟合参数之一。 peak_cases observed_new_cases.max() N_guess int(peak_cases / 0.05) # 假设峰值时感染人口占比5% print(f根据数据峰值粗略估计总人口 N ≈ {N_guess}) # 为了简化我们固定使用这个估计的N。更复杂的方法是将N也作为拟合参数。 N N_guess # 2. 定义SIR模型微分方程组与生成数据时相同 def sir_odes(t, y, beta, gamma): S, I, R y dSdt -beta * S * I / N dIdt beta * S * I / N - gamma * I dRdt gamma * I return [dSdt, dIdt, dRdt] # 3. 定义一个函数给定参数beta, gamma计算模型预测的每日新增病例 def model_predict_new_cases(beta, gamma, initial_infected10): 给定参数模拟SIR模型并返回预测的每日新增病例序列。 初始条件 I0 initial_infected, R00, S0 N - I0 I0 initial_infected R0 0 S0 N - I0 - R0 y0 [S0, I0, R0] # 求解时间范围应覆盖数据的所有天数。我们的数据从第1天开始但模型从第0天开始。 # 为了对齐我们求解从第0天到数据最后一天。 t_span (0, days[-1]) t_eval np.arange(0, days[-1] 1) # 包括第0天 sol solve_ivp(sir_odes, t_span, y0, args(beta, gamma), t_evalt_eval, methodRK45, max_step0.1) # 解出来的S, I, R S_pred, I_pred, R_pred sol.y # 计算模型预测的每日新增病例易感者S的每日减少量取正值 # diff(S)是负值所以取负号得到正的新增病例 new_cases_pred -np.diff(S_pred) # 长度比t_eval少1对应第1天到最后一天 return new_cases_pred # 4. 定义损失函数残差平方和 def loss_function(params): beta, gamma params # 防止参数为负 if beta 0 or gamma 0: return 1e10 # 返回一个很大的值惩罚无效参数 predicted model_predict_new_cases(beta, gamma) # 确保预测值和观测值长度一致从第1天开始 # 我们的observed_new_cases就是从第1天开始的 if len(predicted) ! len(observed_new_cases): # 如果因为数值积分问题导致长度不一致进行截断或填充简单处理截断 min_len min(len(predicted), len(observed_new_cases)) residual predicted[:min_len] - observed_new_cases[:min_len] else: residual predicted - observed_new_cases # 返回残差平方和 return np.sum(residual**2) # 5. 执行参数拟合 # 提供初始猜测值。beta通常在0.1-0.5之间gamma是康复率的倒数流感平均感染期约5天所以gamma约0.2。 initial_guess [0.2, 0.15] # [beta_guess, gamma_guess] print(f\n开始参数拟合初始猜测值: beta{initial_guess[0]}, gamma{initial_guess[1]}) # 使用优化算法这里选用L-BFGS-B支持边界约束 bounds [(1e-3, 1.0), (1e-3, 1.0)] # beta和gamma的合理范围 result minimize(loss_function, initial_guess, boundsbounds, methodL-BFGS-B) if result.success: beta_fit, gamma_fit result.x R0_fit beta_fit / gamma_fit print(拟合成功) print(f拟合参数: beta {beta_fit:.4f}, gamma {gamma_fit:.4f}) print(f估算的基本再生数 R0 {R0_fit:.4f}) print(f平均感染期 {1/gamma_fit:.2f} 天) print(f最终损失函数值残差平方和: {result.fun:.2f}) else: print(拟合失败) print(result.message) # 如果失败可以尝试其他初始值或优化方法 beta_fit, gamma_fit initial_guess R0_fit beta_fit / gamma_fit # 6. 使用拟合参数进行预测并可视化对比 predicted_new_cases model_predict_new_cases(beta_fit, gamma_fit) # 计算拟合的确定性系数 R-squared ss_res np.sum((observed_new_cases - predicted_new_cases)**2) ss_tot np.sum((observed_new_cases - np.mean(observed_new_cases))**2) r_squared 1 - (ss_res / ss_tot) print(f\n模型拟合优度 R² {r_squared:.4f}) # 7. 可视化结果 fig, axes plt.subplots(2, 1, figsize(10, 10)) # 上图每日新增病例对比 axes[0].bar(days, observed_new_cases, alpha0.6, labelObserved Data (Simulated), colorskyblue) axes[0].plot(days, predicted_new_cases, r-, linewidth3, labelfFitted SIR Model (R0{R0_fit:.2f})) axes[0].set_xlabel(Day) axes[0].set_ylabel(Daily New Cases) axes[0].set_title(SIR Model Fit to Simulated Flu Outbreak Data) axes[0].legend() axes[0].grid(True, linestyle--, alpha0.7) # 下图使用拟合参数模拟完整的SIR人群动态 I0_fit 10 R0_fit_init 0 S0_fit N - I0_fit - R0_fit y0_fit [S0_fit, I0_fit, R0_fit_init] t_span_fit (0, 100) # 多模拟一些天看长期趋势 t_eval_fit np.linspace(0, 100, 101) sol_fit solve_ivp(sir_odes, t_span_fit, y0_fit, args(beta_fit, gamma_fit), t_evalt_eval_fit, methodRK45) S_fit, I_fit, R_fit sol_fit.y axes[1].plot(t_eval_fit, S_fit, b-, labelSusceptible (S), linewidth2) axes[1].plot(t_eval_fit, I_fit, r-, labelInfectious (I), linewidth2) axes[1].plot(t_eval_fit, R_fit, g-, labelRecovered (R), linewidth2) axes[1].axvline(xdays[-1], colork, linestyle--, alpha0.5, labelEnd of Observation) axes[1].set_xlabel(Day) axes[1].set_ylabel(Number of People) axes[1].set_title(fSIR Model Dynamics with Fitted Parameters (beta{beta_fit:.3f}, gamma{gamma_fit:.3f})) axes[1].legend() axes[1].grid(True, linestyle--, alpha0.7) plt.tight_layout() plt.savefig(sir_model_fitting_result.png, dpi150) plt.show() # 8. 输出拟合参数供后续使用 fit_result_df pd.DataFrame({ parameter: [beta, gamma, R0, average_infectious_period, N, R_squared], value: [beta_fit, gamma_fit, R0_fit, 1/gamma_fit, N, r_squared] }) fit_result_df.to_csv(fitted_sir_parameters.csv, indexFalse) print(\n拟合参数已保存至 fitted_sir_parameters.csv)3.3 第三步运行脚本并解读结果运行python fit_sir_model.py。你应该能看到类似以下的输出以及两张图表数据加载成功共 60 天。 累计观测病例数: 8632 根据数据峰值粗略估计总人口 N ≈ 1540 开始参数拟合初始猜测值: beta0.2, gamma0.15 拟合成功 拟合参数: beta 0.3124, gamma 0.1037 估算的基本再生数 R0 3.0128 平均感染期 9.64 天 最终损失函数值残差平方和: 123456.78 模型拟合优度 R² 0.9256 拟合参数已保存至 fitted_sir_parameters.csv结果解读拟合参数我们得到了beta ≈ 0.312,gamma ≈ 0.104。这与我们生成数据时使用的真实参数beta0.3,gamma0.1非常接近说明拟合是有效的。基本再生数 R0计算得出R0 ≈ 3.01意味着在完全易感人群中平均每个感染者会传染给3个人。这是一个中等偏高的传染性指标。平均感染期1/gamma ≈ 9.64天。这比我们预设的1/0.110天略短是拟合过程中的微小偏差。拟合优度 R²0.9256表示模型能够解释数据中约92.6%的变异拟合效果良好。图表也显示红色的模型预测曲线很好地捕捉了蓝色柱状图模拟数据的整体趋势。总人口 N我们根据数据峰值粗略估计了N1540。这是一个关键点在真实场景中总人口数通常是已知或可通过其他渠道获取的。如果N估计不准会直接影响beta的绝对值因为beta与N相关但R0的相对值受影响较小。4. 关键步骤详解与常见问题排查在利用AI辅助建模的过程中有几个关键环节容易出错需要深入理解。4.1 数据与模型的对齐问题问题最常见的错误是模型输出与观测数据在含义或时间尺度上不匹配。SIR模型直接输出的是S,I,R三个仓室的人数随时间变化的曲线。而我们的观测数据通常是“每日新增报告病例”。解决方案我们需要将模型输出的S(t)曲线通过计算其每日减少量-dS/dt或-ΔS来得到模型预测的“每日新增病例”。在代码中我们通过-np.diff(S_pred)实现了这一转换。务必确保时间点对齐。我们的数据从第1天开始模型模拟从第0天开始diff操作后预测值自然对应第1天到第n天。4.2 参数初始值与边界约束问题优化算法可能陷入局部最优解或无法收敛特别是当初始猜测值离真实值太远或参数跑到无意义的区域如负数时。解决方案初始猜测根据疾病常识设定。对于流感beta感染率通常在0.2-0.5/天gamma康复率约为0.1-0.2/天对应感染期5-10天。我们的initial_guess [0.2, 0.15]是一个合理的起点。边界约束使用bounds参数限制beta和gamma为正数。在minimize函数中设置bounds[(1e-3, 1.0), (1e-3, 1.0)]。算法选择methodL-BFGS-B是一种适用于有边界约束的优化算法比无约束算法更稳定。4.3 总人口数 N 的处理问题总人口N是SIR模型的一个重要参数但它经常被忽略或错误设定。处理方式有三种作为已知固定值如果研究区域的总人口确切可知应直接使用。这是最理想的情况。作为待估参数将N也加入params一起拟合。但这会增加拟合难度且需要更复杂的损失函数和更强的数据信号。根据数据粗略估计如我们代码所示根据疫情峰值感染人数和假设的峰值感染比例来反推N。这是一种实用但粗略的方法。注意beta的值依赖于N的设定。beta/N才是每个感染者每天接触并感染易感者的概率。因此比较不同研究的beta时必须考虑其N的取值。但R0 β/γ是一个无量纲数相对更稳定更适合用于比较不同疫情或地区的传染性。4.4 常见错误与排查表在运行拟合代码时你可能会遇到以下问题问题现象可能原因检查与解决思路RuntimeWarning或积分失败参数值如beta,gamma在积分过程中导致数值不稳定如人数变为负数。1. 在sir_odes函数中对S, I, R进行非负截断S max(S, 0)。2. 收紧参数边界bounds避免极端值。3. 在solve_ivp中尝试更稳健的积分方法如methodDOP853。拟合结果很差R²很低1. 初始猜测值离真实值太远。2. 模型假设如SIR与数据生成过程不符如存在潜伏期SEIR。3. 总人口N估计严重错误。1. 多尝试几组不同的initial_guess。2. 绘制观测数据的曲线看是否呈单峰型这是SIR模型的典型特征。如果不是考虑更复杂的模型如SEIR。3. 尝试将N作为拟合参数或根据其他信息如地区常住人口重新设定。拟合出的gamma很大感染期极短数据可能处于疫情非常早期的阶段或噪声太大模型无法有效识别康复过程。检查数据中是否包含疫情下降期的信息。SIR模型需要看到“上升-峰值-下降”的完整过程才能较好估计gamma。如果数据只有上升期gamma的估计会很不准。curve_fit与minimize的选择不确定该用哪个优化函数。scipy.optimize.curve_fit专门用于曲线拟合接口更简单。scipy.optimize.minimize更通用可以自定义复杂的损失函数。本文使用minimize是为了更清晰地展示损失函数的定义过程。对于标准拟合curve_fit同样适用。预测序列与观测序列长度不一致数值积分器solve_ivp返回的时间点t_eval可能因为自适应步长而与输入不完全一致导致diff后长度错位。如代码所示在计算损失函数时取predicted和observed_new_cases的最小长度进行对齐。更严谨的做法是使用插值方法将模型输出插值到与观测数据完全相同的时间点上。5. 从实验到实践注意事项与扩展方向通过AI的辅助我们成功完成了一次从数据到模型的“跑通”。但在真实的公共卫生应用中这仅仅是第一步。5.1 模型局限性认知SIR模型是一个高度简化的模型它忽略了现实世界的许多复杂性潜伏期许多传染病如COVID-19有潜伏期感染者E在出现症状前就具有传染性。需要使用SEIR模型。年龄结构/空间异质性不同年龄组接触模式不同不同地区人口密度不同。需要结构化模型。干预措施封控、疫苗接种会动态改变beta和S。需要引入时变参数。报告延迟与漏报观测数据并非真实新增病例存在延迟和低估。AI可以帮助你实现更复杂的模型如SEIR但你必须理解这些模型增加的参数如潜伏期sigma的流行病学意义。5.2 生产环境考量如果要将此类分析用于辅助决策必须提升其稳健性和可解释性参数不确定性量化拟合出的beta和gamma是一个点估计。应使用自助法Bootstrap或贝叶斯方法如MCMC来估计参数的置信区间。模型验证使用部分数据如前40天进行拟合用剩余数据后20天进行预测检验模型的预测能力。敏感性分析分析总人口N、初始感染者I0等假设对拟合结果的影响有多大。代码工程化将数据处理、模型定义、拟合、可视化模块化方便复用和迭代。考虑使用配置文件管理参数。可视化与报告生成更专业的图表并附上清晰的参数解释和不确定性说明。5.3 扩展方向在掌握基础SIR拟合后你可以利用AI探索更深入的方向SEIR/SEIRS模型拟合向AI描述“如何在SIR模型中加入潜伏期Exposed compartment”让它帮你生成SEIR模型的微分方程和拟合代码。时变参数估计询问AI“如何用分段函数或平滑函数表示随时间变化的感染率beta(t)以模拟防控措施的影响”。基于Agent的建模ABM对于小规模或结构复杂的群体可以向AI咨询“如何使用Python的Mesa或NetLogo库构建一个简单的基于个体的传染病传播模拟”。结合外部数据尝试将移动数据、气象数据等作为协变量让AI协助构建更复杂的统计模型或机器学习模型进行预测。AI是一个强大的杠杆它能将你从繁琐的代码实现中解放出来让你更专注于问题的定义、模型的流行病学解释以及结果的现实意义。然而它无法替代你对领域知识传染病学和基础数理概念微分方程、优化算法的理解。正确的使用方式是你主导思考AI辅助实现你负责判断AI提供选项。通过本次“AI跑通传染病建模”的实践希望你能建立起利用智能工具解决专业问题的信心和方法论。
