加权最小二乘法推导

### 加权最小二乘法的概念 加权最小二乘法(Weighted Least Squares, WLS)是一种改进的最小二乘法,用于处理观测数据存在异方差的情况。在标准最小二乘法中,假设误差项具有相同的方差(同方差性),但在实际应用中,这一假设往往无法满足。为了克服这一局限性,WLS通过引入权重矩阵来调整不同观测点的重要性,从而提高模型拟合精度。 具体而言,在加权最小二乘法中,目标是最小化带权重的目标函数: \[ S(x) = (b' - A'x)^T W (b' - A'x), \] 其中 \(W\) 是一个对角矩阵,其对角元素代表各观测点对应的权重[^3]。 --- ### 加权最小二乘法的原理 #### 1. 权重的选择 权重通常基于观测值的可靠性或噪声水平选取。如果某个观测值受到较大的噪声影响,则赋予较小的权重;反之亦然。这种机制允许模型更加关注可靠的数据点,减少异常值的影响。 #### 2. 迭代优化过程 对于某些复杂场景,初始权重可能未知或者难以精确指定。此时可以通过迭代方式逐步更新权重和参数估计值。例如,在直线拟合问题中,可以先忽略权重完成初步拟合,再依据残差大小动态调整权重并重复计算直到收敛[^1]。 #### 3. 数学推导 设原始线性回归模型为 \(y_i = a_0 + a_1 x_i + e_i\),其中 \(e_i\) 表示误差项。若已知每条记录对应的标准偏差 \(\sigma_i\),那么定义第 i 个样本点上的相对重要程度 w_i=\frac{1}{σ²\_i} 。于是新的损失函数变为如下形式: \[ L(a)=∑w\_i(y\_i-(a\_0+a\_1*x\_i))^2. \] 通过对该表达式求偏导数等于零可得最优解公式。 --- ### 加权最小二乘法实现方法 以下是 Python 中使用 NumPy 库实现简单的一维加权最小二乘的例子: ```python import numpy as np def weighted_least_squares(X, Y, weights): """ Perform Weighted Least Squares regression. Parameters: X : array-like of shape (n_samples,) Independent variable values. Y : array-like of shape (n_samples,) Dependent variable values. weights : array-like of shape (n_samples,) Weights associated with each sample point. Returns: coefficients : ndarray of shape (2,) Coefficients [intercept, slope]. """ # Convert inputs to column vectors and ensure proper shapes X = np.array(X).reshape(-1, 1) ones_column = np.ones((X.shape[0], 1)) design_matrix = np.hstack([ones_column, X]) sqrt_weights = np.sqrt(weights).reshape(-1, 1) W = np.diag(sqrt_weights.flatten()) Wy = W @ Y.reshape(-1, 1) WX = W @ design_matrix beta_hat = np.linalg.inv(WX.T @ WX) @ WX.T @ Wy return beta_hat.flatten() # Example usage if __name__ == "__main__": X = np.array([1, 2, 3, 4]) # Predictor variables Y = np.array([2.1, 4.9, 7.8, 9.5]) # Response variables weights = np.array([1, 2, 3, 4]) # Associated weights result = weighted_least_squares(X, Y, weights) print(f"Coefficients: intercept={result[0]}, slope={result[1]}") ``` 此代码片段展示了如何构建设计矩阵以及利用平方根变换后的权重构造最终解决方案的过程。 ---
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值