回归分析中的评价方法

以你之姓@ 2022-06-03 07:48 902阅读 0赞

回归（Regression）不同于分类问题，在回归方法中我们预测一系列连续的值，在预测完后有个问题是如何评价预测的结果好坏，关于这个问题目前学术界也没有统一的标准。下面是我在论文中的看到的一些常用方法，希望对有缘人有用。

**1 MAE（Mean Absolute Error）平均绝对差值**

In [statistics][], the **mean absolute error (MAE)** is a quantity used to measure how close forecasts or predictions are to the eventual outcomes. The mean absolute error is given by

![\\mathrm\{MAE\} = \\frac\{1\}\{n\}\\sum\_\{i=1\}^n \\left| f\_i-y\_i\\right| =\\frac\{1\}\{n\}\\sum\_\{i=1\}^n \\left| e\_i \\right|.][mathrm_MAE_ _ _frac_1_n_sum_i_1_n _left_ f_i-y_i_right_ _frac_1_n_sum_i_1_n _left_ e_i _right_.]

As the name suggests, the mean absolute error is an average of the absolute errors ![e\_i = |f\_i - y\_i|][e_i _ _f_i - y_i], where ![f\_i][f_i] is the prediction and ![y\_i][y_i] the true value. Note that alternative formulations may include relative frequencies as weight factors.

注意：MAE与 MAD（Mean Absolute Difference）等价。此外，MAE很容易跟absolute deviation(绝对偏差)混淆，它们的定义很类似，但使用的场景完全不同；绝对偏差针对的就是一组数据，而MAD针对的两组数据（预测值一组，真实值一组）。

**2 MSE(Mean Square Error)均方误差**

If ![\\hat\{Y\}][hat_Y] is a vector of n predictions, and ![Y][] is the vector of the true values, then the (estimated) MSE of the predictor is: ![\\operatorname\{MSE\}=\\frac\{1\}\{n\}\\sum\_\{i=1\}^n(\\hat\{Y\_i\} - Y\_i)^2.][operatorname_MSE_frac_1_n_sum_i_1_n_hat_Y_i_ - Y_i_2.]

**3 RMSE(Root Mean Square error)均方根误差**

RMSE跟RMSD(Root-mean-square deviation)均方根偏差的定义等价，RMSE实际上就是MSE的平方根。

The RMSD of an [estimator][] ![\\hat\{\\theta\}][hat_theta] with respect to an estimated parameter ![\\theta][theta] is defined as the square root of the [mean square error][]:

![\\operatorname\{RMSD\}(\\hat\{\\theta\}) = \\sqrt\{\\operatorname\{MSE\}(\\hat\{\\theta\})\} = \\sqrt\{\\operatorname\{E\}((\\hat\{\\theta\}-\\theta)^2)\}.][operatorname_RMSD_hat_theta_ _ _sqrt_operatorname_MSE_hat_theta_ _ _sqrt_operatorname_E_hat_theta_-_theta_2_.]

For an [unbiased estimator][], the RMSD is the square root of the variance, known as the [standard error][].

