<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta>
      <journal-title-group>
        <journal-title>May</journal-title>
      </journal-title-group>
    </journal-meta>
    <article-meta>
      <title-group>
        <article-title>On linear regression for fuzzy data of different quality </article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Serhii Mashchenko</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Oleksandr Marchenko</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Taras Shevchenko National University of Kyiv</institution>
          ,
          <addr-line>64/13, Volodymyrska Street, City of Kyiv, 01601</addr-line>
          ,
          <country country="UA">Ukraine</country>
        </aff>
      </contrib-group>
      <pub-date>
        <year>2023</year>
      </pub-date>
      <volume>3</volume>
      <issue>2023</issue>
      <fpage>0000</fpage>
      <lpage>0003</lpage>
      <abstract>
        <p>   The present paper is devoted to the linear regression for a fuzzy set of fuzzy data samples. This model allows one to take into account the data of different quality. It is shown that regression parameters are type-2 fuzzy sets. Furthermore, the corresponding type-2 membership functions are given. The decomposition approach is used to investigate the T2FSs of linear regression parameters. It is shown that each T2FS of regression parameter can be decomposed according to secondary membership grades into a finite collection of fuzzy numbers. Each of these fuzzy number is the corresponding fuzzy regression parameter for a set of data numbers. This set is the corresponding α-cut of the original fuzzy set of fuzzy data samples. The illustrative example is given.</p>
      </abstract>
      <kwd-group>
        <kwd>1  Linear regression</kwd>
        <kwd>fuzzy least squares estimator</kwd>
        <kwd>fuzzy number</kwd>
        <kwd>type-2 fuzzy set</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>1. Introduction </title>
      <p>The classical regression analysis is based on crisp data and a crisp relationship between the
dependent variable and the independent variables. In practice, there are many situations in which
observations cannot be measured as crisp quantities, because the information is often fuzzy,
incomplete, linguistic or noisy. Fuzzy regression analysis is a non-statistical method based on a fuzzy
set theory rather than probability theory (see [1]). In a general model of a fuzzy regression both input
and output are fuzzy. In this regard, the fuzzy regression model contains fuzzy parameters instead of
error terms.</p>
      <p>The three main fields can be distinguished in fuzzy regression analysis. These are possibilistic
regression analysis, fuzzy least squares methods and machine learning techniques. The probabilistic
approach in fuzzy regression analysis was first proposed by Tanaka et al. [2]. Unlike conventional
regression analysis, where deviations between observed and predicted values reflect a measurement
error, deviations in fuzzy regression reflect the uncertainty of the system structure expressed by fuzzy
parameters of the regression model. Fuzzy parameters of the model are considered to be distributions
of possibilities and determined by solving a linear programming problem that allows one to minimize
fuzzy deviations subject to membership degrees constraints. Since the membership functions (MFs) of
fuzzy sets (FSs) can be viewed as probability distributions, this approach was called ‘possibilistic
regression analysis.’ The possibilistic approach was explored and improved by many authors.
Reviews of possibilistic regression analysis can be found in D’Urso [3].</p>
      <p>Fuzzy regression analysis was also considered from the viewpoint of generalizing the classical
least squares method to the case of fuzzy data. The idea of this approach is to minimize in some sense
the different distance measures between the predicted fuzzy values and the given fuzzy data. Fuzzy
least squares methods were first proposed by Celmins [4]. Later, this approach was significantly
developed by many researchers. In [1] one can find qualitative reviews on the fuzzy least squares and
fuzzy least absolute methods.</p>
      <p>Machine learning techniques made it possible to generalize fuzzy regression analysis through the
use of genetic algorithms, neural networks, and support vector machines. Relevant references can be
found in Chukhrova and Johannssen [1] and Hastie et al. [5].</p>
      <p>Often, when solving applied problems, data of different quality can be used [6]. For example,
according to [7], wind tunnel experiments provide high simulation accuracy (a source of high-fidelity
data). Also, experiments based on computational physical models have a higher error (source of low
fidelity data). In some applications, a significantly more accurate regression model can be built if
lowprecision data are also used. In this case, the problem arises of constructing a regression based on data
of different quality.</p>
      <p>Using data of different quality with the aim of improving model accuracy is not a new concept.
For instance, Hevesi et al. [8] predict average annual precipitation values near a potential nuclear
waste disposal site using a set of precipitation measurements from the region along with more easily
obtainable elevation map of the area. Kennedy and O’Hagan [9] approach the subject from the
perspective of model construction using data resulting from computational simulations of varying
fidelities and costs. In [7], to construct a Gaussian regression model, the problem of planning an
experiment is solved with the choice of the ratio between the sizes of samples of low-precision and
high-precision data. Also, to process data with a variable degree of certainty the methods of transfer
learning [10], space mapping [11], and others [12] are used.</p>
      <p>Data quality (degree of usability of data) is a complex concept. It is characterized by objectivity,
integrity, relevance, measurability, controllability, etc. In some cases, the quality of data may not be
crisp defined [13]. According to [14], where using data without expert knowledge, the choice of a
representative sample becomes an NP-complete problem. Therefore, samples have to be found within
a reasonable time, and this justifies the use of fuzzy methods that formalize expert knowledge
expressed in natural language words. For instance, in the framework of quality control fuzzy expert
assessments are used in [15] to construct acceptance sampling plans (how many units can be selected
from a consignment and how many defective units are allowed in this sample).</p>
      <p>In this article, we intend to investigate the method of constructing a fuzzy linear regression in the
case when fuzzy data samples are of different qualities. Furthermore, the degrees of membership to a
FS are known for these samples. This leads to a possibility that the regression takes into account
fuzziness of the quality assessments of different samples, rather than just uncertainty of data.
Examples of such FSs data samples could be: ‘High quality data samples’, ‘Questionable data
samples’, ‘Actual data samples’, etc.</p>
      <p>The main result of the article justifies the fact that a FS of different quality samples of type-1 fuzzy
data generates a type-2 fuzzy regression model. In this model, the regression parameters are T2FSs on
the real line with constant secondary grades. Although, in general, a T2FS is a rather complicated
mathematical object, T2FSs with constant secondary grades are simple enough for practical use. This
feature allows us to decompose this set by secondary grades into a collection of corresponding fuzzy
numbers. Each of them represents the corresponding fuzzy regression parameter for a crisp set of data
samples. The set in the focus is the α-cut of the original fuzzy set of data samples. We note that the
well-known type 2 fuzzy regression models use crisp collections of type 2 data sets, while the model
proposed in the article uses a fuzzy collection of type 1 fuzzy data samples. This is a principal
difference between them. It should also be added that this article continues the line of research in the
field of mathematical operations with a fuzzy set of operands, first introduced in the context of
intersections and unions of fuzzy sets [16, 17].</p>
    </sec>
    <sec id="sec-2">
      <title>2. Materials and Methods </title>
      <p>In this section, we briefly review some existing theories and definitions.</p>
    </sec>
    <sec id="sec-3">
      <title>2.1. Linear  regression  analysis  for  fuzzy  input  and  output  data  using  the  extension principle </title>
      <p>The article focuses on a linear regression in the case when data samples form a FS. We stress that
