我正试图计算一个有540个观测值的模型中预测因子的方差。从预测器矩阵(X_blocked(中，我想取第二列(第一列用于截距，第二列用于第一个预测器(，计算其方差，并将值存储在变量中：
blocked_pred1_var = np.sum(X_blocked[:, 1] - np.mean(X_blocked[:, 1]) ** 2) / 539
但这给了我一个错误的数字。
括号位置错误。目前，你正在对平均值求平方，然后从你的预测值中减去它。你需要从预测值中减去平均值，然后将其平方
blocked_pred1_var = np.sum((X_blocked[:, 1] - np.mean(X_blocked[:, 1])) ** 2) / 539