參考資料:
機器學習實戰
'''
@version: 0.0.1
@Author: tqrs
@dev: python3 vscode
@Date: 2019-11-12 09:29:53
@LastEditTime: 2019-11-12 12:24:30
@FilePath: \\機器學習實戰\\12-FP-growth算法\\FPGrowth.py
@Descripttion: 只需對數據庫進行兩次掃描,第一次對所有元素項出現次數進行統計;第二次只考慮頻繁元素,它能夠更爲高效的挖掘數據
'''
class treeNode:
def __init__(self, nameValue, numOccur, parentNode):
self.name = nameValue
self.count = numOccur
self.nodeLink = None
self.parent = parentNode
self.children = {}
def inc(self, numOccur):
self.count += numOccur
def disp(self, ind=1):
print(' ' * ind, self.name, ' ', self.count)
for child in self.children.values():
child.disp(ind + 1)
def createTree(dataSet, minSup=1):
"""
[summary]:將數據集轉化爲FP樹
1. 遍歷數據集,統計各元素項出現次數,創建頭指針表
2. 移除頭指針表中不滿足最小值尺度的元素項
3. 第二次遍歷數據集,創建FP樹。對每個數據集中的項集:
3.1 初始化空FP樹
3.2 對每個項集進行過濾和重排序
3.3 使用這個項集更新FP樹,從FP樹的根節點開始:
3.3.1 如果當前項集的第一個元素項存在於FP樹當前節點的子節點中,則更新這個子節點的計數值
3.3.2 否則,創建新的子節點,更新頭指針表
3.3.3 對當前項集的其餘元素項和當前元素項的對應子節點遞歸3.3的過程
Arguments:
dataSet -- 數據集
Keyword Arguments:
minSup {int} -- 最小支持度 (default: {1})
Returns:
[type] -- [description]
"""
headerTable = {}
for trans in dataSet:
for item in trans:
headerTable[item] = headerTable.get(item, 0) + dataSet[trans]
for k in list(headerTable.keys()):
if headerTable[k] < minSup:
del (headerTable[k])
freqItemSet = set(headerTable.keys())
if len(freqItemSet) == 0:
return None, None
for k in headerTable:
headerTable[k] = [headerTable[k], None]
retTree = treeNode('Null Set', 1, None)
for tranSet, count in dataSet.items():
localD = {}
for item in tranSet:
if item in freqItemSet:
localD[item] = headerTable[item][0]
if len(localD) > 0:
orderedItems = [
v[0] for v in sorted(
localD.items(), key=lambda p: p[1], reverse=True)
]
updateTree(orderedItems, retTree, headerTable, count)
return retTree, headerTable
def updateTree(items, inTree, headerTabel, count):
if items[0] in inTree.children:
inTree.children[items[0]].inc(count)
else:
inTree.children[items[0]] = treeNode(items[0], count, inTree)
if headerTabel[items[0]][1] == None:
headerTabel[items[0]][1] = inTree.children[items[0]]
else:
updateHeader(headerTabel[items[0]][1], inTree.children[items[0]])
if len(items) > 1:
updateTree(items[1::], inTree.children[items[0]], headerTabel, count)
def updateHeader(nodeToTest, targetNode):
while (nodeToTest.nodeLink != None):
nodeToTest = nodeToTest.nodeLink
nodeToTest.nodeLink = targetNode
def loadSimpDat():
simpDat = [['r', 'z', 'h', 'j', 'p'],
['z', 'y', 'x', 'w', 'v', 'u', 't', 's'], ['z'],
['r', 'x', 'n', 'o', 's'], ['y', 'r', 'x', 'z', 'q', 't', 'p'],
['y', 'z', 'x', 'e', 'q', 's', 't', 'm']]
return simpDat
def createInitSet(dataSet):
retDict = {}
for trans in dataSet:
retDict[frozenset(trans)] = 1
return retDict
def test_fp():
simpDat = loadSimpDat()
initSet = createInitSet(simpDat)
myFPtree, myHeaderTab = createTree(initSet, 3)
myFPtree.disp()
def ascendTree(leafNode, prefixPath):
if leafNode.parent != None:
prefixPath.append(leafNode.name)
ascendTree(leafNode.parent, prefixPath)
def findPrefixPath(basePat, treeNode):
condPats = {}
while treeNode != None:
prefixPath = []
ascendTree(treeNode, prefixPath)
if len(prefixPath) > 1:
condPats[frozenset(prefixPath[1:])] = treeNode.count
treeNode = treeNode.nodeLink
return condPats
def test_pre():
simpDat = loadSimpDat()
initSet = createInitSet(simpDat)
myFPtree, myHeaderTab = createTree(initSet, 3)
condPats = findPrefixPath('r', myHeaderTab['r'][1])
print(condPats)
def mineTree(inTree, headerTable, minSup, preFix, freqItemList):
"""
[summary]:遞歸查找頻繁項集
Arguments:
inTree {[type]} -- [description]
headerTable {[type]} -- [description]
minSup {[type]} -- [description]
preFix {[type]} -- [description]
freqItemList -- 頻繁項集列表
"""
bigL = [v[0] for v in sorted(headerTable.items(), key=lambda p: p[1][0])]
for basePat in bigL:
newFreqSet = preFix.copy()
newFreqSet.add(basePat)
freqItemList.append(newFreqSet)
condPattBases = findPrefixPath(basePat, headerTable[basePat][1])
myCondTree, myHead = createTree(condPattBases, minSup)
if myHead != None:
print('conditional tree for: ', newFreqSet)
myCondTree.disp(1)
mineTree(myCondTree, myHead, minSup, newFreqSet, freqItemList)
def test_mineTree():
simpDat = loadSimpDat()
initSet = createInitSet(simpDat)
myFPtree, myHeaderTab = createTree(initSet, 3)
freqItems = []
mineTree(myFPtree, myHeaderTab, 3, set([]), freqItems)
if __name__ == '__main__':
paresdDat = [
line.split()
for line in open(r'.\12-FP-growth算法\kosarak.dat').readlines()
]
initSet = createInitSet(paresdDat)
myFPtree, myHeaderTab = createTree(initSet, 100000)
myFreaList = []
mineTree(myFPtree, myHeaderTab, 100000, set([]), myFreaList)
print('len:', len(myFreaList), 'myFreaList', myFreaList)