# pyFTS.models.sadaei — pyFTS 1.4 documentation

# Source code for pyFTS.models.sadaei

    """
    First Order Exponentialy Weighted Fuzzy Time Series by Sadaei et al. (2013)

    H. J. Sadaei, R. Enayatifar, A. H. Abdullah, and A. Gani, “Short-term load forecasting using a hybrid model with a
    refined exponentially weighted fuzzy time series and an improved harmony search,” Int. J. Electr. Power Energy Syst., vol. 62, no. from 2005, pp. 118–129, 2014.
    """

    import numpy as np
    from pyFTS.common import FuzzySet,FLR,fts, flrg

    default_c = 1.1

    class ExponentialyWeightedFLRG(flrg.FLRG):
        """First Order Exponentialy Weighted Fuzzy Logical Relationship Group"""
        def __init__(self, LHS, **kwargs):
            super(ExponentialyWeightedFLRG, self).__init__(1, **kwargs)
            self.LHS = LHS
            self.RHS = []
            self.count = 0.0
            self.c = kwargs.get("c",default_c)
            self.w = None

        def append_rhs(self, c, **kwargs):
            count = kwargs.get('count', 1.0)
            self.RHS.append(c)
            self.count += count

        def weights(self):
            if self.w is None:
                wei = [self.c ** k for k in np.arange(0.0, self.count, 1.0)]
                tot = sum(wei)
                self.w = np.array([k / tot for k in wei])
            return self.w

        def __str__(self):
            tmp = self.LHS + " -> "
            tmp2 = ""
            cc = 0
            wei = [self.c ** k for k in np.arange(0.0, self.count, 1.0)]
            tot = sum(wei)
            for c in sorted(self.RHS):
                if len(tmp2) > 0:
                    tmp2 = tmp2 + ","
                tmp2 = tmp2 + c + "(" + str(wei[cc] / tot) + ")"
                cc = cc + 1
            return tmp + tmp2

        def __len__(self):
            return len(self.RHS)

    class ExponentialyWeightedFTS(fts.FTS):
        """First Order Exponentialy Weighted Fuzzy Time Series"""
        def __init__(self, **kwargs):
            super(ExponentialyWeightedFTS, self).__init__(order=1, name="EWFTS", **kwargs)
            self.name = "Exponentialy Weighted FTS"
            self.detail = "Sadaei"
            self.c = kwargs.get('c', default_c)

        def generate_flrg(self, flrs, c):
            for flr in flrs:
                if flr.LHS in self.flrgs:
                    self.flrgs[flr.LHS].append_rhs(flr.RHS)
                else:
                    self.flrgs[flr.LHS] = ExponentialyWeightedFLRG(flr.LHS, c=c);
                    self.flrgs[flr.LHS].append_rhs(flr.RHS)

        def train(self, data, **kwargs):
            tmpdata = FuzzySet.fuzzyfy(data, partitioner=self.partitioner, method='maximum', mode='sets')
            flrs = FLR.generate_recurrent_flrs(tmpdata)
            self.generate_flrg(flrs, self.c)

        def forecast(self, ndata, **kwargs):

            explain = kwargs.get('explain', False)

            if self.partitioner is not None:
                ordered_sets = self.partitioner.ordered_sets
            else:
                ordered_sets = FuzzySet.set_ordered(self.sets)

            data = np.array(ndata)

            l = len(ndata)

            ret = []

            for k in np.arange(0, l):

                actual = FuzzySet.get_maximum_membership_fuzzyset(ndata[k], self.sets, ordered_sets)

                if explain:
                    print("Fuzzyfication:\n\n {} -> {} \n".format(ndata[k], actual.name))

                if actual.name not in self.flrgs:
                    ret.append(actual.centroid)

                    if explain:
                        print("Rules:\n\n {} -> {} (Naïve)\t Midpoint: {}  \n\n".format(actual.name, actual.name,actual.centroid))

                else:
                    flrg = self.flrgs[actual.name]
                    mp = flrg.get_midpoints(self.sets)

                    final = mp.dot(flrg.weights())

                    ret.append(final)

                    if explain:
                        print("Rules:\n\n {} \n\n ".format(str(flrg)))
                        print("Midpoints: \n\n {}\n\n".format(mp))

                        print("Deffuzyfied value: {} \n".format(final))

            return ret