项目需要,用随机森林和决策树对已有50个事件做预测

import time
import pandas as pd
from sklearn import metrics
from sklearn.model_selection import train_test_split


# Random Forest Classifier
def random_forest_classifier(train_X, train_y):
    from sklearn.ensemble import RandomForestClassifier
    model = RandomForestClassifier(n_estimators=10)
    model.fit(train_X, train_y)  #拟合模型 由X_train, y_train训练数据集建模;
    return model


# Decision Tree Classifier
def decision_tree_classifier(train_X, train_y):
    from sklearn import tree
    model = tree.DecisionTreeClassifier()
    model.fit(train_X, train_y)
    return model


def read_data(data_file):
    col = ['ric']
    col += [str(i) for i in range(50)]
    col += ['label']
    data = pd.read_csv(data_file, names=col)
    # data.head()

    y = data['label'].apply(lambda x : 0 if x is -1 else 1)
    X_train, X_test, y_train, y_test = train_test_split(data.drop(['ric', 'label'], axis=1), y, test_size=0.2, random_state=0)
    return X_train, X_test, y_train, y_test


if __name__ == '__main__':
    data_file = 'C:\\Python36\\TestCode\\RandomForest\\part-00000'
    test_classifiers = ['Random Forest', 'Decision Tree']
    classifiers = {'Random Forest': random_forest_classifier,
                    'Decision Tree': decision_tree_classifier}

    print('****Reading training and testing data...****')

    X_train, X_test, y_train, y_test = read_data(data_file)
    num_train, num_feat = X_train.shape
    num_test, num_feat = X_test.shape
    print('******************** Data Info *********************')
    print('#training data: %d, #testing_data: %d, dimension: %d' % (num_train, num_test, num_feat))

    for classifier in test_classifiers:
        print('******************* %s ********************' % classifier)
        start_time = time.time()
        model = classifiers[classifier](X_train, y_train)
        print('training took %fs!' % (time.time() - start_time))
        predict = model.predict(X_test)  #模型预测, X_test测试数据集预测;对训练数据集测试得分(因为有时根本不知道测试数据集对应的真实y值)
        
        #score = model.score(X_test, y_test)
        precision = metrics.precision_score(y_test, predict) #
        accuracy = metrics.accuracy_score(y_test, predict)
        recall = metrics.recall_score(y_test, predict)
        F_Measure = metrics.f1_score(y_test, predict)

        #print('Model score is: %.2f%%' % (100 * score))  #评估模型准确率
        print('precision ratio: %.2f%%' % (100 * precision)) #准确率
        print('Accuracy ratio: %.2f%%' % (100 * accuracy))
        print('Recall ratio: %.2f%%' % (100 * recall))
        print('F-Measure ratio: %.2f%%' % (100 * F_Measure))

 

精确率和准确率大概在54%,召回率只有30%+,还需要继续调

 

posted on 2018-07-10 16:04  Fia  阅读(1737)  评论(0编辑  收藏  举报