本文于2020-10-10更新。 如发现问题或者有建议,欢迎提交 Issue
1 为什么要做归一化?
tbl1 <- tibble(
person = 1:4,
age_yr = c(35,40,35,40),
height_cm = c(190,190,160,160)
)
tbl1
## # A tibble: 4 x 3
## person age_yr height_cm
## <int> <dbl> <dbl>
## 1 1 35 190
## 2 2 40 190
## 3 3 35 160
## 4 4 40 160
tbl1 %>%
ggplot(aes(x = age_yr, y = height_cm, col = as.factor(person))) +
geom_point() +
xlim(0,60) +
ylim(150,200)
看这张图我们会发现1
和2
比较相似,3
和4
比较相似。
但是当我们把身高单位变换后,
tbl2 <- tbl1 %>%
mutate(height_ft = height_cm / 30.48)
tbl2 %>%
ggplot(aes(x = age_yr, y = height_ft, col = as.factor(person))) +
geom_point() +
xlim(35,40) +
ylim(4,7)
我们会发现1
和3
比较相似,2
和4
比较相似,因此我们需要对数据进行标准化。
tbl1 %>%
mutate(age_yr_scaled = scale(age_yr), height_cm = scale(height_cm))
## # A tibble: 4 x 4
## person age_yr height_cm[,1] age_yr_scaled[,1]
## <int> <dbl> <dbl> <dbl>
## 1 1 35 0.866 -0.866
## 2 2 40 0.866 0.866
## 3 3 35 -0.866 -0.866
## 4 4 40 -0.866 0.866
我们会发现他们四个都很相似。
例子参考 Kalisch (2012) 和 Yu (2015), 但是当数据单位一致时,从实际出发,没必要标准化 (Kalisch 2012; Yu 2015),比如宏观数据,都是统一的货币单位时。
2 梯度下降的方式理解
忆臻 (2017) 给出了到位的解释,感觉蛮到位。
就是说,我们假设一个二元回归,
y=β1x1+β2x2+u
损失函数,
J=min∑(y−ˆβ1x1+ˆβ2x2)2
做了归一化的ˆβ1和ˆβ2接近,因此在以ˆβ1和ˆβ2为坐标系的直线坐标系中,图像更接近于圆。
因此梯度下降的走位就不会那么风骚了。
3 max min 标准化
归一化是一种简化计算的方式,即将有量纲的表达式,经过变换,化为无量纲的表达式,成为标量。在多种计算中都经常用到这种方法。
举个例子,将X={x1,x2,...,xn}转化为{˜x1,˜x2,...,˜xn},公式变化为(一种情况)
˜xi=x−min(x)max(X)−min
4 z标准化在python中
可以在python中自行构建代码 (Aruliah 2017)。
˜xi=x−ˉxσ(X)
Python的函数
def zscore(series):
return (series - series.mean()) / series.std()
这里的series
是未知数,可以写成l
、x
等等。
5 StandardScalar
在python中
StandardScalar就是标准化,所谓的标准化就是μ=0和σ=1。
˜x=x−ˉxσx
E(˜x)=E(x−ˉxσx)=E(xσx)−E(ˉxσx)=0
Var(˜x)=Var(x−ˉxσx)=1σ2x⋅Var(x−ˉx)=Var(x)σ2x=1
from sklearn.preprocessing import StandardScalar
引入。
6 聚类需要使用归一化
比如点x(0.1,10000)和y(0.9,9999)表示样本内的两个点,其中有两个特征, 第一个特征的range=[0.1,0.9], 第二个特征的range=[9999,10000], 显然不归一化的话,第一个特征完全被弱化了。 (张腾 2017)
7 标准化与归一化的区别
简单来说,标准化是依照特征矩阵的列处理数据,其通过求z-score的方法,将样本的特征值转换到同一量纲下。归一化是依照特征矩阵的行处理数据,其目的在于样本向量在点乘运算或其他核函数计算相似性时,拥有统一的标准,也就是说都转化为“单位向量”1。规则为l2的归一化公式如下:
˜x=x|x|
|x|2=m∑i=1xi
其中x表示一个行向量,即一个用户的数据。 m表示有m个特征。
参考文献
Aruliah, Dhavide. 2017. “Using Vectorized Functions.” 2017. https://campus.datacamp.com/courses/manipulating-dataframes-with-pandas/extracting-and-transforming-data?ex=16.
Kalisch, Markus. 2012. “Measuring Distances: Applied Multivariate Statistics – Spring 2012.” Eidgenössische Technische Hochschule Zürich. 2012. https://stat.ethz.ch/education/semesters/ss2012/ams/slides/v4.2.pdf.
Yu, Lijia. 2015. “度量距离时对数据做不做比例归一化.” 2015. http://yulijia.net/cn/%E6%89%80%E6%80%9D%E6%89%80%E6%82%9F/2015/11/06/To-scale-or-not-to-scale.html.
张腾. 2017. “聚类时,特征要做归一化/标准化/正则化吗?.” 2017. https://www.zhihu.com/question/58302000/answer/265061045.
忆臻. 2017. “为什么要对数据进行归一化处理?.” 2017. https://mp.weixin.qq.com/s/-Pr7e9tfVI_I2MyX8AXaUA.
如果x2+y2+z2=1,则向量[x,y,z]称为单位向量。 只要模为1的向量,就称为单位向量,单位向量有无穷多个,在任何一个方向上都有一个单位向量。↩