one could have exploited different known methods of constructing a linear regression for fuzzy data.
As an alternative, we intend to modify the method [18] based on the extension principle.</p>
      <p>Let K  {1,..., K } be the set of indices of data samples {yi , xi1,..., xip}, i  K , where K is the
cardinality of K . A crisp statistical linear regression has the form</p>
      <p>
        yi  xi ( (K ))T   i ,  i  K ,  (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ) 
where for each i  K , yi is the dependent variable; xi  (1, xi1,..., xip ) is the vector of independent
 (K )  (0 (K ),..., p (K ))  (l (K ))l0,..., p
symbol
denotes
the
      </p>
      <p>transpose,
variables (factors, regressors) xil , l  1,..., p ;  i is the independent normal random variable. The
T
p
is
is
the
the number
vector of</p>
      <p>
        of independent variables,
regression parameters. Let
y(K )  X (K )( (K ))T   (K ) be the matrix notation of equations (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ) with X (K )  {xis}iK ,s0,...p and
xi0  1 for all i  K , y(K )  ( yi )TiK and  (K )  ( i )TiK . For the convenience of presentation, the set
K of sample indices is indicated hereinafter as a parameter in these formulae. According to the least
squares method, the estimate ˆ(K )  (ˆl (K ))l0,..., p of the parameter vector  (K )  (l (K ))l0,..., p has
the form
      </p>
      <p>
        ˆ(K )  [ X T (K ) X (K )]1 X T (K ) y(K ) . 
Assume that the data is fuzzy. To generalize formula (
        <xref ref-type="bibr" rid="ref2">2</xref>
        ) we denote by
      </p>
      <p>fl ( X (K ), y(K )) ˆl (K )  
the l -th element, l  0,..., p of the estimate ˆ(K ) . We also denote by</p>
      <p>X (K )  {xis}iK,s0,...p ,  y(K )  ( yi )TiK  </p>
      <p>
        1, r [l (K )] ;
1[l (K )]u (r)  0, r [l (K )]uu ;  
r 
, u [0,1] . According to [18] and Remark 1, formula (
        <xref ref-type="bibr" rid="ref5">5</xref>
        ) implies that u -cut [l (K )]u of the FN
l (K ) (r)  X (mK)a,yx(K ){iKm,si0n,..., p
      </p>
      <p>{ xis (xis ), yi ( yi )}: r  fl ( X (K ), y(K )),
X (K )  {xis}iK;s0,..., p 
p K , y(K )  {yi}iK </p>
      <p>
        K }
 
by Zadeh’s extension principle [20]. Here, fl ( X (K ), y(K )) is the l-th element of the vector
ˆ(K )  [ X T (K ) X (K )]1 X T (K ) y(K ) by (
        <xref ref-type="bibr" rid="ref2">2</xref>
        ) and (
        <xref ref-type="bibr" rid="ref3">3</xref>
        ). According to Remark 1, the maximum in (
        <xref ref-type="bibr" rid="ref5">5</xref>
        )
exists. As shown in [21], the representation of FNs by u -cuts is simpler for calculations than the
functional approach. Therefore, for each l  0,..., p , we represent the MF of the FN l (K ) in the
form
l (K ) (r)  max u1[l (K )]u (r) ,  (
        <xref ref-type="bibr" rid="ref6">6</xref>
        ) 
      </p>
      <p>u[0,1]
where [l (K )]u is the u -cut of the FN l (K ) . This u -cut is the set [l (K )]u  {r  : l (K ) (r)  u}
(K )  [ X T (K ) X (K )]1 X T (K ) y (K ) according to the least squares method [18]. For each l  0,..., p ,
l (K )  fl ( X (K ), y(K ))  {(r,l (K ) (r)) : r  } is the FN with the MF
with the MF
l (K ) has the form
the matrix of independent variables and the vector of dependent variables, respectively, where for
each i  K , xi0  1 is the crisp number which is equal to 1; xis , s  1,..., p and yi are fuzzy numbers
(FNs) with the MFs  xis (xis ) , xis  , s  1,..., p and  yi ( yi ) , yi  , i  K , respectively. Here,
is the real line.</p>
      <p>Remark 1. Recall that a FN is a normal FS on
concave MF (for example, see [19]).</p>
      <p>
        The vector (K )  (l (K ))l0,..., p
with the upper semicontinuous and
quasiof fuzzy parameters of the regression has the form
(
        <xref ref-type="bibr" rid="ref2">2</xref>
        ) 
(
        <xref ref-type="bibr" rid="ref3">3</xref>
        ) 
(
        <xref ref-type="bibr" rid="ref4">4</xref>
        ) 
(
        <xref ref-type="bibr" rid="ref5">5</xref>
        ) 
(
        <xref ref-type="bibr" rid="ref7">7</xref>
        ) 
respectively, of the real line
[l (K )]u  [[l (K )]uD ,[l (K )]uH ] 
      </p>
      <p>too. Equality (8) entails
