Matlab统计建模三步法：数据清洗、可视化与假设验证实战

1. 项目概述：统计建模的“三步醒肤法”

如果你在Matlab里搞统计建模，结果总感觉模型“睡不醒”——预测不准、解释力弱、结果不稳定，那问题很可能出在最开始的几步。很多人拿到数据就急着跑回归、拟合分布，却忽略了建模前那些看似基础、实则决定性的准备工作。这就好比不洗脸、不补水就直接上妆，妆容怎么可能服帖持久？所谓“醒肤”，就是让数据“清醒”过来，呈现出它最真实、最可被模型理解的状态。根据我多年的数据分析经验，无论模型多复杂，有三步基础操作是绝对不能跳过的：数据质量诊断与清洗、变量关系的可视化探索、以及模型假设的初步验证。这三步是连接原始数据与统计模型的桥梁，缺了任何一步，你的建模工作都可能建立在流沙之上。接下来，我就结合Matlab这个强大的工具，把这“三步醒肤法”的实操细节、背后的统计原理，以及我踩过的坑，给你掰开揉碎了讲清楚。

2. 核心思路：为什么是这三步？

在深入代码之前，我们得先想明白为什么是这三步，而不是别的。统计建模的本质，是从充满噪声的现实数据中，提炼出稳定、可靠的关系模式。Matlab提供了琳琅满目的工具箱，从最基本的fitlm做线性回归，到复杂的fitrgp做高斯过程回归，函数调用可能只需要一行代码。但**“Garbage in, garbage out”（垃圾进，垃圾出）** 这条计算领域的铁律，在统计建模中体现得尤为深刻。

第一步，数据质量诊断与清洗，目标是解决“输入是不是垃圾”的问题。原始数据几乎总是存在缺失值、异常值、记录错误或量纲不一的情况。不处理缺失值，Matlab的很多建模函数会直接报错或剔除整行，导致信息损失；不处理异常值，一个离群点就可能把回归线“拉偏”，让模型结论完全失真。这一步是确保数据“可用”。

第二步，变量关系的可视化探索，目标是回答“数据大概长什么样”的问题。统计模型有很多种（线性、非线性、参数、非参数），选错模型类型就像用螺丝刀去敲钉子。通过散点图矩阵、箱线图、分布直方图等可视化手段，我们可以直观地看到变量间是线性关系还是曲线关系，是否存在明显的分组效应，方差是否均匀。这一步是为模型选择提供“视觉证据”，避免盲目试错。

第三步，模型假设的初步验证，目标是提前检验“我想用的模型基本假设成立吗”。例如，普通最小二乘线性回归的核心假设包括：线性关系、误差项独立同分布、同方差性、正态性等。如果数据明显违背这些假设，那么模型参数的估计值、显著性检验的p值都将不可信。在Matlab里，我们可以在正式拟合复杂模型前，用一些简单方法（如残差图、相关性检验）进行快速筛查。这一步是给建模思路“踩刹车”，防止方向性错误。

这三步环环相扣，构成了一个严谨的数据分析工作流。它们消耗的时间可能占整个项目的30%-50%，但能避免后续80%的返工和错误结论。

3. 第一步：数据质量诊断与清洗实战

拿到数据集后，千万别急着plot或fit。我的习惯是，先把它当成一个需要体检的病人。

3.1 缺失值诊断与处理

在Matlab中，缺失值通常用NaN表示。首先是用summary函数或自编脚本来个快速普查。

% 假设你的数据在一个表T中 disp('各变量缺失值数量：') missing_summary = sum(ismissing(T)); disp(missing_summary) % 更直观的，查看缺失值比例 missing_percentage = (sum(ismissing(T)) / height(T)) * 100; disp('各变量缺失值比例(%)：') disp(missing_percentage)

如果某个变量的缺失比例超过30%，你就需要慎重考虑是否还要保留它。因为无论用什么方法填补，其信息可靠性都存疑。

对于需要处理的缺失值，常见方法有：

1. 直接删除：如果缺失值很少（如<5%），且是随机缺失，用rmmissing函数删除含有缺失值的行是最简单的。但要注意，这可能减少样本量，或引入偏差（如果缺失不是完全随机的）。

   T_clean = rmmissing(T); % 删除任何列包含NaN的行

2. 均值/中位数/众数填补：适用于数值型变量。用fillmissing函数。

   % 用列均值填补 T_filled_mean = fillmissing(T, 'constant', mean(T.VarName, 'omitnan')); % 用列中位数填补（对异常值更稳健） T_filled_median = fillmissing(T, 'constant', median(T.VarName, 'omitnan'));

