匿名化のためのPCA
Pythonで学ぶデータプライバシーと匿名化
Rebeca Gonzalez
Data engineer
主成分分析（PCA）
$$ $$ 大規模データの次元を削減するためによく使われる次元削減法。
PCAによるデータマスキング
PCAは元の特徴量の線形変換から新しい「主成分」を作る。
ビールのデータセットでは、PCAは新しい特徴量を構成する。例:
$$ 2\times AlcoholicVolume - BitternessLevel$$
PCAによるデータマスキング
データの新しい射影
PCAによるデータマスキング
次元削減なしのPCA
$$
- 元の空間の回転にすぎない。
- 距離は保たれる。
- 予測タスクやアルゴリズムに有利。
PCAによるデータマスキング
- 結果の値を説明なしで公開しても、アルゴリズムはそれで学習し高精度で予測できる。
- 攻撃者はそのマスク済み値の解釈ができない。
PCAによるデータマスキング
# Explore the dataset
heart_df.head()
age sex cp trestbps chol fbs restecg thalach exang oldpeak slope ca thal target
0 63 1 3 145 233 1 0 150 0 2.3 0 0 1 1
1 37 1 2 130 250 0 1 187 0 3.5 0 0 2 1
2 41 0 1 130 204 0 0 172 0 1.4 2 0 2 1
3 56 1 1 120 236 0 1 178 0 0.8 2 0 2 1
4 57 0 0 120 354 0 1 163 1 0.6 2 0 2 1
PCAとScikit-learnによるデータマスキング
# Obtain the data without the target column x_data = df.drop(['target'], axis = 1)# Target column as array of values y = df.target.values
PCAとScikit-learnによるデータマスキング
# Import PCA from Scikit-learn from sklearn.decomposition import PCA# Initialize PCA with number of components to be the same as the number of columns pca = PCA(n_components=len(x_data.columns))# Apply PCA to the data x_data_pca = pca.fit_transform(x_data)
PCAとScikit-learnによるデータマスキング
# See the data
x_data_pca
array([[-1.22673448e+01, 2.87383781e+00, 1.49698788e+01, ...,
7.31102828e-01, -2.90393586e-01, 5.12575925e-01],
[ 2.69013712e+00, -3.98713736e+01, 8.77882303e-01, ...,
4.04206943e-01, -4.25920179e-01, -1.48124511e-01],
[-4.29502141e+01, -2.36368199e+01, 1.75944589e+00, ...,
-9.15397287e-01, 2.17828257e-01, 7.97593843e-02],
...,
PCAとScikit-learnによるデータマスキング
# Create a DataFrame from the resulting PCA transformed data df_x_data_pca = pd.DataFrame(x_data_pca)# Inspect the shape of the dataset df_x_data_pca.shape
(1213, 13)
PCAマスキング後のデータ有用性
- ロジスティック回帰で分類し、元データと変換後データの精度低下を比較する。
$$
- ロジスティック回帰は分類アルゴリズムで、独立変数から二値の結果を予測する。
PCAマスキング後のデータ有用性
ロジスティック回帰で分類し、精度低下を確認する。
# Split the resulting dataset into training and test data x_train, x_test, y_train, y_test = train_test_split(x_data_pca, y, test_size=0.2)# Create the model lr = LogisticRegression(max_iter=200)# Fit train the model lr.fit(x_train,y_train)# Run the model and perform predictions to obtain accuracy score acc = lr.score(x_test, y_test) * 100 print("Test Accuracy is ", acc)
Test Accuracy is 85.24590163934425
PCAマスキング前のデータ有用性
ロジスティック回帰で分類し、元データでのスコアを確認する。
# Split the resulting dataset into training and test data
x_train, x_test, y_train, y_test = train_test_split(x_data.to_numpy(),y,test_size = 0.2)
# Create the model
lr = LogisticRegression(max_iter=200)
# Fit train the model
lr.fit(x_train,y_train)
# Run the model and perform predictions to obtain accuracy score
acc = lr.score(x_test,y_test) * 100
print("Test Accuracy is ", acc)
Test Accuracy is 85.24590163934425
Ayo berlatih!
Pythonで学ぶデータプライバシーと匿名化
