numpyの練習を兼ねてソフトマージンSVMを書きました. 行列演算に関してはRと比べても大きな不自由は感じないですね. 最適化部分は二次計画法のソルバー(openopt, cvxopt)に丸投げ. openoptから呼べるソルバーのうちデフォルトで含まれているqlcpではデータによっては解けなかったのでcvxoptを利用しました. ソフトマージンSVM - 人工知能に関する断創録の記事が非常にわかりやすくておすすめです. ソースコード 適当にガウス分布を2つ作って学習してます. 必要なライブラリ,numpy,openopt,cvxopt. import numpy as np import openopt as oo #settings C = 1 EPS = 0.00001 #data generation x1 = np.random.multivariate_normal([0,