3. 插值法：对于时间序列数据，fillmissing也支持线性、样条插值。

   T_filled_linear = fillmissing(T, 'linear');

实操心得：不要盲目使用均值填补！如果数据存在明显的组别或趋势，分组填补或模型预测填补（如用回归预测缺失值）会更合理。对于类别变量，可以单独增加一个“缺失”类别。

3.2 异常值检测与处理

异常值不一定是错误，但会极大影响模型。我常用的检测方法是“可视化初筛+统计量确认”。

可视化初筛：箱线图是利器。

figure; boxplot(T.VarName); title('变量VarName的箱线图'); ylabel('值');

箱线图能直观显示中位数、四分位数和“须”的范围，落在须外的点通常被视为候选异常值。

统计量确认：常用的是基于标准差或四分位距的方法。

data = T.VarName; data = data(~isnan(data)); % 先去除NaN % 方法1：3σ原则（假设数据近似正态） mean_val = mean(data); std_val = std(data); outliers_3sigma = data(data < mean_val - 3*std_val | data > mean_val + 3*std_val); % 方法2：基于四分位距（IQR），更稳健 Q = quantile(data, [0.25 0.75]); IQR = Q(2) - Q(1); lower_bound = Q(1) - 1.5 * IQR; upper_bound = Q(2) + 1.5 * IQR; outliers_iqr = data(data < lower_bound | data > upper_bound); disp(['基于IQR发现的异常值数量：', num2str(length(outliers_iqr))]);

处理异常值需要结合业务知识：

• 修正或删除：如果是明显的记录错误（如身高2.5米），应修正或删除。
• 保留但处理：如果是真实但极端的值，可以考虑在建模时使用对异常值不敏感的模型（如决策树、分位数回归），或在回归中引入鲁棒性方法（如fitlm中的RobustOpts选项）。
• 转换：对偏态严重的变量进行对数转换、平方根转换，可以压缩极端值的影响。

  T.VarName_log = log(T.VarName + 1); % +1防止有0值

3.3 数据一致性检查与转换

检查分类变量的取值是否一致，比如“男”、“Male”、“M”应该统一。对于数值变量，检查量纲，考虑是否需要标准化（zscore）或归一化（缩放到[0,1]），特别是当变量单位差异巨大时（如收入以元计，年龄以年计）。

% 标准化 (Z-score) T.VarName_Z = (T.VarName - mean(T.VarName, 'omitnan')) / std(T.VarName, 'omitnan'); % 最小-最大归一化 min_val = min(T.VarName, [], 'omitnan'); max_val = max(T.VarName, [], 'omitnan'); T.VarName_Norm = (T.VarName - min_val) / (max_val - min_val);

注意事项：数据转换（特别是标准化）需要在训练集上计算参数（均值、标准差），然后将其应用于测试集，以避免数据泄露。可以用cvpartition划分数据后分别处理。

4. 第二步：变量关系的可视化探索

清洗后的数据，我们要像侦探一样观察它。Matlab的绘图功能在这里大放异彩。

4.1 单变量分布探查

了解每个变量的分布形态，是后续选择模型和检验假设的基础。

figure; subplot(1,2,1); histogram(T.VarName, 'Normalization', 'pdf'); hold on; % 叠加核密度估计，更平滑 [f, xi] = ksdensity(T.VarName); plot(xi, f, 'LineWidth', 2); title('直方图与核密度估计'); xlabel('VarName'); ylabel('密度'); subplot(1,2,2); boxplot(T.VarName); title('箱线图');

通过分布图，你可以判断变量是否接近正态分布（钟形曲线）。严重的偏态或双峰分布会违背许多参数模型的假设。

4.2 双变量关系探查

这是探索的核心，目标是看预测变量（X）和响应变量（Y）之间是什么关系。

散点图：最基本，也最有效。

figure; scatter(T.X1, T.Y, 20, 'filled'); xlabel('预测变量 X1'); ylabel('响应变量 Y'); title('X1与Y的散点图'); grid on; % 可以加个简单的趋势线看看 hold on; p = polyfit(T.X1(~isnan(T.X1) & ~isnan(T.Y)), T.Y(~isnan(T.X1) & ~isnan(T.Y)), 1); x_fit = linspace(min(T.X1), max(T.X1), 100); y_fit = polyval(p, x_fit); plot(x_fit, y_fit, 'r-', 'LineWidth', 2); legend('数据点', '线性趋势线', 'Location', 'best');

