Question

我正在学习线性回归并在apache commons数学库3.5中使用OLSMultipleLinearRegression。我运行了以下用于计算回归参数的示例代码。

// estimate weights from heights and waist
OLSMultipleLinearRegression regression = new OLSMultipleLinearRegression();
// weight
double[] y = new double[] { 50, 60, 65, 65, 70, 75, 80, 85, 90, 95 };
// height, waist
double[][] x = new double[10][];
x[0] = new double[] { 165, 65 };
x[1] = new double[] { 170, 68 };
x[2] = new double[] { 172, 70 };
x[3] = new double[] { 175, 65 };
x[4] = new double[] { 170, 80 };
x[5] = new double[] { 172, 85 };
x[6] = new double[] { 183, 78 };
x[7] = new double[] { 187, 79 };
x[8] = new double[] { 180, 95 };
x[9] = new double[] { 185, 97 };
regression.newSampleData(y, x);

double[] coe = regression.estimateRegressionParameters();
for (double p : coe) {
    System.out.println(p);
}

我想计算每一行的估算值，以便绘制它们。但只要我阅读该文档，该库就不提供这样做的方法。当然，这并不困难，可以实现如下，

double[] z = new double[10];
for (int i = 0; i < z.length; i++) {
    z[i] = coe[0];
    for (int j = 1; j < coe.length; j++) z[i] += coe[j] * x[i][j - 1];
}

但是有更好的方法吗？我找到了calculateResiduals，它似乎与我的目的很接近，但它是一种受保护的方法。

Answer 1

我决定通过继承OLSMultipleLinearRegression类来创建MyOLSMultipleLinearRegression类，并实现以下功能。

public double[] calculateEstimatedValues() {
    RealVector b = calculateBeta();
    return getX().operate(b).toArray();
}

通过OLSMultipleLinearRegression计算估计值

1 个答案: