Soft Methods for Data Science 1st Edition Maria Brigida Ferraro pdf download
Soft Methods for Data Science 1st Edition Maria Brigida Ferraro pdf download
https://textbookfull.com/product/soft-methods-for-data-
science-1st-edition-maria-brigida-ferraro/
https://textbookfull.com/product/java-for-data-science-1st-
edition-reese/
https://textbookfull.com/product/r-for-data-science-1st-edition-
garrett-grolemund/
https://textbookfull.com/product/algorithms-for-data-science-1st-
edition-brian-steele/
https://textbookfull.com/product/data-science-for-
mathematicians-1st-edition-nathan-carter/
Mathematical Problems in Data Science Theoretical and
Practical Methods 1st Edition Li M. Chen
https://textbookfull.com/product/mathematical-problems-in-data-
science-theoretical-and-practical-methods-1st-edition-li-m-chen/
https://textbookfull.com/product/r-programming-for-data-
science-1st-edition-roger-peng/
https://textbookfull.com/product/programming-skills-for-data-
science-1st-edition-michael-freeman/
https://textbookfull.com/product/cultural-anthropology-an-
applied-perspective-gary-ferraro/
https://textbookfull.com/product/data-science-and-social-
research-ii-methods-technologies-and-applications-paolo-mariani/
Advances in Intelligent Systems and Computing 456
Maria Brigida Ferraro
Paolo Giordani
Barbara Vantaggi
Marek Gagolewski
María Ángeles Gil
Przemysław Grzegorzewski
Olgierd Hryniewicz Editors
Soft Methods
for Data
Science
Advances in Intelligent Systems and Computing
Volume 456
Series editor
Janusz Kacprzyk, Polish Academy of Sciences, Warsaw, Poland
e-mail: kacprzyk@ibspan.waw.pl
About this Series
The series “Advances in Intelligent Systems and Computing” contains publications on theory,
applications, and design methods of Intelligent Systems and Intelligent Computing. Virtually
all disciplines such as engineering, natural sciences, computer and information science, ICT,
economics, business, e-commerce, environment, healthcare, life science are covered. The list
of topics spans all the areas of modern intelligent systems and computing.
The publications within “Advances in Intelligent Systems and Computing” are primarily
textbooks and proceedings of important conferences, symposia and congresses. They cover
significant recent developments in the field, both of a foundational and applicable character.
An important characteristic feature of the series is the short publication time and world-wide
distribution. This permits a rapid and broad dissemination of research results.
Advisory Board
Chairman
Nikhil R. Pal, Indian Statistical Institute, Kolkata, India
e-mail: nikhil@isical.ac.in
Members
Rafael Bello, Universidad Central “Marta Abreu” de Las Villas, Santa Clara, Cuba
e-mail: rbellop@uclv.edu.cu
Emilio S. Corchado, University of Salamanca, Salamanca, Spain
e-mail: escorchado@usal.es
Hani Hagras, University of Essex, Colchester, UK
e-mail: hani@essex.ac.uk
László T. Kóczy, Széchenyi István University, Győr, Hungary
e-mail: koczy@sze.hu
Vladik Kreinovich, University of Texas at El Paso, El Paso, USA
e-mail: vladik@utep.edu
Chin-Teng Lin, National Chiao Tung University, Hsinchu, Taiwan
e-mail: ctlin@mail.nctu.edu.tw
Jie Lu, University of Technology, Sydney, Australia
e-mail: Jie.Lu@uts.edu.au
Patricia Melin, Tijuana Institute of Technology, Tijuana, Mexico
e-mail: epmelin@hafsamx.org
Nadia Nedjah, State University of Rio de Janeiro, Rio de Janeiro, Brazil
e-mail: nadia@eng.uerj.br
Ngoc Thanh Nguyen, Wroclaw University of Technology, Wroclaw, Poland
e-mail: Ngoc-Thanh.Nguyen@pwr.edu.pl
Jun Wang, The Chinese University of Hong Kong, Shatin, Hong Kong
e-mail: jwang@mae.cuhk.edu.hk
Olgierd Hryniewicz
Editors
123
Editors
Maria Brigida Ferraro María Ángeles Gil
Department of Statistical Sciences Department of Statistics and Operational
Sapienza University of Rome Research and Mathematics Didactics
Rome University of Oviedo
Italy Oviedo
Spain
Paolo Giordani
Department of Statistical Sciences Przemysław Grzegorzewski
Sapienza University of Rome Department of Stochastic Methods, Systems
Rome Research Institute
Italy Polish Academy of Sciences
Warsaw
Barbara Vantaggi Poland
Department of Basic and Applied Sciences
for Engineering Olgierd Hryniewicz
Sapienza University of Rome Department of Stochastic Methods, Systems
Rome Research Institute
Italy Polish Academy of Sciences
Warsaw
Marek Gagolewski Poland
Department of Stochastic Methods, Systems
Research Institute
Polish Academy of Sciences
Warsaw
Poland
v
vi Preface
General Chairs
Maria Brigida Ferraro, Rome, Italy
Paolo Giordani, Rome, Italy
Barbara Vantaggi, Rome, Italy
Executive Board
María Ángeles Gil, Oviedo, Spain
Przemysław Grzegorzewski, Warsaw, Poland
Olgierd Hryniewicz, Warsaw, Poland
Program Committee
Bernard de Baets, Ghent, Belgium
Christian Borgelt, Mieres, Spain
Giulianella Coletti, Perugia, Italy
Ana Colubi, Oviedo, Spain
Ines Couso, Oviedo, Spain
Gert De Cooman, Ghent, Belgium
Thierry Denoeux, Compiégne, France
Didier Dubois, Toulouse, France
Fabrizio Durante, Bolzano, Italy
Pierpaolo D’Urso, Rome, Italy
Gisella Facchinetti, Lecce, Italy
Peter Filzmoser, Vienna, Austria
Lluís Godo, Barcelona, Spain
Gil González Rodríguez, Oviedo, Spain
Mario Guarracino, Naples, Italy
vii
viii Organization
Publication Chair
Marek Gagolewski, Warsaw, Poland
Additional Reviewers
Serena Arima, Salem Benferhat, Patrizia Berti, Angela Blanco-Fernández, Christian
Borgelt, Andrea Capotorti, Marco Cattaneo, Jasper de Bock, Alessandro De
Gregorio, Sara De la Rosa de Sáa, Marco Di Zio, Marta Disegna, Scott Ferson, Gianna
Figá-Talamanca, Enrico Foscolo, Marek Gagolewski, Luis Angel García-Escudero,
José Luis García-Lapresta, Angelo Gilio, Rita Giuliano, Stefania Gubbiotti,
Maria Letizia Guerra, Radim Jiroušek, Tomas Kroupa, Brunero Liseo, Massimo
Marinacci, Agustin Mayo-Iscar, Enrique Miranda, Piotr Nowak, David Over,
Sonia Pérez-Fernández, Davide Petturiti, Niki Pfeifer, Ana B. Ramos-Guajardo,
Organization ix
Pietro Rigo, Giuseppe Sanfilippo, Fabian Schmidt, Rudolf Seising, Beatriz Sinova,
Damjan Skulj, Yann Soullard, Andrea Tancredi, Luca Tardella, Pedro Terán,
Matthias Troffaes, Tiziana Tuoto, Arthur Van Camp, Paolo Vicig.
Local Organization
Marco Brandi
Massimiliano Coppola—conference website chair
Francesco Dotto
Liana Francescangeli—administrative staff
Tullia Padellini
Alessandra Pelorosso—administrative chair
Marco Stefanucci
Contents
xi
xii Contents
Abstract We propose a definition of mean value and variance for fuzzy numbers
whose membership functions are upper-semicontinuous but are not necessarily con-
tinuous. Our proposal uses the total variation of bounded variation functions.
1 Introduction
In this work we face the problem of defining the mean value and variance of fuzzy
numbers whose membership functions are upper-semicontinuous but are not neces-
sarily continuous. The literature presents a high number of definitions in the con-
tinuous case [3–6, 9] also because the average has often been counted among the
ranking modes of fuzzy numbers. The starting point is that a fuzzy number is defined
by a membership function that is of bounded variation. Even in this more general
context, it is possible to introduce a weighted average by means of a classical vari-
ational formulation and by α-cuts, as well as many authors do in the continuous
case. In the non-continuous case, we introduce the lower and upper weighted mean
values, but the generality of the weights doesn’t let to obtain the weighted mean
value as the middle point of the previous ones. This property will be obtained either
in the continuous case or for particular weight functions. An interesting property of
this new version is connected to the view of the weighted mean value in a possi-
bilistic framework, as in the continuous case Carlsson and Fullér [3] and Fullér and
Majlender [5] do. This property is interesting and harbinger of future developments
also in the non-continuous case. Following the same line, we pass to introduce the
concept of variance and we suggest two different definitions as happens in the case
of continuous membership functions.
In this section we recall some basic properties of functions of bounded variation. For
more details see, e.g., [7].
Let I ⊂ R be an interval and u : I → R be a function. The total variation of u on
n
I is defined as V I [u] = sup |u(xi ) − u(xi−1 )|, where the supremum is taken over
i=1
all finite partitions {x0 , x1 , . . . , xn } ⊂ I with x0 < x1 < · · · < xn . We say that u is a
bounded variation function on I if V I [u] < +∞. We denote BV (I ) the space of all
bounded variation functions on I . The following properties hold: if u is monotone
and bounded on I then u ∈ BV (I ) and V I [u] = sup I u − inf I u; if u 1 , u 2 ∈ BV (I )
and k is a constant then ku 1 , u 1 + u 2 , u 1 − u 2 , u 1 u 2 belong to BV (I ).
Let u : R → R be a bounded variation function. The total variation function of u is
the increasing function vu defined by vu (x) = V−∞ x
[u], where V−∞
x
[u] denotes the
total variation of u on ] − ∞, x].
The total variation measure |Du| of u is defined as the Lebesgue-Stieltjes measure
dvu associated to vu . The positive and negative variations of u are defined, respec-
tively, by the increasing functions
In this section we propose a definition of f -weighted mean value for a fuzzy number
whose membership function is upper-semicontinuous but not necessarily continuous.
These properties produce that the fuzzy number has membership function of bounded
variation and so we realize to introduce its weighted mean value by its total variation
measure.
A fuzzy number A is a fuzzy set of R with a normal, (fuzzy) convex and upper-
semicontinuous membership function μ : R → [0, 1] of bounded support [2]. From
Mean Value and Variance of Fuzzy Numbers … 3
the previous definition there exist two functions μ L , μ R : R → [0, 1], where μ L is
nondecreasing and right-continuous and μ R is nonincreasing and left-continuous,
such that ⎧
⎪
⎪ 0 x < a1 or x > a4
⎨
μ L (x) a1 ≤ x < a2
μ(x) =
⎪
⎪ 1 a2 ≤ x ≤ a3
⎩
μ R (x) a3 < x ≤ a4
In order to define the lower and upper mean values of A, we consider the positive
and negative variations of μ, as defined in (1), given by, respectively,
⎧ ⎧
⎪
⎨0 x < a1 ⎪
⎨0 x ≤ a3
+ −
μ (x) = μ L (x) a1 ≤ x < a2 μ (x) = 1 − μ R (x) a3 < x ≤ a4 (4)
⎪
⎩ ⎪
⎩
1 x ≥ a2 , 1 x > a4 .
Following the line present in the classical case of continuous fuzzy numbers, we face
the problem to give an expression of lower and upper f -weighted mean values using
α-cuts. This different version lets the possibility to write the previous definitions by
the left and right extreme points of fuzzy number α-cuts, exactly as happens in the
continuous case.
Proposition 2
1 1
a L (α) f (μ(a L (α))) dα a R (α) f (μ(a R (α))) dα
E ∗ (A; f ) = 0 1 , E ∗ (A; f ) = 0 1 . (6)
0 f (μ(a L (α))) dα 0 f (μ(a R (α))) dα
Since a L (α) = min{x ≥ a1 ; μ L (x) ≥ α}, α ∈ [0, 1], applying the change of vari-
+∞ a2
able formula [7, Theorem 5.42] we obtain g(x) f (μ(x)) dμ+ (x) = g(x)
−∞
1 a1
f (μ(x)) dμ L (x) = 0 g(a L (α)) f (μ(a L (α))) dα and thus (7). We now prove (8).
+∞ a3 −a3
We have g(x) f (μ(x)) dμ− (x) = g(x) f (μ(x)) dμ R (x) = g(−y)
−∞ −a3 a4 −a4
f (μ(−y)) dμ R (−y) = −a4 g(−y) f (μ(−y)) dk(y)
where we have used the change of variable y = −x and let k(x) = μ R (−x).
Since a R (α) = max{x; x ≤ a4 , μ R (x) ≥ α}, α ∈ [0, 1], letting h(α) = −a R (α) =
min{y; y ≥ −a4 , k(y) ≥ α} and using the change of variable formula [7, Theorem
5.42] (observing that k is increasing), we have, continuing the above chain of equal-
ities,
k(−a3 ) μ R (a3 )
= g(−h(α)) f (μ(−h(α))) dα = g(a R (α)) f (μ(a R (α))) dα =
k(−a4 ) μ R (a4 )
1
0 g(a R (α)) f (μ(a R (α))) dα. Then (8) is proved. Substituting (7) and (8) in (5)
with g(x) = x and g(x) = 1 we obtain for E ∗ (A; f ) and E ∗ (A; f ) the expressions
given in (6).
We now consider the special case f (α) = 1. We denote by E(A), E ∗ (A) and E ∗ (A),
respectively, the mean value and lower and upper mean values of a fuzzy number A
computed using f (α) = 1. This particular case is interesting as we may show that
a sufficient condition so that the weight w is equal to 1/2 is that f (α) = 1. This
condition is not necessary as we will show later in an example.
+∞ 1
Proposition 3 We have E(A) = 21 −∞ x |Dμ| = 21 0 (a L (α) + a R (α)) dα,
+∞ 1 +∞ 1
E ∗ (A) = −∞ x dμ+ (x) = 0 a L (α) dα and E ∗ (A) = −∞ x dμ− (x) = 0 a R (α)
dα. In particular E(A) = (E ∗ (A) + E ∗ (A))/2.
+∞ +∞ +
Proof Since μ+ (x) is increasing, we have (see Sect. 2) −∞ dμ+ (x) = V−∞ [μ ] =
+ +
+∞
sup μ − inf μ = 1, where last equality follows from (4). Similarly we get −∞
+∞ +∞ +∞
dμ− (x) = 1. Then, from (2), −∞ |Dμ| = −∞ dμ+ (x) + −∞ dμ− (x) = 2. Sub-
stituting in (3), with f = 1, we have the first equality of the thesis for E(A). The
second equality follows observing that w, as defined in Proposition 1, is equal to 1/2
and using (6) in Proposition 1 with f (α) = 1.
3.3 Example
We now present an example to show that, when f (α) = 1, we may have w = 1/2.
We compute the mean values of two fuzzy numbers shown in Fig. 1 for f (α) = 2α.
+∞
First, we consider the fuzzy number (a). We have −∞ x f (μ(x)) dμ+ (x) = 25/12,
+∞ +
+∞
−∞ f (μ(x)) dμ (x) = 5/4 and thus E ∗ (A; f ) = 5/3. Furthermore, −∞ x
+∞
f (μ(x)) dμ− (x) = 59/16, −∞ f (μ(x)) dμ− (x) = 17/16 and then E ∗ (A; f ) =
59/17. Moreover E(A; f ) = 277/111 ≈ 2.5 and w = 17/37 +∞ = 1/2.
We now consider the fuzzy number (b). We have −∞ x f (μ(x)) dμ+ (x) =
+∞ +∞
25/12, −∞ f (μ(x)) dμ+ (x) = 5/4 and thus E ∗ (A; f ) = 5/3. Furthermore, −∞
+∞
x f (μ(x)) dμ− (x) = 55/12, −∞ f (μ(x)) dμ− (x) = 5/4 and then E ∗ (A; f ) =
11/3. Moreover E(A; f ) = 8/3 ≈ 2.7 and w = 1/2.
We observe that our proposal, which starts from a variational measure, may be viewed
in a possibilistic framework. Indeed, since Poss(A ≤ a L (α)) = supx≤aL (α) μ(x) =
μ(a L (α)) and Poss(A ≥ a R (α)) = supx≥a R (α) μ(x) = μ(a R (α)), from (6) we get
1 1
E ∗ (A; f ) = 0 a L (α) f (Poss(A ≤ a L (α))) dα 0 f (Poss(A ≤ a L (α))) dα ,
6 L. Anzilli and G. Facchinetti
(a) (b)
1 1
1/2 1/2
1/4
1 1 3 4
3 4
1 1
E ∗ (A; f ) = 0 a R (α) f (Poss(A ≥ a R (α))) dα 0 f (Poss(A ≥ a R (α))) dα.
Note that A does not need to be continuous. In the special case when A is a continuous
fuzzy number we retrieve the lower and upper possibilistic mean values proposed
in [3] and [5] (see Proposition 4). Thus our approach offers an extension of possi-
bilistic mean values to the case of fuzzy numbers whose membership functions are
upper-semicontinuous but are not necessarily continuous.
If we use the variational approach for continuous fuzzy numbers we obtain that the
weighted mean value is the simple average of upper and lower f -weighted mean
values. This means that for every weighting function f the weight w = 1/2.
Proof The assertions easily follow from previous results observing that, since μ
is continuous, we have μ(a L (α)) = α and μ(a R (α)) = α. In particular, we have
+∞ 1 1
f (μ(x)) dμ+ (x) = 0 f (μ(a L (α))) dα = 0 f (α) dα = 1 and, similarly,
−∞
+∞ −
−∞ f (μ(x)) dμ (x) = 1.
Mean Value and Variance of Fuzzy Numbers … 7
Remark 1 If A is a continuous fuzzy number and f (α) = 2α we obtain that V ar1 (A)
matches the variance proposed by Zhang and Nie [9].
5 Conclusion
In this paper we have presented the weighted mean value and variance definitions for
fuzzy numbers with upper semicontinuous membership functions that are of bounded
variation. In the same context in a previous paper we have looked for a definition of
evaluation and ranking for fuzzy numbers [1]. The idea to use the space of bounded
variation functions has allowed us to view the mean value and the variance of fuzzy
numbers either in a variational context or in a classical way by α-cuts. This choice,
thanks to the freedom that we have left to the weights, has shown interesting results
that are able to generalize the previous ones but even to recover others present in
literature for the continuous case. We have intention to continue in this direction to
see if the space of bounded variation functions offers a new research field for fuzzy
sets.
References
1. Anzilli L, Facchinetti G (2013) The total variation of bounded variation functions to evaluate
and rank fuzzy quantities. Int J Intell Syst 28(10):927–956
2. Bede B (2013) Mathematics of fuzzy sets and fuzzy logic. Springer
3. Carlsson C, Fullér R (2001) On possibilistic mean value and variance of fuzzy numbers. Fuzzy
Sets Syst 122(2):315–326
4. Dubois D, Prade H (1987) The mean value of a fuzzy number. Fuzzy Sets Syst 24(3):279–300
5. Fullér R, Majlender P (2003) On weighted possibilistic mean and variance of fuzzy numbers.
Fuzzy Sets Syst 136(3):363–374
6. Georgescu I (2009) Possibilistic risk aversion. Fuzzy Sets Syst 160(18):2608–2619
7. Leoni G (2009) A first course in Sobolev spaces, vol 105. American Mathematical Society
Providence, RI
8. Li X, Guo S, Yu L (2015) Skewness of fuzzy numbers and its applications in portfolio selection.
Fuzzy Syst IEEE Trans 23(6):2135–2143
9. Zhang WG, Nie ZK (2003) On possibilistic variance of fuzzy numbers. In: Rough sets, fuzzy
sets, data mining, and granular computing.Springer, pp 398–402
On the Construction of Radially Symmetric
Trivariate Copulas
1 Introduction
n
n
Ĉn (x1 , . . . , xn ) = x j − (n − 1) + Cn (1, . . . , 1 − xi , . . . , 1 − x j , . . . , 1)
j=1 i< j
n
− Cn (1, . . . , 1 − xi , . . . , 1 − x j , . . . 1 − xk , . . . , 1) + . . .
i< j<k
+(−1) Cn (1 − x1 , 1 − x2 , . . . , 1 − xn ) .
n
(1)
C2 (x, y) + (1 − C2 (1 − x, 1 − y)) = x + y .
3 The Construction
In [8], the authors study several transformations of bivariate copulas. They show that
every radially symmetric 2-copula has the following form
1
SC2 ,D2 (x, y, z) = [1 − x − y − z + C2 (x, y) + C2 (x, z) + C2 (y, z)]
2
1
+ [HC2 ,D2 (x, y, z) − HC2 ,D2 (1 − x, 1 − y, 1 − z)] , (3)
2
where
HC2 ,D2 (x, y, z) = D2 (x, C2 (y, z)) + D2 (y, C2 (x, z)) + D2 (z, C2 (x, y))
2
− [D2 (x, D2 (y, z)) + D2 (y, D2 (x, z)) + D2 (z, D2 (x, y))] .
3
Note that the function HC2 ,D2 satisfies the boundary conditions of a 3-copula,
and from this is easy to prove that SC2 ,D2 also satisfies the boundary conditions of a
copula, and that the bivariate marginals of SC2 ,D2 are all equal to C2 .
On the Construction of Radially Symmetric Trivariate Copulas 13
Proof It can be shown after some tedious computations that we can rewrite SC2 ,D2
as
i.e., from the definition of survival copula given by Eq. (1), it follows that if SC2 ,D2
is a 3-copula, then SC2 ,D2 is a radially symmetric 3-copula because it coincides with
its associated survival copula.
However, SC2 ,D2 is not necessarily a 3-copula, as it may even not be an increas-
ing function. For example, if C2 = D2 = F (−2) , we can easily see that S F (−2) ,F (−2)
( 21 , 10
1
, 10
1
) < 0 = S F (−2) ,F (−2) (0, 10
1
, 10
1
). We now provide some examples where the
construction effectively yields a 3-copula.
Example 1 Consider the Frank family of 2-copulas. From the results in [11], we
know that for α ≥ − ln(2), the 3-dimensional version of the Frank 2-copula, given by
F3(α) (x, y, z) = F (α) (x, F (α) (y, z)), is a 3-copula. From this, it follows immediately
that if α ≥ − ln(2) then S F (α) ,F (α) is a 3-copula. However, with some computational
help, it can be shown that S F (α) ,F (α) is also a 3-copula for α ≥ − ln(3).
Example 3 For any 2-copula C2 , SC2 ,Π2 is a 3-copula if and only if for any
x1 , x2 , y1 , y2 , z 1 , z 2 ∈ [0, 1] such that x1 ≤ x2 , y1 ≤ y2 , z 1 ≤ z 2 it holds that
∂ 2 C2 ∂ 2 C2 ∂ 2 C2
(x, y) + (x, z) + (y, z) ≥ 2 . (4)
∂ x∂ y ∂ x∂z ∂ y∂z
Random documents with unrelated
content Scribd suggests to you:
Gutenberg” appears, or with which the phrase “Project
Gutenberg” is associated) is accessed, displayed, performed,
viewed, copied or distributed:
1.E.6. You may convert to and distribute this work in any binary,
compressed, marked up, nonproprietary or proprietary form,
including any word processing or hypertext form. However, if
you provide access to or distribute copies of a Project
Gutenberg™ work in a format other than “Plain Vanilla ASCII” or
other format used in the official version posted on the official
Project Gutenberg™ website (www.gutenberg.org), you must,
at no additional cost, fee or expense to the user, provide a copy,
a means of exporting a copy, or a means of obtaining a copy
upon request, of the work in its original “Plain Vanilla ASCII” or
other form. Any alternate format must include the full Project
Gutenberg™ License as specified in paragraph 1.E.1.
• You pay a royalty fee of 20% of the gross profits you derive
from the use of Project Gutenberg™ works calculated using the
method you already use to calculate your applicable taxes. The
fee is owed to the owner of the Project Gutenberg™ trademark,
but he has agreed to donate royalties under this paragraph to
the Project Gutenberg Literary Archive Foundation. Royalty
payments must be paid within 60 days following each date on
which you prepare (or are legally required to prepare) your
periodic tax returns. Royalty payments should be clearly marked
as such and sent to the Project Gutenberg Literary Archive
Foundation at the address specified in Section 4, “Information
about donations to the Project Gutenberg Literary Archive
Foundation.”
• You comply with all other terms of this agreement for free
distribution of Project Gutenberg™ works.
1.F.
Most people start at our website which has the main PG search
facility: www.gutenberg.org.
Our website is not just a platform for buying books, but a bridge
connecting readers to the timeless values of culture and wisdom. With
an elegant, user-friendly interface and an intelligent search system,
we are committed to providing a quick and convenient shopping
experience. Additionally, our special promotions and home delivery
services ensure that you save time and fully enjoy the joy of reading.
textbookfull.com