我正在尝试使用UCI机器学习团队发布的数据集来实现朴素贝叶斯分类器。我是机器学习的新手,并试图理解用于我的工作相关问题的技术,所以我认为最好先理解理论。
我正在使用pima数据集(Link to Data - UCI-ML),我的目标是为K类问题构建朴素贝叶斯单变量高斯分类器(数据仅用于K = 2)。我已经完成了分割数据,并计算每个类的平均值,标准偏差,每个类的先验,但在此之后我有点卡住,因为我不确定在此之后我应该做什么以及如何做。我觉得我应该计算后验概率,
这是我的代码,我使用百分比作为向量,因为我希望看到行为,因为我将训练数据大小从80:20分割增加。基本上如果你通过[10 20 30 40],它将从80:20分割中获得该百分比,并使用80%的10%作为训练。
function[classMean] = naivebayes(file, iter, percent)
dm = load(file);
for i=1:iter
idx = randperm(size(dm.data,1))
%Using same idx for data and labels
shuffledMatrix_data = dm.data(idx,:);
shuffledMatrix_label = dm.labels(idx,:);
percent_data_80 = round((0.8) * length(shuffledMatrix_data));
%Doing 80-20 split
train = shuffledMatrix_data(1:percent_data_80,:);
test = shuffledMatrix_data(percent_data_80+1:length(shuffledMatrix_data),:);
train_labels = shuffledMatrix_label(1:percent_data_80,:)
test_labels = shuffledMatrix_data(percent_data_80+1:length(shuffledMatrix_data),:);
%Getting the array of percents
for pRows = 1:length(percent)
percentOfRows = round((percent(pRows)/100) * length(train));
new_train = train(1:percentOfRows,:)
new_trin_label = shuffledMatrix_label(1:percentOfRows)
%get unique labels in training
numClasses = size(unique(new_trin_label),1)
classMean = zeros(numClasses,size(new_train,2));
for kclass=1:numClasses
classMean(kclass,:) = mean(new_train(new_trin_label == kclass,:))
std(new_train(new_trin_label == kclass,:))
priorClassforK = length(new_train(new_trin_label == kclass))/length(new_train)
priorClassforK_1 = 1 - priorClassforK
end
end
end
end
答案 0 :(得分:2)
首先,根据频率计算计算evey类标签的概率。对于给定的数据样本和数据集中的给定类,可以计算evey特征的概率。之后,将样本中所有特征的条件概率相乘,并将所考虑的类标签的概率相乘。最后,比较所有类标签的值,并选择具有最大概率的类的标签(贝叶斯分类规则)。
对于计算条件概率,您可以简单地使用正态分布函数。