[l (K )]u  { fl ( X (K ), y(K )) : xis [xis ]u , s  0,..., p; yi [ yi ]u , i  K} , 
where for each i  K , the u -cuts of the FNs xis , s  0,... p and yi are closed intervals
]D ,[xis ]uH ], s  0,..., p  and  [ yi ]u  [[ yi ]uD ,[ yi ]uH ] , 
 [xil ]u  [[xis u</p>
      <p>
        . Since l (K ) is a FN, then its u -cut [l (K )]u is the interval
[l (K )]uD  min{ fl ( X (K ), y(K )) : xis [xis ]u , s  0,..., p; yi [ yi ]u , i  K} , 
[l (K )]uH  max{ fl ( X (K ), y(K )) : xis [xis ]u , , s  0,..., p; yi [ yi ]u , i  K} . 
Thus, (
        <xref ref-type="bibr" rid="ref6">6</xref>
        ) ensures that formula (
        <xref ref-type="bibr" rid="ref5">5</xref>
        ) and the FN l (K ) have the forms
      </p>
      <p>l (K ) (r)  max{u [0,1]: [l (K )]uD  r  [l (K )]uH } ,  r 
l (K )  {(r,u) : r [[l (K )]uL ,[l (K )]uH ],u [0,1]}  {([[l (K )]uL ,[l (K )]uH ],u) : u [0,1]} ,
respectively. Problems (10) and (11) are rather complicated. In view of this, it is suggested in [18] to
use, for l  0,..., p , the approximate value of the l -th fuzzy parameter</p>
      <p>l (K )  {([[l (K )]uL ,[l (K )]uH ],u) : u [0,1]} , 
l (K ) (r)  max{u [0,1] : [l (K )]uD  r  [l (K )]uH } ,  r 
fl ([ X (K )]uH ,[ y(K )]uD ), fl ([ X (K )]uH ,[ y(K )]uH )},
[l (K )]uD  min{ fl ([ X (K )]uD ,[ y(K )]uD ), fl ([ X (K )]uD ,[ y(K )]H ,
u  
fl ([ X (K )]uH ,[ y(K )]uD ), fl ([ X (K )]uH ,[ y(K )]uH )}.
[l (K )]uH  max{ fl ([ X (K )]uD ,[ y(K )]uD ), fl ([ X (K )]uD ,[ y(K )]H ,
u  
where
Here, the matrices [ X (K )]uD , [ X (K )]uH and the vectors [ y(K )]uD , [ y (K )]uH are comprised in of the
elements [xis ]uD ,[xis ]uH , s  0,..., p; [ yi ]uD ,[ yi ]uH , i  K , respectively (see (9)). It is clear that
[l (K )]uD  [l (K )]uD , [l (K )]uH  [l (K )]uH . Therefore, the inclusion [l (K )]u  [l (K )]u is hold for
the closed interval [l (K )]u  [[l (K )]uD ,[l (K )]uH ] . This entails l (K ) (r)  l (K ) (r) and thereupon
the value l (K ) (r) is a lower bound of l (K ) (r) for r 
.
2.2.</p>
    </sec>
    <sec id="sec-4">
      <title>Type‐2 fuzzy sets </title>
      <p>Zadeh [22] introduced the notion of type-2 fuzzy set (T2FS) as a generalization of type-1 fuzzy set
(T1FS) (that is, an ordinary FS). Unlike a T1FS, the membership degree of elements in a T2FS is a FS
on closed interval [0, 1]. Based on the ideas of Karnik and Mendel [23], Mendel and John [24] gave a
different definition. The T2FS C on a set X is the collection</p>
      <p>C  {((x,u),C (x,u)) : x  X , u  J x  [0,1]} , 
where  C (x,u) is a type-2 membership function (T2MF), J x is the set of primary membership
degrees u of x  X to C . The value  C (x,u) is a crisp number from the closed interval [0, 1] which
is called the secondary grade of the pair (x,u) .</p>
      <p>Remark 2. According to the comments of Harding et al. [25] and Aisbett et al. [26], one has to
define the T2MF  C (x,u) for all x  X and u [0,1] . To this end, one should put  C (x,u)  0 for all
u  J x , x  X .</p>
      <p>Remark 3. The primary membership degree u is deemed as the degree of manifestation of some
property (which determines the given FS) of x  X . According to [27], we interpret the secondary
grade  C (x,u) as the degree of truth of the corresponding primary degree u of this property for x.</p>
      <p>Following
[24],
we
define
embedded
, 
, 
(8) 
(9) 
(10) 
(11) 
(12) 
(13) 
(14) 
(15) 
(16) 
(17) 
C  {((x,u),C (x,u)) : x  X ,u [0,1]} . Assume that ux  Ce1 (x) [0,1] is a unique primary degree of
membership, for
each
x  X ,
where
Ce1 (x) ,
x  X
is the</p>
      <p>MF
of the</p>
      <p>T1FS
Ce1  {(x,Ce1 (x)) : x  X }. This T1FS is called embedded in the T2FS C . We define the embedded
T2FS C e2 in C in the form C e2  {((x,ux ),Ce2 (x,ux )) : x  X } with Ce2 (x,ux )) C (x,Ce1 (x))) for
all x  X .</p>
      <p>Remark 4. According to [24], each element of the type-2 fuzzy collection
C  {((x,u),C (x,u)) : x  X ,u [0,1]} is interpreted as a subset. Thus, the collection is represented as
the classical union of its elements in the sense of T1FSs. In [24], Mendel and John stated that each
T2FS can be represented as a collection of all possible embedded T2FSs.
2.3.</p>
    </sec>
    <sec id="sec-5">
      <title>Collections of T2FSs with constant secondary grades </title>
      <p>We shall need one special case of a T2FS to be defined according to [28, 29, 30]. Let
A  { C (x,u) : C (x,u)  0, x  X ,u [0,1]} be the set of all possible positive values of secondary
grades for the T2FS C  {((x,u),C (x,u)) : x  X , u [0,1]} . Assume that the set A is finite.</p>
      <p>Definition 1. We say that an embedded T2FS C e2 ( )  {((x,ux ),Ce2 ( ) (x,ux )) : x  X } in the T2FS
 has a constant secondary grade   A if, for each x  X , the unique primary degree
C
ux  Ce1( ) (x) [0,1] exists for which Ce2 ( ) (x,ux )  , i.e. C e2 ( )  {((x,Ce1( ) (x)), ) : x  X } .</p>
      <p>Here Ce1( ) (x) , x  X is the MF of the embedded T1FS C e1( )  {(x,Ce1( ) (x)) : x  X } in the
T2FS C .</p>
      <p>Remark 5. Obviously, for the T2FS C and each   A , there is the unique embedded T1FS
C e1( )  {(x,Ce1( ) (x)) : x  X } which is corresponding to the embedded T2FS C e2 ( ) with
constant secondary grade  . Hence,
 {((x,Ce1( ) (x)), ) : x  X } .</p>
      <p>C e2 ( )  {(Ce1( ), )}  {({(x,Ce1( ) (x)) : x  X }, )} </p>
    </sec>
    <sec id="sec-6">
      <title>3. Formulation of the problem and the main idea </title>
      <p>
        Let N  {1,..., n} be the set of indices of fuzzy data samples {yi , xi1,..., xip}, i  N in the form of
FNs with the MFs  xis (xis ) , xis  , s  1,..., p ,  yi ( yi ) , yi  , i  N , respectively. The matrix
X (N ) of independent variables and the vector y(N ) of dependent variables are given by formula
(
        <xref ref-type="bibr" rid="ref4">4</xref>
        ). Assume that qualities of data samples {yi , xi1,..., xip}, i  N are different. Furthermore, the
degrees  I ( j), j  N of membership to the FS I  {( j, I ( j)) : j  N} of data samples indices are
known. The following question arises: ‘what are linear regression parameters in the case when fuzzy
data samples are involved in the calculation with the corresponding degrees  I ( j) , j  N of
membership?’ In other words: ‘what are the fuzzy parameters of the regression for the FS I of the
data samples indices? We investigate this problem for the l -th regression parameter. First, we
generalize formula (12) for the case of an arbitrary subset K  N of sample indices and represent it
in a convenient form for us. For each r  , we consider the mapping Ulr : 2N  [0,1] given by
Ulr (K )  max{u [0,1] : [l (K )]uD  r  [l (K )]uH }, K  N . 
(18) 
According to (12), the mapping Ulr associates each subset K  N of data sample indices with the
value of the MF l (K ) (r) of the fuzzy parameter l (K ) (see (13)). The latter is the FN
(19) 
(20) 
(21) 
(22) 
      </p>
      <p>l (K )  {(r,l (K ) (r)) : r  }, 
l (K ) (r)  Ulr (K ) ,   r  supp(l (K ))  {r  :Ulr (K )  0} , 
where supp(l (K )) is the support of the FN l (K ) . Next, we generalize formulae (19) and (20) to
the case of the FS I of sample indices. We denote by Bl the l -th regression parameter, and by
MBl (r) , r  its MF for the FS I of data samples indices. In this case, the value of the MF MBl (r)
for each fixed r  r * coincides with the image U lr* (I ) of the FS I under Ulr* , i.e. M Bl (r*)  Ulr* (I )
. According to Zadeh’s extension principle [20], the image of the FS I under the mapping
Ulr* : 2N  [0,1] (see (18)) is the FS Ulr*(I )  {(u,Ulr*(I ) (u)) : u [0,1]} with the MF
Ulr*(I ) (u)  max{ [0,1] : u  Ulr* (I ( ))} ,  u [0,1] . 
Here, I ( )  { j  N :  I ( j)  } is the  -cut,  [0,1] of the FS I  {( j, I ( j)) : j  N} of sample
indices;</p>
      <p>Ulr* (I ( ))  l (I ( )) (r*) , 
is the image of the  -cut I ( ) ,  [0,1] of the FS I of the samples indices in the mapping Ulr* (see
(18)). The value Ulr* (I ( )) is equal to the MF value l (I ( )) (r*) of the l -th fuzzy parameter
l (I ( )) for the set I ( ) of sample indices.</p>
      <p>Remark 6. Let A  { I ( j) : j  N} be the set of membership degrees values of the fuzzy set
I  {( j, I ( j)) : j  N} of sample indices. Note that the cardinality of the set A is |A| ≤ n. The
situation |A| &lt; n may occur if the degrees of membership  I ( j) coincide for different indices j  N
of samples. It is clear that while obtaining  -cuts I ( )  { j  N :  I ( j)  }   of the fuzzy set I
we can assume that   A rather than  [0,1] .</p>
      <p>Thus, according to (20) and (21), for fixed r  r * , values of M Bl (r*) form the T1FS
{(u, MBl (r*) (u)) : u [0,1]} on [0, 1] with the MF  MBl (r*) (u)  Ulr*(I ) (u)  max{  A : u  Ulr* (I ( ))},
u  supp(M Bl (r*))  {u [0,1] : u  Ulr* (I ( )),   A} . Then (22) entails
 MBl (r*) (u)  max{  A : u  l (I ( )) (r*)},  
(23) 
where u supp(M Bl (r*))  {u [0,1]: u  l (I ( )) (r*),   A}. Therefore, we conclude that Bl is a
FS on with the MF whose values form T1FS on [0,1] . Then, according to [22], Bl is the T2FS on
. In the manner of vertical slices [27] the T2FS Bl on has the form:</p>
      <p>Bl  {(r, M Bl (r)) : r  R}  {(r,{(u, MBl (r) (u)) : u  Jr }) : r  } ,  (24) 
where  MBl (r) (u) , u [0,1] is the MF of the T1FS M Bl (r)  {{(u, MBl (r) (u)) : u [0,1]} of values of
fuzzy degree of membership of the element r  to the T2FS Bl and Jr  supp(M Bl (r)) is the set of
primary membership degrees, where supp(M Bl (r)) is the support of the T1FS MBl (r) for r  .
According to Section 2.2, we can also characterize the T2FS Bl of the l -th regression parameter by
means of the T2MF  Bl (r,u)   MBl (r) (u), u  Jr and  Bl (r,u)  0, u  Jr . This conclusion allows us
to introduce the following notion.</p>
      <p>Definition 2. By the regression parameter with index l  0,..., p for the FS I of sample indices is
meant the T2FS</p>
      <p>Bl  {((r,u), Bl (r,u)) : u [0,1], r  } , 
(25) 
with the T2MF</p>
      <p>Here,</p>
      <p>max{  A : u  l (I ( )) (r)}, u  Jr ;  
 Bl (r,u)   0, u  Jr .</p>
      <p>Jr  {u [0,1] : u  l (I ( )) (r),   A}, 
4. Regression for a fuzzy set of sample indices 
4.1. Decomposition of T2FSs of regression parameters 
is the set of primary membership degrees u [0,1] with strictly positive secondary grades  Bl (r,u)
which coincides with the support supp(M Bl (r)) (see (23)) of the T1FS MBl (r) of fuzzy membership
degrees of the element r </p>
      <p>;
is the MF of the l -th fuzzy parameter
l (I ( )) (r)  max{u [0,1] : [l (I ( ))]uD  r  [l (I ( ))]uH } , </p>
      <p>l (I ( ))  {(r,l (I ( )) (r)) : r  } , 
for the set I ( ) of sample indices (see (19)-(20) for K  I ( ) );
fl ([ X (I ( ))]uH ,[ y(I ( ))]uD ), fl ([ X (I ( ))]uH ,[ y(I ( ))]uH )}
[l (I ( ))]uD  min{ fl ([ X (I ( ))]uD ,[ y(I ( ))]uD ), fl ([ X (I ( ))]uD ,[ y(I ( ))]uH ,  
and
[l (I ( ))]uH  max{ fl ([ X (I ( ))]uD ,[ y(I ( ))]uD ), fl ([ X (I ( ))]uD ,[ y(I ( ))]uH ,  
fl ([ X (I ( ))]uH ,[ y(I ( ))]uD ), fl ([ X (I ( ))]uH ,[ y(I ( ))]uH )}
are the estimates of the lower and upper bounds (see (15)-(16) for K  I ( ) ) of the u -cut
[l (I ( ))]u of the FN l (I ( )) ;</p>
      <p>I ( )  { j  N :  I ( j)  }   (32) 
is the  -cut of the FS I of sample indices;</p>
      <p>A is the set of the membership degrees values  I ( j), j  N of the FS I  {( j, I ( j)) : j  N} of
sample indices (see Remark 6). According to (26), the set A includes all possible positive values of
secondary grades for the T2FS Bl of the l -th regression parameter.
(26) 
(27) 
(28) 
(29) 
(30) 
(31) 
(33) </p>
      <p>For each l -th regression parameter, l  0,..., p , we apply a decomposition approach to represent
the T2FS Bl in a more convenient form. Theorem 1 justifies the representation of the T2FS Bl in the
form of a collection of the embedded T2FSs with constant secondary grades.</p>
      <p>Theorem 1. The T2FS Bl of the l -th regression parameter is represented in the form of the
collection Bl  {Ble2 (I ( )) :   A} of embedded T2FSs</p>
      <p>Ble2 (I ( ))  {(l (I ( )), ))} ,  
with the constant secondary grades   A . For each
  A , the embedded</p>
      <p>T1FS
l (I ( ))  {(r,l (I ( )) (r)) : r  } is the FN which is the l -th fuzzy parameter for the set
I ( )  { j  N :  I ( j)  } of sample indices, with the MF l (I ( )) (r) in form (28).</p>
      <p>Proof. According to (25), the T2FS of the regression parameter with the index l  0,..., p has the
form Bl  {((r,u), Bl (r,u)) : u [0,1], r  } . Hence,</p>
      <p>Bl  {{((r,u), max{  A : u  l (I ( )) (r)}) : u  Jr }  {((r,u), 0) : u  Jr }: r  } ,  
by (26). Remark 4 allows us to ignore the pairs (r,u) that have secondary grades equal to 0. Thus,
(34) </p>
      <p>Bl  {((r,u), max{  A : u  l (I ( )) (r)}) : u  J x , r  } ,  
and thereupon</p>
      <p>Bl  {((r,l (I ( )) (r)) :   A) : r  }
by (27). Note that the collection
{(l (I ( )) (r), ) :   A} is the T1FS which is formed by the unique value u  l (I ( )) (r) of the fuzzy
degree of membership of r  . The different values   A may correspond to u  l (I ( )) (r) .
Therefore, the equality {(l (I ( )) (r), ) :  A}  {(u, max{ : u  l (I ( )) (r)})} holds true. Further,
A
regrouping elements yields
  Bl  {(r, (l (I ( )) (r), )) :  A, r  }  {{((r,l (I ( )) (r)), ) : r  }:  A} . 
(35) 
(36) 
Finally, by virtue of formula (29) we conclude that Bl  {(l (I ( )), ) :   A} .</p>
      <p>The proof of Theorem 1 is complete.
4.2.</p>
    </sec>
    <sec id="sec-7">
      <title>Calculation of T2FSs of regression parameters </title>
      <p>First, we construct the set A  { I ( j) : j  N} of membership degrees values of the FS
I  {( j, I ( j)) : j  N} of sample indices. For each   A , according to (32), we construct the  -cut
I ( )  { j  N :  I ( j)  } of the FS I . Further, for each l  0,..., p , we use the representation of the
T2FS Bl in the form of a collection of embedded T2FSs with constant secondary grades (see
Theorem 1). This leads to the following sequence of calculations for each   A .</p>
      <p>We construct the embedded T1FS l (I ( ))  {(r,l (I ( )) (r)) : r  } . This is the FN, which is the
l -th fuzzy parameter of the regression for the set I ( ) of sample indices. To construct the FN
l (I ( )) one can use any known methods. An application of the method worked out in Section 2.1 of
[18] yields l (I ( ))  {([[l (I ( ))]uD ,[l (I ( ))]uH ],u) : u [0,1]} (see (13), (30), (31) for K  I ( ) )
with the MF
l (I ( )) (rl )  max{u [0,1] : [l (I ( ))]uD  rl  [l (I ( ))]uH } ,  rl 
. 
(37) 
Formula (37) is justified by (14) with K  I ( ) . Then, the corresponding embedded T2FS with the
constant secondary grade  has the form Ble2 (I ( ))  {(l (I ( )), ))} according to (33).</p>
      <p>Once all embedded T2FSs Ble2 (I ( )) with constant secondary grades   A have been obtained,
the
resulting</p>
      <p>T2FS
of
the
l -th
regression
parameter
has
the
form
Bl  {Ble2 (I ( )) :   A}  {(l (I ( )), )) :   A} by Theorem 1. According to Remark 3, the T2FSs
Bl can be interpreted as follows. For fuzzy data of different quality, the l -th regression parameter Bl
is equal to the l -th parameter l (I ( )) of the regression for the corresponding crisp set I ( ) of data
samples with the degree of truth being equal to  ,   A .</p>
    </sec>
    <sec id="sec-8">
      <title>5. Illustration and discussion </title>
    </sec>
    <sec id="sec-9">
      <title>5.1. Example </title>
      <p>This example is devised to illustrate our approach. We stress that this example does not use real
data. All data are given to the second decimal place.</p>
      <p>To conduct the historical research, we need to find out the dependence of weight of a middle-aged
person on he or her tall in the 10th century. The data borrowed from four authentic historical
documents are located in lines 1-4 of Table 1. There is a reason to believe that the values in line 3 are
slightly less reliable than the rest. We assume that in those days the height and the weight were
measured with an accuracy to ±1Kg and ±5%, respectively. We denote by a  (aL , aC , aU ) a
‘triangular’ FN with the MF
(r  a L ) / (aC  a L ), r [a L  aC ];

 a (r)  (a R  r) / (a R  aC ), r [aC  a R ];  
0, otherwise.
(38) 
yi1  ( yiL , yiC , yiU ) , i {1,..., 4} . The fuzzy data xi1 and yi be denoted as xi1  appr( xiC1 ) and
yi  appr( yiC ) , and can be interpreted as ‘approximately xiC1 and yiC ’, respectively. The dependence
of the weight on the height for majority of data in rows 1-4 of Table 1 may seem strange, since weight
gain decreases as the height increases. We assume that this can be explained by the small sample size
and its non-representativeness.</p>
      <p>Table 1 
Weight depending on height 
i 
1 
2 
3 
4 
5 
6 
7 
8 
 
 
 
 
 
 
 
 
 </p>
      <p>
        Height  xi1  
 
To solve this problem, we supplement the data sample. We know from the medical sources [31] that
the Broca’s formula for determining the height and the optimal weight has the form W = H − 100,
where H is the height in cm and W is the optimal weight in Kg. According to this formula, we
calculate the optimal weights for four different heights and place these data in rows 5-8 of Table 1.
Thus, we have three sources of data of different quality. The high-quality source is historical, which
we fully trust. We evaluate the degree of reliability of data from this source. The results which are
located in rows 1,2,4 of Table 1 have the values of the degree of reliability are equal to 1. The degree
of less reliable data from the second source of historical data which are located in line 3 of Table 1 is
estimated at 0,9. The third source of data is medical. It is known that modern persons are taller than
persons living in the 10th century and having the same weight. As a consequence, the degree of
reliability of data from a medical source is not very high (for example, we estimate it as 0,7), since
this source only provides information on the ratio of the weight and the height for modern people.
Thus, we construct the FS ‘Reliable data samples’
I  {(1;1), (2;1), (3;0,9), (4;1), (5;0, 7), (6;0, 7), (7;0, 7), (8;0, 7)} on the set N  {1,...,8} of data sample
indexes with the MF values  I (
        <xref ref-type="bibr" rid="ref1">1</xref>
        )   I (
        <xref ref-type="bibr" rid="ref2">2</xref>
        )   I (
        <xref ref-type="bibr" rid="ref4">4</xref>
        )  1 ,  I (
        <xref ref-type="bibr" rid="ref3">3</xref>
        )  0,9 и  I (
        <xref ref-type="bibr" rid="ref5">5</xref>
        )  ...   I (8)  0, 7 .
      </p>
      <p>
        In view of Remark 6, the set of membership degrees values of the FS I takes the form
A  {0, 7;0,9;1}. For   1, 0 ;   0, 9 and   0, 7 , according to (32), we construct the
corresponding α-cuts I (
        <xref ref-type="bibr" rid="ref1">1</xref>
        )  {1, 2, 4} , I (0,9)  {1,..., 4} and I (
        <xref ref-type="bibr" rid="ref7">0, 7</xref>
        )  {1,...,8} of the FS I of data
samples indices. Next, we construct FNs of parameters of linear least-squares regressions (see Section
2.1). In Example 1, these FNs are of the ‘triangular’ type (in general case, this is not necessary). Here,
these FNs are of the ‘triangular’ type (in general case, this is not necessary).
      </p>
      <p>
        For   1, we get embedded T1FSs 0 (I (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ))  (86, 51; 82, 2; 79, 6)  appr(82, 2) and
1(I (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ))  (0,88; 0, 91; 0, 94)  appr(0, 91)
in
the
      </p>
      <p>T2FSs</p>
      <p>
        B0e2 (I (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ))  {(0 (I (
        <xref ref-type="bibr" rid="ref1">1</xref>
        )),1))}
and
B1e2 (I (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ))  {(1(I (
        <xref ref-type="bibr" rid="ref1">1</xref>
        )),1))}, respectively, with the constant secondary grade (the degree of truth) being
