0

scikitのLOGO(1つのグループを除外)をクロスバリデーションメソッドとして、学習曲線と組み合わせて使用​​するのが好きです。これは、私が扱うほとんどのケースではうまく機能しますが、私は経験から得た最も重要な2つのパラメータ(最大の特徴と推定数)を使用するだけです(効率的に)。以下の私のコードの例:Scikit-learnでRandomizedSearchCV(またはGridSearcCV)とLeaveOneGroupOutのクロスバリデーションを組み合わせる

Fscorer = make_scorer(f1_score, average = 'micro') 
    gp = training_data["GP"].values 
    logo = LeaveOneGroupOut() 
    from sklearn.ensemble import RandomForestClassifier 
    RF_clf100 = RandomForestClassifier (n_estimators=100, n_jobs=-1, random_state = 49) 
    RF_clf200 = RandomForestClassifier (n_estimators=200, n_jobs=-1, random_state = 49) 
    RF_clf300 = RandomForestClassifier (n_estimators=300, n_jobs=-1, random_state = 49) 
    RF_clf400 = RandomForestClassifier (n_estimators=400, n_jobs=-1, random_state = 49) 
    RF_clf500 = RandomForestClassifier (n_estimators=500, n_jobs=-1, random_state = 49) 
    RF_clf600 = RandomForestClassifier (n_estimators=600, n_jobs=-1, random_state = 49) 

    param_name = "max_features" 
    param_range = param_range = [5, 10, 15, 20, 25, 30] 


    plt.figure() 
    plt.suptitle('n_estimators = 100', fontsize=14, fontweight='bold') 
    _, test_scores = validation_curve(RF_clf100, X, y, cv=logo.split(X, y, groups=gp), 
             param_name=param_name, param_range=param_range, 
             scoring=Fscorer, n_jobs=-1) 
    test_scores_mean = np.mean(test_scores, axis=1) 
    plt.plot(param_range, test_scores_mean) 
    plt.xlabel(param_name) 
    plt.xlim(min(param_range), max(param_range)) 
    plt.ylabel("F1") 
    plt.ylim(0.47, 0.57) 
    plt.legend(loc="best") 
    plt.show() 


    plt.figure() 
    plt.suptitle('n_estimators = 200', fontsize=14, fontweight='bold') 
    _, test_scores = validation_curve(RF_clf200, X, y, cv=logo.split(X, y, groups=gp), 
             param_name=param_name, param_range=param_range, 
             scoring=Fscorer, n_jobs=-1) 
    test_scores_mean = np.mean(test_scores, axis=1) 
    plt.plot(param_range, test_scores_mean) 
    plt.xlabel(param_name) 
    plt.xlim(min(param_range), max(param_range)) 
    plt.ylabel("F1") 
    plt.ylim(0.47, 0.57) 
    plt.legend(loc="best") 
    plt.show() 
    ... 
    ... 

私が本当にかかわらたいことは、より徹底したパラメータ空間検索のため、グリッドサーチ、またはランダム化された検索でロゴを組み合わせることです。

私のコードは次のようになります。今のところ:

param_dist = {"n_estimators": [100, 200, 300, 400, 500, 600], 
       "max_features": sp_randint(5, 30), 
       "max_depth": sp_randint(2, 18), 
       "criterion": ['entropy', 'gini'], 
       "min_samples_leaf": sp_randint(2, 17)} 

clf = RandomForestClassifier(random_state = 49) 

n_iter_search = 45 
random_search = RandomizedSearchCV(clf, param_distributions=param_dist, 
            n_iter=n_iter_search, 
            scoring=Fscorer, cv=8, 
            n_jobs=-1) 
random_search.fit(X, y) 

私はcv=logo.split(X, y, groups=gp)cv = 8を交換するとき、私はこのエラーメッセージを得る:(1)何が起こっているのかについての

--------------------------------------------------------------------------- 
TypeError         Traceback (most recent call last) 
<ipython-input-10-0092e11ffbf4> in <module>() 
---> 35 random_search.fit(X, y) 


/Applications/anaconda/lib/python2.7/site-packages/sklearn/model_selection/_search.pyc in fit(self, X, y, groups) 
    1183           self.n_iter, 
    1184           random_state=self.random_state) 
-> 1185   return self._fit(X, y, groups, sampled_params) 

/Applications/anaconda/lib/python2.7/site-packages/sklearn/model_selection/_search.pyc in _fit(self, X, y, groups, parameter_iterable) 
    540 
    541   X, y, groups = indexable(X, y, groups) 
--> 542   n_splits = cv.get_n_splits(X, y, groups) 
    543   if self.verbose > 0 and isinstance(parameter_iterable, Sized): 
    544    n_candidates = len(parameter_iterable) 

/Applications/anaconda/lib/python2.7/site-packages/sklearn/model_selection/_split.pyc in get_n_splits(self, X, y, groups) 
    1489    Returns the number of splitting iterations in the cross-validator. 
    1490   """ 
-> 1491   return len(self.cv) # Both iterables and old-cv objects support len 
    1492 
    1493  def split(self, X=None, y=None, groups=None): 

TypeError: object of type 'generator' has no len() 

任意の提案をさらに重要なのは、(2)RandomizedSearchCVとLeaveOneGroupOutを組み合わせて動作させる方法ですか?

* UPDATE 2017年2月8日*

それはあなただけにlogoようcvオブジェクトを渡し、RandomizedSearchCVにlogo.split()を渡すべきではありませんrandom_search.fit(X, y, wells)

答えて

1

の@Vivekクマー」の提案とcv=logoを使用して働いていましたそれ。 RandomizedSearchCVは列車のテスト指数を生成するためにsplit()を内部的に呼び出します。 gpグループを、RandomizedSearchCVまたはGridSearchCVオブジェクトのfit()コールに渡すことができます。代わりにこれを行うの

random_search.fit(X, y) 

これを行う:

random_search.fit(X, y, gp) 

EDIT:あなたはまた、辞書のようにパラメータfit_paramsにGridSearchCVまたはRandomizedSearchCVのコンストラクタにGPを渡すことができます。

+0

わかりません。 'cv.get_n_splits'をどこに渡しますか? – MyCarta

+0

@MyCarta申し訳ありません、私は 'cv.get_n_splits'ではなく、' logo.split() 'について話していました。私は混乱を取り除くために私の答えを編集しました。 –

+0

@ Vivek Kumar okこれは少しはっきりしています。回避策がないと言っていますか? – MyCarta

関連する問題