Source code for PAMI.periodicFrequentPattern.basic.PFPGrowth

# PFPGrowth is one of the fundamental algorithm to discover periodic-frequent patterns in a transactional database.
#
# **Importing this algorithm into a python program**
#
#             from PAMI.periodicFrequentPattern.basic import PFPGrowth as alg
#
#             iFile = 'sampleDB.txt'
#
#             minSup = 10  # can also be specified between 0 and 1
#
#             maxPer = 20 # can also be specified between 0 and 1
#
#             obj = alg.PFPGrowth(iFile, minSup, maxPer)
#
#             obj.mine()
#
#             periodicFrequentPatterns = obj.getPatterns()
#
#             print("Total number of Periodic Frequent Patterns:", len(periodicFrequentPatterns))
#
#             obj.save("periodicFrequentPatterns")
#
#             Df = obj.getPatternsAsDataFrame()
#
#             memUSS = obj.getMemoryUSS()
#
#             print("Total Memory in USS:", memUSS)
#
#             memRSS = obj.getMemoryRSS()
#
#             print("Total Memory in RSS", memRSS)
#
#             run = obj.getRuntime()
#
#             print("Total ExecutionTime in seconds:", run)
#


__copyright__ = """
 Copyright (C)  2021 Rage Uday Kiran

     This program is free software: you can redistribute it and/or modify
     it under the terms of the GNU General Public License as published by
     the Free Software Foundation, either version 3 of the License, or
     (at your option) any later version.

     This program is distributed in the hope that it will be useful,
     but WITHOUT ANY WARRANTY; without even the implied warranty of
     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
     GNU General Public License for more details.

     You should have received a copy of the GNU General Public License
     along with this program.  If not, see <https://www.gnu.org/licenses/>.
     Copyright (C)  2021 Rage Uday Kiran

"""

from PAMI.periodicFrequentPattern.basic import abstract as _ab
from typing import Dict, Tuple
import pandas as pd
from deprecated import deprecated
import numpy as np

_maxPer = float()
_minSup = float()
_lno = int()


class _Node(object):
    """
    A class used to represent the node of frequentPatternTree

    :**Attributes**:    - **item** (*int or None*) -- *Storing item of a node.*
                        - **timeStamps** (*list*) -- *To maintain the timestamps of a database at the end of the branch.*
                        - **parent** (*list*) -- *To maintain the parent of every node.*
                        - **children** (*list*) -- *To maintain the children of a node.*

    :**Methods**:    -**addChild(itemName)** -- *Storing the children to their respective parent nodes.*
    """

    def __init__(self, item, locations, parent=None):
        self.item = item
        self.locations = locations
        self.parent = parent
        self.children = {}

    def addChild(self, item, locations):
        """
        This method takes an item and locations as input, adds a new child node
        if the item does not already exist among the current node's children, or
        updates the locations of the existing child node if the item is already present.

        :param item: Represents the distinct item to be added as a child node.
        :type item: Any
        :param locations: Represents the locations associated with the item.
        :type locations: list
        :return: The child node associated with the item.
        :rtype: _Node
        """
        if item not in self.children:
            self.children[item] = _Node(item, locations, self)
        else:
            self.children[item].locations = locations + self.children[item].locations
            
        return self.children[item]

    def traverse(self):
        """
        This method constructs a transaction by traversing from the current node to the root node, collecting items along the way.

        :return: A tuple containing the transaction and the locations associated with the current node.
        :rtype: tuple(list, Any)
        """
        transaction = []
        locs = self.locations
        node = self.parent
        while node.parent is not None:
            transaction.append(node.item)
            node = node.parent
        return transaction[::-1], locs