观察点是随机散布，还是有明显的线性、曲线趋势？点云是均匀的漏斗形，还是其他形状？这直接关系到你该用线性模型还是非线性模型。

散点图矩阵：当有多个预测变量时，plotmatrix或gplotmatrix（统计和机器学习工具箱）可以一次性查看所有两两关系。

figure; [S, AX, BigAx, H, HAx] = plotmatrix(table2array(T(:, {'X1', 'X2', 'X3', 'Y'}))); title('变量间散点图矩阵');

对角线位置通常是每个变量的直方图，非对角线位置是两两散点图。这能帮你快速发现强相关的预测变量（共线性问题），以及每个X与Y的关系模式。

分组可视化：如果数据有分类变量（如性别、地区），用gscatter可以按组着色。

figure; gscatter(T.X1, T.Y, T.Group); xlabel('X1'); ylabel('Y'); legend('Location', 'best'); title('按组别着色的散点图');

如果不同组别的数据点呈现出不同的趋势或截距，那你在建模时就需要考虑加入交互项或使用分层模型。

4.3 相关性分析

可视化之后，用数值量化关系强度。corr函数计算相关系数矩阵。

corr_matrix = corr(T{:,:}, 'Rows', 'pairwise'); % 'pairwise'成对删除缺失值 figure; imagesc(corr_matrix); colorbar; title('变量相关系数矩阵热图'); set(gca, 'XTick', 1:width(T), 'XTickLabel', T.Properties.VariableNames); set(gca, 'YTick', 1:width(T), 'YTickLabel', T.Properties.VariableNames);

注意：corr计算的是线性相关系数（默认皮尔逊相关）。如果散点图显示非线性关系，皮尔逊相关可能很低，但实际存在强关联。此时应考虑秩相关（'type', 'Spearman'）。

5. 第三步：模型假设的初步验证

假设我们初步判断可以用线性模型。在调用fitlm正式拟合之前，我们可以做几项快速的“体检”。

5.1 线性与独立性初步判断

线性关系我们在散点图里已经看了。独立性通常指误差项独立，这往往与数据采集过程有关（如时间序列数据常不独立）。一个快速的检查是看残差与拟合值或预测变量的散点图是否有规律模式。虽然正式残差要等模型拟合后才有，但我们可以用“局部回归”或“平滑散点图”来预判。

% 使用平滑散点图（lowess）观察趋势 figure; scatter(T.X1, T.Y, 'filled'); hold on; % 使用 smoothdata 进行局部加权散点平滑 ysmooth = smoothdata(T.Y, 'loess', T.X1); % 需要 R2017a 或更高版本 plot(T.X1, ysmooth, 'r-', 'LineWidth', 3); xlabel('X1'); ylabel('Y'); title('数据点与Lowess平滑趋势线'); legend('数据', '平滑趋势', 'Location', 'best');

如果平滑曲线明显不是一条直线，而是曲线，那么强行用线性模型拟合效果就不会好。

5.2 同方差性（方差齐性）快速观察

同方差性是指误差的方差在不同X水平上恒定。在散点图中，如果随着X增大，Y值的波动范围（垂直方向的分散程度）明显变宽或变窄，就提示可能存在异方差性。

% 绘制“残差风格”的图：将数据按X分箱，查看每个箱内Y的波动 edges = linspace(min(T.X1), max(T.X1), 10); % 将X1分成10个区间 [N, edges, bin] = histcounts(T.X1, edges); figure; for i = 1:length(edges)-1 y_in_bin = T.Y(bin == i); if ~isempty(y_in_bin) % 画出每个箱的中位数和四分位距（类似箱线图的简化） median_y = median(y_in_bin); iqr_y = iqr(y_in_bin); x_center = (edges(i) + edges(i+1)) / 2; plot([x_center, x_center], [median_y - 0.5*iqr_y, median_y + 0.5*iqr_y], 'k-', 'LineWidth', 2); hold on; plot(x_center, median_y, 'ro', 'MarkerSize', 8, 'MarkerFaceColor', 'r'); end end xlabel('X1 (分箱后)'); ylabel('Y (箱内中位数及范围)'); title('按X1分箱后Y的分布范围观察'); grid on;

如果这些垂直线段的高度（代表波动范围）随着X中心的变化而有系统地增大或减小，则异方差性可能存在。

5.3 正态性初步筛查

