• CART决策树算法:使用Gini指数选择最优划分特征。

    • 递归构建二叉树结构(但代码实际生成多叉树,与标准CART不同)。

    • 终止条件为节点样本全属同一类别(返回该类)或者无剩余特征可用(返回多数类)。

程序代码:

import math
import json
import numpy as np
import pandas as pd
from sklearn.model_selection import train_test_split
import random

file_path = '西瓜数据集.xlsx'
label_list = []
data = pd.read_excel(file_path, skiprows=1,nrows=17)
data_dict = data.to_dict(orient='list')
data_list = list(zip(data['色泽'], data['根蒂'], data['敲声'], data['纹理'], data['脐部'], data['触感'], data['好瓜']))
train_set, test_set = train_test_split(data_list, test_size=2, random_state=random.randint(1, 1000))

print("data_list:",len(data_list),data_list)
print("train_list:",len(train_set),train_set)
print("test_list:",len(test_set),test_set)

keys = []
#keys.append(None)
for index in data.keys():
    keys.append(index)
keys.remove('编号')
keys.remove('好瓜')
print('keys:',keys)

def predict(tree, sample):
    if isinstance(tree, dict):
        feature, subtree = next(iter(tree.items()))
        value = sample[feature]
        if value in subtree:
            return predict(subtree[value], sample)
    else:
        return tree

def evaluate(tree, test_set, label_index):
    correct_predictions = 0
    total_samples = len(test_set)

    for sample in test_set:
        prediction = predict(tree, sample)
        if prediction == sample[label_index]:
            correct_predictions += 1

    accuracy = correct_predictions / total_samples
    return accuracy


def calculate_gini(data, label_index):
    labels = [entry[label_index] for entry in data]
    unique_labels = set(labels)
    total_samples = len(data)
    gini = 1.0

    for label in unique_labels:
        prob = labels.count(label) / total_samples
        gini -= prob ** 2

    return gini


def calculate_feature_gini(data, feature_index, label_index):
    feature_values = set([entry[feature_index] for entry in data])
    gini = 0

    for value in feature_values:
        subset_data = [entry for entry in data if entry[feature_index] == value]
        prob = len(subset_data) / len(data)
        gini += prob * calculate_gini(subset_data, label_index)

    return gini


def choose_best_feature_cart(data, features, label_index):
    best_gini = float('inf')
    best_feature = None

    for feature in features:
        feature_gini = calculate_feature_gini(data, features.index(feature), label_index)

        print(f"Feature: {feature}, Gini: {feature_gini}")

        if feature_gini < best_gini:
            best_gini = feature_gini
            best_feature = feature

    return best_feature


def cart_decision_tree(data, features, label_index):
    if len(set(data[i][label_index] for i in range(len(data)))) == 1:
        return data[0][label_index]

    if len(features) == 0:
        max_class = max(set(entry[label_index] for entry in data), key=lambda x: data.count(x))
        return max_class

    root_feature = choose_best_feature_cart(data, features, label_index)
    if root_feature is None:
        max_class = max(set(entry[label_index] for entry in data), key=lambda x: data.count(x))
        return max_class

    print(root_feature)
    tree = {root_feature: {}}
    for value in set(data[i][features.index(root_feature)] for i in range(len(data))):
        subset_data = [data[i] for i in range(len(data)) if data[i][features.index(root_feature)] == value]
        subset_features = [feat for feat in features if feat != root_feature]
        tree[root_feature][value] = cart_decision_tree(subset_data, subset_features, label_index)

    return tree

# 测试
#'色泽','根蒂','敲声','纹理','脐部','触感'
selected_features = [0,1,2,3,4,5]
label_index = 6
cart_decision_tree = cart_decision_tree(train_set, selected_features, label_index)
print(json.dumps(cart_decision_tree, indent=2, ensure_ascii=False))

test_accuracy_cart = evaluate(cart_decision_tree, test_set, label_index)
print("CART Test Accuracy: {:.2%}".format(test_accuracy_cart))

运行结果:

