Python 機(jī)器學(xué)習(xí)之 SVM 預(yù)測(cè)買(mǎi)賣(mài)（標(biāo)的物：比特幣）

andong777 發(fā)布于2019-07-30 14:55 / 1605人閱讀

摘要：入門(mén)簡(jiǎn)單策略機(jī)器學(xué)習(xí)庫(kù)的使用什么是支持向量機(jī)。就是試圖把棍放在最佳位置，好讓在棍的兩邊有盡可能大的間隙。現(xiàn)在即使魔鬼放了更多的球，棍仍然是一個(gè)好的分界線(xiàn)。魔鬼看到大俠已經(jīng)學(xué)會(huì)了一個(gè)，于是魔鬼給了大俠一個(gè)新的挑戰(zhàn)。

Python入門(mén)簡(jiǎn)單策略 sklearn 機(jī)器學(xué)習(xí)庫(kù)的使用

什么是 SVM ？

支持向量機(jī)/support vector machine (SVM)。

當(dāng)然首先看一下wiki.

Support Vector Machines are learning models used for classification: which individuals in a population belong where? So… how do SVM and the mysterious “kernel” work?

好吧，故事是這樣子的：

在很久以前的情人節(jié)，大俠要去救他的愛(ài)人，但魔鬼和他玩了一個(gè)游戲。

魔鬼在桌子上似乎有規(guī)律放了兩種顏色的球，說(shuō)：“你用一根棍分開(kāi)它們？要求：盡量在放更多球之后，仍然適用。”

于是大俠這樣放，干的不錯(cuò)？

然后魔鬼，又在桌上放了更多的球，似乎有一個(gè)球站錯(cuò)了陣營(yíng)。

SVM就是試圖把棍放在最佳位置，好讓在棍的兩邊有盡可能大的間隙。

現(xiàn)在即使魔鬼放了更多的球，棍仍然是一個(gè)好的分界線(xiàn)。

然后，在SVM 工具箱中有另一個(gè)更加重要的 trick。魔鬼看到大俠已經(jīng)學(xué)會(huì)了一個(gè)trick，于是魔鬼給了大俠一個(gè)新的挑戰(zhàn)。

現(xiàn)在，大俠沒(méi)有棍可以很好幫他分開(kāi)兩種球了，現(xiàn)在怎么辦呢？當(dāng)然像所有武俠片中一樣大俠桌子一拍，球飛到空中。然后，憑借大俠的輕功，大俠抓起一張紙，插到了兩種球的中間。

現(xiàn)在，從魔鬼的角度看這些球，這些球看起來(lái)像是被一條曲線(xiàn)分開(kāi)了。

再之后，無(wú)聊的大人們，把這些球叫做「data」，把棍子叫做「classifier」, 最大間隙trick 叫做「optimization」，拍桌子叫做「kernelling」, 那張紙叫做「hyperplane」。

參考：

Please explain Support Vector Machines (SVM) like I am a 5 year old. : MachineLearning

Support Vector Machines explained well

https://www.youtube.com/watch...

[以上例子引用自網(wǎng)絡(luò)，侵刪]

一個(gè) SVM 應(yīng)用于 bitcoin trading 的 DEMO

回測(cè)系統(tǒng)自帶的庫(kù)有
實(shí)現(xiàn)語(yǔ)言 Python 2

numpy pandas TA-Lib scipy statsmodels sklearn cvxopt hmmlearn pykalman arch matplotlib

實(shí)盤(pán)需要在托管者所在機(jī)器安裝策略需要的庫(kù)

OK , Talk is cheap, Show code to you!

from sklearn import svm
import numpy as np

def main():
    preTime = 0
    n = 0
    success = 0
    predict = None
    pTime = None
    marketPosition = 0
    initAccount = exchange.GetAccount()
    Log("Running...")
    while True:
        r = exchange.GetRecords()
        if len(r) < 60:
            continue
        bar = r[len(r)-1]
        if bar.Time > preTime:
            preTime = bar.Time
            if pTime is not None and r[len(r)-2].Time == pTime:
                diff = r[len(r)-2].Close - r[len(r)-3].Close
                if diff > SpreadVal:
                    success += 1 if predict == 0 else 0
                elif diff < -SpreadVal:
                    success += 1 if predict == 1 else 0
                else:
                    success += 1 if predict == 2 else 0
                pTime = None
                LogStatus("預(yù)測(cè)次數(shù)", n, "成功次數(shù)", success, "準(zhǔn)確率:", "%.3f %%" % round(float(success) * 100 / n, 2))
        else:
            Sleep(1000)
            continue
        inputs_X, output_Y = [], []
        sets = [None, None, None]
        for i in xrange(1, len(r)-2, 1):
            inputs_X.append([r[i].Open, r[i].Close])
            Y = 0
            diff = r[i+1].Close - r[i].Close
            if diff > SpreadVal:
                Y = 0
                sets[0] = True
            elif diff < -SpreadVal:
                Y = 1
                sets[1] = True
            else:
                Y = 2
                sets[2] = True
            output_Y.append(Y)
        if None in sets:
            Log("樣本不足, 無(wú)法預(yù)測(cè) ...")
            continue
        n += 1
        clf = svm.LinearSVC()
        clf.fit(inputs_X, output_Y)
        predict = clf.predict(np.array([bar.Open, bar.Close]).reshape((1, -1)))
        pTime = bar.Time
        
        Log("預(yù)測(cè)當(dāng)前Bar結(jié)束:", bar.Time, ["漲", "跌", "橫"][predict])
        if marketPosition == 0:
            if predict == 0:
                exchange.Buy(initAccount.Balance/2)
                marketPosition = 1
            elif predict == 1:
                exchange.Sell(initAccount.Stocks/2)
                marketPosition = -1
        else:
            nowAccount = exchange.GetAccount()
            if marketPosition > 0 and predict != 0:
                exchange.Sell(nowAccount.Stocks - initAccount.Stocks)
                nowAccount = exchange.GetAccount()
                marketPosition = 0
            elif marketPosition < 0 and predict != 1:
                while True:
                    dif = initAccount.Stocks - nowAccount.Stocks
                    if dif < 0.01:
                        break
                    ticker = exchange.GetTicker()
                    exchange.Buy(ticker.Sell + (ticker.Sell-ticker.Buy)*2, dif)
                    while True:
                        Sleep(1000)
                        orders = exchange.GetOrders()
                        for order in orders:
                            exchange.CancelOrder(order.Id)
                        if len(orders) == 0:
                            break
                    nowAccount = exchange.GetAccount()
                marketPosition = 0
            if marketPosition == 0:
                LogProfit(_N(nowAccount.Balance - initAccount.Balance, 4), nowAccount)

小樣本測(cè)試預(yù)測(cè)對(duì)的概率是三分之一，是不是很有趣！（預(yù)測(cè)情況分三種漲、跌、橫盤(pán)）
原文鏈接： https://www.botvs.com/strateg...