许多统计检验（如t检验）和模型（如线性回归的系数推断）依赖于误差项的正态性假设。我们可以对响应变量Y（或后续的残差）做快速检查。

figure; subplot(1,2,1); qqplot(T.Y); title('Y的Q-Q图'); ylabel('Y的分位数'); xlabel('标准正态分位数'); subplot(1,2,2); probplot('normal', T.Y); title('Y的概率图');

Q-Q图将数据分位数与理论正态分布分位数对比。如果点大致落在一条直线上，则正态性假设大致满足。如果呈现明显的曲线，则偏离正态。概率图是类似的工具。

常见问题：我的数据看起来不正态怎么办？首先，许多模型（特别是线性回归）对于中度偏离正态性是稳健的，尤其是样本量较大时（中心极限定理）。其次，可以考虑对Y进行变换（如Box-Cox变换）。最后，也可以选择不依赖于正态假设的模型，如广义线性模型或非参数方法。

6. 整合三步：一个完整的建模前工作流示例

让我们用一个模拟数据集，把这三步串起来走一遍。假设我们研究学习时间(StudyHours)和考前睡眠(SleepHours)对考试成绩(ExamScore)的影响。

%% 1. 模拟生成一些“脏”数据 rng(42); % 设定随机种子，确保可重复性 n = 100; StudyHours = 10 + 3*randn(n,1); % 正态分布，均值10h，标准差3h SleepHours = 7 + 1.5*randn(n,1); % 正态分布，均值7h，标准差1.5h % 真实的模型关系：分数与学习时间正相关，与睡眠时间在一定范围内正相关 true_score = 50 + 3*StudyHours + 5*SleepHours - 0.3*(SleepHours-7).^2 + 10*randn(n,1); ExamScore = true_score; % 人为制造一些“问题” % a) 加入缺失值 missing_idx = randperm(n, 5); ExamScore(missing_idx(1:3)) = NaN; StudyHours(missing_idx(4:5)) = NaN; % b) 加入异常值 outlier_idx = n+1; StudyHours(outlier_idx) = 30; % 一个异常高的学习时间 SleepHours(outlier_idx) = 7; ExamScore(outlier_idx) = 40; % 对应的分数却不高 n = n + 1; % c) 加入记录错误 error_idx = randi([1, n], 2,1); SleepHours(error_idx(1)) = 20; % 不可能睡20小时 ExamScore(error_idx(2)) = -5; % 分数不可能为负 % 创建数据表 T_raw = table(StudyHours, SleepHours, ExamScore, 'VariableNames', {'StudyHours', 'SleepHours', 'ExamScore'}); disp('原始数据前几行：'); disp(head(T_raw)); %% 2. 第一步：数据质量诊断与清洗 fprintf('\n--- 第一步：数据质量诊断 ---\n'); % 2.1 缺失值诊断 missing_summary = sum(ismissing(T_raw)); disp('缺失值统计：'); disp(missing_summary); % 2.2 处理明显错误（业务逻辑清洗） % 睡眠时间应在合理范围，分数应为正 T_clean = T_raw; T_clean.SleepHours(T_clean.SleepHours > 15 | T_clean.SleepHours < 0) = NaN; T_clean.ExamScore(T_clean.ExamScore < 0 | T_clean.ExamScore > 100) = NaN; % 2.3 异常值检测（IQR方法） vars = {'StudyHours', 'SleepHours', 'ExamScore'}; for i = 1:length(vars) data = T_clean.(vars{i}); data_valid = data(~isnan(data)); Q = quantile(data_valid, [0.25 0.75]); IQR = Q(2) - Q(1); lb = Q(1) - 3 * IQR; % 使用更严格的3倍IQR，因为模拟数据中异常值明显 ub = Q(2) + 3 * IQR; outlier_mask = (data < lb) | (data > ub); fprintf('变量 %s 发现 %d 个异常值（基于3*IQR）.\n', vars{i}, sum(outlier_mask)); % 这里选择将异常值标记为NaN，后续统一处理 T_clean.(vars{i})(outlier_mask) = NaN; end % 2.4 缺失值填补（使用中位数，更稳健） for i = 1:length(vars) col_data = T_clean.(vars{i}); if any(isnan(col_data)) median_val = median(col_data, 'omitnan'); T_clean.(vars{i})(isnan(col_data)) = median_val; fprintf('变量 %s 用中位数 %.2f 填补了缺失值。\n', vars{i}, median_val); end end % 2.5 最终检查 fprintf('\n清洗后数据维度：%d 行 x %d 列\n', height(T_clean), width(T_clean)); disp('清洗后数据摘要：'); disp(summary(T_clean)); %% 3. 第二步：变量关系可视化探索 fprintf('\n--- 第二步：可视化探索 ---\n'); figure('Position', [100, 100, 1200, 800]); % 3.1 单变量分布 subplot(2,3,1); histogram(T_clean.StudyHours, 'FaceColor', [0.2 0.6 0.8]); title('学习时间分布'); xlabel('小时'); ylabel('频数'); subplot(2,3,2); histogram(T_clean.SleepHours, 'FaceColor', [0.8 0.4 0.2]); title('睡眠时间分布'); xlabel('小时'); ylabel('频数'); subplot(2,3,3); histogram(T_clean.ExamScore, 'FaceColor', [0.4 0.8 0.4]); title('考试成绩分布'); xlabel('分数'); ylabel('频数'); % 3.2 双变量关系（散点图与趋势） subplot(2,3,4); scatter(T_clean.StudyHours, T_clean.ExamScore, 30, 'filled'); xlabel('学习时间 (小时)'); ylabel('考试成绩'); title('学习时间 vs 成绩'); grid on; lsline; % 添加最小二乘线 subplot(2,3,5); scatter(T_clean.SleepHours, T_clean.ExamScore, 30, 'filled'); xlabel('睡眠时间 (小时)'); ylabel('考试成绩'); title('睡眠时间 vs 成绩'); grid on; lsline; % 3.3 双变量关系（按学习时间分组看睡眠与成绩） subplot(2,3,6); % 将学习时间离散化为高/低两组 study_median = median(T_clean.StudyHours); high_study = T_clean.StudyHours > study_median; gscatter(T_clean.SleepHours, T_clean.ExamScore, high_study, 'br', 'o*'); xlabel('睡眠时间 (小时)'); ylabel('考试成绩'); legend('低学习时间', '高学习时间', 'Location', 'best'); title('睡眠vs成绩（按学习时间分组）'); grid on; %% 4. 第三步：模型假设初步思考 fprintf('\n--- 第三步：模型假设初步验证 ---\n'); % 从散点图观察： % - StudyHours 与 ExamScore 呈现明显的正相关线性趋势。 % - SleepHours 与 ExamScore 似乎存在曲线关系（倒U型？），线性趋势不明显。 % - 在“睡眠vs成绩”图中，两组（高/低学习时间）的点似乎沿着不同的“曲线”分布，提示可能存在交互作用。 % 计算相关系数 corr_matrix = corr(table2array(T_clean), 'Rows', 'all'); fprintf('\n变量间相关系数矩阵：\n'); disp(array2table(corr_matrix, 'VariableNames', vars, 'RowNames', vars)); % 初步结论： % 1. 学习时间与成绩强相关（r=0.85）。 % 2. 睡眠时间与成绩相关较弱（r=0.21），但散点图提示可能是非线性。 % 3. 因此，建模时不应简单使用 ExamScore ~ StudyHours + SleepHours 的线性模型。 % 应考虑加入 SleepHours 的二次项，以及 StudyHours 与 SleepHours 的交互项。 % 例如：ExamScore ~ StudyHours + SleepHours + SleepHours^2 + StudyHours*SleepHours fprintf('\n“三步醒肤”完成。基于探索性分析，建议的模型形式已初步明确。\n'); fprintf('接下来可以基于 T_clean 数据集，使用 fitlm 函数拟合包含非线性项和交互项的线性模型。\n');

