teratail header banner
teratail header banner
質問するログイン新規登録

質問編集履歴

1

やってみたことについての記載がなかったため、追加しました。他にも資料やウェブページなどを見ましたが理解して実装することができなかったので、訂正もしくはサンプルコードを教えていただけると嬉しいです。

2018/09/24 06:49

投稿

su142ie
su142ie

スコア12

title CHANGED
File without changes
body CHANGED
@@ -1,4 +1,54 @@
1
+ ### やりたいこと
1
- 説明変数が複数あるSVRの実装をしたいのですが、
2
+ - 説明変数が複数あるSVRの実装
3
+
4
+ ### 起きている問題
5
+ ウェブページを参考にやってみましたが、
6
+ 計算結果が出てきませんでした。
7
+
2
- まったくどこから取り組めばいいのかわかりません
8
+ 結果が表示されない原因があれ教えてください。
3
- どなたかデータセット(ex.[wine quality](http://archive.ics.uci.edu/ml/datasets/Wine+Quality), boston)を使って
4
- サンプルコードを作っていただけませんでょうか
9
+ また、以下についも教えていただければ嬉いです
10
+ - パラメータの表示
11
+ - 予測値をグラフに表す
12
+
13
+ 使用したデータセットはUCIのBoston house-prices(米国ボストン市郊外における地域別の住宅価格のデータセット)です。
14
+
15
+
16
+ ```python
17
+ import numpy as np
18
+ from sklearn.svm import SVR
19
+ from matplotlib import pyplot as plt
20
+ from sklearn.datasets import load_boston
21
+ boston = load_boston()
22
+ import pandas as pd
23
+ pd.DataFrame(boston.data, columns=boston.feature_names)
24
+
25
+ #簡易的にトレーニングデータとテストデータを分けた
26
+ X = boston.data[:354, :]
27
+ y = boston.target[:354]
28
+
29
+ svr_rbf = SVR(kernel='rbf', C=1e3, gamma=0.1)
30
+ svr_lin = SVR(kernel='linear', C=1e3)
31
+ svr_poly = SVR(kernel='poly', C=1e3, degree=3)
32
+ y_rbf = svr_rbf.fit(X, y).predict(X)
33
+ y_lin = svr_lin.fit(X, y).predict(X)
34
+ y_poly = svr_poly.fit(X, y).predict(X)
35
+
36
+ test_X = boston.data[354:, :]
37
+ test_y = boston.target[354:]
38
+
39
+ test_rbf = svr_rbf.predict(test_X)
40
+ test_lin = svr_lin.predict(test_X)
41
+ test_poly = svr_poly.predict(test_X)
42
+
43
+ from sklearn.metrics import mean_squared_error
44
+ from math import sqrt
45
+
46
+ # RMSEを計算
47
+ rbf_rmse = sqrt(mean_squared_error(test_y, test_rbf))
48
+ lin_rmse = sqrt(mean_squared_error(test_y, test_lin))
49
+ poly_rmse = sqrt(mean_squared_error(test_y, test_poly))
50
+
51
+ ```
52
+
53
+ 参考URL
54
+ [scikit-learnで多変数の回帰モデル - SVRの比較検証やってみた](https://qiita.com/Leonhalt2714/items/dc704f841d21627e988a)