equal to 1.
      </p>
      <p>For   0, 9 , we get embedded T1FSs</p>
      <p>0 (I (0, 9))  (68, 47; 64, 4; 60, 41)  appr(64, 4) and
1(I (0, 9))  (0, 77; 0,81; 0,85)  appr(0,81) in the T2FSs
B0e2 (I (0,9))  {(0 (I (0, 9)); 0, 9))}
and
B1e2 (I (0,9))  {(1(I (0,9));0,9))} , respectively, with the degree of truth being equal to 0,9.</p>
      <p>
        For   0, 7 , we get embedded T1FSs 0 (I (
        <xref ref-type="bibr" rid="ref7">0, 7</xref>
        ))  (81, 79; 77; 72,3)  appr(77) and
1(I (
        <xref ref-type="bibr" rid="ref7">0, 7</xref>
        ))  (0,85;0,9;0,95)  appr(0,9) in the T2FSs
B0e2 (I (
        <xref ref-type="bibr" rid="ref7">0, 7</xref>
        ))  {(0 (I (
        <xref ref-type="bibr" rid="ref7">0, 7</xref>
        ));0, 7))} and
B1e2 (I (
        <xref ref-type="bibr" rid="ref7">0, 7</xref>
        ))  {(1(I (
        <xref ref-type="bibr" rid="ref7">0, 7</xref>
        ));0, 7))} , respectively, with the degree of truth) being equal to 0,7.
      </p>
      <p>According to Theorem 1, the resulting T2FSs of regression parameters have the forms