data_list: 17 [('青绿', '蜷缩', '浊响', '清晰', '凹陷', '硬滑', '是'), ('乌黑', '蜷缩', '沉闷', '清晰', '凹陷', '硬滑', '是'), ('乌黑', '蜷缩', '浊响', '清晰', '凹陷', '硬滑', '是'), ('青绿', '蜷缩', '沉闷', '清晰', '凹陷', '硬滑', '是'), ('浅白', '蜷缩', '浊响', '清晰', '凹陷', '硬滑', '是'), ('青绿', '稍蜷', '浊响', '清晰', '稍凹', '软粘', '是'), ('乌黑', '稍蜷', '浊响', '稍糊', '稍凹', '软粘', '是'), ('乌黑', '稍蜷', '浊响', '清晰', '稍凹', '硬滑', '是'), ('乌黑', '稍蜷', '沉闷', '稍糊', '稍凹', '硬滑', '否'), ('青绿', '硬挺', '清脆', '清晰', '平坦', '软粘', '否'), ('浅白', '硬挺', '清脆', '模糊', '平坦', '硬滑', '否'), ('浅白', '蜷缩', '浊响', '模糊', '平坦', '软粘', '否'), ('青绿', '稍蜷', '浊响', '稍糊', '凹陷', '硬滑', '否'), ('浅白', '稍蜷', '沉闷', '稍糊', '凹陷', '硬滑', '否'), ('乌黑', '稍蜷', '浊响', '清晰', '稍凹', '软粘', '否'), ('浅白', '蜷缩', '浊响', '模糊', '平坦', '硬滑', '否'), ('青绿', '蜷缩', '沉闷', '稍糊', '稍凹', '硬滑', '否')]
train_list: 15 [('浅白', '硬挺', '清脆', '模糊', '平坦', '硬滑', '否'), ('青绿', '蜷缩', '沉闷', '清晰', '凹陷', '硬滑', '是'), ('青绿', '硬挺', '清脆', '清晰', '平坦', '软粘', '否'), ('乌黑', '稍蜷', '浊响', '清晰', '稍凹', '硬滑', '是'), ('浅白', '蜷缩', '浊响', '模糊', '平坦', '硬滑', '否'), ('浅白', '稍蜷', '沉闷', '稍糊', '凹陷', '硬滑', '否'), ('青绿', '稍蜷', '浊响', '清晰', '稍凹', '软粘', '是'), ('浅白', '蜷缩', '浊响', '清晰', '凹陷', '硬滑', '是'), ('青绿', '蜷缩', '浊响', '清晰', '凹陷', '硬滑', '是'), ('浅白', '蜷缩', '浊响', '模糊', '平坦', '软粘', '否'), ('青绿', '稍蜷', '浊响', '稍糊', '凹陷', '硬滑', '否'), ('乌黑', '稍蜷', '沉闷', '稍糊', '稍凹', '硬滑', '否'), ('乌黑', '蜷缩', '浊响', '清晰', '凹陷', '硬滑', '是'), ('乌黑', '稍蜷', '浊响', '稍糊', '稍凹', '软粘', '是'), ('乌黑', '稍蜷', '浊响', '清晰', '稍凹', '软粘', '否')]
test_list: 2 [('乌黑', '蜷缩', '沉闷', '清晰', '凹陷', '硬滑', '是'), ('青绿', '蜷缩', '沉闷', '稍糊', '稍凹', '硬滑', '否')]
keys: ['色泽', '根蒂', '敲声', '纹理', '脐部', '触感']
Feature: 0, Gini: 0.4266666666666665
Feature: 1, Gini: 0.4063492063492064
Feature: 2, Gini: 0.40888888888888886
Feature: 3, Gini: 0.30000000000000004
Feature: 4, Gini: 0.3377777777777778
Feature: 5, Gini: 0.4933333333333333
3
Feature: 0, Gini: 0.25
Feature: 1, Gini: 0.375
Feature: 2, Gini: 0.25
Feature: 4, Gini: 0.375
Feature: 5, Gini: 0.25
0
Feature: 1, Gini: 0.5
Feature: 2, Gini: 0.5
Feature: 4, Gini: 0.0
Feature: 5, Gini: 0.5
4
Feature: 0, Gini: 0.3541666666666667
Feature: 1, Gini: 0.16666666666666669
Feature: 2, Gini: 0.20833333333333326
Feature: 4, Gini: 0.375
Feature: 5, Gini: 0.16666666666666669
1
Feature: 0, Gini: 0.3333333333333333
Feature: 2, Gini: 0.4444444444444445
Feature: 4, Gini: 0.4444444444444445
Feature: 5, Gini: 0.4444444444444445
0
Feature: 2, Gini: 0.5
Feature: 4, Gini: 0.5
Feature: 5, Gini: 0.5
2
Feature: 4, Gini: 0.5
Feature: 5, Gini: 0.5
4
Feature: 5, Gini: 0.5
5
{
  "3": {
    "稍糊": {
      "0": {
        "乌黑": {
          "4": {
            "浊响": "是",
            "沉闷": "否"
          }
        },
        "青绿": "否",
        "浅白": "否"
      }
    },
    "模糊": "否",
    "清晰": {
      "1": {
        "硬挺": "否",
        "稍蜷": {
          "0": {
            "乌黑": {
              "2": {
                "乌黑": {
                  "4": {
                    "乌黑": {
                      "5": {
                        "乌黑": "是"
                      }
                    }
                  }
                }
              }
            },
            "青绿": "是"
          }
        },
        "蜷缩": "是"
      }
    }
  }
}
CART Test Accuracy: 100.00%

Logo

北京人形旗下天工造物具身智能开源社区,聚焦具身天工与慧思开物两大平台

更多推荐