运行这段代码，你会得到一个完整的报告。从输出和图形中，你可以清晰地看到数据如何被清洗，变量间的关系如何，以及应该如何调整你的建模策略。这就是“三步醒肤法”的价值：它让数据自己说话，指导你选择正确的模型，而不是让你拍脑袋决定。

7. 常见问题与排查技巧实录

即使按照上述流程，实践中还是会遇到各种问题。下面是我总结的一些高频问题和解决思路。

7.1 数据清洗中的两难选择

问题1：缺失值太多，是删除变量还是填补？

• 排查：计算每个变量的缺失比例。如果超过40%-50%，通常考虑删除该变量，因为填补会引入过多噪声。在20%-40%之间，需要谨慎，尝试使用多重插补等高级方法（如fillmissing的'movmean'或第三方工具箱），并评估填补后模型的稳定性。低于20%，可以用中位数或模型预测填补。
• 技巧：创建一个“缺失指示器”变量（1表示缺失，0表示未缺失），作为一个新的预测变量加入模型。有时，数据缺失本身（如用户不愿填写收入）可能就是有意义的模式。

问题2：异常值删不掉，一删样本量就锐减？

• 排查：检查异常值是否集中在某个分组或条件下。如果是，可能代表一个特殊的子群体，不应简单删除，而应考虑分层建模或使用混合模型。
• 技巧：使用对异常值不敏感的统计量（如中位数、MAD）或模型（如分位数回归、Huber回归）。在Matlab中，fitlm的'RobustOpts'选项可以启用稳健回归。

