如何使用 sklearn 優(yōu)雅地進(jìn)行數(shù)據(jù)挖掘!
共 14848字,需瀏覽 30分鐘
·
2024-06-26 07:53
??我的小冊(cè) 45章教程:(小白零基礎(chǔ)用Python量化股票分析小冊(cè)) ,原價(jià)299,限時(shí)特價(jià)2杯咖啡,滿100人漲10元。
轉(zhuǎn)自:網(wǎng)絡(luò)
一、使用sklearn數(shù)據(jù)挖掘
def fit(self, X, y=None):
"""Do nothing and return the estimator unchanged
This method is just there to implement the usual API and hence
work in pipelines."""
X = check_array(X, accept_sparse='csr')
return self
在此,我們?nèi)匀皇褂肐RIS數(shù)據(jù)集來(lái)進(jìn)行說(shuō)明。為了適應(yīng)提出的場(chǎng)景,對(duì)原數(shù)據(jù)集需要稍微加工:
ffrom numpy import hstack, vstack, array, median, nan
from numpy.random import choice
from sklearn.datasets import load_iris
iris = load_iris()
#特征矩陣加工
#使用vstack增加一行含缺失值的樣本(nan, nan, nan, nan)
#使用hstack增加一列表示花的顏色(0-白、1-黃、2-紅),花的顏色是隨機(jī)的,意味著顏色并不影響花的分類
iris.data = hstack((choice([0, 1, 2], size=iris.data.shape[0]+1).reshape(-1,1), vstack((iris.data, array([nan, nan, nan, nan]).reshape(1,-1)))))
#目標(biāo)值向量加工
#增加一個(gè)目標(biāo)值,對(duì)應(yīng)含缺失值的樣本,值為眾數(shù)
iris.target = hstack((iris.target, array([median(iris.target)])))
并行處理,流水線處理,自動(dòng)化調(diào)參,持久化是使用sklearn優(yōu)雅地進(jìn)行數(shù)據(jù)挖掘的核心。并行處理和流水線處理將多個(gè)特征處理工作,甚至包括模型訓(xùn)練工作組合成一個(gè)工作(從代碼的角度來(lái)說(shuō),即將多個(gè)對(duì)象組合成了一個(gè)對(duì)象)。
在組合的前提下,自動(dòng)化調(diào)參技術(shù)幫我們省去了人工調(diào)參的反鎖。訓(xùn)練好的模型是貯存在內(nèi)存中的數(shù)據(jù),持久化能夠?qū)⑦@些數(shù)據(jù)保存在文件系統(tǒng)中,之后使用時(shí)無(wú)需再進(jìn)行訓(xùn)練,直接從文件系統(tǒng)中加載即可。
二、并行處理
pipeline包提供了FeatureUnion類來(lái)進(jìn)行整體并行處理:
from numpy import log1p
from sklearn.preprocessing import FunctionTransformer
from sklearn.preprocessing import Binarizer
from sklearn.pipeline import FeatureUnion
#新建將整體特征矩陣進(jìn)行對(duì)數(shù)函數(shù)轉(zhuǎn)換的對(duì)象
step2_1 = ('ToLog', FunctionTransformer(log1p))
#新建將整體特征矩陣進(jìn)行二值化類的對(duì)象
step2_2 = ('ToBinary', Binarizer())
#新建整體并行處理對(duì)象
#該對(duì)象也有fit和transform方法,fit和transform方法均是并行地調(diào)用需要并行處理的對(duì)象的fit和transform方法
#參數(shù)transformer_list為需要并行處理的對(duì)象列表,該列表為二元組列表,第一元為對(duì)象的名稱,第二元為對(duì)象
step2 = ('FeatureUnion', FeatureUnion(transformer_list=[step2_1, step2_2]))
from sklearn.pipeline import FeatureUnion, _fit_one_transformer, _fit_transform_one, _transform_one
from sklearn.externals.joblib import Parallel, delayed
from scipy import sparse
import numpy as np
#部分并行處理,繼承FeatureUnion
class FeatureUnionExt(FeatureUnion):
#相比FeatureUnion,多了idx_list參數(shù),其表示每個(gè)并行工作需要讀取的特征矩陣的列
def __init__(self, transformer_list, idx_list, n_jobs=1, transformer_weights=None):
self.idx_list = idx_list
FeatureUnion.__init__(self, transformer_list=map(lambda trans:(trans[0], trans[1]), transformer_list), n_jobs=n_jobs, transformer_weights=transformer_weights)
#由于只部分讀取特征矩陣,方法fit需要重構(gòu)
def fit(self, X, y=None):
transformer_idx_list = map(lambda trans, idx:(trans[0], trans[1], idx), self.transformer_list, self.idx_list)
transformers = Parallel(n_jobs=self.n_jobs)(
#從特征矩陣中提取部分輸入fit方法
delayed(_fit_one_transformer)(trans, X[:,idx], y)
for name, trans, idx in transformer_idx_list)
self._update_transformer_list(transformers)
return self
#由于只部分讀取特征矩陣,方法fit_transform需要重構(gòu)
def fit_transform(self, X, y=None, **fit_params):
transformer_idx_list = map(lambda trans, idx:(trans[0], trans[1], idx), self.transformer_list, self.idx_list)
result = Parallel(n_jobs=self.n_jobs)(
#從特征矩陣中提取部分輸入fit_transform方法
delayed(_fit_transform_one)(trans, name, X[:,idx], y,
self.transformer_weights, **fit_params)
for name, trans, idx in transformer_idx_list)
Xs, transformers = zip(*result)
self._update_transformer_list(transformers)
if any(sparse.issparse(f) for f in Xs):
Xs = sparse.hstack(Xs).tocsr()
else:
Xs = np.hstack(Xs)
return Xs
#由于只部分讀取特征矩陣,方法transform需要重構(gòu)
def transform(self, X):
transformer_idx_list = map(lambda trans, idx:(trans[0], trans[1], idx), self.transformer_list, self.idx_list)
Xs = Parallel(n_jobs=self.n_jobs)(
#從特征矩陣中提取部分輸入transform方法
delayed(_transform_one)(trans, name, X[:,idx], self.transformer_weights)
for name, trans, idx in transformer_idx_list)
if any(sparse.issparse(f) for f in Xs):
Xs = sparse.hstack(Xs).tocsr()
else:
Xs = np.hstack(Xs)
return Xs
from numpy import log1p
from sklearn.preprocessing import OneHotEncoder
from sklearn.preprocessing import FunctionTransformer
from sklearn.preprocessing import Binarizer
#新建將部分特征矩陣進(jìn)行獨(dú)熱編碼的對(duì)象
step2_1 = ('OneHotEncoder', OneHotEncoder(sparse=False))
#新建將部分特征矩陣進(jìn)行對(duì)數(shù)函數(shù)轉(zhuǎn)換的對(duì)象
step2_2 = ('ToLog', FunctionTransformer(log1p))
#新建將部分特征矩陣進(jìn)行二值化類的對(duì)象
step2_3 = ('ToBinary', Binarizer())
#新建部分并行處理對(duì)象
#參數(shù)transformer_list為需要并行處理的對(duì)象列表,該列表為二元組列表,第一元為對(duì)象的名稱,第二元為對(duì)象
#參數(shù)idx_list為相應(yīng)的需要讀取的特征矩陣的列
step2 = ('FeatureUnionExt', FeatureUnionExt(transformer_list=[step2_1, step2_2, step2_3], idx_list=[[0], [1, 2, 3], [4]]))
三、流水線處理
from numpy import log1p
from sklearn.preprocessing import Imputer
from sklearn.preprocessing import OneHotEncoder
from sklearn.preprocessing import FunctionTransformer
from sklearn.preprocessing import Binarizer
from sklearn.preprocessing import MinMaxScaler
from sklearn.feature_selection import SelectKBest
from sklearn.feature_selection import chi2
from sklearn.decomposition import PCA
from sklearn.linear_model import LogisticRegression
from sklearn.pipeline import Pipeline
#新建計(jì)算缺失值的對(duì)象
step1 = ('Imputer', Imputer())
#新建將部分特征矩陣進(jìn)行定性特征編碼的對(duì)象
step2_1 = ('OneHotEncoder', OneHotEncoder(sparse=False))
#新建將部分特征矩陣進(jìn)行對(duì)數(shù)函數(shù)轉(zhuǎn)換的對(duì)象
step2_2 = ('ToLog', FunctionTransformer(log1p))
#新建將部分特征矩陣進(jìn)行二值化類的對(duì)象
step2_3 = ('ToBinary', Binarizer())
#新建部分并行處理對(duì)象,返回值為每個(gè)并行工作的輸出的合并
step2 = ('FeatureUnionExt', FeatureUnionExt(transformer_list=[step2_1, step2_2, step2_3], idx_list=[[0], [1, 2, 3], [4]]))
#新建無(wú)量綱化對(duì)象
step3 = ('MinMaxScaler', MinMaxScaler())
#新建卡方校驗(yàn)選擇特征的對(duì)象
step4 = ('SelectKBest', SelectKBest(chi2, k=3))
#新建PCA降維的對(duì)象
step5 = ('PCA', PCA(n_components=2))
#新建邏輯回歸的對(duì)象,其為待訓(xùn)練的模型作為流水線的最后一步
step6 = ('LogisticRegression', LogisticRegression(penalty='l2'))
#新建流水線處理對(duì)象
#參數(shù)steps為需要流水線處理的對(duì)象列表,該列表為二元組列表,第一元為對(duì)象的名稱,第二元為對(duì)象
pipeline = Pipeline(steps=[step1, step2, step3, step4, step5, step6])
四、自動(dòng)化調(diào)參
from sklearn.grid_search import GridSearchCV
iris = load_iris()
#新建網(wǎng)格搜索對(duì)象
#第一參數(shù)為待訓(xùn)練的模型
#param_grid為待調(diào)參數(shù)組成的網(wǎng)格,字典格式,鍵為參數(shù)名稱(格式“對(duì)象名稱__子對(duì)象名稱__參數(shù)名稱”),值為可取的參數(shù)值列表
grid_search = GridSearchCV(pipeline, param_grid={'FeatureUnionExt__ToBinary__threshold':[1.0, 2.0, 3.0, 4.0], 'LogisticRegression__C':[0.1, 0.2, 0.4, 0.8]})
#訓(xùn)練以及調(diào)參
grid_search.fit(iris.data, iris.target)
五、持久化
#持久化數(shù)據(jù)
#第一個(gè)參數(shù)為內(nèi)存中的對(duì)象
#第二個(gè)參數(shù)為保存在文件系統(tǒng)中的名稱
#第三個(gè)參數(shù)為壓縮級(jí)別,0為不壓縮,3為合適的壓縮級(jí)別
dump(grid_search, 'grid_search.dmp', compress=3)
#從文件系統(tǒng)中加載數(shù)據(jù)到內(nèi)存中
grid_search = load('grid_search.dmp')
?回顧
最后推薦一下我們團(tuán)隊(duì)寫(xiě)的量化小冊(cè)的內(nèi)容,45篇內(nèi)容!從Python安裝,入門(mén),數(shù)據(jù)分析,爬取股票基金的歷史+實(shí)時(shí)數(shù)據(jù),以及如何寫(xiě)一個(gè)簡(jiǎn)單量化策略,策略回測(cè),如何看資金曲線統(tǒng)統(tǒng)都有介紹!非常超值!
歡迎訂閱:原價(jià)299 早鳥(niǎo)價(jià)2杯咖啡錢(qián),即可永久閱讀。滿400人又要漲價(jià)了,現(xiàn)在的價(jià)格非常非常低,只要2杯奶茶,就可以終身訂閱+課程源碼,還有永久陪伴群。48小時(shí)無(wú)理由退款,放心食用!
往期推薦
量化: 如何用Python爬取創(chuàng)業(yè)板歷史+實(shí)時(shí)股票數(shù)據(jù)!|實(shí)戰(zhàn)股票分析篇利用Pandas 9招挖掘五糧液股價(jià)!|實(shí)戰(zhàn)股票數(shù)據(jù)分析篇 Pandas滾動(dòng)操作 |量化股票第一步,用Python畫(huà)股票K線,雙均線圖,可視化你的股票數(shù)據(jù)!|如何用Python爬取全部800多只ETF基金數(shù)據(jù)!|如何用Python寫(xiě)一個(gè)雙均線策略 |如何用Python開(kāi)發(fā)一個(gè)多策略機(jī)器人!上篇!|Python量化系列-用布林策略買五糧液能賺多少錢(qián)?|只要4秒鐘!用Python 獲取上證指數(shù)34年的歷史日線數(shù)據(jù)!
入門(mén): 最全的零基礎(chǔ)學(xué)Python的問(wèn)題 | 零基礎(chǔ)學(xué)了8個(gè)月的Python | 實(shí)戰(zhàn)項(xiàng)目 |學(xué)Python就是這條捷徑
干貨:爬取豆瓣短評(píng),電影《后來(lái)的我們》 | 38年NBA最佳球員分析 | 從萬(wàn)眾期待到口碑撲街!唐探3令人失望 | 笑看新倚天屠龍記 | 燈謎答題王 |用Python做個(gè)海量小姐姐素描圖 |碟中諜這么火,我用機(jī)器學(xué)習(xí)做個(gè)迷你推薦系統(tǒng)電影
趣味:彈球游戲 | 九宮格 | 漂亮的花 | 兩百行Python《天天酷跑》游戲!
AI: 會(huì)做詩(shī)的機(jī)器人 | 給圖片上色 | 預(yù)測(cè)收入 | 碟中諜這么火,我用機(jī)器學(xué)習(xí)做個(gè)迷你推薦系統(tǒng)電影
小工具: Pdf轉(zhuǎn)Word,輕松搞定表格和水印! | 一鍵把html網(wǎng)頁(yè)保存為pdf!| 再見(jiàn)PDF提取收費(fèi)! | 用90行代碼打造最強(qiáng)PDF轉(zhuǎn)換器,word、PPT、excel、markdown、html一鍵轉(zhuǎn)換 | 制作一款釘釘?shù)蛢r(jià)機(jī)票提示器! |60行代碼做了一個(gè)語(yǔ)音壁紙切換器天天看小姐姐!