[docs] class PFPGrowth(_ab._periodicFrequentPatterns): """ **About this algorithm** :**Description**: PFPGrowth is one of the fundamental algorithm to discover periodic-frequent patterns in a transactional database. :**Reference**: Syed Khairuzzaman Tanbeer, Chowdhury Farhan, Byeong-Soo Jeong, and Young-Koo Lee, "Discovering Periodic-Frequent Patterns in Transactional Databases", PAKDD 2009, https://doi.org/10.1007/978-3-642-01307-2_24 :**Parameters**: - **iFile** (*str or URL or dataFrame*) -- *Name of the Input file to mine complete set of frequent patterns.* - **oFile** (*str*) -- *Name of the output file to store complete set of frequent patterns.* - **minSup** (*int or float or str*) -- *The user can specify minSup either in count or proportion of database size. If the program detects the data type of minSup is integer, then it treats minSup is expressed in count. Otherwise, it will be treated as float.* - **maxPer** (*int or float or str*) -- *The user can specify maxPer either in count or proportion of database size. It controls the maximum number of transactions in which any two items within a pattern can reappear.* - **sep** (*str*) -- *This variable is used to distinguish items from one another in a transaction. The default seperator is tab space. However, the users can override their default separator.* :**Attributes**: - **startTime** (*float*) -- *To record the start time of the mining process.* - **endTime** (*float*) -- *To record the completion time of the mining process.* - **finalPatterns** (*dict*) -- *Storing the complete set of patterns in a dictionary variable.* - **memoryUSS** (*float*) -- *To store the total amount of USS memory consumed by the program.* - **memoryRSS** (*float*) -- *To store the total amount of RSS memory consumed by the program.* - **Database** (*list*) -- *To store the transactions of a database in list.* - **mapSupport** (*Dictionary*) -- *To maintain the information of item and their frequency.* - **lno** (*int*) -- *It represents the total no of transactions* - **tree** (*class*) -- *it represents the Tree class.* - **itemSetCount** (*int*) -- *it represents the total no of patterns.* :**Methods**: - **mine()** -- *Mining process will start from here.* - **getPatterns()** -- *Complete set of patterns will be retrieved with this function.* - **save(oFile)** -- *Complete set of periodic-frequent patterns will be loaded in to a output file.* - **getPatternsAsDataFrame()** -- *Complete set of periodic-frequent patterns will be loaded in to a dataframe.* - **getMemoryUSS()** -- *Total amount of USS memory consumed by the mining process will be retrieved from this function.* - **getMemoryRSS()** -- *Total amount of RSS memory consumed by the mining process will be retrieved from this function.* - **getRuntime()** -- *Total amount of runtime taken by the mining process will be retrieved from this function.* - **creatingItemSets(fileName)** -- *Scans the dataset and stores in a list format.* - **PeriodicFrequentOneItem()** -- *Extracts the one-periodic-frequent patterns from database.* - **updateDatabases()** -- *Update the database by removing aperiodic items and sort the Database by item decreased support.* - **buildTree()** -- *After updating the Database, remaining items will be added into the tree by setting root node as null.* - **convert()** -- *This methos is used to convert the user specified value.* **Execution methods** **Terminal command** .. code-block:: console Format: (.venv) $ python3 PFPGrowth.py <inputFile> <outputFile> <minSup> <maxPer> Example usage: (.venv) $ python3 PFPGrowth.py sampleDB.txt patterns.txt 10.0 20.0 .. note:: minSup will be considered in percentage of database transactions **Calling from a python program** .. code-block:: python from PAMI.periodicFrequentPattern.basic import PFPGrowth as alg iFile = 'sampleDB.txt' minSup = 10 # can also be specified between 0 and 1 maxPer = 20 # can also be specified between 0 and 1 obj = alg.PFPGrowth(iFile, minSup, maxPer) obj.mine() periodicFrequentPatterns = obj.getPatterns() print("Total number of Periodic Frequent Patterns:", len(periodicFrequentPatterns)) obj.save("periodicFrequentPatterns") Df = obj.getPatternsAsDataFrame() memUSS = obj.getMemoryUSS() print("Total Memory in USS:", memUSS) memRSS = obj.getMemoryRSS() print("Total Memory in RSS", memRSS) run = obj.getRuntime() print("Total ExecutionTime in seconds:", run) **Credits** The complete program was written by P. Likhitha and revised by Tarun Sreepada under the supervision of Professor Rage Uday Kiran. """ _startTime = float() _endTime = float() _minSup = str() _maxPer = float() _finalPatterns = {} _iFile = " " _oFile = " " _sep = " " _memoryUSS = float() _memoryRSS = float() _Database = [] _rank = {} _rankedUp = {} _lno = 0 def _creatingItemSets(self) -> None: """ Storing the complete transactions of the database/input file in a database variable :return: None """ self._Database = [] if isinstance(self._iFile, _ab._pd.DataFrame): data, ts = [], [] if self._iFile.empty: print("its empty..") i = self._iFile.columns.values.tolist() if 'TS' in i: ts = self._iFile['TS'].tolist() if 'Transactions' in i: data = self._iFile['Transactions'].tolist() for i in range(len(data)): if data[i]: tr = [str(ts[i])] + [x for x in data[i].split(self._sep)] self._Database.append(tr) else: self._Database.append([str(ts[i])]) if isinstance(self._iFile, str): if _ab._validators.url(self._iFile): data = _ab._urlopen(self._iFile) for line in data: line.strip() line = line.decode("utf-8") temp = [i.rstrip() for i in line.split(self._sep)] temp = [x for x in temp if x] self._Database.append(temp) else: try: with open(self._iFile, 'r', encoding='utf-8') as f: for line in f: line.strip() temp = [i.rstrip() for i in line.split(self._sep)] temp = [x for x in temp if x] self._Database.append(temp) except IOError: print("File Not Found") quit() def _convert(self, value) -> int: """ To convert the given user specified value :param value: user specified value :return: converted value """ if type(value) is int: value = int(value) if type(value) is float: value = (len(self._Database) * value) if type(value) is str: if '.' in value: value = float(value) value = (len(self._Database) * value) else: value = int(value) return value
[docs] @deprecated("It is recommended to use 'mine()' instead of 'mine()' for mining process. Starting from January 2025, 'mine()' will be completely terminated.") def startMine(self) -> None: self.mine()
def _getMaxPer(self, arr, maxTS): """ This method appends `0` and `maxTS` to the input array, sorts the array, and then computes the differences between consecutive elements. :param arr: The input array of elements. :type arr: numpy.ndarray :param maxTS: The maximum timestamp to be appended to the array. :type maxTS: int or float :return: None """ arr = np.append(arr, [0, maxTS]) arr = np.sort(arr) arr = np.diff(arr) return np.max(arr) def _construct(self, items, data, minSup, maxPer, maxTS, patterns): """ This method filters the items based on the minimum support (minSup) and maximum period (maxPer). It then constructs a tree structure from the filtered items and data. :param items: A dictionary where keys are items and values are lists of timestamps. :type items: dict :param data: The dataset used to construct the tree, where each entry is a list with an index followed by items. :type data: list of lists :param minSup: The minimum support threshold. :type minSup: int :param maxPer: The maximum period threshold. :type maxPer: int or float :param maxTS: The maximum timestamp. :type maxTS: int or float :param patterns: A dictionary to store the patterns discovered during the construction. :type patterns: dict :return: A tuple containing the root node of the constructed tree and a dictionary of item nodes. :rtype: tuple(_Node, dict) """ # maxPerItems = {k: self.getMaxPer(v, maxTS) for k, v in items.items() if len(v) >= minSup} items = {k: v for k, v in items.items() if len(v) >= minSup and self._getMaxPer(v, maxTS) <= maxPer} #tested ok for item, ts in items.items(): # pat = "\t".join(item) # self.patCount += 1 # patterns[pat] = (len(ts), self.getMaxPer(ts, maxTS)) patterns[tuple([item])] = [len(ts), self._getMaxPer(ts, maxTS)] root = _Node([], None, None) itemNodes = {} for line in data: currNode = root index = int(line[0]) line = line[1:] line = sorted([item for item in line if item in items], key = lambda x: len(items[x]), reverse = True) for item in line: currNode = currNode.addChild(item, [index]) # heavy if item in itemNodes: itemNodes[item].add(currNode) else: itemNodes[item] = set([currNode]) return root, itemNodes def _recursive(self, root, itemNode, minSup, maxPer, patterns, maxTS): """ This method recursively constructs a pattern tree from the given root node, filtering items based on the minimum support (minSup) and maximum period (maxPer). It updates the patterns dictionary with the discovered patterns. :param root: The current root node of the pattern tree. :type root: _Node :param itemNode: A dictionary where keys are items and values are sets of nodes associated with those items. :type itemNode: dict :param minSup: The minimum support threshold. :type minSup: int :param maxPer: The maximum period threshold. :type maxPer: int or float :param patterns: A dictionary to store the patterns discovered during the recursion. :type patterns: dict :param maxTS: The maximum timestamp. :type maxTS: int or float """ for item in itemNode: newRoot = _Node(root.item + [item], None, None) itemLocs = {} transactions = {} for node in itemNode[item]: transaction, locs = node.traverse() if len(transaction) < 1: continue # transactions.append((transaction, locs)) if tuple(transaction) in transactions: transactions[tuple(transaction)].extend(locs) else: transactions[tuple(transaction)] = locs for item in transaction: if item in itemLocs: itemLocs[item] += locs else: itemLocs[item] = list(locs) # Precompute getMaxPer results for itemLocs maxPerResults = {item: self._getMaxPer(itemLocs[item], maxTS) for item in itemLocs if len(itemLocs[item]) >= minSup} # Filter itemLocs based on minSup and maxPer itemLocs = {k: len(v) for k, v in itemLocs.items() if k in maxPerResults and maxPerResults[k] <= maxPer} # Iterate over filtered itemLocs for item in itemLocs: # pat = "\t".join([str(x) for x in newRoot.item + [item]]) # self.patCount += 1 # patterns[pat] = [itemLocs[item], maxPerResults[item]] patterns[tuple(newRoot.item + [item])] = [itemLocs[item], maxPerResults[item]] if not itemLocs: continue newItemNodes = {} for transaction, locs in transactions.items(): transaction = sorted([item for item in transaction if item in itemLocs], key = lambda x: itemLocs[x], reverse = True) if len(transaction) < 1: continue currNode = newRoot for item in transaction: currNode = currNode.addChild(item, locs) if item in newItemNodes: newItemNodes[item].add(currNode) else: newItemNodes[item] = set([currNode]) self._recursive(newRoot, newItemNodes, minSup, maxPer, patterns, _lno)
[docs] def mine(self) -> None: """ Mining process will start from this function :return: None """ global _minSup, _maxPer, _lno self._startTime = _ab._time.time() if self._iFile is None: raise Exception("Please enter the file path or file name:") if self._minSup is None: raise Exception("Please enter the Minimum Support") if self._maxPer is None: raise Exception("Please enter the Maximum Periodicity") if self._sep is None: raise Exception("Default separator is tab space, please enter the separator if you have different separator in the input file") self._creatingItemSets() self._minSup = self._convert(self._minSup) self._maxPer = self._convert(self._maxPer) #tested ok _minSup, _maxPer, _lno = self._minSup, self._maxPer, len(self._Database) if self._minSup > len(self._Database): raise Exception("Please enter the minSup in range between 0 to 1") items = {} # tested ok for line in self._Database: index = int(line[0]) for item in line[1:]: if item not in items: items[item] = [] items[item].append(index) root, itemNodes = self._construct(items, self._Database, _minSup, _maxPer, _lno, self._finalPatterns) self._recursive(root, itemNodes, _minSup, _maxPer, self._finalPatterns, _lno) newPattern = {} for k, v in self._finalPatterns.items(): newPattern["\t".join([str(x) for x in k])] = v self._finalPatterns = newPattern self._endTime = _ab._time.time() process = _ab._psutil.Process(_ab._os.getpid()) self._memoryUSS = float() self._memoryRSS = float() self._memoryUSS = process.memory_full_info().uss self._memoryRSS = process.memory_info().rss print("Periodic Frequent patterns were generated successfully using PFPGrowth algorithm ")
[docs] def getMemoryUSS(self) -> float: """ Total amount of USS memory consumed by the mining process will be retrieved from this function :return: returning USS memory consumed by the mining process :rtype: float """ return self._memoryUSS
[docs] def getMemoryRSS(self) -> float: """ Total amount of RSS memory consumed by the mining process will be retrieved from this function :return: returning RSS memory consumed by the mining process :rtype: float """ return self._memoryRSS
[docs] def getRuntime(self) -> float: """ Calculating the total amount of runtime taken by the mining process :return: returning total amount of runtime taken by the mining process :rtype: float """ return self._endTime - self._startTime
[docs] def getPatternsAsDataFrame(self) -> _ab._pd.DataFrame: """ Storing final periodic-frequent patterns in a dataframe :return: returning periodic-frequent patterns in a dataframe :rtype: pd.DataFrame """ dataFrame = {} data = [] for a, b in self._finalPatterns.items(): data.append([a, b[0], b[1]]) dataFrame = _ab._pd.DataFrame(data, columns=['Patterns', 'Support', 'Periodicity']) return dataFrame
[docs] def save(self, outFile: str) -> None: """ Complete set of periodic-frequent patterns will be loaded in to an output file :param outFile: name of the output file :type outFile: csv file :return: None """ self._oFile = outFile writer = open(self._oFile, 'w+') for x, y in self._finalPatterns.items(): s1 = x + ":" + str(y[0]) + ":" + str(y[1]) #s1 = x.replace(' ', '\t').strip() + ":" + str(y[0]) + ":" + str(y[1]) writer.write("%s \n" % s1)
[docs] def getPatterns(self) -> Dict[str, Tuple[int, int]]: """ Function to send the set of periodic-frequent patterns after completion of the mining process :return: returning periodic-frequent patterns :rtype: dict """ return self._finalPatterns
[docs] def printResults(self) -> None: """ This function is used to print the results :return: None """ print("Total number of Periodic Frequent Patterns:", len(self.getPatterns())) print("Total Memory in USS:", self.getMemoryUSS()) print("Total Memory in RSS", self.getMemoryRSS()) print("Total ExecutionTime in ms:", self.getRuntime())
if __name__ == "__main__": _ap = str() if len(_ab._sys.argv) == 5 or len(_ab._sys.argv) == 6: if len(_ab._sys.argv) == 6: _ap = PFPGrowth(_ab._sys.argv[1], _ab._sys.argv[3], _ab._sys.argv[4], _ab._sys.argv[5]) if len(_ab._sys.argv) == 5: _ap = PFPGrowth(_ab._sys.argv[1], _ab._sys.argv[3], _ab._sys.argv[4]) _ap.mine() _ap.mine() print("Total number of Periodic-Frequent Patterns:", len(_ap.getPatterns())) _ap.save(_ab._sys.argv[2]) print("Total Memory in USS:", _ap.getMemoryUSS()) print("Total Memory in RSS", _ap.getMemoryRSS()) print("Total ExecutionTime in ms:", _ap.getRuntime()) else: print("Error! The number of input parameters do not match the total number of parameters provided")