7.2 可视化中的陷阱

问题3：散点图点太密，什么都看不出来？

• 排查：数据量过大（>10000点）会导致点重叠严重。
• 技巧：
  1. 抽样：随机抽取一部分数据（如10%）画图。
  2. 透明度：使用scatter(X,Y, 'filled', 'MarkerFaceAlpha', 0.2)降低点的不透明度。
  3. 二维直方图：使用histogram2或scatterhist查看点的密度分布。
  4. 数据抖动：对离散数据或存在大量重复值时，加入微小的随机噪声jitter。

问题4：如何判断非线性关系是二次、指数还是其他？

• 排查：散点图上的趋势线（lsline）是线性的。尝试绘制局部加权散点平滑（smoothdata的'loess'方法），它能更好地揭示局部趋势。
• 技巧：在散点图上叠加不同模型的拟合线。例如，先用polyfit拟合一条二次曲线画上去，与线性趋势对比。如果图形上难以判断，可以后续在模型中分别加入二次项、对数项等，并通过anova函数比较模型拟合优度的提升是否显著。

7.3 模型假设验证的后续步骤

问题5：正式建模后，残差图显示异方差（漏斗形），怎么办？

• 排查：这是线性回归常见问题。使用plotResiduals(mdl, 'fitted')绘制残差与拟合值图。
• 解决思路：
  1. 变换响应变量Y：尝试对Y做对数变换（log(Y)）、平方根变换。这常能稳定方差。
  2. 使用加权最小二乘法：在fitlm中指定'Weights'参数，给方差较小的观测以更高权重。
  3. 改用广义线性模型：如果Y是计数数据（如泊松分布）或比例数据（如二项分布），使用fitglm并指定合适的分布族和连接函数。
  4. 接受并使用稳健标准误：如果不便变换，在报告结果时使用异方差稳健的标准误（Huber-White标准误）来计算置信区间和p值，这可以在后续使用coefTest或其他统计检验函数时指定相关选项。

问题6：Q-Q图显示残差尾部偏离直线，正态性假设不满足，影响大吗？

• 排查：主要关注偏离的严重程度和样本量。轻微的偏态或厚尾，在大样本下（n>100）对参数估计和显著性检验的影响有限（中心极限定理）。严重的非正态性可能影响预测区间。
• 解决思路：
  1. 变换Y或X：Box-Cox变换是专门用于使数据更接近正态分布的方法。
  2. 自助法：使用bootstrp函数进行参数估计，不依赖于正态性假设。
  3. 转向非参数方法：如使用决策树、随机森林等算法，它们对分布没有要求。

7.4 Matlab函数使用小贴士

• fitlm与stepwiselm：fitlm用于拟合指定形式的线性模型。如果你不确定该加入哪些变量，可以使用stepwiselm进行逐步回归，让Matlab基于信息准则（如AIC）自动选择变量。但慎用，因为这是数据驱动的，可能产生过拟合或误导性模型，最好结合业务知识。
• ttest与ttest2的区别：这是热词中提到的问题。ttest用于单样本或配对样本t检验（比较一组数据的均值与某个常数，或比较两组配对数据的差异）。ttest2用于独立双样本t检验（比较两组独立数据的均值）。用错会导致完全错误的结论。
• 模型诊断图：拟合模型mdl后，直接用plot(mdl)会生成四个诊断图（残差vs拟合值、Q-Q图、尺度-位置图、残差vs杠杆值），这是检验线性模型假设的标准流程，务必会看。

记住，统计建模没有一成不变的“正确”流程，只有“更合适”的流程。这“三步醒肤法”提供的是一个稳健的起点，能帮你避开大多数低级错误，建立对数据的直觉，从而构建出更可靠、更有解释力的模型。在Matlab的强大算力背后，真正值钱的，始终是你对数据本身的理解和思考。