Title: | Fuzzy Rule-Based Systems for Classification and Regression Tasks |
---|---|
Description: | An implementation of various learning algorithms based on fuzzy rule-based systems (FRBSs) for dealing with classification and regression tasks. Moreover, it allows to construct an FRBS model defined by human experts. FRBSs are based on the concept of fuzzy sets, proposed by Zadeh in 1965, which aims at representing the reasoning of human experts in a set of IF-THEN rules, to handle real-life problems in, e.g., control, prediction and inference, data mining, bioinformatics data processing, and robotics. FRBSs are also known as fuzzy inference systems and fuzzy models. During the modeling of an FRBS, there are two important steps that need to be conducted: structure identification and parameter estimation. Nowadays, there exists a wide variety of algorithms to generate fuzzy IF-THEN rules automatically from numerical data, covering both steps. Approaches that have been used in the past are, e.g., heuristic procedures, neuro-fuzzy techniques, clustering methods, genetic algorithms, squares methods, etc. Furthermore, in this version we provide a universal framework named 'frbsPMML', which is adopted from the Predictive Model Markup Language (PMML), for representing FRBS models. PMML is an XML-based language to provide a standard for describing models produced by data mining and machine learning algorithms. Therefore, we are allowed to export and import an FRBS model to/from 'frbsPMML'. Finally, this package aims to implement the most widely used standard procedures, thus offering a standard package for FRBS modeling to the R community. |
Authors: | Lala Septem Riza, Christoph Bergmeir, Francisco Herrera, and Jose Manuel Benitez |
Maintainer: | Christoph Bergmeir <[email protected]> |
License: | GPL (>= 2) | file LICENSE |
Version: | 3.2-0 |
Built: | 2024-10-31 22:13:49 UTC |
Source: | CRAN |
Fuzzy rule-based systems (FRBSs) are based on the fuzzy concept proposed by Zadeh in 1965, which represents the reasoning of human experts in production rules (a set of IF-THEN rules) to handle real-life problems from domains such as control, prediction and inference, data mining, bioinformatics data processing, robotics, and speech recognition. FRBSs are also known as fuzzy inference systems and fuzzy models. When applied to specific tasks, they may also be known under specific names such as fuzzy associative memories or fuzzy controllers. In this package, we consider systems with multi-inputs and single-output (MISO), with real-valued data.
FRBSs are a competitive alternative to other classic models and algorithms in order to solve classification and regression problems. Generally, an FRBS consists of four functional parts:
a fuzzification interface which transforms the crisp inputs into degrees
of membership functions of the linguistic term of each variable.
See fuzzifier
.
a knowledge base consisting of a database (DB) and a rulebase (RB). While the database includes the fuzzy set definitions, the rulebase contains fuzzy IF-THEN rules. We will represent the knowledge as a set of rules. Each one has the following structure.
IF premise (antecedent) THEN conclusion (consequent)
See rulebase
.
an inference engine which performs the inference operations on the fuzzy IF-THEN rules.
There are two kinds of inference for fuzzy systems based on linguistic rules:
The Mamdani and the Takagi Sugeno Kang model. See inference
.
a defuzzification process to obtain the crisp values from linguistic values. There are several methods for
defuzzification such as the weighted average, centroid, etc.
See defuzzifier
.
Since it may be difficult to obtain information from human experts in the form required, an alternative and effective way to acquire the knowledge is to generate the fuzzy IF-THEN rules automatically from the numerical training data. In general, when modeling an FRBS, there are two important processes which should be conducted, namely structure identification and parameter estimation. Structure identification is a process to find appropriate fuzzy IF-THEN rules and to determine the overall number of rules. Parameter estimation is applied to tune parameters of membership functions. Many approaches have been proposed in order to perform this modeling such as a table-lookup scheme, heuristic procedures, neuro-fuzzy techniques, clustering methods, genetic algorithms, least squares methods, gradient descent, etc. In this package, the following approaches to generate fuzzy IF-THEN rules have been implemented:
FRBS based on space partition
FRBS based on neural networks
FRBS based on clustering approach
FRBS based on genetic algorithms
The Thrift's method (GFS.THRIFT
): It is used to solve regression tasks. See GFS.Thrift
.
The Genetic fuzzy systems for fuzzy rule learning based on the MOGUL methodology (GFS.FR.MOGUL
):
It is used to solve regression tasks. See GFS.FR.MOGUL
.
The Ishibuchi's method based on genetic cooperative-competitive learning (GFS.GCCL
):
It is used to solve classification tasks. See GFS.GCCL
.
The Ishibuchi's method based on hybridization of genetic cooperative-competitive learning (GCCL) and Pittsburgh (FH.GBML
):
It is used to solve classification tasks. See FH.GBML
.
The structural learning algorithm on vague environtment (SLAVE
):
It is used to solve classification tasks. See SLAVE
.
The genetic for lateral tuning and rule selection of linguistic fuzzy system (GFS.LT.RS
):
It is used to solve regression tasks. See GFS.LT.RS
.
FRBS based on the gradient descent method
The functions documented in the manual for the single methods are all called internally
by frbs.learn
, which is the central function of the package.
However, in the documentation of each of the internal learning functions,
we give some theoretical background and references to the original literature.
Usage of the package:
First of all, if you have problems using the package, find a bug, or have suggestions, please contact the package maintainer by email, instead of writing to the general R lists or to other internet forums and mailing lists.
The main functions of the package are the following:
The function frbs.learn
allows to generate the model by
creating fuzzy IF-THEN rules or cluster centers from training data.
In other words, users just need to call this function to generate an FRBS model from
training data. The different algorithms mentioned above are all accessible through this function.
The outcome of the function is an frbs-object
.
Even though the main purpose of this package is to generate
the FRBS models from training data automatically, we provide the function frbs.gen
,
which can be used to build a model manually without using a learning method.
Moreover, we provide the following features: linguistic hedges, the "and" and "or" operators,
and the "dont_care" value for representing the degree of 1. The higher degree of interpretability
can also be achieved by using the "dont_care" value. If we want to define various length of rules
in the rulebase, we can also define the "dont_care" value. See rulebase
.
The purpose of the function predict
is to obtain predicted values
according to the testing data and the model (analogous to the predict
function
that is implemented in many other R packages).
There exist functions summary.frbs
and plotMF
to
show a summary about an frbs-object
, and to plot the shapes of
the membership functions.
Exporting an FRBS model to the frbsPMML format can be done by executing frbsPMML
and write.frbsPMML
.
The frbsPMML format is a universal framework adopted from the Predictive Model Markup Language (PMML) format. Then,
in order to consume/import the frbsPMML format to an FRBS model, we call read.frbsPMML
.
To get started with the package, the user can have a look at some examples included in
the documentation of the functions frbs.learn
and frbs.gen
for generating models and
predict
for the prediction phase.
Also, there are many demos that ship with the package. To get a list of them, type:
demo()
Then, to start a demo, type demo(<demo_name_here>)
. All the demos are present as
R scripts in the package sources in the "demo"
subdirectory. Note that
some of them may take quite a long time which depends on specification hardwares.
Currently, there are the following demos available:
Regression using the Gas Furnance dataset:
demo(WM.GasFur)
,
demo(SBC.GasFur)
,
demo(ANFIS.GasFur)
,
demo(FS.HGD.GasFur)
,
demo(DENFIS.GasFur)
,
demo(HyFIS.GasFur)
,
demo(FIR.DM.GasFur)
,
demo(GFS.FR.MOGUL.GasFur)
,
demo(GFS.THRIFT.GasFur)
,
demo(GFS.LT.RS.GasFur)
.
Regression using the Mackey-Glass dataset:
demo(WM.MG1000)
,
demo(SBC.MG1000)
,
demo(ANFIS.MG1000)
,
demo(FS.HGD.MG1000)
,
demo(DENFIS.MG1000)
,
demo(HyFIS.MG1000)
,
demo(GFS.THRIFT.MG1000)
,
demo(FIR.DM.MG1000)
,
demo(GFS.FR.MOGUL.MG1000)
,
demo(GFS.LT.RS.MG1000)
.
Classification using the Iris dataset:
demo(FRBCS.W.Iris)
,
demo(FRBCS.CHI.Iris)
,
demo(GFS.GCCL.Iris)
,
demo(FH.GBML.Iris)
,
demo(SLAVE.Iris)
.
Generating FRBS model without learning process:
demo(FRBS.Mamdani.Manual)
,
demo(FRBS.TSK.Manual)
,
demo(FRBS.Manual)
.
Exporting/importing to/from frbsPMML:
demo(WM.GasFur.PMML)
,
demo(ANFIS.GasFur.PMML)
,
demo(GFS.GCCL.Iris.PMML)
.
The Gas Furnance data and Mackey-Glass data are included in the package,
please see frbsData
. The Iris data is the standard Iris dataset that
ships with R.
Also have a look at the package webpage http://sci2s.ugr.es/dicits/software/FRBS, where we provide a more extensive introduction as well as additional explanations of the procedures.
Lala Septem Riza [email protected],
Christoph Bergmeir [email protected],
Francisco Herrera [email protected],
and Jose Manuel Benitez [email protected]
DiCITS Lab, SCI2S group, DECSAI, University of Granada.
http://sci2s.ugr.es/dicits/, http://sci2s.ugr.es
A. Guazzelli, M. Zeller, W.C. Lin, and G. Williams., "pmml: An open standard for sharing models", The R Journal, Vol. 1, No. 1, pp. 60-65 (2009).
C.C. Lee, "Fuzzy Logic in control systems: Fuzzy logic controller part I", IEEE Trans. Syst., Man, Cybern., vol. 20, no. 2, pp. 404 - 418 (1990).
C.C. Lee, "Fuzzy Logic in control systems: Fuzzy logic controller part II", IEEE Trans. Syst., Man, Cybern., vol. 20, no. 2, pp. 419 - 435 (1990).
E.H. Mamdani and S. Assilian, "An experiment in linguistic synthesis with a fuzzy logic controller," International Journal of Man Machine Studies, vol. 7, no. 1, pp. 1 - 13 (1975).
W. Pedrycz, "Fuzzy Control and Fuzzy Systems," New York: Wiley (1989).
L.S. Riza, C. Bergmeir, F. Herrera, and J.M. Benitez, "frbs: Fuzzy Rule-Based Systems for Classification and Regression in R," Journal of Statistical Software, vol. 65, no. 6, pp. 1 - 30 (2015).
M. Sugeno and G.T. Kang, "Structure identification of fuzzy model," Fuzzy Sets Syst., vol. 28, pp. 15 - 33 (1988).
T. Takagi and M. Sugeno, "Fuzzy identification of systems and its application to modelling and control", IEEE Transactions on Systems, Man and Cybernetics, vol. 15, no. 1, pp. 116 - 132 (1985).
L.A. Zadeh, "Fuzzy sets", Information and Control, vol. 8, pp. 338 - 353 (1965).
frbs.learn
, frbs.gen
, frbsPMML
, and predict
.
################################## ## I. Regression Problem ## In this example, we are using the gas furnace dataset that ## contains two input and one output variables. ################################## ## Input data: Using the Gas Furnace dataset ## then split the data to be training and testing datasets data(frbsData) data.train <- frbsData$GasFurnance.dt[1 : 204, ] data.tst <- frbsData$GasFurnance.dt[205 : 292, 1 : 2] real.val <- matrix(frbsData$GasFurnance.dt[205 : 292, 3], ncol = 1) ## Define interval of data range.data <-apply(data.train, 2, range) ## Set the method and its parameters, ## for example, we use Wang and Mendel's algorithm method.type <- "WM" control <- list(num.labels = 15, type.mf = "GAUSSIAN", type.defuz = "WAM", type.tnorm = "MIN", type.snorm = "MAX", type.implication.func = "ZADEH", name = "sim-0") ## Learning step: Generate an FRBS model object.reg <- frbs.learn(data.train, range.data, method.type, control) ## Predicting step: Predict for newdata res.test <- predict(object.reg, data.tst) ## Display the FRBS model summary(object.reg) ## Plot the membership functions plotMF(object.reg) ################################## ## II. Classification Problem ## In this example, we are using the iris dataset that ## contains four input and one output variables. ################################## ## Input data: Using the Iris dataset data(iris) set.seed(2) ## Shuffle the data ## then split the data to be training and testing datasets irisShuffled <- iris[sample(nrow(iris)), ] irisShuffled[, 5] <- unclass(irisShuffled[, 5]) tra.iris <- irisShuffled[1 : 105, ] tst.iris <- irisShuffled[106 : nrow(irisShuffled), 1 : 4] real.iris <- matrix(irisShuffled[106 : nrow(irisShuffled), 5], ncol = 1) ## Define range of input data. Note that it is only for the input variables. range.data.input <- apply(iris[, -ncol(iris)], 2, range) ## Set the method and its parameters. In this case we use FRBCS.W algorithm method.type <- "FRBCS.W" control <- list(num.labels = 7, type.mf = "GAUSSIAN", type.tnorm = "MIN", type.snorm = "MAX", type.implication.func = "ZADEH") ## Learning step: Generate fuzzy model object.cls <- frbs.learn(tra.iris, range.data.input, method.type, control) ## Predicting step: Predict newdata res.test <- predict(object.cls, tst.iris) ## Display the FRBS model summary(object.cls) ## Plot the membership functions plotMF(object.cls) ################################################# ## III. Constructing an FRBS model from human expert. ## In this example, we only consider the Mamdani model for regression. However, ## other models can be done in the same way. ## Note: ## In the examples, let us consider four input and one output variables. ################################################# ## Define a matrix representing shape and parameters of membership functions of input variables. ## The matrix has 5 rows where the first row represent the type of the membership function whereas ## others are values of its parameters. ## Detailed explanation can be seen in the fuzzifier function to construct the matrix. varinp.mf <- matrix(c(2, 0, 20, 40, NA, 4, 20, 40, 60, 80, 3, 60, 80, 100, NA, 2, 0, 35, 75, NA, 3, 35, 75, 100, NA, 2, 0, 20, 40, NA, 1, 20, 50, 80, NA, 3, 60, 80, 100, NA, 2, 0, 20, 40, NA, 4, 20, 40, 60, 80, 3, 60, 80, 100, NA), nrow = 5, byrow = FALSE) ## Define number of linguistic terms of input variables. ## Suppose, we have 3, 2, 3, and 3 numbers of linguistic terms ## for the first, second, third and fourth variables, respectively. num.fvalinput <- matrix(c(3, 2, 3, 3), nrow=1) ## Give the names of the linguistic terms of each input variables. varinput.1 <- c("low", "medium", "high") varinput.2 <- c("yes", "no") varinput.3 <- c("bad", "neutral", "good") varinput.4 <- c("low", "medium", "high") names.varinput <- c(varinput.1, varinput.2, varinput.3, varinput.4) ## Set interval of data. range.data <- matrix(c(0, 100, 0, 100, 0, 100, 0, 100, 0, 100), nrow = 2) ## Define inference parameters. ## Detailed information about values can be seen in the inference function. type.defuz <- "WAM" type.tnorm <- "MIN" type.snorm <- "MAX" type.implication.func <- "ZADEH" ## Give the name of simulation. name <- "Sim-0" ## Provide new data for testing. newdata<- matrix(c(25, 40, 35, 15, 45, 75, 78, 70), nrow = 2, byrow = TRUE) ## the names of variables colnames.var <- c("input1", "input2", "input3", "input4", "output1") ## Define number of linguistic terms of output variable. ## In this case, we set the number of linguistic terms to 3. num.fvaloutput <- matrix(c(3), nrow = 1) ## Give the names of the linguistic terms of the output variable. varoutput.1 <- c("bad", "neutral", "good") names.varoutput <- c(varoutput.1) ## Define the shapes and parameters of the membership functions of the output variables. varout.mf <- matrix(c(2, 0, 20, 40, NA, 4, 20, 40, 60, 80, 3, 60, 80, 100, NA), nrow = 5, byrow = FALSE) ## Set type of model which is "MAMDANI". type.model <- "MAMDANI" ## Define the fuzzy IF-THEN rules; ## In this example we are using the Mamdani model ## Note: e.g., ## "a1", "and", "b1, "->", "e1" means that ## "IF inputvar.1 is a1 and inputvar.2 is b1 THEN outputvar.1 is e1" ## Make sure that each rule has a "->" sign. rule <- matrix( c("low", "and", "yes", "and", "bad", "and", "low", "->", "bad", "medium", "and", "no", "and", "neutral", "and", "medium", "->", "neutral", "high", "and", "no", "and", "neutral", "and", "low", "->", "good"), nrow = 3, byrow = TRUE) ## Generate a fuzzy model with frbs.gen. object <- frbs.gen(range.data, num.fvalinput, names.varinput, num.fvaloutput, varout.mf, names.varoutput, rule, varinp.mf, type.model, type.defuz, type.tnorm, type.snorm, func.tsk = NULL, colnames.var, type.implication.func, name) ## Plot the membership function. plotMF(object) ## Predicting using new data. res <- predict(object, newdata)$predicted.val ################################################# ## IV. Specifying an FRBS model in the frbsPMML format. ## other examples can be seen in the frbsPMML function. ################################################# ## Input data data(frbsData) data.train <- frbsData$GasFurnance.dt[1 : 204, ] data.fit <- data.train[, 1 : 2] data.tst <- frbsData$GasFurnance.dt[205 : 292, 1 : 2] real.val <- matrix(frbsData$GasFurnance.dt[205 : 292, 3], ncol = 1) range.data<-matrix(c(-2.716, 2.834, 45.6, 60.5, 45.6, 60.5), nrow = 2) ## Set the method and its parameters method.type <- "WM" control <- list(num.labels = 3, type.mf = "GAUSSIAN", type.defuz = "WAM", type.tnorm = "MIN", type.snorm = "MAX", type.implication.func = "ZADEH", name="sim-0") ## Generate fuzzy model object <- frbs.learn(data.train, range.data, method.type, control) ## 2. Constructing the frbsPMML format frbsPMML(object)
################################## ## I. Regression Problem ## In this example, we are using the gas furnace dataset that ## contains two input and one output variables. ################################## ## Input data: Using the Gas Furnace dataset ## then split the data to be training and testing datasets data(frbsData) data.train <- frbsData$GasFurnance.dt[1 : 204, ] data.tst <- frbsData$GasFurnance.dt[205 : 292, 1 : 2] real.val <- matrix(frbsData$GasFurnance.dt[205 : 292, 3], ncol = 1) ## Define interval of data range.data <-apply(data.train, 2, range) ## Set the method and its parameters, ## for example, we use Wang and Mendel's algorithm method.type <- "WM" control <- list(num.labels = 15, type.mf = "GAUSSIAN", type.defuz = "WAM", type.tnorm = "MIN", type.snorm = "MAX", type.implication.func = "ZADEH", name = "sim-0") ## Learning step: Generate an FRBS model object.reg <- frbs.learn(data.train, range.data, method.type, control) ## Predicting step: Predict for newdata res.test <- predict(object.reg, data.tst) ## Display the FRBS model summary(object.reg) ## Plot the membership functions plotMF(object.reg) ################################## ## II. Classification Problem ## In this example, we are using the iris dataset that ## contains four input and one output variables. ################################## ## Input data: Using the Iris dataset data(iris) set.seed(2) ## Shuffle the data ## then split the data to be training and testing datasets irisShuffled <- iris[sample(nrow(iris)), ] irisShuffled[, 5] <- unclass(irisShuffled[, 5]) tra.iris <- irisShuffled[1 : 105, ] tst.iris <- irisShuffled[106 : nrow(irisShuffled), 1 : 4] real.iris <- matrix(irisShuffled[106 : nrow(irisShuffled), 5], ncol = 1) ## Define range of input data. Note that it is only for the input variables. range.data.input <- apply(iris[, -ncol(iris)], 2, range) ## Set the method and its parameters. In this case we use FRBCS.W algorithm method.type <- "FRBCS.W" control <- list(num.labels = 7, type.mf = "GAUSSIAN", type.tnorm = "MIN", type.snorm = "MAX", type.implication.func = "ZADEH") ## Learning step: Generate fuzzy model object.cls <- frbs.learn(tra.iris, range.data.input, method.type, control) ## Predicting step: Predict newdata res.test <- predict(object.cls, tst.iris) ## Display the FRBS model summary(object.cls) ## Plot the membership functions plotMF(object.cls) ################################################# ## III. Constructing an FRBS model from human expert. ## In this example, we only consider the Mamdani model for regression. However, ## other models can be done in the same way. ## Note: ## In the examples, let us consider four input and one output variables. ################################################# ## Define a matrix representing shape and parameters of membership functions of input variables. ## The matrix has 5 rows where the first row represent the type of the membership function whereas ## others are values of its parameters. ## Detailed explanation can be seen in the fuzzifier function to construct the matrix. varinp.mf <- matrix(c(2, 0, 20, 40, NA, 4, 20, 40, 60, 80, 3, 60, 80, 100, NA, 2, 0, 35, 75, NA, 3, 35, 75, 100, NA, 2, 0, 20, 40, NA, 1, 20, 50, 80, NA, 3, 60, 80, 100, NA, 2, 0, 20, 40, NA, 4, 20, 40, 60, 80, 3, 60, 80, 100, NA), nrow = 5, byrow = FALSE) ## Define number of linguistic terms of input variables. ## Suppose, we have 3, 2, 3, and 3 numbers of linguistic terms ## for the first, second, third and fourth variables, respectively. num.fvalinput <- matrix(c(3, 2, 3, 3), nrow=1) ## Give the names of the linguistic terms of each input variables. varinput.1 <- c("low", "medium", "high") varinput.2 <- c("yes", "no") varinput.3 <- c("bad", "neutral", "good") varinput.4 <- c("low", "medium", "high") names.varinput <- c(varinput.1, varinput.2, varinput.3, varinput.4) ## Set interval of data. range.data <- matrix(c(0, 100, 0, 100, 0, 100, 0, 100, 0, 100), nrow = 2) ## Define inference parameters. ## Detailed information about values can be seen in the inference function. type.defuz <- "WAM" type.tnorm <- "MIN" type.snorm <- "MAX" type.implication.func <- "ZADEH" ## Give the name of simulation. name <- "Sim-0" ## Provide new data for testing. newdata<- matrix(c(25, 40, 35, 15, 45, 75, 78, 70), nrow = 2, byrow = TRUE) ## the names of variables colnames.var <- c("input1", "input2", "input3", "input4", "output1") ## Define number of linguistic terms of output variable. ## In this case, we set the number of linguistic terms to 3. num.fvaloutput <- matrix(c(3), nrow = 1) ## Give the names of the linguistic terms of the output variable. varoutput.1 <- c("bad", "neutral", "good") names.varoutput <- c(varoutput.1) ## Define the shapes and parameters of the membership functions of the output variables. varout.mf <- matrix(c(2, 0, 20, 40, NA, 4, 20, 40, 60, 80, 3, 60, 80, 100, NA), nrow = 5, byrow = FALSE) ## Set type of model which is "MAMDANI". type.model <- "MAMDANI" ## Define the fuzzy IF-THEN rules; ## In this example we are using the Mamdani model ## Note: e.g., ## "a1", "and", "b1, "->", "e1" means that ## "IF inputvar.1 is a1 and inputvar.2 is b1 THEN outputvar.1 is e1" ## Make sure that each rule has a "->" sign. rule <- matrix( c("low", "and", "yes", "and", "bad", "and", "low", "->", "bad", "medium", "and", "no", "and", "neutral", "and", "medium", "->", "neutral", "high", "and", "no", "and", "neutral", "and", "low", "->", "good"), nrow = 3, byrow = TRUE) ## Generate a fuzzy model with frbs.gen. object <- frbs.gen(range.data, num.fvalinput, names.varinput, num.fvaloutput, varout.mf, names.varoutput, rule, varinp.mf, type.model, type.defuz, type.tnorm, type.snorm, func.tsk = NULL, colnames.var, type.implication.func, name) ## Plot the membership function. plotMF(object) ## Predicting using new data. res <- predict(object, newdata)$predicted.val ################################################# ## IV. Specifying an FRBS model in the frbsPMML format. ## other examples can be seen in the frbsPMML function. ################################################# ## Input data data(frbsData) data.train <- frbsData$GasFurnance.dt[1 : 204, ] data.fit <- data.train[, 1 : 2] data.tst <- frbsData$GasFurnance.dt[205 : 292, 1 : 2] real.val <- matrix(frbsData$GasFurnance.dt[205 : 292, 3], ncol = 1) range.data<-matrix(c(-2.716, 2.834, 45.6, 60.5, 45.6, 60.5), nrow = 2) ## Set the method and its parameters method.type <- "WM" control <- list(num.labels = 3, type.mf = "GAUSSIAN", type.defuz = "WAM", type.tnorm = "MIN", type.snorm = "MAX", type.implication.func = "ZADEH", name="sim-0") ## Generate fuzzy model object <- frbs.learn(data.train, range.data, method.type, control) ## 2. Constructing the frbsPMML format frbsPMML(object)
This is the internal function that implements the adaptive-network-based
fuzzy inference system (ANFIS). It is used to solve regression tasks.
Users do not need to call it directly,
but just use frbs.learn
and predict
.
ANFIS(data.train, num.labels, max.iter = 10, step.size = 0.01, type.tnorm = "MIN", type.snorm = "MAX", type.implication.func = "ZADEH")
ANFIS(data.train, num.labels, max.iter = 10, step.size = 0.01, type.tnorm = "MIN", type.snorm = "MAX", type.implication.func = "ZADEH")
data.train |
a matrix ( |
num.labels |
a matrix ( |
max.iter |
the maximal number of iterations. |
step.size |
a real number between 0 and 1 representing the step size of the gradient descent. |
type.tnorm |
the type of t-norm. For more detail, please have a look at |
type.snorm |
the type of s-norm. For more detail, please have a look at |
type.implication.func |
a value representing the type of implication functions.
For more detail, please have a look at |
This method was proposed by J. S. R. Jang. It uses the Takagi Sugeno Kang model on the consequent part of the fuzzy IF-THEN rules. The ANFIS architecture consists of two processes, the forward and the backward stage. The forward stage has five layers as follows:
Layer 1: The fuzzification process which transforms crisp values into linguistic terms using the Gaussian function as the shape of the membership function.
Layer 2: The inference stage using the t-norm operator (the AND operator).
Layer 3: Calculating the ratio of the strengths of the rules.
Layer 4: Calculating the consequent parameters.
Layer 5: Calculating the overall output as the sum of all incoming signals.
The backward stage is a process of parameter learning. In this step, the least squares method is used in order to obtain the parameters, which are coefficients of linear equations on the consequent part, and mean and variance on the antecedent part.
J.S.R. Jang, "ANFIS: Adaptive-network-based fuzzy inference system", IEEE Transactions on Systems, Man, and Cybernetics, vol. 23, no. 3, pp. 665 - 685 (1993).
J.S.R. Jang, C.T. Sun, and E. Mizutani., "Neuro-fuzzy and soft computing: a computational approach to learning and machine intelligence", Prentice-Hall, Inc (1997).
ANFIS.update
, frbs.learn
, and predict
The role of this function is to update parameters in the ANFIS method.
This function is called by the main function of the ANFIS method, ANFIS
.
ANFIS.update(data.train, def, rule.data.num, miu.rule, func.tsk, varinp.mf, step.size = 0.01)
ANFIS.update(data.train, def, rule.data.num, miu.rule, func.tsk, varinp.mf, step.size = 0.01)
data.train |
a matrix ( |
def |
a predicted value |
rule.data.num |
a matrix containing the rule base in integer form. |
miu.rule |
a matrix with the degrees of rules. See |
func.tsk |
a matrix of parameters of the function on the consequent part using the Takagi Sugeno Kang model. |
varinp.mf |
a matrix of parameters of membership functions of the input variables. |
step.size |
a real number between 0 and 1 representing the step size of the gradient descent. |
The purpose of this function is to generate data, which contains two input variables and one output variable, automatically for all values on a plane.
data.gen3d(range.input, num.grid = 10)
data.gen3d(range.input, num.grid = 10)
range.input |
the range of the input variables, as a matrix ( |
num.grid |
a number representing the size of the grid on the plane. |
the data
range.input <- matrix(c(0, 100, 0, 100), nrow=2) num.grid <- 10 data.test <- data.gen3d(range.input, num.grid)
range.input <- matrix(c(0, 100, 0, 100), nrow=2) num.grid <- 10 data.test <- data.gen3d(range.input, num.grid)
Defuzzification is a transformation that extracts the crisp values from the linguistic terms.
defuzzifier(data, rule = NULL, range.output = NULL, names.varoutput = NULL, varout.mf = NULL, miu.rule, type.defuz = NULL, type.model = "TSK", func.tsk = NULL)
defuzzifier(data, rule = NULL, range.output = NULL, names.varoutput = NULL, varout.mf = NULL, miu.rule, type.defuz = NULL, type.model = "TSK", func.tsk = NULL)
data |
a matrix ( |
rule |
a list or matrix of fuzzy IF-THEN rules, as discussed in |
range.output |
a matrix ( |
names.varoutput |
a list for giving names to the linguistic terms. See |
varout.mf |
a matrix constructing the membership function of the output variable.
See |
miu.rule |
the results of the inference module. See |
type.defuz |
the type of defuzzification to be used as follows.
|
type.model |
the type of the model that will be used in the simulation.
Here, |
func.tsk |
a matrix used to build the linear equation for the consequent part
if we are using Takagi Sugeno Kang. See also |
In this function, there exist two kinds of models which are based on the Mamdani and
Takagi Sugeno Kang model.
For the Mamdani model there are five methods for defuzzifying a linguistic term of a universe
of discourse
.
They are as follows:
weighted average method (WAM
).
first of maxima (FIRST.MAX
).
last of maxima (LAST.MAX
)
mean of maxima (MEAN.MAX
).
modified center of gravity (COG
).
A matrix of crisp values
fuzzifier
, rulebase
, and inference
This is the internal function that implements the dynamic evolving neural-fuzzy inference system (DENFIS).
It is used to handle regression tasks. Users do not need to call it directly,
but just use frbs.learn
and predict
.
DENFIS(data.train, range.data.ori, Dthr = 0.1, max.iter = 100, step.size = 0.01, d = 2)
DENFIS(data.train, range.data.ori, Dthr = 0.1, max.iter = 100, step.size = 0.01, d = 2)
data.train |
a matrix ( |
range.data.ori |
a matrix ( |
Dthr |
the threshold value for the evolving clustering method (ECM), between 0 and 1. |
max.iter |
the maximal number of iterations. |
step.size |
the step size of the least squares method, between 0 and 1. |
d |
a parameter for the width of the triangular membership function. |
This method was proposed by Nikola K. Kasabov and Q. Song. There are several steps in this method that are to determine the cluster centers using the evolving clustering method (ECM), to partition the input space and to find optimal parameters on the consequent part (Takagi Sugeno Kang model) for the IF-THEN rule using a least squares estimator.
ECM is a distance-based clustering method which is determined by a threshold value, Dthr
. This parameter
influences how many clusters are created. In the beginning of the clustering process, the first instance from the
training data is chosen to be a cluster center, and the determining radius is set to zero. Afterwards, using the
next instance, cluster centers and radius are changed based on certain mechanisms of ECM (please see ECM
).
All of the cluster centers are then obtained after evaluating all the training data.
The next step is to update the parameters on the consequent part with the assumption that the antecedent part which we got from ECM is fixed.
Actually, ECM can perform well as an online clustering method, but in this package it is used in an offline mode.
N.K. Kasabov and Q. Song, "DENFIS: Dynamic evolving neural-fuzzy inference system and its Application for time-series prediction", IEEE Transactions on Fuzzy Systems, vol. 10, no. 2, pp. 144 - 154 (2002).
DENFIS.eng
, frbs.learn
, and predict
This function is an internal function for the prediction phase using the DENFIS method.
The user should use this function not directly, but with calling predict
.
DENFIS.eng(object, newdata)
DENFIS.eng(object, newdata)
object |
the frbs model. See |
newdata |
a matrix ( |
a matrix of predicted values
This function is to transform from normalized data into real-valued data.
denorm.data(dt.norm, range.data, min.scale = 0, max.scale = 1)
denorm.data(dt.norm, range.data, min.scale = 0, max.scale = 1)
dt.norm |
a matrix ( |
range.data |
a matrix ( |
min.scale |
the minimum value within normalization. |
max.scale |
the maximum value within normalization. |
the real-valued data
The role of this function is to update the parameters of the fuzzy inference rules by descent method (FIR.DM).
This function is called by the main function of the FIR.DM method, FIR.DM
.
DM.update(data.train, rule.data.num, miu.rule, func.tsk, varinp.mf, step.size = 0.01, def)
DM.update(data.train, rule.data.num, miu.rule, func.tsk, varinp.mf, step.size = 0.01, def)
data.train |
a matrix ( |
rule.data.num |
a matrix containing the rulebase. Its elements are integers, see |
miu.rule |
a matrix with the degrees of rules which is a result of the |
func.tsk |
a matrix of parameters of the functions on the consequent part of the Takagi Sugeno Kang model. |
varinp.mf |
a matrix of parameters of the membership functions of the input variables. |
step.size |
the step size of the descent method, between 0 and 1. |
def |
a matrix which is obtained from the defuzzification. Please have a look at |
frbs.learn
, predict
, and FIR.DM
.
This function is a part of the DENFIS method to generate cluster centers.
ECM(data.train, Dthr)
ECM(data.train, Dthr)
data.train |
a matrix ( |
Dthr |
the threshold value for the evolving clustering method (ECM), between 0 and 1. |
a matrix of cluster centers
DENFIS
and DENFIS.eng
This is the internal function that implements the Ishibuchi's method based on
hybridization of genetic cooperative-competitive learning (GCCL) and Pittsburgh (FH.GBML). It is used to solve classification tasks.
Users do not need to call it directly,
but just use frbs.learn
and predict
.
FH.GBML(data.train, popu.size = 10, max.num.rule = 5, persen_cross = 0.6, persen_mutant = 0.3, max.gen = 10, num.class, range.data.input, p.dcare = 0.5, p.gccl = 0.5)
FH.GBML(data.train, popu.size = 10, max.num.rule = 5, persen_cross = 0.6, persen_mutant = 0.3, max.gen = 10, num.class, range.data.input, p.dcare = 0.5, p.gccl = 0.5)
data.train |
a matrix ( |
popu.size |
the size of the population which is generated in each generation. |
max.num.rule |
the maximum number of rules. |
persen_cross |
a real number between 0 and 1 determining the probability of crossover. |
persen_mutant |
a real number between 0 and 1 determining the probability of mutation. |
max.gen |
the maximal number of generations for the genetic algorithms. |
num.class |
a number of the classes. |
range.data.input |
a matrix containing the ranges of the normalized input data. |
p.dcare |
a probability of "don't care" attributes occurred. |
p.gccl |
a probability of GCCL process occurred. |
This method is based on Ishibuchi's method using the hybridization of GCCL and the Pittsburgh approach for genetic fuzzy systems. The algorithm of this method is as follows:
Step 1: Generate population where each individual in the population is a fuzzy rule set.
Step 2: Calculate the fitness value of each rule set in the current population.
Step 3: Generate new rule sets by the selection, crossover, and mutation in the same manner as the Pittsburgh-style algorithm. Then, apply iterations of the GCCL to each of the generated rule sets with a probability.
Step 4: Add the best rule set in the current population to newly generated rule sets to form the next population.
Step 5: Return to Step 2 if the prespecified stopping condition is not satisfied.
H. Ishibuchi, T. Yamamoto, and T. Nakashima, "Hybridization of fuzzy GBML approaches for pattern classification problems," IEEE Trans. on Systems, Man, and Cybernetics-Part B: Cybernetics, vol. 35, no. 2, pp. 359 - 365 (2005).
This is the internal function that implements the fuzzy inference rules by descent method (FIR.DM).
It is used to solve regression tasks. Users do not need to call it directly,
but just use frbs.learn
and predict
.
FIR.DM(data.train, num.labels, max.iter, step.size, type.tnorm = "MIN", type.snorm = "MAX", type.implication.func = "ZADEH")
FIR.DM(data.train, num.labels, max.iter, step.size, type.tnorm = "MIN", type.snorm = "MAX", type.implication.func = "ZADEH")
data.train |
a matrix ( |
num.labels |
a matrix ( |
max.iter |
the maximal number of iterations. |
step.size |
the step size of the descent method, between 0 and 1. |
type.tnorm |
the type of t-norm. For more detail, please have a look at |
type.snorm |
the type of s-norm. For more detail, please have a look at |
type.implication.func |
a value representing type of implication function. For more detail, please have a look at |
This method was proposed by H. Nomura, I. Hayashi, and N. Wakami. FIR.DM uses simplified fuzzy reasoning where the consequent part is a real number (a particular case within the Takagi Sugeno Kang model), while the membership function on the antecedent part is expressed by an isosceles triangle. So, in the learning phase, FIR.DM updates three parameters which are center and width of the triangular and a real number on the consequent part using a descent method.
H. Nomura, I. Hayashi and N. Wakami, "A learning method of fuzzy inference rules by descent method", IEEE International Conference on Fuzzy Systems, pp. 203 - 210 (1992).
DM.update
, frbs.learn
, and predict
.
This is the internal function that implements the fuzzy rule-based classification
system using Chi's technique (FRBCS.CHI). It is used to solve classification tasks.
Users do not need to call it directly,
but just use frbs.learn
and predict
. This method is
suitable only for classification problems.
FRBCS.CHI(range.data, data.train, num.labels, num.class, type.mf = "TRIANGLE", type.tnorm = "MIN", type.snorm = "MAX", type.implication.func = "ZADEH")
FRBCS.CHI(range.data, data.train, num.labels, num.class, type.mf = "TRIANGLE", type.tnorm = "MIN", type.snorm = "MAX", type.implication.func = "ZADEH")
range.data |
a matrix ( |
data.train |
a matrix ( |
num.labels |
a matrix ( |
num.class |
an integer number representing the number of labels (linguistic terms). |
type.mf |
the type of the shape of the membership functions. See |
type.tnorm |
the type of t-norm. See |
type.snorm |
the type of s-norm. See |
type.implication.func |
the type of implication function. See |
This method was proposed by Z. Chi, H. Yan, and T. Pham that extends
Wang and Mendel's method for tackling classification problems.
Basically, the algorithm is quite similar as Wang and Mendel's technique.
However, since it is based on the FRBCS model, Chi's method only takes class labels on each data
to be consequent parts of fuzzy IF-THEN rules. In other words, we generate rules as in
Wang and Mendel's technique (WM
) and then we replace consequent parts with their classes.
Regarding calculating degress of each rule, they are determined by antecedent parts of the rules.
Redudant rules can be deleted by considering their degrees. Lastly, we obtain fuzzy IF-THEN rules
based on the FRBCS model.
Z. Chi, H. Yan, T. Pham, "Fuzzy algorithms with applications to image processing and pattern recognition", World Scientific, Singapore (1996).
FRBCS.eng
, frbs.learn
, and predict
This function is the internal function of the fuzzy rule-based classification systems (FRBCS) to compute the predicted values.
FRBCS.eng(object, newdata)
FRBCS.eng(object, newdata)
object |
the |
newdata |
a matrix ( |
A matrix of predicted values.
This is the internal function that implements the fuzzy rule-based classification
system with weight factor (FRBCS.W). It is used to solve classification tasks.
Users do not need to call it directly,
but just use frbs.learn
and predict
. This method is
suitable only for classification problems.
FRBCS.W(range.data, data.train, num.labels, num.class, type.mf, type.tnorm = "MIN", type.snorm = "MAX", type.implication.func = "ZADEH")
FRBCS.W(range.data, data.train, num.labels, num.class, type.mf, type.tnorm = "MIN", type.snorm = "MAX", type.implication.func = "ZADEH")
range.data |
a matrix ( |
data.train |
a matrix ( |
num.labels |
a matrix ( |
num.class |
an integer number representing the number of labels (linguistic terms). |
type.mf |
the type of the shape of the membership functions. |
type.tnorm |
the type of t-norm. See |
type.snorm |
the type of s-norm. See |
type.implication.func |
the type of implication function. See |
This method is adopted from Ishibuchi and Nakashima's paper. Each fuzzy IF-THEN rule consists of antecedent linguistic values and a single consequent class with certainty grades (weights). The antecedent part is determined by a grid-type fuzzy partition from the training data. The consequent class is defined as the dominant class in the fuzzy subspace corresponding to the antecedent part of each fuzzy IF-THEN rule and the certainty grade is calculated from the ratio among the consequent class. A class of the new instance is determined by the consequent class of the rule with the maximal product of the compatibility grade and the certainty grade.
H. Ishibuchi and T. Nakashima, "Effect of rule weights in fuzzy rule-based classification systems", IEEE Transactions on Fuzzy Systems, vol. 1, pp. 59 - 64 (2001).
FRBCS.eng
, frbs.learn
, and predict
This function is one of the main internal functions of the package. It determines the values within the prediction phase.
frbs.eng(object, newdata)
frbs.eng(object, newdata)
object |
the |
newdata |
a matrix ( |
This function involves four different processing steps on fuzzy rule-based systems.
Firstly, the rulebase (see rulebase
) validates
the consistency of the fuzzy IF-THEN rules form. Then, the fuzzification
(see fuzzifier
) transforms crisp values
into linguistic terms. Next, the inference calculates the degree of rule strengths using
the t-norm and the s-norm.
Finally, the defuzzification process calculates the results of the model using the Mamdani
or the Takagi Sugeno Kang model.
A list with the following items:
rule |
the fuzzy IF-THEN rules |
varinp.mf |
a matrix to generate the shapes of the membership functions for the input variables |
MF |
a matrix of the degrees of the membership functions |
miu.rule |
a matrix of the degrees of the rules |
func.tsk |
a matrix of the Takagi Sugeno Kang model for the consequent part of the fuzzy IF-THEN rules |
predicted.val |
a matrix of the predicted values |
fuzzifier
, rulebase
, inference
and defuzzifier
.
The purpose of this function is to generate a FRBS model from user-given input without a learning process.
frbs.gen(range.data, num.fvalinput, names.varinput, num.fvaloutput = NULL, varout.mf = NULL, names.varoutput = NULL, rule, varinp.mf, type.model = "MAMDANI", type.defuz = "WAM", type.tnorm = "MIN", type.snorm = "MAX", func.tsk = NULL, colnames.var = NULL, type.implication.func = "ZADEH", name = "Sim-0")
frbs.gen(range.data, num.fvalinput, names.varinput, num.fvaloutput = NULL, varout.mf = NULL, names.varoutput = NULL, rule, varinp.mf, type.model = "MAMDANI", type.defuz = "WAM", type.tnorm = "MIN", type.snorm = "MAX", func.tsk = NULL, colnames.var = NULL, type.implication.func = "ZADEH", name = "Sim-0")
range.data |
a matrix ( |
num.fvalinput |
a matrix representing the number of linguistic terms of each input variables. For example: means that there are two variables where the first variable has three linguistic terms and the second one has two linguistic terms. |
names.varinput |
a list containing names to the linguistic terms for input variables. See |
num.fvaloutput |
the number of linguistic terms of the output variable. This parameter is required for the Mamdani model only. For example: means there are 3 linguistic terms for the output variable. |
varout.mf |
a matrix for constructing the membership functions of the output variable.
The form is the same as for the |
names.varoutput |
a list giving names of the linguistic terms for the output variable. The form is the same as
for the |
rule |
a list of fuzzy IF-THEN rules. There are some types of rule structures, for example: Mamdani, Takagi Sugeno Kang,
and fuzzy rule-based classification systems (FRBCS). If we use the Mamdani model then the consequent part is a linguistic term,
but if we use Takagi Sugeno Kang then we build a matrix representing linear equations in the consequent part.
e.g., "a1", "and", "b1, "->", "e1" means that
"IF inputvar.1 is a1 and inputvar.2 is b1 THEN outputvar.1 is e1".
Make sure that each rule has a "->" sign.
Furthermore, we are allowed to use linguistic hedges (e.g., "extremely", "slightly", etc), negation (i.e., "not"),
and the "dont_care" value representing degree of membership is always 1.
For more detail, see |
varinp.mf |
a matrix for constructing the shapes of the membership functions. See how to construct it in |
type.model |
the type of the model. There are three types available as follows.
|
type.defuz |
the type of the defuzzification method. It is used in the Mamdani model only.
See |
type.tnorm |
the type of the t-norm method. See |
type.snorm |
the type of the s-norm method. See |
func.tsk |
a matrix of parameters of the function on the consequent part using the Takagi Sugeno Kang model.
This parameter must be defined when we are using Takagi Sugeno Kang. See |
colnames.var |
a list of names of input and output variables. |
type.implication.func |
a type of implication function. See |
name |
a name of the simulation. |
It can be used if rules have already been obtained manually, without employing the
learning process.
In the examples shown, we generate a fuzzy model using frbs.gen
and generate the
fuzzy rule-based systems step by step manually. Additionally, the examples show several scenarios as follows.
Using frbs.gen
for constructing the Mamdani model on a regression task.
Using frbs.gen
for constructing the Takagi Sugeno Kang model on a regression task.
Constructing the Mamdani model by executing internal functions such as rulebase
, fuzzifier
,
inference
, and defuzzifier
for the Mamdani model.
Using frbs.gen
for constructing fuzzy rule-based classification systems (FRBCS) model.
The frbs-object
.
################################################# ## 1. The following codes show how to generate a fuzzy model ## using the frbs.gen function for regression tasks. ## The following are three scenarios: ## 1a. Using the Mamdani model ## 1b. Using the Takagi Sugeno Kang model ## 1c. Using the Mamdani model and internal functions: fuzzifier, etc. ## Note: ## In the examples, let us consider four input variabels and one output variable. ## Some variables could be shared together for other examples. ################################################# ## Define shape and parameters of membership functions of input variables. ## Please see the fuzzifier function to construct the matrix. ## It can be seen that in this case we employ TRAPEZOID as the membership functions. varinp.mf <- matrix(c(2, 0, 20, 40, NA, 4, 20, 40, 60, 80, 3, 60, 80, 100, NA, 2, 0, 35, 75, NA, 3, 35, 75, 100, NA, 2, 0, 20, 40, NA, 1, 20, 50, 80, NA, 3, 60, 80, 100, NA, 2, 0, 20, 40, NA, 4, 20, 40, 60, 80, 3, 60, 80, 100, NA), nrow = 5, byrow = FALSE) ## Define number of linguistic terms of the input variables. ## Suppose, we have 3, 2, 3, and 3 numbers of linguistic terms ## for the first, second, third and fourth variables, respectively. num.fvalinput <- matrix(c(3, 2, 3, 3), nrow=1) ## Give the names of the linguistic terms of each input variables. varinput.1 <- c("a1", "a2", "a3") varinput.2 <- c("b1", "b2") varinput.3 <- c("c1", "c2", "c3") varinput.4 <- c("d1", "d2", "d3") names.varinput <- c(varinput.1, varinput.2, varinput.3, varinput.4) ## Set interval of data. range.data <- matrix(c(0,100, 0, 100, 0, 100, 0, 100, 0, 100), nrow=2) ## Define inference parameters. type.defuz <- "WAM" type.tnorm <- "MIN" type.snorm <- "MAX" type.implication.func <- "ZADEH" ## Give the name of simulation. name <- "Sim-0" ## Provide new data for testing. newdata <- matrix(c(15, 80, 85, 85, 45, 75, 78, 70), nrow = 2, byrow = TRUE) ## the names of variables colnames.var <- c("input1", "input2", "input3", "input4", "output1") ################################################################### ## 1a. Using the Mamdani Model #################################################################### ## Define number of linguistic terms of output variable. ## In this case, we set the number of linguistic terms to 3. num.fvaloutput <- matrix(c(3), nrow = 1) ## Give the names of the linguistic terms of the output variable. varoutput.1 <- c("e1", "e2", "e3") names.varoutput <- c(varoutput.1) ## Define the shapes and parameters of the membership functions of the output variables. varout.mf <- matrix(c(2, 0, 20, 40, NA, 4, 20, 40, 60, 80, 3, 60, 80, 100, NA), nrow = 5, byrow = FALSE) ## Set type of model which is "MAMDANI" or "TSK" for Mamdani or ## Takagi Sugeno Kang models, respectively. ## In this case, we choose the Mamdani model. type.model <- "MAMDANI" ## Define the fuzzy IF-THEN rules; In this case, we provide two scenarios using different operators: rule.or <- matrix(c("a1", "or", "b1", "or", "c1", "or", "d1", "->", "e1", "a2", "and", "b2", "and", "c2", "and", "d2", "->", "e2", "a3", "and", "b2", "and", "c2", "and", "d1", "->", "e3"), nrow = 3, byrow = TRUE) ## Define the fuzzy IF-THEN rules; rule.and <- matrix(c("a1", "and", "b1", "and", "c1", "and", "d1", "->", "e1", "a2", "and", "b2", "and", "c2", "and", "d2", "->", "e2", "a3", "and", "b2", "and", "c2", "and", "d1", "->", "e3"), nrow = 3, byrow = TRUE) ## Generate a fuzzy model with frbs.gen. object.or <- frbs.gen(range.data, num.fvalinput, names.varinput, num.fvaloutput, varout.mf, names.varoutput, rule.or, varinp.mf, type.model, type.defuz, type.tnorm, type.snorm, func.tsk = NULL, colnames.var, type.implication.func, name) object.and <- frbs.gen(range.data, num.fvalinput, names.varinput, num.fvaloutput, varout.mf, names.varoutput, rule.and, varinp.mf, type.model, type.defuz, type.tnorm, type.snorm, func.tsk = NULL, colnames.var, type.implication.func, name) ## Plot the membership function. plotMF(object.and) ## Predicting using new data. res.or <- predict(object.or, newdata)$predicted.val res.and <- predict(object.and, newdata)$predicted.val ##################################################################### ## 1b. Using the Takagi Sugeno Kang (TSK) Model ##################################################################### ## Define "TSK" for the Takagi Sugeno Kang model type.model <- "TSK" ## Define linear equations for consequent parts. ## The following command means that we have three equation related to the rules we have. ## e.g., the first equation is 1*inputvar.1 + 1*inputvar.2 + 5*inputvar.3 + 2*inputvar.4 + 1, ## where inputvar.i is a value of the i-th input variable. func.tsk <- matrix(c(1, 1, 5, 2, 1, 3, 1, 0.5, 0.1, 2, 1, 3, 2, 2, 2), nrow = 3, byrow = TRUE) ## Define the fuzzy IF-THEN rules; ## For TSK model, it isn't necessary to put linguistic term in consequent parts. ## Make sure that each rule has a "->" sign. rule <- matrix(c("a1", "and", "b1", "and", "c1", "and", "d1", "->", "a2", "and", "b2", "and", "c2", "and", "d2", "->", "a3", "and", "b2", "and", "c2", "and", "d1", "->"), nrow = 3, byrow = TRUE) ## Generate a fuzzy model with frbs.gen. ## It should be noted that for TSK model, we do not need to input: ## num.fvaloutput, varout.mf, names.varoutput, type.defuz. object <- frbs.gen(range.data, num.fvalinput, names.varinput, num.fvaloutput = NULL, varout.mf = NULL, names.varoutput = NULL, rule, varinp.mf, type.model, type.defuz = NULL, type.tnorm, type.snorm, func.tsk, colnames.var, type.implication.func, name) ## Plot the membership function. plotMF(object) ## Predicting using new data. res <- predict(object, newdata)$predicted.val ###################### ## 1c. Using the same data as in the previous example, this example performs ## step by step of the generation of a fuzzy rule-based system ###################### ## Using the Mamdani model. type.model <- "MAMDANI" ## Construct rules. rule <- matrix(c("a1", "and", "b1", "and", "c1", "and", "d1", "->", "e1", "a2", "and", "b2", "and", "c2", "and", "d2", "->", "e2", "a3", "and", "b2", "and", "c2", "and", "d1", "->", "e3"), nrow = 3, byrow = TRUE) ## Check input data given by user. rule <- rulebase(type.model, rule, func.tsk = NULL) ## Fuzzification Module: ## In this function, we convert crisp into linguistic values/terms ## based on the data and the parameters of the membership function. ## The output: a matrix representing the degree of the membership of the data num.varinput <- ncol(num.fvalinput) MF <- fuzzifier(newdata, num.varinput, num.fvalinput, varinp.mf) ## Inference Module: ## In this function, we will calculate the confidence factor on the antecedent for each rule ## considering t-norm and s-norm. miu.rule <- inference(MF, rule, names.varinput, type.tnorm, type.snorm) ## Defuzzification Module. ## In this function, we calculate and convert the linguistic values back into crisp values. range.output <- range.data[, ncol(range.data), drop = FALSE] result <- defuzzifier(newdata, rule, range.output, names.varoutput, varout.mf, miu.rule, type.defuz, type.model, func.tsk = NULL) ################################################# ## 2. The following codes show how to generate a fuzzy model ## using the frbs.gen function for classification tasks using the Mamdani model. ################################################# ## define range of data. ## Note. we only define range of input data. range.data.input <- matrix(c(0, 1, 0, 1, 0, 1, 0, 1), nrow=2) ## Define shape and parameters of membership functions of input variables. ## Please see fuzzifier function to construct the matrix. ## In this case, we are using TRIANGLE for membership functions. varinp.mf <- matrix(c(1, 0, 0, 0.5, NA, 1, 0, 0.5, 1, NA, 1, 0.5, 1, 1, NA, 1, 0, 0, 0.5, NA, 1, 0, 0.5, 1, NA, 1, 0.5, 1, 1, NA, 1, 0, 0, 0.5, NA, 1, 0, 0.5, 1, NA, 1, 0.5, 1, 1, NA, 1, 0, 0, 0.5, NA, 1, 0, 0.5, 1, NA, 1, 0.5, 1, 1, NA), nrow = 5, byrow = FALSE) ## Define number of linguistic terms of input variables. ## Suppose, we have 3, 3, 3, and 3 numbers of linguistic terms ## for first up to fourth variables, respectively. num.fvalinput <- matrix(c(3, 3, 3, 3), nrow=1) ## Give the names of the linguistic terms of each input variable. varinput.1 <- c("v.1_a.1", "v.1_a.2", "v.1_a.3") varinput.2 <- c("v.2_a.1", "v.2_a.2", "v.2_a.3") varinput.3 <- c("v.3_a.1", "v.3_a.2", "v.3_a.3") varinput.4 <- c("v.4_a.1", "v.4_a.2", "v.4_a.3") names.varinput <- c(varinput.1, varinput.2, varinput.3, varinput.4) ## Provide inference parameters. type.tnorm <- "MIN" type.snorm <- "MAX" type.implication.func <- "ZADEH" type.model <- "FRBCS" ## Give the name of simulation. name <- "Sim-0" ## Provide new data for testing. newdata<- matrix(c(0.45, 0.5, 0.89, 0.44, 0.51, 0.99, 0.1, 0.98, 0.51, 0.56, 0.55, 0.5), nrow = 3, byrow = TRUE) ## the names of variables colnames.var <- c("input1", "input2", "input3", "input4", "output1") ## Construct rules. ## It should be noted that on consequent parts we define categorical values instead of ## linguistic terms. rule <- matrix( c("v.1_a.2", "and", "v.2_a.2", "and", "v.3_a.3", "and", "v.4_a.2", "->", "3", "v.1_a.2", "and", "v.2_a.3", "and", "v.3_a.1", "and", "v.4_a.3", "->", "1", "v.1_a.2", "and", "v.2_a.2", "and", "v.3_a.2", "and", "v.4_a.2", "->", "2"), nrow = 3, byrow = TRUE) ## Generate frbs object. object <- frbs.gen(range.data = range.data.input, num.fvalinput, names.varinput, num.fvaloutput = NULL, varout.mf = NULL, names.varoutput = NULL, rule, varinp.mf, type.model, type.defuz = NULL, type.tnorm, type.snorm, func.tsk = NULL, colnames.var, type.implication.func, name) ## Plot the shape of membership functions. plotMF(object) ## Predicting using new data. res <- predict(object, newdata) #################################################### ## 3. The following example shows how to convert ## the frbs model into frbsPMML #################################################### ## In this example, we are using the last object of FRBS. ## Display frbsPMML in R objPMML <- frbsPMML(object) ## Write into a file with .frbsPMML extention ## Not run: write.frbsPMML(objPMML, fileName="obj_frbsPMML") ## Read the frbsPMML file into an R object of FRBS obj <- read.frbsPMML("obj_frbsPMML.frbsPMML") ## End(Not run)
################################################# ## 1. The following codes show how to generate a fuzzy model ## using the frbs.gen function for regression tasks. ## The following are three scenarios: ## 1a. Using the Mamdani model ## 1b. Using the Takagi Sugeno Kang model ## 1c. Using the Mamdani model and internal functions: fuzzifier, etc. ## Note: ## In the examples, let us consider four input variabels and one output variable. ## Some variables could be shared together for other examples. ################################################# ## Define shape and parameters of membership functions of input variables. ## Please see the fuzzifier function to construct the matrix. ## It can be seen that in this case we employ TRAPEZOID as the membership functions. varinp.mf <- matrix(c(2, 0, 20, 40, NA, 4, 20, 40, 60, 80, 3, 60, 80, 100, NA, 2, 0, 35, 75, NA, 3, 35, 75, 100, NA, 2, 0, 20, 40, NA, 1, 20, 50, 80, NA, 3, 60, 80, 100, NA, 2, 0, 20, 40, NA, 4, 20, 40, 60, 80, 3, 60, 80, 100, NA), nrow = 5, byrow = FALSE) ## Define number of linguistic terms of the input variables. ## Suppose, we have 3, 2, 3, and 3 numbers of linguistic terms ## for the first, second, third and fourth variables, respectively. num.fvalinput <- matrix(c(3, 2, 3, 3), nrow=1) ## Give the names of the linguistic terms of each input variables. varinput.1 <- c("a1", "a2", "a3") varinput.2 <- c("b1", "b2") varinput.3 <- c("c1", "c2", "c3") varinput.4 <- c("d1", "d2", "d3") names.varinput <- c(varinput.1, varinput.2, varinput.3, varinput.4) ## Set interval of data. range.data <- matrix(c(0,100, 0, 100, 0, 100, 0, 100, 0, 100), nrow=2) ## Define inference parameters. type.defuz <- "WAM" type.tnorm <- "MIN" type.snorm <- "MAX" type.implication.func <- "ZADEH" ## Give the name of simulation. name <- "Sim-0" ## Provide new data for testing. newdata <- matrix(c(15, 80, 85, 85, 45, 75, 78, 70), nrow = 2, byrow = TRUE) ## the names of variables colnames.var <- c("input1", "input2", "input3", "input4", "output1") ################################################################### ## 1a. Using the Mamdani Model #################################################################### ## Define number of linguistic terms of output variable. ## In this case, we set the number of linguistic terms to 3. num.fvaloutput <- matrix(c(3), nrow = 1) ## Give the names of the linguistic terms of the output variable. varoutput.1 <- c("e1", "e2", "e3") names.varoutput <- c(varoutput.1) ## Define the shapes and parameters of the membership functions of the output variables. varout.mf <- matrix(c(2, 0, 20, 40, NA, 4, 20, 40, 60, 80, 3, 60, 80, 100, NA), nrow = 5, byrow = FALSE) ## Set type of model which is "MAMDANI" or "TSK" for Mamdani or ## Takagi Sugeno Kang models, respectively. ## In this case, we choose the Mamdani model. type.model <- "MAMDANI" ## Define the fuzzy IF-THEN rules; In this case, we provide two scenarios using different operators: rule.or <- matrix(c("a1", "or", "b1", "or", "c1", "or", "d1", "->", "e1", "a2", "and", "b2", "and", "c2", "and", "d2", "->", "e2", "a3", "and", "b2", "and", "c2", "and", "d1", "->", "e3"), nrow = 3, byrow = TRUE) ## Define the fuzzy IF-THEN rules; rule.and <- matrix(c("a1", "and", "b1", "and", "c1", "and", "d1", "->", "e1", "a2", "and", "b2", "and", "c2", "and", "d2", "->", "e2", "a3", "and", "b2", "and", "c2", "and", "d1", "->", "e3"), nrow = 3, byrow = TRUE) ## Generate a fuzzy model with frbs.gen. object.or <- frbs.gen(range.data, num.fvalinput, names.varinput, num.fvaloutput, varout.mf, names.varoutput, rule.or, varinp.mf, type.model, type.defuz, type.tnorm, type.snorm, func.tsk = NULL, colnames.var, type.implication.func, name) object.and <- frbs.gen(range.data, num.fvalinput, names.varinput, num.fvaloutput, varout.mf, names.varoutput, rule.and, varinp.mf, type.model, type.defuz, type.tnorm, type.snorm, func.tsk = NULL, colnames.var, type.implication.func, name) ## Plot the membership function. plotMF(object.and) ## Predicting using new data. res.or <- predict(object.or, newdata)$predicted.val res.and <- predict(object.and, newdata)$predicted.val ##################################################################### ## 1b. Using the Takagi Sugeno Kang (TSK) Model ##################################################################### ## Define "TSK" for the Takagi Sugeno Kang model type.model <- "TSK" ## Define linear equations for consequent parts. ## The following command means that we have three equation related to the rules we have. ## e.g., the first equation is 1*inputvar.1 + 1*inputvar.2 + 5*inputvar.3 + 2*inputvar.4 + 1, ## where inputvar.i is a value of the i-th input variable. func.tsk <- matrix(c(1, 1, 5, 2, 1, 3, 1, 0.5, 0.1, 2, 1, 3, 2, 2, 2), nrow = 3, byrow = TRUE) ## Define the fuzzy IF-THEN rules; ## For TSK model, it isn't necessary to put linguistic term in consequent parts. ## Make sure that each rule has a "->" sign. rule <- matrix(c("a1", "and", "b1", "and", "c1", "and", "d1", "->", "a2", "and", "b2", "and", "c2", "and", "d2", "->", "a3", "and", "b2", "and", "c2", "and", "d1", "->"), nrow = 3, byrow = TRUE) ## Generate a fuzzy model with frbs.gen. ## It should be noted that for TSK model, we do not need to input: ## num.fvaloutput, varout.mf, names.varoutput, type.defuz. object <- frbs.gen(range.data, num.fvalinput, names.varinput, num.fvaloutput = NULL, varout.mf = NULL, names.varoutput = NULL, rule, varinp.mf, type.model, type.defuz = NULL, type.tnorm, type.snorm, func.tsk, colnames.var, type.implication.func, name) ## Plot the membership function. plotMF(object) ## Predicting using new data. res <- predict(object, newdata)$predicted.val ###################### ## 1c. Using the same data as in the previous example, this example performs ## step by step of the generation of a fuzzy rule-based system ###################### ## Using the Mamdani model. type.model <- "MAMDANI" ## Construct rules. rule <- matrix(c("a1", "and", "b1", "and", "c1", "and", "d1", "->", "e1", "a2", "and", "b2", "and", "c2", "and", "d2", "->", "e2", "a3", "and", "b2", "and", "c2", "and", "d1", "->", "e3"), nrow = 3, byrow = TRUE) ## Check input data given by user. rule <- rulebase(type.model, rule, func.tsk = NULL) ## Fuzzification Module: ## In this function, we convert crisp into linguistic values/terms ## based on the data and the parameters of the membership function. ## The output: a matrix representing the degree of the membership of the data num.varinput <- ncol(num.fvalinput) MF <- fuzzifier(newdata, num.varinput, num.fvalinput, varinp.mf) ## Inference Module: ## In this function, we will calculate the confidence factor on the antecedent for each rule ## considering t-norm and s-norm. miu.rule <- inference(MF, rule, names.varinput, type.tnorm, type.snorm) ## Defuzzification Module. ## In this function, we calculate and convert the linguistic values back into crisp values. range.output <- range.data[, ncol(range.data), drop = FALSE] result <- defuzzifier(newdata, rule, range.output, names.varoutput, varout.mf, miu.rule, type.defuz, type.model, func.tsk = NULL) ################################################# ## 2. The following codes show how to generate a fuzzy model ## using the frbs.gen function for classification tasks using the Mamdani model. ################################################# ## define range of data. ## Note. we only define range of input data. range.data.input <- matrix(c(0, 1, 0, 1, 0, 1, 0, 1), nrow=2) ## Define shape and parameters of membership functions of input variables. ## Please see fuzzifier function to construct the matrix. ## In this case, we are using TRIANGLE for membership functions. varinp.mf <- matrix(c(1, 0, 0, 0.5, NA, 1, 0, 0.5, 1, NA, 1, 0.5, 1, 1, NA, 1, 0, 0, 0.5, NA, 1, 0, 0.5, 1, NA, 1, 0.5, 1, 1, NA, 1, 0, 0, 0.5, NA, 1, 0, 0.5, 1, NA, 1, 0.5, 1, 1, NA, 1, 0, 0, 0.5, NA, 1, 0, 0.5, 1, NA, 1, 0.5, 1, 1, NA), nrow = 5, byrow = FALSE) ## Define number of linguistic terms of input variables. ## Suppose, we have 3, 3, 3, and 3 numbers of linguistic terms ## for first up to fourth variables, respectively. num.fvalinput <- matrix(c(3, 3, 3, 3), nrow=1) ## Give the names of the linguistic terms of each input variable. varinput.1 <- c("v.1_a.1", "v.1_a.2", "v.1_a.3") varinput.2 <- c("v.2_a.1", "v.2_a.2", "v.2_a.3") varinput.3 <- c("v.3_a.1", "v.3_a.2", "v.3_a.3") varinput.4 <- c("v.4_a.1", "v.4_a.2", "v.4_a.3") names.varinput <- c(varinput.1, varinput.2, varinput.3, varinput.4) ## Provide inference parameters. type.tnorm <- "MIN" type.snorm <- "MAX" type.implication.func <- "ZADEH" type.model <- "FRBCS" ## Give the name of simulation. name <- "Sim-0" ## Provide new data for testing. newdata<- matrix(c(0.45, 0.5, 0.89, 0.44, 0.51, 0.99, 0.1, 0.98, 0.51, 0.56, 0.55, 0.5), nrow = 3, byrow = TRUE) ## the names of variables colnames.var <- c("input1", "input2", "input3", "input4", "output1") ## Construct rules. ## It should be noted that on consequent parts we define categorical values instead of ## linguistic terms. rule <- matrix( c("v.1_a.2", "and", "v.2_a.2", "and", "v.3_a.3", "and", "v.4_a.2", "->", "3", "v.1_a.2", "and", "v.2_a.3", "and", "v.3_a.1", "and", "v.4_a.3", "->", "1", "v.1_a.2", "and", "v.2_a.2", "and", "v.3_a.2", "and", "v.4_a.2", "->", "2"), nrow = 3, byrow = TRUE) ## Generate frbs object. object <- frbs.gen(range.data = range.data.input, num.fvalinput, names.varinput, num.fvaloutput = NULL, varout.mf = NULL, names.varoutput = NULL, rule, varinp.mf, type.model, type.defuz = NULL, type.tnorm, type.snorm, func.tsk = NULL, colnames.var, type.implication.func, name) ## Plot the shape of membership functions. plotMF(object) ## Predicting using new data. res <- predict(object, newdata) #################################################### ## 3. The following example shows how to convert ## the frbs model into frbsPMML #################################################### ## In this example, we are using the last object of FRBS. ## Display frbsPMML in R objPMML <- frbsPMML(object) ## Write into a file with .frbsPMML extention ## Not run: write.frbsPMML(objPMML, fileName="obj_frbsPMML") ## Read the frbsPMML file into an R object of FRBS obj <- read.frbsPMML("obj_frbsPMML.frbsPMML") ## End(Not run)
This is one of the central functions of the package. This function is used to generate/learn the model from numerical data using fuzzy rule-based systems.
frbs.learn(data.train, range.data = NULL, method.type = c("WM"), control = list())
frbs.learn(data.train, range.data = NULL, method.type = c("WM"), control = list())
data.train |
a data frame or matrix ( |
range.data |
a matrix ( |
method.type |
this parameter determines the learning algorithm to be used. The following methods are implemented:
|
control |
a list containing all arguments, depending on the learning algorithm to use. The following list are parameters required for each methods, whereas their descriptions will be explained later on.
Description of the
|
This function makes accessible all learning methods that are implemented
in this package. All of the methods use this function as interface for the learning
stage, so users do not need to call other functions in the learning phase.
In order to obtain good results, users need to adjust some parameters such as the
number of labels, the type of the shape of the membership function, the maximal number of iterations,
the step size of the gradient descent, or other method-dependent parameters which are collected in the control
parameter. After creating the model using this function, it can be used to predict new data with predict
.
The frbs-object
.
predict
for the prediction phase, and
the following main functions of each of the methods for theoretical background and references: WM
, SBC
,
HyFIS
, ANFIS
, FIR.DM
, DENFIS
,
FS.HGD
, FRBCS.W
, FRBCS.CHI
, GFS.FR.MOGUL
,
GFS.Thrift
, GFS.GCCL
, FH.GBML
, GFS.LT.RS
, and SLAVE
.
################################## ## I. Regression Problem ## Suppose data have two input variables and one output variable. ## We separate them into training, fitting, and testing data. ## data.train, data.fit, data.test, and range.data are inputs ## for all regression methods. ################################### ## Take into account that the simulation might take a long time ## depending on the hardware you are using. The chosen parameters ## may not be optimal. ## Data must be in data.frame or matrix form and the last column ## is the output variable/attribute. ## The training data must be expressed in numbers (numerical data). data.train <- matrix(c(5.2, -8.1, 4.8, 8.8, -16.1, 4.1, 10.6, -7.8, 5.5, 10.4, -29.0, 5.0, 1.8, -19.2, 3.4, 12.7, -18.9, 3.4, 15.6, -10.6, 4.9, 1.9, -25.0, 3.7, 2.2, -3.1, 3.9, 4.8, -7.8, 4.5, 7.9, -13.9, 4.8, 5.2, -4.5, 4.9, 0.9, -11.6, 3.0, 11.8, -2.1, 4.6, 7.9, -2.0, 4.8, 11.5, -9.0, 5.5, 10.6, -11.2, 4.5, 11.1, -6.1, 4.7, 12.8, -1.0, 6.6, 11.3, -3.6, 5.1, 1.0, -8.2, 3.9, 14.5, -0.5, 5.7, 11.9, -2.0, 5.1, 8.1, -1.6, 5.2, 15.5, -0.7, 4.9, 12.4, -0.8, 5.2, 11.1, -16.8, 5.1, 5.1, -5.1, 4.6, 4.8, -9.5, 3.9, 13.2, -0.7, 6.0, 9.9, -3.3, 4.9, 12.5, -13.6, 4.1, 8.9, -10.0, 4.9, 10.8, -13.5, 5.1), ncol = 3, byrow = TRUE) colnames(data.train) <- c("inp.1", "inp.2", "out.1") data.fit <- data.train[, -ncol(data.train)] data.test <- matrix(c(10.5, -0.9, 5.8, -2.8, 8.5, -0.6, 13.8, -11.9, 9.8, -1.2, 11.0, -14.3, 4.2, -17.0, 6.9, -3.3, 13.2, -1.9), ncol = 2, byrow = TRUE) range.data <- matrix(apply(data.train, 2, range), nrow = 2) ############################################################# ## I.1 Example: Constructing an FRBS model using Wang & Mendel ############################################################# method.type <- "WM" ## collect control parameters into a list ## num.labels = 3 means we define 3 as the number of linguistic terms control.WM <- list(num.labels = 3, type.mf = "GAUSSIAN", type.tnorm = "MIN", type.defuz = "WAM", type.implication.func = "ZADEH", name = "Sim-0") ## generate the model and save it as object.WM object.WM <- frbs.learn(data.train, range.data, method.type, control.WM) ############################################################# ## I.2 Example: Constructing an FRBS model using SBC ############################################################# ## Not run: method.type <- "SBC" control.SBC <- list(r.a = 0.5, eps.high = 0.5, eps.low = 0.15, name = "Sim-0") object.SBC <- frbs.learn(data.train, range.data, method.type, control.SBC) ## End(Not run) ############################################################# ## I.3 Example: Constructing an FRBS model using HYFIS ############################################################# ## Not run: method.type <- "HYFIS" control.HYFIS <- list(num.labels = 5, max.iter = 50, step.size = 0.01, type.tnorm = "MIN", type.defuz = "COG", type.implication.func = "ZADEH", name = "Sim-0") object.HYFIS <- frbs.learn(data.train, range.data, method.type, control.HYFIS) ## End(Not run) ############################################################# ## I.4 Example: Constructing an FRBS model using ANFIS ############################################################# ## Not run: method.type <- "ANFIS" control.ANFIS <- list(num.labels = 5, max.iter = 10, step.size = 0.01, type.tnorm = "MIN", type.implication.func = "ZADEH", name = "Sim-0") object.ANFIS <- frbs.learn(data.train, range.data, method.type, control.ANFIS) ## End(Not run) ############################################################# ## I.5 Example: Constructing an FRBS model using DENFIS ############################################################# ## Not run: control.DENFIS <- list(Dthr = 0.1, max.iter = 10, step.size = 0.001, d = 2, name = "Sim-0") method.type <- "DENFIS" object.DENFIS <- frbs.learn(data.train, range.data, method.type, control.DENFIS) ## End(Not run) ############################################################# ## I.6 Example: Constructing an FRBS model using FIR.DM ############################################################# ## Not run: method.type <- "FIR.DM" control.DM <- list(num.labels = 5, max.iter = 10, step.size = 0.01, type.tnorm = "MIN", type.implication.func = "ZADEH", name = "Sim-0") object.DM <- frbs.learn(data.train, range.data, method.type, control.DM) ## End(Not run) ############################################################# ## I.7 Example: Constructing an FRBS model using FS.HGD ############################################################# ## Not run: method.type <- "FS.HGD" control.HGD <- list(num.labels = 5, max.iter = 10, step.size = 0.01, alpha.heuristic = 1, type.tnorm = "MIN", type.implication.func = "ZADEH", name = "Sim-0") object.HGD <- frbs.learn(data.train, range.data, method.type, control.HGD) ## End(Not run) ############################################################# ## I.8 Example: Constructing an FRBS model using GFS.FR.MOGUL ############################################################# ## Not run: method.type <- "GFS.FR.MOGUL" control.GFS.FR.MOGUL <- list(persen_cross = 0.6, max.iter = 5, max.gen = 2, max.tune = 2, persen_mutant = 0.3, epsilon = 0.8, name="sim-0") object.GFS.FR.MOGUL <- frbs.learn(data.train, range.data, method.type, control.GFS.FR.MOGUL) ## End(Not run) ############################################################# ## I.9 Example: Constructing an FRBS model using Thrift's method (GFS.THRIFT) ############################################################# ## Not run: method.type <- "GFS.THRIFT" control.Thrift <- list(popu.size = 6, num.labels = 3, persen_cross = 1, max.gen = 5, persen_mutant = 1, type.tnorm = "MIN", type.defuz = "COG", type.implication.func = "ZADEH", name="sim-0") object.Thrift <- frbs.learn(data.train, range.data, method.type, control.Thrift) ## End(Not run) ############################################################## ## I.10 Example: Constructing an FRBS model using ## genetic for lateral tuning and rule selection (GFS.LT.RS) ############################################################# ## Set the method and its parameters ## Not run: method.type <- "GFS.LT.RS" control.lt.rs <- list(popu.size = 5, num.labels = 5, persen_mutant = 0.3, max.gen = 10, mode.tuning = "LOCAL", type.tnorm = "MIN", type.implication.func = "ZADEH", type.defuz = "WAM", rule.selection = TRUE, name="sim-0") ## Generate fuzzy model object.lt.rs <- frbs.learn(data.train, range.data, method.type, control.lt.rs) ## End(Not run) ############################################################# ## II. Classification Problems ############################################################# ## The iris dataset is shuffled and divided into training and ## testing data. Bad results in the predicted values may result ## from casual imbalanced classes in the training data. ## Take into account that the simulation may take a long time ## depending on the hardware you use. ## One may get better results with other parameters. ## Data are in data.frame or matrix form and the last column is ## the output variable/attribute ## The data must be expressed in numbers (numerical data). data(iris) irisShuffled <- iris[sample(nrow(iris)),] irisShuffled[,5] <- unclass(irisShuffled[,5]) tra.iris <- irisShuffled[1:105,] tst.iris <- irisShuffled[106:nrow(irisShuffled),1:4] real.iris <- matrix(irisShuffled[106:nrow(irisShuffled),5], ncol = 1) ## Please take into account that the interval needed is the range of input data only. range.data.input <- matrix(apply(iris[, -ncol(iris)], 2, range), nrow = 2) ######################################################### ## II.1 Example: Constructing an FRBS model using ## FRBCS with weighted factor based on Ishibuchi's method ############################################################### ## generate the model ## Not run: method.type <- "FRBCS.W" control <- list(num.labels = 3, type.mf = "TRIANGLE", type.tnorm = "MIN", type.implication.func = "ZADEH", name = "sim-0") object <- frbs.learn(tra.iris, range.data.input, method.type, control) ## conduct the prediction process res.test <- predict(object, tst.iris) ## End(Not run) ######################################################### ## II.2 Example: Constructing an FRBS model using ## FRBCS based on Chi's method ############################################################### ## generate the model ## Not run: method.type <- "FRBCS.CHI" control <- list(num.labels = 7, type.mf = "TRIANGLE", type.tnorm = "MIN", type.implication.func = "ZADEH", name = "sim-0") object <- frbs.learn(tra.iris, range.data.input, method.type, control) ## conduct the prediction process res.test <- predict(object, tst.iris) ## End(Not run) ######################################################### ## II.3 The example: Constructing an FRBS model using GFS.GCCL ############################################################### ## Not run: method.type <- "GFS.GCCL" control <- list(popu.size = 5, num.class = 3, num.labels = 5, persen_cross = 0.9, max.gen = 2, persen_mutant = 0.3, name="sim-0") ## Training process ## The main result of the training is a rule database which is used later for prediction. object <- frbs.learn(tra.iris, range.data.input, method.type, control) ## Prediction process res.test <- predict(object, tst.iris) ## End(Not run) ######################################################### ## II.4 Example: Constructing an FRBS model using FH.GBML ############################################################### ## Not run: method.type <- "FH.GBML" control <- list(popu.size = 5, max.num.rule = 5, num.class = 3, persen_cross = 0.9, max.gen = 2, persen_mutant = 0.3, p.dcare = 0.5, p.gccl = 1, name="sim-0") ## Training process ## The main result of the training is a rule database which is used later for prediction. object <- frbs.learn(tra.iris, range.data.input, method.type, control) ## Prediction process res.test <- predict(object, tst.iris) ## End(Not run) ######################################################### ## II.5 The example: Constructing an FRBS model using SLAVE ############################################################### ## Not run: method.type <- "SLAVE" control <- list(num.class = 3, num.labels = 5, persen_cross = 0.9, max.iter = 5, max.gen = 3, persen_mutant = 0.3, k.lower = 0.25, k.upper = 0.75, epsilon = 0.1, name="sim-0") ## Training process ## The main result of the training is a rule database which is used later for prediction. object <- frbs.learn(tra.iris, range.data.input, method.type, control) ## Prediction process res.test <- predict(object, tst.iris) ## End(Not run)
################################## ## I. Regression Problem ## Suppose data have two input variables and one output variable. ## We separate them into training, fitting, and testing data. ## data.train, data.fit, data.test, and range.data are inputs ## for all regression methods. ################################### ## Take into account that the simulation might take a long time ## depending on the hardware you are using. The chosen parameters ## may not be optimal. ## Data must be in data.frame or matrix form and the last column ## is the output variable/attribute. ## The training data must be expressed in numbers (numerical data). data.train <- matrix(c(5.2, -8.1, 4.8, 8.8, -16.1, 4.1, 10.6, -7.8, 5.5, 10.4, -29.0, 5.0, 1.8, -19.2, 3.4, 12.7, -18.9, 3.4, 15.6, -10.6, 4.9, 1.9, -25.0, 3.7, 2.2, -3.1, 3.9, 4.8, -7.8, 4.5, 7.9, -13.9, 4.8, 5.2, -4.5, 4.9, 0.9, -11.6, 3.0, 11.8, -2.1, 4.6, 7.9, -2.0, 4.8, 11.5, -9.0, 5.5, 10.6, -11.2, 4.5, 11.1, -6.1, 4.7, 12.8, -1.0, 6.6, 11.3, -3.6, 5.1, 1.0, -8.2, 3.9, 14.5, -0.5, 5.7, 11.9, -2.0, 5.1, 8.1, -1.6, 5.2, 15.5, -0.7, 4.9, 12.4, -0.8, 5.2, 11.1, -16.8, 5.1, 5.1, -5.1, 4.6, 4.8, -9.5, 3.9, 13.2, -0.7, 6.0, 9.9, -3.3, 4.9, 12.5, -13.6, 4.1, 8.9, -10.0, 4.9, 10.8, -13.5, 5.1), ncol = 3, byrow = TRUE) colnames(data.train) <- c("inp.1", "inp.2", "out.1") data.fit <- data.train[, -ncol(data.train)] data.test <- matrix(c(10.5, -0.9, 5.8, -2.8, 8.5, -0.6, 13.8, -11.9, 9.8, -1.2, 11.0, -14.3, 4.2, -17.0, 6.9, -3.3, 13.2, -1.9), ncol = 2, byrow = TRUE) range.data <- matrix(apply(data.train, 2, range), nrow = 2) ############################################################# ## I.1 Example: Constructing an FRBS model using Wang & Mendel ############################################################# method.type <- "WM" ## collect control parameters into a list ## num.labels = 3 means we define 3 as the number of linguistic terms control.WM <- list(num.labels = 3, type.mf = "GAUSSIAN", type.tnorm = "MIN", type.defuz = "WAM", type.implication.func = "ZADEH", name = "Sim-0") ## generate the model and save it as object.WM object.WM <- frbs.learn(data.train, range.data, method.type, control.WM) ############################################################# ## I.2 Example: Constructing an FRBS model using SBC ############################################################# ## Not run: method.type <- "SBC" control.SBC <- list(r.a = 0.5, eps.high = 0.5, eps.low = 0.15, name = "Sim-0") object.SBC <- frbs.learn(data.train, range.data, method.type, control.SBC) ## End(Not run) ############################################################# ## I.3 Example: Constructing an FRBS model using HYFIS ############################################################# ## Not run: method.type <- "HYFIS" control.HYFIS <- list(num.labels = 5, max.iter = 50, step.size = 0.01, type.tnorm = "MIN", type.defuz = "COG", type.implication.func = "ZADEH", name = "Sim-0") object.HYFIS <- frbs.learn(data.train, range.data, method.type, control.HYFIS) ## End(Not run) ############################################################# ## I.4 Example: Constructing an FRBS model using ANFIS ############################################################# ## Not run: method.type <- "ANFIS" control.ANFIS <- list(num.labels = 5, max.iter = 10, step.size = 0.01, type.tnorm = "MIN", type.implication.func = "ZADEH", name = "Sim-0") object.ANFIS <- frbs.learn(data.train, range.data, method.type, control.ANFIS) ## End(Not run) ############################################################# ## I.5 Example: Constructing an FRBS model using DENFIS ############################################################# ## Not run: control.DENFIS <- list(Dthr = 0.1, max.iter = 10, step.size = 0.001, d = 2, name = "Sim-0") method.type <- "DENFIS" object.DENFIS <- frbs.learn(data.train, range.data, method.type, control.DENFIS) ## End(Not run) ############################################################# ## I.6 Example: Constructing an FRBS model using FIR.DM ############################################################# ## Not run: method.type <- "FIR.DM" control.DM <- list(num.labels = 5, max.iter = 10, step.size = 0.01, type.tnorm = "MIN", type.implication.func = "ZADEH", name = "Sim-0") object.DM <- frbs.learn(data.train, range.data, method.type, control.DM) ## End(Not run) ############################################################# ## I.7 Example: Constructing an FRBS model using FS.HGD ############################################################# ## Not run: method.type <- "FS.HGD" control.HGD <- list(num.labels = 5, max.iter = 10, step.size = 0.01, alpha.heuristic = 1, type.tnorm = "MIN", type.implication.func = "ZADEH", name = "Sim-0") object.HGD <- frbs.learn(data.train, range.data, method.type, control.HGD) ## End(Not run) ############################################################# ## I.8 Example: Constructing an FRBS model using GFS.FR.MOGUL ############################################################# ## Not run: method.type <- "GFS.FR.MOGUL" control.GFS.FR.MOGUL <- list(persen_cross = 0.6, max.iter = 5, max.gen = 2, max.tune = 2, persen_mutant = 0.3, epsilon = 0.8, name="sim-0") object.GFS.FR.MOGUL <- frbs.learn(data.train, range.data, method.type, control.GFS.FR.MOGUL) ## End(Not run) ############################################################# ## I.9 Example: Constructing an FRBS model using Thrift's method (GFS.THRIFT) ############################################################# ## Not run: method.type <- "GFS.THRIFT" control.Thrift <- list(popu.size = 6, num.labels = 3, persen_cross = 1, max.gen = 5, persen_mutant = 1, type.tnorm = "MIN", type.defuz = "COG", type.implication.func = "ZADEH", name="sim-0") object.Thrift <- frbs.learn(data.train, range.data, method.type, control.Thrift) ## End(Not run) ############################################################## ## I.10 Example: Constructing an FRBS model using ## genetic for lateral tuning and rule selection (GFS.LT.RS) ############################################################# ## Set the method and its parameters ## Not run: method.type <- "GFS.LT.RS" control.lt.rs <- list(popu.size = 5, num.labels = 5, persen_mutant = 0.3, max.gen = 10, mode.tuning = "LOCAL", type.tnorm = "MIN", type.implication.func = "ZADEH", type.defuz = "WAM", rule.selection = TRUE, name="sim-0") ## Generate fuzzy model object.lt.rs <- frbs.learn(data.train, range.data, method.type, control.lt.rs) ## End(Not run) ############################################################# ## II. Classification Problems ############################################################# ## The iris dataset is shuffled and divided into training and ## testing data. Bad results in the predicted values may result ## from casual imbalanced classes in the training data. ## Take into account that the simulation may take a long time ## depending on the hardware you use. ## One may get better results with other parameters. ## Data are in data.frame or matrix form and the last column is ## the output variable/attribute ## The data must be expressed in numbers (numerical data). data(iris) irisShuffled <- iris[sample(nrow(iris)),] irisShuffled[,5] <- unclass(irisShuffled[,5]) tra.iris <- irisShuffled[1:105,] tst.iris <- irisShuffled[106:nrow(irisShuffled),1:4] real.iris <- matrix(irisShuffled[106:nrow(irisShuffled),5], ncol = 1) ## Please take into account that the interval needed is the range of input data only. range.data.input <- matrix(apply(iris[, -ncol(iris)], 2, range), nrow = 2) ######################################################### ## II.1 Example: Constructing an FRBS model using ## FRBCS with weighted factor based on Ishibuchi's method ############################################################### ## generate the model ## Not run: method.type <- "FRBCS.W" control <- list(num.labels = 3, type.mf = "TRIANGLE", type.tnorm = "MIN", type.implication.func = "ZADEH", name = "sim-0") object <- frbs.learn(tra.iris, range.data.input, method.type, control) ## conduct the prediction process res.test <- predict(object, tst.iris) ## End(Not run) ######################################################### ## II.2 Example: Constructing an FRBS model using ## FRBCS based on Chi's method ############################################################### ## generate the model ## Not run: method.type <- "FRBCS.CHI" control <- list(num.labels = 7, type.mf = "TRIANGLE", type.tnorm = "MIN", type.implication.func = "ZADEH", name = "sim-0") object <- frbs.learn(tra.iris, range.data.input, method.type, control) ## conduct the prediction process res.test <- predict(object, tst.iris) ## End(Not run) ######################################################### ## II.3 The example: Constructing an FRBS model using GFS.GCCL ############################################################### ## Not run: method.type <- "GFS.GCCL" control <- list(popu.size = 5, num.class = 3, num.labels = 5, persen_cross = 0.9, max.gen = 2, persen_mutant = 0.3, name="sim-0") ## Training process ## The main result of the training is a rule database which is used later for prediction. object <- frbs.learn(tra.iris, range.data.input, method.type, control) ## Prediction process res.test <- predict(object, tst.iris) ## End(Not run) ######################################################### ## II.4 Example: Constructing an FRBS model using FH.GBML ############################################################### ## Not run: method.type <- "FH.GBML" control <- list(popu.size = 5, max.num.rule = 5, num.class = 3, persen_cross = 0.9, max.gen = 2, persen_mutant = 0.3, p.dcare = 0.5, p.gccl = 1, name="sim-0") ## Training process ## The main result of the training is a rule database which is used later for prediction. object <- frbs.learn(tra.iris, range.data.input, method.type, control) ## Prediction process res.test <- predict(object, tst.iris) ## End(Not run) ######################################################### ## II.5 The example: Constructing an FRBS model using SLAVE ############################################################### ## Not run: method.type <- "SLAVE" control <- list(num.class = 3, num.labels = 5, persen_cross = 0.9, max.iter = 5, max.gen = 3, persen_mutant = 0.3, k.lower = 0.25, k.upper = 0.75, epsilon = 0.1, name="sim-0") ## Training process ## The main result of the training is a rule database which is used later for prediction. object <- frbs.learn(tra.iris, range.data.input, method.type, control) ## Prediction process res.test <- predict(object, tst.iris) ## End(Not run)
The package includes embedded versions of the Mackey-Glass chaotic time series and the Gas Furnance dataset.
Mackey-Glass chaotic time series
The Mackey-Glass chaotic time series is defined by the following delayed differential equation:
For this dataset, we generated 1000 samples, with input parameters as follows:
The dataset is embedded in the following way:
input variables: ,
,
,
output variable:
Gas Furnance dataset
The Gas Furnance dataset is taken from Box and Jenkins. It consists of 292 consecutive
values of methane at time , and the CO2 produced in a furnance at time
as input
variables, with the produced CO2 at time
as an output variable. So, each training data
point consists of
, where
is methane and
is CO2.
G. E. P. Box and G. M. Jenkins, "Time series analysis, forecasting and control", San Fransisco, CA: Holden Day (1970).
M. Mackey and L. Glass, "Oscillation and chaos in physiological control systems", Science, vol. 197, pp. 287 - 289 (1977).
This function creates objects of type frbs
. Currently, its
implementation is very basic and does no argument checking, as
it is only used internally.
frbsObjectFactory(mod)
frbsObjectFactory(mod)
mod |
a list containing all the attributes for the object |
The members of the frbs
object depend on the used learning method. The following list describes all of the members that can be present.
num.labels
the number of linguistic terms for the variables
varout.mf
a matrix to generate the shapes of the membership functions for the output variable. The first row represents the shape of the membership functions, the other rows contain the parameters that have been generated. Whether the values of parameters within the matrix are normalized to lie between 0 and 1 or not depends on the selected method.
rule
the fuzzy IF-THEN rules; In the GFS.FR.MOGUL
case, a rule refers to the parameter values of the membership function
which represents the rule.
rule.data.num
the fuzzy IF-THEN rules in integer format.
varinp.mf
a matrix to generate the shapes of the membership functions for the input variables.
The first row represents the shape of the membership functions,
the other rows contain the non NA
values representing the parameters related with their type of membership function.
For example, TRAPEZOID
, TRIANGLE
, and GAUSSIAN
have four, three, and two values as their parameters, respectively.
Whether the values of parameters within the matrix are normalized to lie between 0 and 1 or not depends on the selected method.
type.model
the type of model. Here, MAMDANI
refers to the Mamdani model, and TSK
refers to the Takagi Sugeno Kang model on the consequence part.
func.tsk
a matrix of the Takagi Sugeno Kang model consequent part of the fuzzy IF-THEN rules.
class
a matrix representing classes of FRBCS
model
num.labels
a number of linguistic terms on each variables/attributes.
type.defuz
the type of the defuzzification method.
type.tnorm
the type of the t-norm method.
type.snorm
the type of the s-norm method.
type.mf
the type of shapes of membership functions.
type.implication.func
the type of the implication function.
method.type
the type of the selected method.
name
the name given to the model.
range.data.ori
range of the original data (before normalization).
cls
cluster centers.
Dthr
the boundary parameter of the DENFIS
method.
d
the multiplier parameters of the DENFIS
method.
r.a
the neighborhood factor of SBC
.
degree.rule
certainty degree of rules.
rule.data.num
a matrix representing the rules in integer form.
grade.cert
grade of certainty for classification problems.
alpha.heuristic
a parameter for the heuristic of the FS.HGD
method.
var.mf.tune
a matrix of parameters of membership function for lateral tuning.
mode.tuning
a type of lateral tuning.
rule.selection
a boolean of rule selection.
colnames.var
the names of variables.
an object of type frbs
It is the main function used for generating the frbsPMML format. In this package, we provide interfaces for writing and reading frbsPMML to/from a text file.
See write.frbsPMML
and read.frbsPMML
.
frbsPMML(model, model.name = "frbs_model", app.name = "frbs", description = NULL, copyright = NULL, algorithm.name = model$method.type, ...)
frbsPMML(model, model.name = "frbs_model", app.name = "frbs", description = NULL, copyright = NULL, algorithm.name = model$method.type, ...)
model |
an frbs model. |
model.name |
a string representing the model name. |
app.name |
a string representing an application name. |
description |
a string representing the simulation description. |
copyright |
a copyright of simulation. |
algorithm.name |
a string representing the algorithm name. |
... |
other parameters |
frbsPMML is a universal framework for representing FRBS models, which is a format adopted from the Predictive Model Markup Language (PMML). PMML is a format constructed by an XML-based language to provide a standard for describing models produced by data mining and machine learning algorithms. A main contribution of PMML is to provide interoperable schemata of predictive models. Using PMML, we can easily perform these tasks as our models are documented in an XML-based language. Human experts can also update and modify the model on the files directly.
Since PMML is an XML-based language, the specification is defined by an XML Schema as recommended by the World Wide Web Consortium (W3C). The PMML format is specified by the main tag PMML that contains some components. In the following, we describe the main components:
Header: It contains general information about the PMML document, such as copyright information for the model, its description, application, and timestamp of generation.
DataDictionary: It contains information related to fields or variables, such as number, names, types, and value ranges of variables.
MODEL-ELEMENT: It is a main part of the PMML document that consists of models supported by PMML. In each model, there are several components embedded in the element, such as MiningSchema and Output. MiningSchema specifies outlier treatment, a missing value replacement policy, and missing value treatment, whereas Output shows a description of the output variable. For example, in a clustering model, we define a schema representing the cluster centers that are included in the ClusteringModel element.
Besides these components, there are some optional elements, such as MiningBuildTask, TransformationDictionary, and Extension. More detailed information about PMML can be found in (Guazzelli et al., 2009).
Three models, which can be used for handling regression and classification tasks, are specified by the proposed representations: Mamdani, Takagi Sugeno Kang, and fuzzy rule-based classification systems. There are the following benefits offered by frbsPMML, as follows:
Interoperability: It is a standard format for representing many models without depending on any programming languages (e.g., Java, Python, and C++) and platforms (e.g., Windows, Linux, and Mac).
Tranparency: Since it is formed based on XML Schema containing formal definitions of the available elements, we can understand FRBS models as written in frbsPMML.
Interpretability: frbsPMML expresses rulebase, database, and inference schema in simple ways. For example,
rulebase is constructed recursively, so that besides it meets to the mathematical logic (predicate), we can
define different operators (i.e., and
and or
) in one rule.
Flexibility: Since frbsPMML is based XML, human experts can easily modify and improve a model in the text file directly.
Reproducibility: Sicen frbsPMML is a universal representation, it allows us to store, share, execute, and reproduce an FRBS model.
FRBS model in frbsPMML format
A. Guazzelli, M. Zeller, W.C. Lin, and G. Williams., "pmml: An open standard for sharing models", The R Journal, Vol. 1, No. 1, pp. 60-65 (2009).
Data Mining Group, http://www.dmg.org/.
## This example shows how to construct a frbsPMML file of the frbs model ## Even though we are using MAMDANI model, other models have the same way ## ## 1. Produce frbs model, for example: we perform Wang & Mendel's technique (WM) ## Input data ## Not run: data(frbsData) data.train <- frbsData$GasFurnance.dt[1 : 204, ] data.fit <- data.train[, 1 : 2] data.tst <- frbsData$GasFurnance.dt[205 : 292, 1 : 2] real.val <- matrix(frbsData$GasFurnance.dt[205 : 292, 3], ncol = 1) range.data <- matrix(c(-2.716, 2.834, 45.6, 60.5, 45.6, 60.5), nrow = 2) ## Set the method and its parameters method.type <- "WM" control <- list(num.labels = 3, type.mf = "GAUSSIAN", type.defuz = "WAM", type.tnorm = "MIN", type.implication.func = "ZADEH", name = "sim-0") ## Generate fuzzy model object <- frbs.learn(data.train, range.data, method.type, control) ## 2. Write frbsPMML file ## by calling frbsPMML(), the frbsPMML format will be displayed in R console frbsPMML(object) ## End(Not run)
## This example shows how to construct a frbsPMML file of the frbs model ## Even though we are using MAMDANI model, other models have the same way ## ## 1. Produce frbs model, for example: we perform Wang & Mendel's technique (WM) ## Input data ## Not run: data(frbsData) data.train <- frbsData$GasFurnance.dt[1 : 204, ] data.fit <- data.train[, 1 : 2] data.tst <- frbsData$GasFurnance.dt[205 : 292, 1 : 2] real.val <- matrix(frbsData$GasFurnance.dt[205 : 292, 3], ncol = 1) range.data <- matrix(c(-2.716, 2.834, 45.6, 60.5, 45.6, 60.5), nrow = 2) ## Set the method and its parameters method.type <- "WM" control <- list(num.labels = 3, type.mf = "GAUSSIAN", type.defuz = "WAM", type.tnorm = "MIN", type.implication.func = "ZADEH", name = "sim-0") ## Generate fuzzy model object <- frbs.learn(data.train, range.data, method.type, control) ## 2. Write frbsPMML file ## by calling frbsPMML(), the frbsPMML format will be displayed in R console frbsPMML(object) ## End(Not run)
This is the internal function that implements the simplified TSK fuzzy rule generation method
using heuristics and gradient descent method (FS.HGD). It is used to solve regression tasks.
Users do not need to call it directly,
but just use frbs.learn
and predict
.
FS.HGD(data.train, num.labels, max.iter = 100, step.size = 0.01, alpha.heuristic = 1, type.tnorm = "MIN", type.snorm = "MAX", type.implication.func = "ZADEH")
FS.HGD(data.train, num.labels, max.iter = 100, step.size = 0.01, alpha.heuristic = 1, type.tnorm = "MIN", type.snorm = "MAX", type.implication.func = "ZADEH")
data.train |
a matrix ( |
num.labels |
a matrix ( |
max.iter |
maximal number of iterations. |
step.size |
step size of the descent method. |
alpha.heuristic |
a positive real number which is the heuristic parameter. |
type.tnorm |
the type of t-norm. For more detail, please have a look at |
type.snorm |
the type of s-norm. For more detail, please have a look at |
type.implication.func |
a value representing type of implication function.
For more detail, please have a look at |
This method was proposed by K. Nozaki, H. Ishibuchi, and H. Tanaka. It uses fuzzy IF-THEN rules with nonfuzzy singletons (i.e. real numbers) in the consequent parts. The techniques of space partition are implemented to generate the antecedent part, while the initial consequent part of each rule is determined by the weighted mean value of the given training data. Then, the gradient descent method updates the value of the consequent part. Futhermore, the heuristic value given by the user affects the value of weight of each data.
H. Ishibuchi, K. Nozaki, H. Tanaka, Y. Hosaka, and M. Matsuda, "Empirical study on learning in fuzzy systems by rice taste analysis", Fuzzy Set and Systems, vol. 64, no. 2, pp. 129 - 144 (1994).
frbs.learn
, predict
, and HGD.update
Fuzzification refers to the process of transforming a crisp set into linguistic terms.
fuzzifier(data, num.varinput, num.labels.input, varinp.mf)
fuzzifier(data, num.varinput, num.labels.input, varinp.mf)
data |
a matrix of data containing numerical elements. |
num.varinput |
number of input variables. |
num.labels.input |
the number of labels of the input variables. |
varinp.mf |
a matrix containing the parameters to form the membership functions. See the Detail section. |
In this function, there are five shapes of membership functions implemented,
namely TRIANGLE
, TRAPEZOID
, GAUSSIAN
, SIGMOID
, and BELL
.
They are represented by a matrix that the dimension is () where
is
a multiplication the number of linguistic terms/labels and the number of input variables.
The rows of the matrix represent:
The first row is the type of membership function, where 1 means
TRIANGLE
,
2 means TRAPEZOID
in left side,
3 means TRAPEZOID
in right side, 4 means TRAPEZOID
in the middle,
5 means GAUSSIAN
,
6 means SIGMOID
, and 7 means BELL
. And, the second up to fifth row indicate
the corner points to construct the functions.
TRIANGLE
has three parameters (), where
is the center point of the
TRIANGLE
,
and and
are the left and right points, respectively.
TRAPEZOID
has four parameters ().
GAUSSIAN
has two parameters ( and
).
SIGMOID
has two parameters ( and
) for representing steepness of the function and distance from the origin, respectively.
BELL
has three parameters ().
For example:
varinp.mf <- matrix(c(2,1,3,2,3,0,30,60,0,40,20,50,80,
30,80,40,70,100,60,100,0,0,100,0,100), nrow=5, byrow=TRUE)
A matrix of the degree of each linguistic terms based on the shape of the membership functions
defuzzifier
, rulebase
, and inference
This is the internal function that implements genetic fuzzy systems for fuzzy rule learning
based on the MOGUL methodology (GFS.FR.MOGUL). It is used to solve regression tasks.
Users do not need to call it directly, but just use frbs.learn
and predict
.
GFS.FR.MOGUL(data.train, persen_cross = 0.6, persen_mutant = 0.3, max.iter = 10, max.gen = 10, max.tune = 10, range.data.ori, epsilon = 0.4)
GFS.FR.MOGUL(data.train, persen_cross = 0.6, persen_mutant = 0.3, max.iter = 10, max.gen = 10, max.tune = 10, range.data.ori, epsilon = 0.4)
data.train |
a matrix ( |
persen_cross |
a real number between 0 and 1 determining the probability of crossover. |
persen_mutant |
a real number between 0 and 1 determining the probability of mutation. |
max.iter |
the maximal number of iterations. |
max.gen |
the maximal number of generations of the genetic algorithm. |
max.tune |
the maximal number of tuning iterations. |
range.data.ori |
a matrix containing the ranges of the original data. |
epsilon |
a real number between 0 and 1 determining the boundary of covering factor. |
This method was proposed by Herrera et al. GFS.FR.MOGUL implements a genetic algorithm determining the structure of the fuzzy IF-THEN rules and the membership function parameters. There are two general types of fuzzy IF-THEN rules, namely the descriptive and the approximative/free semantic approaches. A descriptive approach means that the linguistic labels represent a real-world semantic; the linguistic labels are uniformly defined for all rules. In contrast, in the approximative approach there isn't any associated linguistic label. This method is based on the latter one. We model a fuzzy IF-THEN rule on a chromosome which consists of the parameter values of the membership function. So, every rule has its own membership function values. A population contains many such generated chromosomes, based on the iterative rule learning approach (IRL). IRL means that the chromosomes will be generated one by one, taking into account the fitness value and covering factor, until there are sufficient chromosomes in the population. After having obtained the population, the genetic algorithm is started, using the genetic operators selection, mutation, and crossover.
F. Herrera, M. Lozano, and J.L. Verdegay, "A learning process for fuzzy control rules using genetic algorithms", Fuzzy Sets and Systems, vol. 100, pp. 143 - 158 (1998).
O. Cordon, M.J. del Jesus, F. Herrera, and M. Lozano, "MOGUL: A methodology to obtain genetic fuzzy rule-based systems under the iterative rule learning approach", International Journal of Intelligent Systems, vol. 14, pp. 1123 - 1153 (1999).
GFS.FR.MOGUL.test
, frbs.learn
, and predict
This function is the internal function of the GFS.FR.MOGUL method to compute the predicted values.
GFS.FR.MOGUL.test(object, newdata)
GFS.FR.MOGUL.test(object, newdata)
object |
the |
newdata |
a matrix ( |
A matrix of predicted values.
This is the internal function that implements the Ishibuchi's method based on
genetic cooperative-competitive learning (GFS.GCCL). It is used to handle classification tasks.
Users do not need to call it directly,
but just use frbs.learn
and predict
.
GFS.GCCL(data.train, popu.size = 10, range.data.input, num.labels, persen_cross = 0.6, persen_mutant = 0.3, max.gen = 10, range.data.ori)
GFS.GCCL(data.train, popu.size = 10, range.data.input, num.labels, persen_cross = 0.6, persen_mutant = 0.3, max.gen = 10, range.data.ori)
data.train |
a matrix ( |
popu.size |
the size of the population which is generated in each generation. |
range.data.input |
a matrix containing the ranges of the normalized input data. |
num.labels |
a matrix describing the number of linguistic terms. |
persen_cross |
a real number between 0 and 1 representing the probability of crossover. |
persen_mutant |
a real number between 0 and 1 representing the probability of mutation. |
max.gen |
the maximal number of generations for the genetic algorithm. |
range.data.ori |
a matrix containing the ranges of the input data. |
This method is based on Ishibuchi's method. In this method, a chromosome describes each linguistic IF-THEN rule using integer as its representation of the antecedent part. In the consequent part of the fuzzy rules, the heuristic method is applied to automatically generate the class. The evaluation is calculated for each rule which means that the performance is not based on the entire rule set. The outline of the method is as follows.
Step 1: Generate an initial population of fuzzy IF-THEN rules.
Step 2: Evaluate each fuzzy IF-THEN rule in the current population.
Step 3: Generate new fuzzy IF-THEN rules by genetic operators.
Step 4: Replace a part of the current population with the newly generated rules.
Step 5: Terminate the algorithm if a stopping condition is satisfied, otherwise return to Step 2.
Additionally, to handle high dimensional data, this method uses "don't care" attributes on the antecedent fuzzy set.
H. Ishibuchi, T. Nakashima, and T. Murata, "Performance evaluation of fuzzy classifier systems for multidimensional pattern classification problems", IEEE trans. on Systems, Man, and Cybernetics - Part B: Sybernetics, vol. 29. no. 5, pp. 601 - 618 (1999).
This function is the internal function of the GFS.GCCL and FH.GBML method to compute the predicted values.
GFS.GCCL.eng(object, newdata)
GFS.GCCL.eng(object, newdata)
object |
the |
newdata |
a matrix ( |
A matrix of predicted values.
This is the internal function that implements genetic lateral tuning and rule selection of linguistic fuzzy systems (GFS.LT.RS).
It is used to solve regression tasks.
Users do not need to call it directly, but just use frbs.learn
and predict
.
GFS.LT.RS(data.train, popu.size = 10, range.data, num.labels, persen_mutant, max.gen = 10, mode.tuning = "GLOBAL", type.tnorm = "MIN", type.snorm = "MAX", type.implication.func = "ZADEH", type.defuz = "WAM", rule.selection = FALSE, range.data.ori)
GFS.LT.RS(data.train, popu.size = 10, range.data, num.labels, persen_mutant, max.gen = 10, mode.tuning = "GLOBAL", type.tnorm = "MIN", type.snorm = "MAX", type.implication.func = "ZADEH", type.defuz = "WAM", rule.selection = FALSE, range.data.ori)
data.train |
a matrix ( |
popu.size |
the size of the population which is generated in each generation. |
range.data |
a matrix representing interval of data. |
num.labels |
a matrix representing the number of linguistic terms in each variables. |
persen_mutant |
a real number between 0 and 1 determining the probability of mutation. |
max.gen |
the maximal number of generations of the genetic algorithm. |
mode.tuning |
a type of tuning which are |
type.tnorm |
a type of t-norm. See |
type.snorm |
a type of s-norm. See |
type.implication.func |
a type of implication function. See |
type.defuz |
a type of defuzzification methods. See |
rule.selection |
a boolean value representing whether performs rule selection or not. |
range.data.ori |
a matrix containing the ranges of the original data. |
This method was proposed by R. Alcala et al. GFS.LT.RS implements a evolutionary algorithm for postprocessing in constructing FRBS model. It uses a new rule representation model based on the linguistic 2-tupples representation that allows the lateral displacement of the labels. This function allows two different tuning which are global and local tuning.
Regarding with evolutionary algorithms, the following are main components:
coding scheme and initial gene pool;
chromosome evalution;
crossover operator;
restarting approach;
evolutionary model;
In first time, population is constructed by Wang & Mendel's technique. Mean square error (MSE) is used to calculate chromosome evaluation. This method performs BLX-a in crossover process. Additionally, rule selection method is performed in order to minimize the number of rules.
R. Alcala, J. Alcala-Fdez, and F. Herrera, "A proposal for the genetic lateral tuning of linguistic fuzzy systems and its interaction with rule selection", IEEE Trans. on Fuzzy Systems, Vol. 15, No. 4, pp. 616 - 635 (2007).
GFS.LT.RS.test
, frbs.learn
, and predict
This function is the internal function of the GFS.LT.RS method to compute the predicted values.
GFS.LT.RS.test(object, newdata)
GFS.LT.RS.test(object, newdata)
object |
the |
newdata |
a matrix ( |
A matrix of predicted values.
This is the internal function that implements the Thrift's technique based
on a genetic algorithm. It is used to tackle regression tasks. Users do not need to call it directly,
but just use frbs.learn
and predict
.
GFS.Thrift(data.train, popu.size = 10, num.labels, persen_cross = 0.6, persen_mutant = 0.3, max.gen = 10, range.data.ori, type.defuz = "WAM", type.tnorm = "MIN", type.snorm = "MAX", type.mf = "TRIANGLE", type.implication.func = "ZADEH")
GFS.Thrift(data.train, popu.size = 10, num.labels, persen_cross = 0.6, persen_mutant = 0.3, max.gen = 10, range.data.ori, type.defuz = "WAM", type.tnorm = "MIN", type.snorm = "MAX", type.mf = "TRIANGLE", type.implication.func = "ZADEH")
data.train |
a matrix ( |
popu.size |
the size of the population which is generated in each generation. |
num.labels |
a matrix describing the number of linguistic terms. |
persen_cross |
a real number between 0 and 1 representing the probability of crossover. |
persen_mutant |
a real number between 0 and 1 representing the probability of mutation. |
max.gen |
the maximal number of generations for the genetic algorithm. |
range.data.ori |
a matrix containing the ranges of the original data. |
type.defuz |
the type of the defuzzification method. For more detail, see |
type.tnorm |
the type of t-norm. For more detail, please have a look at |
type.snorm |
the type of s-norm. For more detail, please have a look at |
type.mf |
the type of shape of membership function. See |
type.implication.func |
the type of implication function. See |
This method was developed by Thrift using Mamdani's model as fuzzy IF-THEN rules. In this method, we consider a table as a genotype with alleles that are fuzzy set indicators over the output domain. The phenotype is produced by the behavior produced by the fuzzification, max-* composition, and defuzzification operations. A chromosome (genotype) is formed from the decision table by going rowwise and producing a string of numbers from the code set. Standard crossover and mutation operators can act on these string.
P. Thrift, "Fuzzy logic synthesis with genetic algorithms", In Proceedings of the Fourth International Conference on Genetic Algorithms (ICGA91), San Diego (United States of America), pp. 509 - 513 (1991).
GFS.Thrift.test
, frbs.learn
, and predict
This function is the internal function of the GFS.Thrift method to compute the predicted values.
GFS.Thrift.test(object, newdata)
GFS.Thrift.test(object, newdata)
object |
the |
newdata |
a matrix ( |
A matrix of predicted values.
The role of this function is to update parameters within the simplified TSK fuzzy rule
generation method using heuristics and the gradient descent method (FS.HGD).
This function is called by the main function
of the FS.HGD method, see FS.HGD
.
HGD.update(data.train, miu.rule, func.tsk, varinp.mf, step.size = 0.01, def)
HGD.update(data.train, miu.rule, func.tsk, varinp.mf, step.size = 0.01, def)
data.train |
a matrix ( |
miu.rule |
a matrix with the degrees of rules which is the result of the |
func.tsk |
a matrix of parameters of the function on the consequent part using the Takagi Sugeno Kang model. See |
varinp.mf |
a matrix of parameters of membership functions of the input variables. |
step.size |
a real number between 0 and 1 representing the step size of the gradient descent. |
def |
a matrix which is obtained by the |
This is the internal function that implements the hybrid neural fuzzy inference
system (HyFIS). It is used to solve regression tasks.
Users do not need to call it directly,
but just use frbs.learn
and predict
HyFIS(data.train, num.labels, max.iter = 10, step.size = 0.01, type.tnorm = "MIN", type.snorm = "MAX", type.defuz = "COG", type.implication.func = "ZADEH")
HyFIS(data.train, num.labels, max.iter = 10, step.size = 0.01, type.tnorm = "MIN", type.snorm = "MAX", type.defuz = "COG", type.implication.func = "ZADEH")
data.train |
a matrix ( |
num.labels |
a matrix ( |
max.iter |
the maximal number of iterations. |
step.size |
step size of the gradient descent method. |
type.tnorm |
the type of t-norm. For more detail, please have a look at |
type.snorm |
the type of s-norm. For more detail, please have a look at |
type.defuz |
the type of aggregation function. For more detail, please have a look at |
type.implication.func |
a value representing type of implication function. For more detail, please have a look at |
This method was proposed by J. Kim and N. Kasabov. There are two phases in this method for learning, namely the knowledge acquisition module and the structure and parameter learning. The knowledge acquition module uses the techniques of Wang and Mendel. The learning of structure and parameters is a supervised learning method using gradient descent-based learning algorithms. This function generates a model which consists of a rule database and parameters of the membership functions. The rules of HyFIS use the Mamdani model on the antecedent and consequent parts. Futhermore, HyFIS uses a Gaussian membership function. So, there are two kinds of parameters that are optimized, mean and variance of the Gaussian function.
J. Kim and N. Kasabov, "HyFIS: Adaptive neuro-fuzzy inference systems and their application to nonlinear dynamical systems", Neural Networks, vol. 12, no. 9, pp. 1301 - 1319 (1999).
HyFIS.update
, frbs.learn
, and predict
.
This function is called by HyFIS
to update the parameters within
the HyFIS method.
HyFIS.update(data.train, def, rule, names.varoutput, var.mf, miu.rule, num.labels, MF, step.size = 0.001, degree.rule)
HyFIS.update(data.train, def, rule, names.varoutput, var.mf, miu.rule, num.labels, MF, step.size = 0.001, degree.rule)
data.train |
a matrix ( |
def |
matrix of defuzzification results. See |
rule |
fuzzy IF-THEN rules. See |
names.varoutput |
a list of names of the output variable. |
var.mf |
a matrix of parameters of the membership functions.
Please see |
miu.rule |
a matrix of degree of rules which is a result of the |
num.labels |
a matrix ( |
MF |
a matrix of parameters of the membership functions
which is a result of the |
step.size |
a real number, the step size of the gradient descent. |
degree.rule |
a matrix of degrees of rules. See |
Inference refers to the process of fuzzy reasoning.
inference(MF, rule, names.varinput, type.tnorm, type.snorm)
inference(MF, rule, names.varinput, type.tnorm, type.snorm)
MF |
a matrix of the degrees of membership functions which is a result of the |
rule |
a matrix or list of fuzzy IF-THEN rules. See |
names.varinput |
a list of names of the input variables. |
type.tnorm |
a value which represents the type of t-norm to be used:
|
type.snorm |
a value which represents the type of s-norm to be used:
|
In this function, fuzzy reasoning is conducted based on Mamdani and Takagi Sugeno Kang model. Furthermore, there are some formula for conjunction and disjunction operators.
The Mamdani model:
A fuzzy system with, e.g., two inputs and
(antecedents) and a single output
(consequent)
is described by the following fuzzy IF-THEN rule:
IF x1 is A1 and x2 is A2 THEN y is B
where and
are the fuzzy sets representing the antecent pairs and
is the fuzzy set representing the consequent.
The Takagi Sugeno Kang model:
Suppose we have two inputs and
and output
, then the fuzzy IF-THEN rule is as follows:
IF x1 is A1 and x2 is A2 THEN y is y = f(x1, x2)
where is a crisp function in the consequent part which is usually a polynomial function,
and
and
are the fuzzy sets representing the antecent pairs.
Futhermore, this function has the following capabilities:
It supports unary operators (not) and binary operators (AND
and OR
).
It provides linguistic hedge (extremely
, very
, somewhat
, and slightly
).
there are several methods for the t-norm and s-norm.
a matrix of the degrees of the rules.
defuzzifier
, rulebase
, and fuzzifier
.
This function is to transform from real-valued data into normalized data.
norm.data(dt.ori, range.data, min.scale = 0, max.scale = 1)
norm.data(dt.ori, range.data, min.scale = 0, max.scale = 1)
dt.ori |
a matrix ( |
range.data |
a matrix ( |
min.scale |
the minimum value within normalization. |
max.scale |
the maximum value within normalization. |
the normalized data
This function can be used to plot the shapes of the membership functions.
plotMF(object)
plotMF(object)
object |
an
|
## The following examples contain two different cases which are ## using an frbs-object and the manual way. ## ## 1. Plotting using frbs object. data(iris) irisShuffled <- iris[sample(nrow(iris)),] irisShuffled[,5] <- unclass(irisShuffled[,5]) tra.iris <- irisShuffled[1:105,] tst.iris <- irisShuffled[106:nrow(irisShuffled),1:4] real.iris <- matrix(irisShuffled[106:nrow(irisShuffled),5], ncol = 1) ## Please take into account that the interval needed is the range of input data only. range.data.input <- matrix(c(4.3, 7.9, 2.0, 4.4, 1.0, 6.9, 0.1, 2.5), nrow=2) ## generate the model method.type <- "FRBCS.W" control <- list(num.labels = 7, type.mf = 1) ## Not run: object <- frbs.learn(tra.iris, range.data.input, method.type, control) ## plot the frbs object ## Not run: plotMF(object) ## 2. Plotting using params. ## Define shape and parameters of membership functions of input variables. ## Please see the fuzzifier function of how to contruct the matrix. varinp.mf <- matrix(c(2, 0, 20, 40, NA, 4, 20, 40, 60, 80, 3, 60, 80, 100, NA, 2, 0, 20, 40, NA, 4, 20, 40, 60, 80, 3, 60, 80, 100, NA, 2, 0, 20, 40, NA, 4, 20, 40, 60, 80, 3, 60, 80, 100, NA, 2, 0, 20, 40, NA, 4, 20, 40, 60, 80, 3, 60, 80, 100, NA), nrow = 5, byrow = FALSE) ## Define the shapes and parameters of the membership functions of the output variables. varout.mf <- matrix(c(2, 0, 20, 40, NA, 4, 20, 40, 60, 80, 3, 60, 80, 100, NA), nrow = 5, byrow = FALSE) var.mf <- cbind(varinp.mf, varout.mf) range.data <- matrix(c(0,100, 0, 100, 0, 100, 0, 100, 0, 100), nrow=2) num.labels <- matrix(c(3,3,3,3,3), nrow = 1) names.variables <- c("input1", "input2", "input3", "input4", "output1") ## plot the membership function. ## Not run: plotMF(object = list(var.mf = var.mf, range.data.ori = range.data, num.labels = num.labels, names.variables = names.variables)) ## End(Not run)
## The following examples contain two different cases which are ## using an frbs-object and the manual way. ## ## 1. Plotting using frbs object. data(iris) irisShuffled <- iris[sample(nrow(iris)),] irisShuffled[,5] <- unclass(irisShuffled[,5]) tra.iris <- irisShuffled[1:105,] tst.iris <- irisShuffled[106:nrow(irisShuffled),1:4] real.iris <- matrix(irisShuffled[106:nrow(irisShuffled),5], ncol = 1) ## Please take into account that the interval needed is the range of input data only. range.data.input <- matrix(c(4.3, 7.9, 2.0, 4.4, 1.0, 6.9, 0.1, 2.5), nrow=2) ## generate the model method.type <- "FRBCS.W" control <- list(num.labels = 7, type.mf = 1) ## Not run: object <- frbs.learn(tra.iris, range.data.input, method.type, control) ## plot the frbs object ## Not run: plotMF(object) ## 2. Plotting using params. ## Define shape and parameters of membership functions of input variables. ## Please see the fuzzifier function of how to contruct the matrix. varinp.mf <- matrix(c(2, 0, 20, 40, NA, 4, 20, 40, 60, 80, 3, 60, 80, 100, NA, 2, 0, 20, 40, NA, 4, 20, 40, 60, 80, 3, 60, 80, 100, NA, 2, 0, 20, 40, NA, 4, 20, 40, 60, 80, 3, 60, 80, 100, NA, 2, 0, 20, 40, NA, 4, 20, 40, 60, 80, 3, 60, 80, 100, NA), nrow = 5, byrow = FALSE) ## Define the shapes and parameters of the membership functions of the output variables. varout.mf <- matrix(c(2, 0, 20, 40, NA, 4, 20, 40, 60, 80, 3, 60, 80, 100, NA), nrow = 5, byrow = FALSE) var.mf <- cbind(varinp.mf, varout.mf) range.data <- matrix(c(0,100, 0, 100, 0, 100, 0, 100, 0, 100), nrow=2) num.labels <- matrix(c(3,3,3,3,3), nrow = 1) names.variables <- c("input1", "input2", "input3", "input4", "output1") ## plot the membership function. ## Not run: plotMF(object = list(var.mf = var.mf, range.data.ori = range.data, num.labels = num.labels, names.variables = names.variables)) ## End(Not run)
This is the main function to obtain a final result as predicted values for all methods in this package.
In order to get predicted values, this function is run using an frbs-object
, which is typically generated using frbs.learn
.
## S3 method for class 'frbs' predict(object, newdata, ...)
## S3 method for class 'frbs' predict(object, newdata, ...)
object |
an |
newdata |
a data frame or matrix ( |
... |
the other parameters (not used) |
The predicted values.
frbs.learn
and frbs.gen
for learning and model generation,
and the internal main functions of each method for the theory:
WM
, SBC
, HyFIS
, ANFIS
,
FIR.DM
, DENFIS
, FS.HGD
, FRBCS.W
,
GFS.FR.MOGUL
, GFS.Thrift
, GFS.GCCL
, FRBCS.CHI
,
FH.GBML
, GFS.LT.RS
, and SLAVE
.
################################## ## I. Regression Problem ################################### ## In this example, we just show how to predict using Wang and Mendel's technique but ## users can do it in the same way for other methods. data.train <- matrix(c(5.2, -8.1, 4.8, 8.8, -16.1, 4.1, 10.6, -7.8, 5.5, 10.4, -29.0, 5.0, 1.8, -19.2, 3.4, 12.7, -18.9, 3.4, 15.6, -10.6, 4.9, 1.9, -25.0, 3.7, 2.2, -3.1, 3.9, 4.8, -7.8, 4.5, 7.9, -13.9, 4.8, 5.2, -4.5, 4.9, 0.9, -11.6, 3.0, 11.8, -2.1, 4.6, 7.9, -2.0, 4.8, 11.5, -9.0, 5.5, 10.6, -11.2, 4.5, 11.1, -6.1, 4.7, 12.8, -1.0, 6.6, 11.3, -3.6, 5.1, 1.0, -8.2, 3.9, 14.5, -0.5, 5.7, 11.9, -2.0, 5.1, 8.1, -1.6, 5.2, 15.5, -0.7, 4.9, 12.4, -0.8, 5.2, 11.1, -16.8, 5.1, 5.1, -5.1, 4.6, 4.8, -9.5, 3.9, 13.2, -0.7, 6.0, 9.9, -3.3, 4.9, 12.5, -13.6, 4.1, 8.9, -10.0, 4.9, 10.8, -13.5, 5.1), ncol = 3, byrow = TRUE) data.fit <- matrix(c(10.5, -0.9, 5.2, 5.8, -2.8, 5.6, 8.5, -0.2, 5.3, 13.8, -11.9, 3.7, 9.8, -1.2, 4.8, 11.0, -14.3, 4.4, 4.2, -17.0, 5.1, 6.9, -3.3, 5.1, 13.2, -1.9, 4.6), ncol = 3, byrow = TRUE) newdata <- matrix(c(10.5, -0.9, 5.8, -2.8, 8.5, -0.2, 13.8, -11.9, 9.8, -1.2, 11.0, -14.3, 4.2, -17.0, 6.9, -3.3, 13.2, -1.9), ncol = 2, byrow = TRUE) range.data<-matrix(c(0.9, 15.6, -29, -0.2, 3, 6.6), ncol=3, byrow = FALSE) ############################################################# ## I.1 Example: Implementation of Wang & Mendel ############################################################# method.type <- "WM" ## collect control parameters into a list ## num.labels = 3 means we define 3 as the number of linguistic terms control.WM <- list(num.labels = 3, type.mf = "GAUSSIAN", type.tnorm = "MIN", type.snorm = "MAX", type.defuz = "WAM", type.implication.func = "ZADEH", name = "Sim-0") ## generate the model and save it as object.WM object.WM <- frbs.learn(data.train, range.data, method.type, control.WM) ## the prediction process ## The following code can be used for all methods res <- predict(object.WM, newdata)
################################## ## I. Regression Problem ################################### ## In this example, we just show how to predict using Wang and Mendel's technique but ## users can do it in the same way for other methods. data.train <- matrix(c(5.2, -8.1, 4.8, 8.8, -16.1, 4.1, 10.6, -7.8, 5.5, 10.4, -29.0, 5.0, 1.8, -19.2, 3.4, 12.7, -18.9, 3.4, 15.6, -10.6, 4.9, 1.9, -25.0, 3.7, 2.2, -3.1, 3.9, 4.8, -7.8, 4.5, 7.9, -13.9, 4.8, 5.2, -4.5, 4.9, 0.9, -11.6, 3.0, 11.8, -2.1, 4.6, 7.9, -2.0, 4.8, 11.5, -9.0, 5.5, 10.6, -11.2, 4.5, 11.1, -6.1, 4.7, 12.8, -1.0, 6.6, 11.3, -3.6, 5.1, 1.0, -8.2, 3.9, 14.5, -0.5, 5.7, 11.9, -2.0, 5.1, 8.1, -1.6, 5.2, 15.5, -0.7, 4.9, 12.4, -0.8, 5.2, 11.1, -16.8, 5.1, 5.1, -5.1, 4.6, 4.8, -9.5, 3.9, 13.2, -0.7, 6.0, 9.9, -3.3, 4.9, 12.5, -13.6, 4.1, 8.9, -10.0, 4.9, 10.8, -13.5, 5.1), ncol = 3, byrow = TRUE) data.fit <- matrix(c(10.5, -0.9, 5.2, 5.8, -2.8, 5.6, 8.5, -0.2, 5.3, 13.8, -11.9, 3.7, 9.8, -1.2, 4.8, 11.0, -14.3, 4.4, 4.2, -17.0, 5.1, 6.9, -3.3, 5.1, 13.2, -1.9, 4.6), ncol = 3, byrow = TRUE) newdata <- matrix(c(10.5, -0.9, 5.8, -2.8, 8.5, -0.2, 13.8, -11.9, 9.8, -1.2, 11.0, -14.3, 4.2, -17.0, 6.9, -3.3, 13.2, -1.9), ncol = 2, byrow = TRUE) range.data<-matrix(c(0.9, 15.6, -29, -0.2, 3, 6.6), ncol=3, byrow = FALSE) ############################################################# ## I.1 Example: Implementation of Wang & Mendel ############################################################# method.type <- "WM" ## collect control parameters into a list ## num.labels = 3 means we define 3 as the number of linguistic terms control.WM <- list(num.labels = 3, type.mf = "GAUSSIAN", type.tnorm = "MIN", type.snorm = "MAX", type.defuz = "WAM", type.implication.func = "ZADEH", name = "Sim-0") ## generate the model and save it as object.WM object.WM <- frbs.learn(data.train, range.data, method.type, control.WM) ## the prediction process ## The following code can be used for all methods res <- predict(object.WM, newdata)
It is used to read the frbsPMML format into an frbs model in R. Detailed information about frbsPMML can be seen in frbsPMML
.
read.frbsPMML(fileName)
read.frbsPMML(fileName)
fileName |
a file name with extension |
an object representing the frbs model.
an frbs object
write.frbsPMML
and frbsPMML
.
## This example shows how to construct and read frbsPMML file of frbs model ## Even though we are using MAMDANI model, other models have the same way ## ## 1. Produce frbs model, for example: we perform Wang & Mendel's technique (WM) ## ## Input data data(frbsData) data.train <- frbsData$GasFurnance.dt[1 : 204, ] data.fit <- data.train[, 1 : 2] data.tst <- frbsData$GasFurnance.dt[205 : 292, 1 : 2] real.val <- matrix(frbsData$GasFurnance.dt[205 : 292, 3], ncol = 1) range.data<-matrix(c(-2.716, 2.834, 45.6, 60.5, 45.6, 60.5), nrow = 2) ## Set the method and its parameters method.type <- "WM" control <- list(num.labels = 15, type.mf = "GAUSSIAN", type.defuz = "WAM", type.tnorm = "MIN", type.implication.func = "ZADEH", name="sim-0") ## Generate fuzzy model ## Not run: object <- frbs.learn(data.train, range.data, method.type, control) ## 2. Write frbsPMML file ## In this step, we provide two ways as follows. ## a. by calling frbsPMML() function directly. ## b. by calling write.frbsPMML() function. ## 2a. by calling frbsPMML(), the format will be displayed in R console ## Not run: frbsPMML(object) ## 2b. by calling write.frbsPMML(), the result will be saved as a file ## in the working directory. ## Not run: write.frbsPMML(object, file = "MAMDANI.GasFur") ## 3. Read frbsPMML file ## Not run: object <- read.frbsPMML("MAMDANI.GasFur.frbsPMML") ## 4. Perform predicting step ## Not run: res.test <- predict(object, data.tst)
## This example shows how to construct and read frbsPMML file of frbs model ## Even though we are using MAMDANI model, other models have the same way ## ## 1. Produce frbs model, for example: we perform Wang & Mendel's technique (WM) ## ## Input data data(frbsData) data.train <- frbsData$GasFurnance.dt[1 : 204, ] data.fit <- data.train[, 1 : 2] data.tst <- frbsData$GasFurnance.dt[205 : 292, 1 : 2] real.val <- matrix(frbsData$GasFurnance.dt[205 : 292, 3], ncol = 1) range.data<-matrix(c(-2.716, 2.834, 45.6, 60.5, 45.6, 60.5), nrow = 2) ## Set the method and its parameters method.type <- "WM" control <- list(num.labels = 15, type.mf = "GAUSSIAN", type.defuz = "WAM", type.tnorm = "MIN", type.implication.func = "ZADEH", name="sim-0") ## Generate fuzzy model ## Not run: object <- frbs.learn(data.train, range.data, method.type, control) ## 2. Write frbsPMML file ## In this step, we provide two ways as follows. ## a. by calling frbsPMML() function directly. ## b. by calling write.frbsPMML() function. ## 2a. by calling frbsPMML(), the format will be displayed in R console ## Not run: frbsPMML(object) ## 2b. by calling write.frbsPMML(), the result will be saved as a file ## in the working directory. ## Not run: write.frbsPMML(object, file = "MAMDANI.GasFur") ## 3. Read frbsPMML file ## Not run: object <- read.frbsPMML("MAMDANI.GasFur.frbsPMML") ## 4. Perform predicting step ## Not run: res.test <- predict(object, data.tst)
This function checks the consistency of a rule definition (given by the user). The rulebase consists of several fuzzy IF-THEN rules. The rules could be in a list or matrix type. Generally, there are three types of rule structures which are rules based on Mamdani, Takagi Sugeno Kang and fuzzy rule-based classification systems (FRBCS).
For rules of the Mamdani model, there are 2 parts in each rule, the antecedent and the consequent part, which are separated by "->".
rulebase(type.model, rule, func.tsk = NULL)
rulebase(type.model, rule, func.tsk = NULL)
type.model |
a value determining the type of model to use.
Here, |
rule |
a matrix or list of rules. |
func.tsk |
a matrix representing the consequent parts of rules in Takagi Sugeno Kang formulation. |
For example: r1 <- c("a1","and","b1","->", "c1")
It means that "IF input.variable1 is a1 and input.variable2 is b1 THEN output.variable is c1"
Here, ("a1", "and", "b1") is the antecedent, with "a1" and "b1" being linguistic terms, and ("c1") is the consequent part.
A fuzzy IF-THEN rule base with several rules is defined in the following way:
r1 <- c("not a1","and","b1", "->", "c1")
r2 <- c("a2","or","b2", "->", "c2")
r3 <- c("a3","or","b2", "->", "c3")
rule <- list(r1,r2,r3)
For rules of the Takagi Sugeno Kang model, the rules are at first defined without the consequent part, e.g.:
r1 <- c("a1",1,"b1","->")
r2 <- c("a2",2,"b2", "->")
r3 <- c("a3","2","b2", "->")
rule <- list(r1,r2,r3)
The consequences are defined then as a matrix fun_tsk
, which contains the linear equations of the consequences of the rules.
The dimension of this matrix is [<number_of_rules>, <number_of_variables> + 1]. The matrix has one extra column for the constants.
If there is no constant, a zero is put.
So, for example, if we have 3 rules and 2 linguistic variables (A, B), the matrix fun_tsk
has dim(3,3), as in:
func.tsk <- matrix(c(1, 1, 5, 2, 1, 3, 1, 2, 2), nrow=3, ncol=3, byrow = TRUE)
Furthermore, we can represent linguistic hedges within the rules. The kinds of hedges used are
"extremely"
reduces the truth value. For example, "extremely a1"
means membership function .
"very"
reduces the truth value. For example, "very a1"
means membership function .
"somewhat"
increases the truth value. For example, "somewhat a1"
means membership function
.
"slightly"
increases the truth value. For example, "slightly a1"
means membership function
An example of fuzzy IF-THEN rules using linguistic hedge is:
r1 <- c("very a1","and","b1","->","c1")
r2 <- c("a2",2,"b2", "->", "c2")
r3 <- c("a3","2","slightly b2", "->", "c3")
rule <- list(r1,r2,r3)
Furthermore, the following is an example in order to give names to the linguistic terms in the input and output variables.
varinput.1 <- c("a1", "a2", "a3")
varinput.2 <- c("b1", "b2")
names.varinput <- c(varinput.1, varinput.2)
names.varoutput <- c("c1", "c2", "c3")
In case of FRBCS model, the structure of rules are quite similar with Takagi Sugeno Kang model. But, instead of using linear equation in consequent part, consequent parts in FRBCS are represented by class. For example, Take into account that consequent parts expresses classes.
rule<-matrix(c("v.1_a.2","and","v.2_a.2","and","v.3_a.3","and","v.4_a.2","->","3",
"v.1_a.2","and","v.2_a.3","and","v.3_a.1","and","v.4_a.3","->","1",
"v.1_a.2","and","v.2_a.2","and","v.3_a.2","and","v.4_a.2","->","2"),
nrow=3, byrow=TRUE)
Where, "1"
, "2"
, "3"
represent class 1
, 2
, and 3
.
Noted that all rules included in rule base must have the same length as many as the number of variables.
However, we can ignore some input variables by defining the "dont_care" value.
For example, the rule ("not a1","and","dont_care", "->", "c1") refers to the rule ("not a1" "->", "c1").
Furthermore, if we are using the learning methods, the fuzzy IF-THEN rules will be generated automatically
as the outputs of frbs.learn
.
fuzzy IF-THEN rule base
defuzzifier
, inference
, and fuzzifier
This is the internal function that implements a combination of the subtractive
clustering method and fuzzy c-means. It is used to solve regression tasks. Users do not need to call it directly,
but just use frbs.learn
and predict
SBC(data.train, range.data.ori, r.a = 0.5, eps.high = 0.5, eps.low = 0.15)
SBC(data.train, range.data.ori, r.a = 0.5, eps.high = 0.5, eps.low = 0.15)
data.train |
a matrix ( |
range.data.ori |
a matrix ( |
r.a |
the radius defining a neighborhood. |
eps.high |
an upper threshold value. |
eps.low |
a lower threshold value. |
This method was proposed by S. Chiu. For generating the rules in the learning phase, the subtractive clustering method is used to obtain the cluster centers. Subtractive clustering (SBC) is an extension of Yager and Filev's mountain method. SBC considers each data point as a potential cluster center by determining the potential of a data point as a function of its distances to all the other data points. A data point has a high potential value if that data point has many nearby neighbors. The highest potential is chosen as the cluster center and then the potential of each data point will be updated. The process of determining new clusters and updating potentials repeats until the remaining potential of all data points falls below some fraction of the potential of the first cluster center. After getting all the cluster centers from subtractive clustering, the cluster centers are optimized by fuzzy c-means.
R. Yager and D. Filev, "Generation of fuzzy rules by mountain clustering," J. of Intelligent and Fuzzy Systems, vol. 2, no. 3, pp. 209 - 219 (1994).
S. Chiu, "Method and software for extracting fuzzy classification rules by subtractive clustering", Fuzzy Information Processing Society, NAFIPS, pp. 461 - 465 (1996).
SBC.test
, frbs.learn
, and predict
This function is the internal function of the SBC method to compute the predicted values.
SBC.test(object, newdata)
SBC.test(object, newdata)
object |
the |
newdata |
a matrix ( |
This is the internal function that implements
the structural learning algorithm on vague environment (SLAVE). It is used to handle classification tasks.
Users do not need to call it directly,
but just use frbs.learn
and predict
.
SLAVE(data.train, persen_cross = 0.6, persen_mutant = 0.3, max.iter = 10, max.gen = 10, num.labels, range.data.input, k.lower = 0.25, k.upper = 0.75, epsilon = 0.1)
SLAVE(data.train, persen_cross = 0.6, persen_mutant = 0.3, max.iter = 10, max.gen = 10, num.labels, range.data.input, k.lower = 0.25, k.upper = 0.75, epsilon = 0.1)
data.train |
a matrix ( |
persen_cross |
a real number between 0 and 1 representing the probability of crossover. |
persen_mutant |
a real number between 0 and 1 representing the probability of mutation. |
max.iter |
the maximal number of iterations. |
max.gen |
the maximal number of generations for the genetic algorithm. |
num.labels |
a number of the linguistic terms. |
range.data.input |
a matrix containing the ranges of the normalized input data. |
k.lower |
a lower bound of the noise threshold. |
k.upper |
an upper bound of the noise threshold. |
epsilon |
a value between 0 and 1 representing the covering factor. |
This method is adopted from A. Gonzalez and R. Perez's paper which is applied for classification problems. SLAVE is based on the iterative rule learning approach which means that we get only one fuzzy rule in each execution of the genetic algorithm. In order to eliminate the irrelevant variables in a rule, SLAVE has a structure composed of two parts: the first part is to represent the relevance of variables and the second one is to define values of the parameters. The following steps are conducted in order to obtain fuzzy rules:
Step 1: Use the genetic algorithm process to obtain ONE RULE for the system.
Step 2: Collect the rule into the final set of rules.
Step 3: Check and penalize this rule.
Step 4: If the stopping criteria is satisfied, the system returns the set of rules as solution. Otherwise, back to Step 1.
This method uses binary codes as representation of the population and applies the basic genetic operators, i.e., selection, crossover, and mutation on it. And, the best rule is obtained by calculating the degree of consistency and completeness.
A. Gonzalez and R. Perez, "Selection of relevant features in a fuzzy genetic learning algorithm", IEEE Transactions on Systems, Man, and Cybernetics, Part B: Cybernetics, vol. 31, no. 3, pp. 417 - 425 (2001).
This function is the internal function of the SLAVE method to compute the predicted values.
SLAVE.test(object, newdata)
SLAVE.test(object, newdata)
object |
the |
newdata |
a matrix ( |
A matrix of predicted values.
This function enables the output of a summary of the frbs-object
.
## S3 method for class 'frbs' summary(object, ...)
## S3 method for class 'frbs' summary(object, ...)
object |
the |
... |
the other parameters (not used) |
This function displays several components of the object. The components of one particular method can be different from components of other methods. The following is a description of all components which might be printed.
The name of the model: A name given by the user representing the name of the simulation or data or model.
Model was trained using: It shows which method we have been used.
The names of attributes: a list of names of training data.
The interval of training data: It is a matrix representing the original interval of data where the first and second rows are minimum and maximum of data, respectively. The number of columns represents the number of variables.
Type of FRBS model: a description expresses one of the following FRBS model available such as "MAMDANI"
, "TSK"
,
"FRBCS"
, "CLUSTERING"
, "APPROXIMATE"
, and "2TUPPLE"
.
Type of membership function: a description expresses one of the following shapes of membership functions:
"GAUSSIAN"
, code"TRIANGLE", "TRAPEZOID"
, "SIGMOID"
, and "BELL"
.
Type of t-norm method: a description expresses one of the following type of t-norm: "MIN"
, "PRODUCT"
, "HAMACHER"
,
"YAGER"
, and "BOUNDED"
.
Type of s-norm method: a description expresses one of the following type of s-norm: "MAX"
, "SUM"
, "HAMACHER"
,
"YAGER"
, and "BOUNDED"
.
Type of defuzzification technique: a description expresses one of the following types: "WAM"
, "FIRST_MAX"
,
"LAST_MAX"
, "MEAN_MAX"
, and "COG"
.
Type of implication function: a description expresses one of the following types:
"DIENES_RESHER"
, "LUKASIEWICZ"
,
"ZADEH"
, "GOGUEN"
, "GODEL"
, "SHARP"
, "MIZUMOTO"
,
"DUBOIS_PRADE"
, and "MIN"
.
The names of linguistic terms of the input variables: These names are generated
automatically by frbs expressing all linguistic terms considered. Generally,
these names are built by two parts which are the name of variables expressed
by "v"
and the name of linguistic terms of each variables represented by "a"
.
For example, "v.1_a.1"
means the linguistic value "a.1"
of the first variable (v.1).
However, we provide different format if we set the number of linguistic terms (num.labels
) to 3, 5, 7.
For example, for the number of label 3, it will be "small"
, "medium"
, and "large"
.
The names of linguistic terms of the output variable: For the Mamdani model, since the frbs package only considers
single output, the names of the linguistic terms for the output variable
are simple and clear and start with "c"
. However, for the Takagi Sugeno Kang model and
fuzzy rule-based classification systems, this component is always NULL
.
The parameter values of membership functions of the input variables (normalized):
It is represented by a matrix () where n depends on the number of
linguistic terms on the input variables and the first row of the matrix describes
a type of membership function, and the rest of rows are
their parameter values.
For example, label
"v.1_a.2"
has value
4.0, 0.23, 0.43, 0.53, 0.73 on its column. It means that the label a.2 of variable v.1
has a parameter as follows.
4.0 on the first row shows TRAPEZOID
shape in the middle position,
while 0.23, 0.43, 0.53, and 0.73 are corner points of a TRAPEZOID
.
Furthermore, the following is the complete list of shapes of membership functions:
TRIANGLE
: 1 on the first row and rows 2, 3, and 4 represent corner points.
TRAPEZOID
: 2, 3, or 4 on the first row means they are TRAPEZOID
in left, right and middle side, respectively,
and rows 2, 3, 4, and 5 represent corner points. But for TRAPEZOID
at left or right side the fifth row is NA
.
GAUSSIAN
: 5 on the first row means it uses GAUSSIAN
and second and third row represent mean and variance.
SIGMOID
: 6 on the first row and two parameters (gamma and c) on second and third rows.
BELL
: 7 on the first row and three parameters (a, b, c) on second, third, and fourth rows.
The fuzzy IF-THEN rules: In this package, there are several models for representing fuzzy IF-THEN rules based on the method used.
the Mamdani model: they are represented as a knowledge base containing two parts: antecedent and consequent parts which are separated by a sign "THEN", as for example in the following rule:
IF var.1 is v.1_a.1 and var.2 is v.2_a.2 THEN var.3 is c.2
the Takagi Sugeno Kang model: In this model, this component only represents the antecedent of rules while the consequent part will be represented by linear equations.
fuzzy rule-based classification systems (FRBCS): This model is quite similar to the Takagi Sugeno Kang model, but the consequent part expresses pre-defined classes instead of a simplify of linear equations.
approximate approach: Especially for GFS.FR.MOGUL
, a matrix of parameters
of membership functions is used to represent the fuzzy IF-THEN rules as well.
The representation of rules and membership functions is a matrix () where
n is the number of rules and m is the number of variables while p is the number of corner points
of the membership function, if we are using
TRIANGLE
or TRAPEZOID
then p = 3 or 4, respectively.
For example, let us consider the triangular membership function and a number of variables of 3.
The representation of rules and membership functions is as follows:
<<a11 a12 a13>> <<b11 b12 b13>> <<c11 c12 c13>>
.
The linear equations on consequent parts of fuzzy IF-THEN rules: It is used in the Takagi Sugeno Kang model.
The weight of the rules or the certainty factor: For the FRBCS.W
method, this shows the weight related to the rules
representing the ratio of dominance among the rules.
The cluster centers: This component is used in clustering methods representing cluster centers.
This is the internal function that implements the model proposed by L. X. Wang and J. M.
Mendel. It is used to solve regression task. Users do not need to call it directly,
but just use frbs.learn
and predict
WM(data.train, num.labels, type.mf = "GAUSSIAN", type.tnorm = "PRODUCT", type.implication.func = "ZADEH", classification = FALSE, range.data = NULL)
WM(data.train, num.labels, type.mf = "GAUSSIAN", type.tnorm = "PRODUCT", type.implication.func = "ZADEH", classification = FALSE, range.data = NULL)
data.train |
a matrix ( |
num.labels |
a matrix ( |
type.mf |
the type of the membership function. See |
type.tnorm |
a value which represents the type of t-norm. See |
type.implication.func |
a value representing type of implication function. Let us consider a rule,
|
classification |
a boolean representing whether it is a classification problem or not. |
range.data |
a matrix representing interval of data. |
The fuzzy rule-based system for learning from L. X. Wang and J. M. Mendel's paper is implemented in this function. For the learning process, there are four stages as follows:
Step 1:
Divide equally the input and output spaces of the given numerical data into
fuzzy regions as the database. In this case, fuzzy regions refers to intervals for each
linguistic term. Therefore, the length of fuzzy regions represents the number of
linguistic terms. For example, the linguistic term "hot" has the fuzzy region .
We can construct a triangular membership function having the corner points
,
, and
where
is a middle point
that its degree of the membership function equals one.
Step 2:
Generate fuzzy IF-THEN rules covering the training data,
using the database from Step 1. First, we calculate degrees of the membership function
for all values in the training data. For each instance in the training data,
we determine a linguistic term having a maximum degree in each variable.
Then, we repeat the process for each instance in the training data to construct
fuzzy rules covering the training data.
Step 3:
Determine a degree for each rule.
Degrees of each rule are determined by aggregating the degree of membership functions in
the antecedent and consequent parts. In this case, we are using the product aggregation operators.
Step 4:
Obtain a final rule base after deleting redundant rules.
Considering degrees of rules, we can delete the redundant rules having lower degrees.
The outcome is a Mamdani model. In the prediction phase, there are four steps: fuzzification, checking the rules, inference, and defuzzification.
L.X. Wang and J.M. Mendel, "Generating fuzzy rule by learning from examples", IEEE Trans. Syst., Man, and Cybern., vol. 22, no. 6, pp. 1414 - 1427 (1992).
frbs.learn
, predict
and frbs.eng
.
It is a function used to save an FRBS model to the .frbsPMML file. Detailed information about frbsPMML can be seen in frbsPMML
.
write.frbsPMML(object, fileName = NULL)
write.frbsPMML(object, fileName = NULL)
object |
a frbsPMML object which is an object produced by |
fileName |
a file name with extension |
a file containing an FRBS model in frbsPMML format
A. Guazzelli, M. Zeller, W.C. Lin, and G. Williams., "pmml: An open standard for sharing models", The R Journal, Vol. 1, No. 1, pp. 60-65 (2009).
Data Mining Group, http://www.dmg.org/.
read.frbsPMML
and frbsPMML
.
## This example shows how to construct frbsPMML file of frbs model ## Even though we are using MAMDANI model, other models have the same way ## ## 1. Produce frbs model, for example: we perform Wang & Mendel's technique (WM) ## ## Input data data(frbsData) data.train <- frbsData$GasFurnance.dt[1 : 204, ] data.fit <- data.train[, 1 : 2] data.tst <- frbsData$GasFurnance.dt[205 : 292, 1 : 2] real.val <- matrix(frbsData$GasFurnance.dt[205 : 292, 3], ncol = 1) range.data<-matrix(c(-2.716, 2.834, 45.6, 60.5, 45.6, 60.5), nrow = 2) ## Set the method and its parameters method.type <- "WM" control <- list(num.labels = 15, type.mf = "GAUSSIAN", type.defuz = "WAM", type.tnorm = "MIN", type.implication.func = "ZADEH", name = "sim-0") ## Generate fuzzy model ## Not run: object <- frbs.learn(data.train, range.data, method.type, control) ## 2. Write frbsPMML file ## In this step, we provide two steps as follows: ## a. by calling frbsPMML() function directly. ## b. by calling write.frbsPMML() function. ## 2a. by calling frbsPMML(), the frbsPMML format will be displayed in R console ## Not run: pmml.obj <- frbsPMML(object) ## 2b. by calling write.frbsPMML(), the result will be saved as a file ## in the working directory. ## Not run: write.frbsPMML(pmml.obj, file = "MAMDANI.GasFur")
## This example shows how to construct frbsPMML file of frbs model ## Even though we are using MAMDANI model, other models have the same way ## ## 1. Produce frbs model, for example: we perform Wang & Mendel's technique (WM) ## ## Input data data(frbsData) data.train <- frbsData$GasFurnance.dt[1 : 204, ] data.fit <- data.train[, 1 : 2] data.tst <- frbsData$GasFurnance.dt[205 : 292, 1 : 2] real.val <- matrix(frbsData$GasFurnance.dt[205 : 292, 3], ncol = 1) range.data<-matrix(c(-2.716, 2.834, 45.6, 60.5, 45.6, 60.5), nrow = 2) ## Set the method and its parameters method.type <- "WM" control <- list(num.labels = 15, type.mf = "GAUSSIAN", type.defuz = "WAM", type.tnorm = "MIN", type.implication.func = "ZADEH", name = "sim-0") ## Generate fuzzy model ## Not run: object <- frbs.learn(data.train, range.data, method.type, control) ## 2. Write frbsPMML file ## In this step, we provide two steps as follows: ## a. by calling frbsPMML() function directly. ## b. by calling write.frbsPMML() function. ## 2a. by calling frbsPMML(), the frbsPMML format will be displayed in R console ## Not run: pmml.obj <- frbsPMML(object) ## 2b. by calling write.frbsPMML(), the result will be saved as a file ## in the working directory. ## Not run: write.frbsPMML(pmml.obj, file = "MAMDANI.GasFur")