The RMSD of predicted values ![\\hat y\_t][hat y_t] for times *t* of a [regression's][regression_s] [dependent variable][] ![y][] is computed for *n* different predictions as the square root of the mean of the squares of the deviations:

![\\operatorname\{RMSD\}=\\sqrt\{\\frac\{\\sum\_\{t=1\}^n (\\hat y\_t - y\_t)^2\}\{n\}\}.][operatorname_RMSD_sqrt_frac_sum_t_1_n _hat y_t - y_t_2_n_.]

In some disciplines, the RMSD is used to compare differences between two things that may vary, neither of which is accepted as the "standard". For example, when measuring the average difference between two time series ![x\_\{1,t\}][x_1_t] and ![x\_\{2,t\}][x_2_t], the formula becomes

![\\operatorname\{RMSD\}= \\sqrt\{\\frac\{\\sum\_\{t=1\}^n (x\_\{1,t\} - x\_\{2,t\})^2\}\{n\}\}.][operatorname_RMSD_ _sqrt_frac_sum_t_1_n _x_1_t_ - x_2_t_2_n_.]  ![38871467][] ![38871467][]

![38871467][]

**4 Normalized root-mean-square deviation归一化均方差跟偏差**

The **normalized root-mean-square deviation** or **error** (**NRMSD** or **NRMSE**) is the RMSD divided by the range of observed values of a variable being predicted,or:

![\\mathrm\{NRMSD\} = \\frac\{\\mathrm\{RMSD\}\}\{y\_\\max -y\_\\min\}][mathrm_NRMSD_ _ _frac_mathrm_RMSD_y_max -y_min]

The value is often expressed as a percentage, where lower values indicate less residual variance.

**5  coefficient of variation of the RMSD**

The **coefficient of variation of the RMSD, CV(RMSD),** or more commonly **CV(RMSE)**, is defined as the RMSD normalized to the [mean][] of the observed values:

![\\mathrm\{CV(RMSD)\} = \\frac \{\\mathrm\{RMSD\}\}\{\\bar y\}.][mathrm_CV_RMSD_ _ _frac _mathrm_RMSD_bar y_.]

It is the same concept as the coefficient of variation except that RMSD replaces the standard deviation.

\-------------------------------------------------------------------------------------------------------------------------------------------

Correlation Coefficient（相关系数）

相关表和相关图可反映两个变量之间的相互关系及其相关方向，但无法确切地表明两个变量之间相关的程度。于是，著名统计学家卡尔·皮尔逊设计了统计指标——相关系数(Correlation coefficient)。相关系数是用以反映变量之间相关关系密切程度的统计指标。相关系数是按积差方法计算，同样以两变量与各自平均值的离差为基础，通过两个离差相乘来反映两变量之间相关程度；着重研究线性的单相关系数。

依据相关现象之间的不同特征，其统计指标的名称有所不同。如将反映两变量间线性相关关系的统计指标称为相关系数（相关系数的平方称为判定系数）；将反映两变量间曲线相关关系的统计指标称为非线性相关系数、非线性判定系数；将反映多元线性相关关系的统计指标称为复相关系数、复判定系数等。

相关关系是一种非确定性的关系，相关系数是研究变量之间线性相关程度的量。由于研究对象的不同，相关系数有如下几种定义方式：

简单相关系数：又叫相关系数或线性相关系数，一般用字母P 表示，用来度量两个变量间的线性关系。

复相关系数：又叫多重相关系数。复相关是指因变量与多个自变量之间的相关关系。例如，某种商品的季节性需求量与其价格水平、职工收入水平等现象之间呈现复相关关系。

典型相关系数：是先对原来各组变量进行主成分分析，得到新的线性关系的综合指标，再通过综合指标之间的线性相关系数来研究原各组变量间相关关系。

下面是两种在回归分析评价中常用的相关系数。

\-------------------------------------------------------------------------------------------------------------------------------------

**6  Pearson's Correlation Coefficient(皮尔逊相关系数)**

**  有的论文里叫COR（相关性）**

相关系数：考察两个事物（在数据里我们称之为变量）之间的相关程度。

如果有两个变量：X、Y，最终计算出的相关系数的含义可以有如下理解：

(1)、当相关系数为0时，X和Y两变量无关系。

(2)、当X的值增大（减小），Y值增大（减小），两个变量为正相关，相关系数在0.00与1.00之间。

(3)、当X的值增大（减小），Y值减小（增大），两个变量为负相关，相关系数在-1.00与0.00之间。

相关系数的绝对值越大，相关性越强，相关系数越接近于1或-1，相关度越强，相关系数越接近于0，相关度越弱。

通常情况下通过以下取值范围判断变量的相关强度：  
相关系数     0.8-1.0     极强相关  
                 0.6-0.8     强相关  
                 0.4-0.6     中等程度相关  
                 0.2-0.4     弱相关  
                 0.0-0.2     极弱相关或无相关

皮尔逊相关也称为积差相关（或积矩相关）是英国统计学家皮尔逊于20世纪提出的一种计算直线相关的方法。

假设有两个变量X、Y，那么两变量间的皮尔逊相关系数可通过以下公式计算：

公式一：

![皮尔逊相关系数计算公式][19961_12788382715I69.gif]

公式二：

![皮尔逊相关系数计算公式][19961_1278838271Lch3.gif]

公式三：

![皮尔逊相关系数计算公式][19961_1278838271OR17.gif]

公式四：

![皮尔逊相关系数计算公式][19961_1278843346L5kz.gif]

以上列出的四个公式等价，其中E是数学期望，cov表示协方差，N表示变量取值的个数。

注意当两个变量的标准差都不为零时，相关系数才有定义，皮尔逊相关系数适用于：

(1)、两个变量之间是线性关系，都是连续数据。

(2)、两个变量的总体是正态分布，或接近正态的单峰分布。

(3)、两个变量的观测值是成对的，每对观测值之间相互独立。

**7 concordance correlation coefficient(一致性相关系数)**

In statistics, the concordance correlation coefficient measures the agreement between two variables, e.g., to evaluate reproducibility or for inter-rater reliability.

**Definition:**

Lawrence Lin has the form of the concordance correlation coefficient ![\\rho\_c][rho_c] as

![\\rho\_c = \\frac\{2\\rho\\sigma\_x\\sigma\_y\}\{\\sigma\_x^2 + \\sigma\_y^2 + (\\mu\_x - \\mu\_y)^2\},][rho_c _ _frac_2_rho_sigma_x_sigma_y_sigma_x_2 _ _sigma_y_2 _ _mu_x - _mu_y_2]

where ![\\mu\_x][mu_x] and ![\\mu\_y][mu_y] are the [means][mean] for the two variables and ![\\sigma^2\_x][sigma_2_x] and ![\\sigma^2\_y][sigma_2_y] are the corresponding [variances][]. ![\\rho][rho] is the correlation coefficient between the two variables.

This follows from its definition[\[1\]][1] as

![\\rho\_c = 1 - \\frac\{\{\\rm Expected\\ orthogonal\\ squared\\ distance\\ from\\ the\\ diagonal\\ \}x=y\}\{\{\\rm Expected\\ orthogonal\\ squared\\ distance\\ from\\ the\\ diagonal\\ \}x=y\{\\rm \\ assuming\\ independence\}\}.][rho_c _ 1 - _frac_rm Expected_ orthogonal_ squared_ distance_ from_ the_ diagonal_ _x_y_rm Expected_ orthogonal_ squared_ distance_ from_ the_ diagonal_ _x_y_rm _ assuming_ independence_.]

When the concordance correlation coefficient is computed on a N-length data set (i.e., two vectors of length *N*) the form is

![\\hat\{\\rho\}\_c = \\frac\{2 s\_\{xy\}\}\{s\_x^2 + s\_y^2 + (\\bar\{x\} - \\bar\{y\})^2\},][hat_rho_c _ _frac_2 s_xy_s_x_2 _ s_y_2 _ _bar_x_ - _bar_y_2]

where the mean is computed as

![\\bar\{x\} = \\frac\{1\}\{N\} \\sum\_\{n=1\}^N x\_n][bar_x_ _ _frac_1_N_ _sum_n_1_N x_n]

and the variance

![s\_x^2 = \\frac\{1\}\{N\} \\sum\_\{n=1\}^N (x\_n - \\bar\{x\})^2][s_x_2 _ _frac_1_N_ _sum_n_1_N _x_n - _bar_x_2]

and the covariance

![s\_\{xy\} = \\frac\{1\}\{N\} \\sum\_\{n=1\}^N (x\_n - \\bar\{x\})(y\_n - \\bar\{y\}) .][s_xy_ _ _frac_1_N_ _sum_n_1_N _x_n - _bar_x_y_n - _bar_y_ .]

Whereas the ordinary [correlation coefficient][] (Pearson's) is immune to whether the biased or unbiased versions for estimation of the variance is used, the concordance correlation coefficient is not. In the original article Lin suggested the 1/N normalization, while in another article Nickerson appears to have used the 1/(N-1), i.e., the concordance correlation coefficient may be computed slightly differently between implementations.

**Relation to other measures of correlation**

The concordance correlation coefficient is nearly identical to some of the measures called [intra-class correlations][], and comparisons of the concordance correlation coefficient with an "ordinary" intraclass correlation on different data sets found only small differences between the two correlations, in one case on the third decimal. It has also been stated that the ideas for concordance correlation coefficient "are quite similar to results already published by Krippendorff in 1970".

In the original article[\[1\]][1] Lin suggested a form for multiple classes (not just 2). Over ten years later a correction to this form was issued.

One example of the use of the concordance correlation coefficient is in a comparison of analysis method for [functional magnetic resonance imaging][] brain scans.

Reference:

http://en.wikipedia.org/wiki/Concordance\_correlation\_coefficient

http://en.wikipedia.org/wiki/Correlation\_coefficient

http://en.wikipedia.org/wiki/Root\_mean\_square\_error

http://blog.csdn.net/wsywl/article/details/5727327

转载自：[/images/20220603/f8b59e54cf984d6cb417471a76d3744e.png][38871467]

#  #

[statistics]: http://en.wikipedia.org/wiki/Statistics
[mathrm_MAE_ _ _frac_1_n_sum_i_1_n _left_ f_i-y_i_right_ _frac_1_n_sum_i_1_n _left_ e_i _right_.]: /images/20220603/874b56214703498aa2f7bd94ee47e5e1.png
[e_i _ _f_i - y_i]: http://upload.wikimedia.org/math/4/7/0/4709dd8136bcaf1abcbc61b8f8a9391d.png
[f_i]: http://upload.wikimedia.org/math/d/b/b/dbb3b84cea7e7015fb1a591964fbd918.png
[y_i]: /images/20220603/e0485d461c094952bdecb5f657398e44.png
[hat_Y]: http://upload.wikimedia.org/math/f/3/e/f3e3835fc00ac33904a281e254e2580f.png
[Y]: http://upload.wikimedia.org/math/5/7/c/57cec4137b614c87cb4e24a3d003a3e0.png
[operatorname_MSE_frac_1_n_sum_i_1_n_hat_Y_i_ - Y_i_2.]: http://upload.wikimedia.org/math/0/9/7/0978fe57785f6ce707efd6bd90550552.png
[estimator]: http://en.wikipedia.org/wiki/Estimator
[hat_theta]: http://upload.wikimedia.org/math/7/4/4/744e7c77ccdce7b7aa6f9a75788ac57d.png
[theta]: http://upload.wikimedia.org/math/5/0/d/50d91f80cbb8feda1d10e167107ad1ff.png
[mean square error]: http://en.wikipedia.org/wiki/Mean_square_error
[operatorname_RMSD_hat_theta_ _ _sqrt_operatorname_MSE_hat_theta_ _ _sqrt_operatorname_E_hat_theta_-_theta_2_.]: /images/20220603/0773c5dc969b40dfadacafe0cee41427.png
[unbiased estimator]: http://en.wikipedia.org/wiki/Unbiased_estimator
[standard error]: http://en.wikipedia.org/wiki/Standard_error_%28statistics%29
[hat y_t]: http://upload.wikimedia.org/math/b/1/8/b1859f53476576f41148b5ef77bafa19.png
[regression_s]: http://en.wikipedia.org/wiki/Regression_analysis
[dependent variable]: http://en.wikipedia.org/wiki/Dependent_variable
[y]: /images/20220603/372391717e904aba916211946c1f7a3a.png
[operatorname_RMSD_sqrt_frac_sum_t_1_n _hat y_t - y_t_2_n_.]: /images/20220603/30a2929476e844a8829ed276a449ef6e.png
[x_1_t]: http://upload.wikimedia.org/math/7/b/2/7b20851421722cf0dc7f6e17867a5786.png
[x_2_t]: /images/20220603/30291813caa543e38ab9d0f51ab7e258.png
[operatorname_RMSD_ _sqrt_frac_sum_t_1_n _x_1_t_ - x_2_t_2_n_.]: http://upload.wikimedia.org/math/3/d/7/3d7a850ae7e65722ba99b04e445df0ff.png
[38871467]: /images/20220603/f8b59e54cf984d6cb417471a76d3744e.png
[mathrm_NRMSD_ _ _frac_mathrm_RMSD_y_max -y_min]: /images/20220603/5c79a01352bd48ba99a0c7e035464e7a.png
[mean]: http://en.wikipedia.org/wiki/Mean
[mathrm_CV_RMSD_ _ _frac _mathrm_RMSD_bar y_.]: /images/20220603/452c2e4730014e119e7592a370b25d31.png
[19961_12788382715I69.gif]: /images/20220603/e4dd7057cca04f778a4bcb364c42adba.png
[19961_1278838271Lch3.gif]: /images/20220603/7037d23d71fa47e9b3f3b6fe7ed119c7.png
[19961_1278838271OR17.gif]: /images/20220603/aa219d83de5946998da49cd875810323.png
[19961_1278843346L5kz.gif]: /images/20220603/fb410bcdfa2e451db978a25ccd9bf6fa.png
[rho_c]: /images/20220603/6538ca34a91a4cf99e566351cd245616.png
[rho_c _ _frac_2_rho_sigma_x_sigma_y_sigma_x_2 _ _sigma_y_2 _ _mu_x - _mu_y_2]: /images/20220603/838cad8eab3448b281746a89bd95ef02.png
[mu_x]: http://upload.wikimedia.org/math/a/4/3/a43b528a736f5a341fdac5ed3bba5c8b.png
[mu_y]: http://upload.wikimedia.org/math/7/f/1/7f1ab53824f9edb858693103e2d07bc6.png
[sigma_2_x]: http://upload.wikimedia.org/math/2/f/6/2f637941b163c57230852e75ab4e5ec0.png
[sigma_2_y]: http://upload.wikimedia.org/math/b/f/0/bf0c9124b79c782f97a5082583727982.png
[variances]: http://en.wikipedia.org/wiki/Variance
[rho]: http://upload.wikimedia.org/math/f/7/f/f7f177957cf064a93e9811df8fe65ed1.png
[1]: http://en.wikipedia.org/wiki/Concordance_correlation_coefficient#cite_note-LinL1989Concordance-1
[rho_c _ 1 - _frac_rm Expected_ orthogonal_ squared_ distance_ from_ the_ diagonal_ _x_y_rm Expected_ orthogonal_ squared_ distance_ from_ the_ diagonal_ _x_y_rm _ assuming_ independence_.]: /images/20220603/da78d7a258444d259b8189752e224ee2.png
[hat_rho_c _ _frac_2 s_xy_s_x_2 _ s_y_2 _ _bar_x_ - _bar_y_2]: /images/20220603/da750e4b02da4d12bf9370cef32961c1.png
[bar_x_ _ _frac_1_N_ _sum_n_1_N x_n]: /images/20220603/7506dd65c3d242dfa2c10777966c5313.png
[s_x_2 _ _frac_1_N_ _sum_n_1_N _x_n - _bar_x_2]: /images/20220603/34d95b0f922047b7966a817bfe60b04e.png
[s_xy_ _ _frac_1_N_ _sum_n_1_N _x_n - _bar_x_y_n - _bar_y_ .]: /images/20220603/f799a4a474dc4092b2975e97b0f4c093.png
[correlation coefficient]: http://en.wikipedia.org/wiki/Pearson_product-moment_correlation_coefficient
[intra-class correlations]: http://en.wikipedia.org/wiki/Intra-class_correlation
[functional magnetic resonance imaging]: http://en.wikipedia.org/wiki/Functional_magnetic_resonance_imaging