B0  {(appr(82, 2);1), (appr(64, 4); 0,9), (appr(77);0, 7)} , </p>
      <p>B1  {(appr(0,91);1), (appr(0,81); 0,9), (appr(0,9); 0, 7)}. 
(39) 
(40) </p>
    </sec>
    <sec id="sec-10">
      <title>The discussion of the results </title>
      <p>
        Let us consider a graphical interpretation obtained regression with T2FSs parameters. On Figure 3
thin solid lines demonstrate the graphs of the regression functions
y  (1(I (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ))0L  (0 (I (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ))0L x  86,51  0,88x , 
(41) 
y  (1(I (
        <xref ref-type="bibr" rid="ref1">1</xref>
        )))0H  (0 (I (
        <xref ref-type="bibr" rid="ref1">1</xref>
        )))0H x  79, 6  0,94x   (42) 
with parameters which correspond to the lower and upper bounds of the 0-cuts
[(0 (I (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ))0L , (0 (I (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ))0H ]  [86,51; 79, 6] and [(1(I (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ))0L , (1(I (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ))1H ]  [0,88;0,94] of the FNs
0 (I (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ))  appr(82, 2) and 1(I (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ))  appr(0,91) , respectively. These FNs are corresponding
(thin lines), and
(thin lines), and
(a thick line). For the crisp set I (
        <xref ref-type="bibr" rid="ref7">0, 7</xref>
        )  {1,...,8} of data sample indices, the dotted lines demonstrate the
graphs of the regression functions y  (1(I (
        <xref ref-type="bibr" rid="ref7">0, 7</xref>
        ))0L  (0 (I (
        <xref ref-type="bibr" rid="ref7">0, 7</xref>
        ))0L x  81, 79  0,85x and
(43) 
(44) 
(45) 
(46) 
y  (1(I (0,9))0H  (0 (I (0,9))0H x  60, 4  0,85x  
y  (1(I (0,9))1L  (0 (I (0,9))1L x  64, 4  0,81x  
y  (1(I (
        <xref ref-type="bibr" rid="ref7">0, 7</xref>
        ))0H  (0 (I (
        <xref ref-type="bibr" rid="ref7">0, 7</xref>
        ))0H x  72,3  0,95x  
      </p>
      <p>
        y  (1(I (
        <xref ref-type="bibr" rid="ref1">1</xref>
        )))0H  (0 (I (
        <xref ref-type="bibr" rid="ref1">1</xref>
        )))0H x  79, 6  0,94x  
Figure 3. The illustration of the regression with the parameter T2FSs 
parameters of the fuzzy linear regression for crisp set I (
        <xref ref-type="bibr" rid="ref1">1</xref>
        )  {1, 2, 4} of data sample indices. For the
same set of sample indices, the thick solid line shows the graph of the regression function
y  (1(I (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ))1L  (0 (I (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ))1L x  82, 2  0,91x
with parameters
which correspond to the 1-cuts
[82, 2; 82, 2] and [0,91;0,91] of the FNs 0 (I (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ))  appr(82, 2) and 1 (I (
        <xref ref-type="bibr" rid="ref1">1</xref>
        ))  appr(0, 91) ,
respectively. For the crisp set I (0,9)  {1,..., 4} of data sample indices, the dashed lines demonstrate the
graphs of the regression functions y  (1 (I (0,9))0L  (0 (I (0,9))0L x  68, 47  0, 77 x and
(a thick line).
      </p>
      <p>Figure 3 allows us to draw the following conclusions. Regressions corresponding to data of
different quality may differ from each other and express a different relationship between the
independent variables and predicted output. Therefore, the resulting regression with T2FSs of
parameters should be taken as a whole as a collection of regressions corresponding to  -cuts
I ( ),  [0,1] of the FS I of data sample indices. Only in this case we get an idea about
dependence between the regression and the quality of the data used. Sometimes understanding this is
important.</p>
      <p>If we need some ‘maximizing’ regression, which has parameters values with the primary
membership degrees equal to 1, then it should be considered as a collection of regressions
corresponding to  -cuts I ( ),  [0,1] of the FS I of data sample indices. For these regressions,
we take the parameters values with the membership degrees equal to 1. In Example 1, these are
regressions with graphs depicted by thick lines in Figure 3.</p>
      <p>The regression for fuzzy data of different quality can be represented as a set of fuzzy regressions
with the degree of truth   A . These fuzzy regressions correspond to data samples of different
quality with indices belonging to  -cuts I ( ),   A of the FS I . Thus, the degree of truth of the
regression is determined by the quality of the data for which it is constructed. The higher the quality
of data samples, the greater the degree of membership of their indices to the FS I and the higher the
degree of truth of the corresponding fuzzy regression in the collection that forms the regression for
fuzzy data of different quality. On the other hand, increasing the amount of low-quality data reduces
the degree of truth of the corresponding fuzzy regression in the collection.</p>
    </sec>
    <sec id="sec-11">
      <title>6. Conclusion </title>
      <p>According to the proposed approach, the linear regression parameters for fuzzy data of different
quality are T2FSs with constant secondary grades. Although, in general, a T2FS is a rather
complicated mathematical object, T2FSs with constant secondary grades are simple enough for
practical use. Therefore, to represent these sets in a form which is convenient for understanding and
applications, we have used a decomposition method. The results obtained have allowed us to
decompose the T2FSs of linear regression parameters according to secondary grades into finite sets of
FNs. These are fuzzy parameters of the regressions which correspond to different  -cuts of FSs of
data sample indices. One direction of future investigation of regressions with a fuzzy set of data
sample indices may be related to the development of a similar approach for possibilistic regression
analysis.</p>
    </sec>
    <sec id="sec-12">
      <title>7. Acknowledgements </title>
      <p>This work has been supported by Ministry of Education and Science of Ukraine: R&amp;D Project
0122U001844 for the period of 2022-2024 at Taras Shevchenko National University of Kyiv.</p>
    </sec>
    <sec id="sec-13">
      <title>8. References </title>
      <p>
        [8] J. Hevesi, A. Flint, J. Istok, Precipitation estimation in mountainous terrain using multivariate
geostatistics. Part II: isohyetal maps. J. Appl. Meteorol. 31 (1992) 677–688.
doi:10.1175/15200450(1992)031!0677:PEIMTUO2.0.CO;2.
[9] M. Kennedy, A. O’Hagan, Predicting the output from a complex computer code when fast
approximations are available, Biometrika 87(
        <xref ref-type="bibr" rid="ref1">1</xref>
        ) (2000) 1–13. doi:10.1093/biomet/87.1.1.
[10] S. Pan, Q. Yang, A survey on transfer learning. IEEE Trans. on Knowledge and Data
      </p>
      <p>Engineering 22(10) (2010) 1345–1359. doi:10.1109/TKDE.2009.191.
[11] J.W. Bandler, Q.S. Cheng, S.A. Dakroury, A.S. Mohamed, M.H. Bakr, K. Madsen, J.</p>
      <p>
        Sondergaard, Space mapping: the state of the art, IEEE Trans. on microwave theory and
techniques 52(
        <xref ref-type="bibr" rid="ref1">1</xref>
        ) (2004) 337–361. https://ieeexplore.ieee.org/document/ 1262727.
[12] C. E. Rasmussen, C. K. I. Williams, Gaussian processes for machine learning, The MIT Press,
2006.
[13] N. Arruda, J. Alcantara, V. Vidal, A. Brayner, M. Casanova, V. Pequeno, W.A. Franco, Fuzzy
Approach for Data Quality Assessment of Linked Datasets, in: Proceedings of the 21th
International Conference on Enterprise Information Syst. (ICEIS 2019), volume 1, 2019, pp.
399–406. doi:10.5220/0007718803990406.
[14] J. Gamez, F. Modave, O. Kosheleva, Selecting the Most Representative Sample is NP-Hard:
Need for Expert (Fuzzy) Knowledge, in: Proceedings of the IEEE International Conference on
Fuzzy Syst. (IEEE World Congress on Computational Intelligence), 2008, pp. 1069–1074.
doi:10.1109/FUZZY.2008. 4630502.
[15] M. Z. Khan, M. F. Khan, M. Aslam, A. R. Mughal, Design of Fuzzy Sampling Plan Using the
      </p>
      <p>
        Birnbaum-Saunders Distribution, Mathematics 7(
        <xref ref-type="bibr" rid="ref1">1</xref>
        ) (2019) 9. doi:10.3390/math7010009.
[16] S. Mashchenko, Intersections and unions of fuzzy sets of operands, Fuzzy Sets and Syst. 352(
        <xref ref-type="bibr" rid="ref1">1</xref>
        )
(2018) 12–25. doi:10.1016/j.fss.2018.04.006.
[17] S.O. Mashchenko, D.O. Kapustian, Decomposition of intersections with fuzzy sets of operands,
in: V.A. Sadovnichiy, M.Z. Zgurovsky (Ed.), Contemporary Approaches and Methods in
Fundamental Mathematics and Mechanics. Understanding Complex Systems, Springer, Cham,
2020, pp. 417–432. https://link.springer.com/book/10.1007/ 978-3-030-50302-4.
[18] H.-C. Wu, Linear regression analysis for fuzzy input and output data using the extension
principle, Computers and Mathematics with Applications 45(12) (2003) 1849–1859.
doi:10.1016/S0898-1221(03)90006-X.
[19] S. Heilpern, Representation and application of fuzzy numbers, Fuzzy Sets and Syst. 91 (1997)
259–268.
[20] L. Zadeh, The concept of a linguistic variable and its application to approximate reasoning – I,
      </p>
      <p>
        Inform. Sci. 8 (1975) 199–249. doi:10.1016/0020-0255(75)90036-5.
[21] D. Dubois, H. Prade, Operations on fuzzy numbers, International J. of Syst. Sci. 9(
        <xref ref-type="bibr" rid="ref6">6</xref>
        ) (1978)
613–626.
[22] L. A. Zadeh, Quantitative fuzzy semantics, Inform. Sci. 3 (1971) 159–176.
doi:10.1016/S00200255(71)80004-X.
[23] N. N. Karnik, J.M. Mendel, Introduction to type-2 fuzzy logic systems, IEEE International
      </p>
      <p>Conference on Fuzzy Syst. 2 (1998) 915–920.
[24] J. M. Mendel, R. I. John, Type-2 fuzzy sets made simple, IEEE Trans. on Fuzzy Syst. 10 (2002)
117–127.
[25] J. Harding, C. Walker, E. Walker, The variety generated by the truth value algebra of T2FSs,</p>
      <p>Fuzzy Sets and Syst. 161 (2010) 735–749.
[26] J. Aisbett, J.T. Rickard, D.G. Morgenthaler, Type-2 fuzzy sets as functions on spaces, IEEE</p>
      <p>
        Trans. on Fuzzy Syst. 18(
        <xref ref-type="bibr" rid="ref4">4</xref>
        ) (2010) 841–844.
[27] J. M. Mendel, Type-2 fuzzy sets: some questions and answers, IEEE Connections. Newsletter of
the IEEE Neural Networks Society 1 (2003) 10–13.
[28] S.O. Mashchenko, Sums of fuzzy set of summands, Fuzzy Sets and Syst. 417 (2020) 140–151.
      </p>
      <p>doi:10.1016/j.fss.2020.10.006.
[29] S.O. Mashchenko, Sum of discrete fuzzy numbers with fuzzy set of summands, Cybernetics and</p>
      <p>
        Syst. Analysis 57(
        <xref ref-type="bibr" rid="ref3">3</xref>
        ) (2021) 374–382. doi:10.1007/s10559-021-00362-w.
[30] S.O. Mashchenko, Minimum of fuzzy numbers with a fuzzy set of operands, Cybernetics and
      </p>
      <p>
        Syst. Analysis 58(
        <xref ref-type="bibr" rid="ref2">2</xref>
        ) (2022) 210–219. doi: 10.1007/s10559-021-00362-w.
[31] J. C. Segen, Concise dictionary of modern medicine, McGraw-Hill, New York, NY,2006.
      </p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <given-names>N.</given-names>
            <surname>Chukhrova</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Johannssen</surname>
          </string-name>
          ,
          <article-title>Fuzzy regression analysis: Systematic review and bibliography</article-title>
          ,
          <source>Appl. Soft Comput. J</source>
          .
          <volume>84</volume>
          (
          <year>2019</year>
          )
          <fpage>1</fpage>
          -
          <lpage>29</lpage>
          . doi:
          <volume>10</volume>
          .1016/j.asoc.
          <year>2019</year>
          .
          <volume>105708</volume>
          .
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <string-name>
            <given-names>H.</given-names>
            <surname>Tanaka</surname>
          </string-name>
          ,
          <string-name>
            <given-names>S.</given-names>
            <surname>Uejima</surname>
          </string-name>
          ,
          <string-name>
            <given-names>K.</given-names>
            <surname>Asai</surname>
          </string-name>
          ,
          <article-title>Linear regression analysis with Fuzzy model</article-title>
          ,
          <source>IEEE Trans. Syst. Man Cybern</source>
          .
          <volume>12</volume>
          (
          <issue>6</issue>
          ) (
          <year>1982</year>
          )
          <fpage>903</fpage>
          -
          <lpage>907</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <surname>P. D'Urso</surname>
          </string-name>
          ,
          <article-title>Exploratory multivariate analysis for empirical information affected by uncertainty and modeled in a fuzzy manner: a review</article-title>
          ,
          <source>Granular Comput</source>
          .
          <volume>2</volume>
          (
          <issue>4</issue>
          ) (
          <year>2017</year>
          )
          <fpage>225</fpage>
          -
          <lpage>247</lpage>
          . doi:
          <volume>10</volume>
          .1007/s41066-017-0040-y.
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [4]
          <string-name>
            <given-names>A.</given-names>
            <surname>Celmins</surname>
          </string-name>
          ,
          <article-title>Least squares model fitting to fuzzy vector data</article-title>
          ,
          <source>Fuzzy Sets and Syst</source>
          .
          <volume>22</volume>
          (
          <issue>3</issue>
          ) (
          <year>1987</year>
          )
          <fpage>245</fpage>
          -
          <lpage>269</lpage>
          . doi:
          <volume>10</volume>
          .1016/
          <fpage>0165</fpage>
          -
          <lpage>0114</lpage>
          (
          <issue>87</issue>
          )
          <fpage>90070</fpage>
          -
          <lpage>4</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [5]
          <string-name>
            <given-names>T.</given-names>
            <surname>Hastie</surname>
          </string-name>
          ,
          <string-name>
            <given-names>R.</given-names>
            <surname>Tibshirani</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J.</given-names>
            <surname>Friedman</surname>
          </string-name>
          ,
          <source>The Elements of Statistical Learning: Data Mining, Inference, and Prediction</source>
          . Springer, New York, NY,
          <year>2009</year>
          . doi:
          <volume>10</volume>
          .1007/978-0-
          <fpage>387</fpage>
          -84858-7.
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          [6]
          <string-name>
            <given-names>A.</given-names>
            <surname>Forrester</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Sobester</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Keane</surname>
          </string-name>
          <article-title>, Multi-fidelity optimization via surrogate modelling</article-title>
          ,
          <source>Proceedings of the Royal society A: Mathematical, Physical and Engineering Sci</source>
          .
          <volume>463</volume>
          (
          <year>2007</year>
          )
          <fpage>3251</fpage>
          -
          <lpage>3269</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          [7]
          <string-name>
            <given-names>A.</given-names>
            <surname>Zaytsev</surname>
          </string-name>
          , E. Burnaev,
          <article-title>Minimax Approach to Variable Fidelity Data Interpolation</article-title>
          ,
          <source>in: Proceedings of the 20th International Conference on Artificial Intelligence and Statistics</source>
          , PMLR, volume
          <volume>54</volume>
          ,
          <string-name>
            <surname>Fort</surname>
            <given-names>Lauderdale</given-names>
          </string-name>
          , FL, USA,
          <year>2017</year>
          , pp.
          <fpage>652</fpage>
          -
          <lpage>661</lpage>
          . https://proceedings.mlr.press/v54/zaytsev 17a.
          <source>html.</source>
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>