上海做网站推荐,成都网站建设网站建设,北京永辉超市有限公司,网站建设前期准备乳腺癌良性恶性预测 1. 特征工程1.1 特征筛选1.2 特征降维 PCA1.3 SVD奇异值分解 2. 代码2.1 逻辑回归、二分类问题2.2 特征降维 PCA2.3 SVD奇异值分解 1. 特征工程
专业上#xff1a;30个人特征来自于临床一线专家#xff0c;每个特征和都有医学内涵#xff1b;数据上30个人特征来自于临床一线专家每个特征和都有医学内涵数据上30个中可能有一些是冗余的。
1.1 特征筛选
从30个选出重要的踢掉不重要的。从数据的角度来讲不建议进行特征筛选做特征筛选可能只是为了解释性给领导看看而已。因为再不重要的特征都有信息筛选掉必然导致信息的丢失。线性回归的系数代表该项特征的重要性系数就代表权重。使用协方差或皮尔逊系数判断两列特征的相关性。
1.2 特征降维 PCA
把原来30个特征中的核心信息中抽取出来融合到新生成的几个特征中新的特征不是原来的任何一个PCAPrincipal Component Analysis主成分分析
1.3 SVD奇异值分解
SVD奇异值分解是一种在信号处理、统计学、机器学习等领域广泛应用的矩阵分解方法。它可以将一个矩阵分解为三个特定的矩阵的乘积这三个矩阵分别是一个正交矩阵、一个对角矩阵其元素为奇异值且按从大到小排列以及另一个正交矩阵的转置。
2. 代码
2.1 逻辑回归、二分类问题 y F ( X ) s i g m o i d ( x 0 w 0 x 1 w 1 x 2 w 2 . . . x 12 w 12 ) yF(X)sigmoid(x_0w_0x_1w_1x_2w_2...x_{12}w_{12}) yF(X)sigmoid(x0w0x1w1x2w2...x12w12)
from sklearn.datasets import load_breast_cancer
X,yload_breast_cancer(return_X_yTrue)
print(X.shape,y.shape) #(569, 30) (569,)from sklearn.model_selection import train_test_split
X_train,X_test,y_train,y_test train_test_split(X,y)
# 逻辑回归实际是二分类就是线性回归加了一个激活函数 sigmoid
from sklearn.linear_model import LogisticRegression
lrLogisticRegression(max_iter10000)
lr.fit(X_train,y_train)
predlr.predict(X_test)
acc(predy_test).mean()
print(acc)2.2 特征降维 PCA
# 特征降维 PCA主成分分析
from sklearn.decomposition import PCA
pca PCA(n_components20)
pca.fit(X_train)
X_train1 pca.transform(X_train)# 原始数据
lr LogisticRegression(max_iter10000)
lr.fit(X_train,y_train)
(lr.predict(X_test)y_test).mean() #0.965034965034965# 10个特征
lr1 LogisticRegression(max_iter10000)
# 特征降维
# 1. 构建
pcaPCA(n_components10)
# 2. 拟合
pca.fit(X_train)
# 3. 转换
X_train1 pca.transform(X_train)
X_test1 pca.transform(X_test)lr1.fit(X_train1,y_train)
(lr1.predict(X_test1)y_test).mean() #0.9580419580419582.3 SVD奇异值分解
import numpy as np
from matplotlib import pyplot as plt
beauty plt.imread(fnamebeauty.jpg)
# [Height, Width, Channel]
print(beauty.shape) #:(627, 481, 3)
rbeauty[:,:,0]
plt.imshow(r,cmapgray)U, S, V np.linalg.svd(ar, full_matricesFalse)
print(U.shape,S.shape, V.shape) #: (627, 481) (481,) (481, 481)
plt.plot(S)降维之后的效果K 20
restore U[:,:K] np.diag(S[:K]) V[:K, :]
plt.imshow(Xrestore, cmapgray)
