35343
35343
https://textbookfull.com/product/uncertainty-modelling-in-data-
science-sebastien-destercke-ed/
https://textbookfull.com/product/beginning-data-science-in-r-
data-analysis-visualization-and-modelling-for-the-data-
scientist-1st-edition-thomas-mailund/
https://textbookfull.com/product/soft-computing-in-data-
science-5th-international-conference-scds-2019-iizuka-japan-
august-28-29-2019-proceedings-michael-w-berry/
https://textbookfull.com/product/advanced-data-analysis-
modelling-in-chemical-engineering-1st-edition-denis-constales/
https://textbookfull.com/product/modelling-based-teaching-in-
science-education-1st-edition-john-k-gilbert/
The Data Science Design Manual 1st Edition Steven S.
Skiena
https://textbookfull.com/product/the-data-science-design-
manual-1st-edition-steven-s-skiena/
https://textbookfull.com/product/research-in-data-science-ellen-
gasparovic/
https://textbookfull.com/product/frontiers-in-data-science-
matthias-dehmer/
https://textbookfull.com/product/uncertainty-analysis-of-
experimental-data-with-r-1st-edition-benjamin-david-shaw/
https://textbookfull.com/product/subsurface-environmental-
modelling-between-science-and-policy-dirk-scheer/
Advances in Intelligent Systems and Computing 832
Sébastien Destercke
Thierry Denoeux · María Ángeles Gil
Przemyslaw Grzegorzewski
Olgierd Hryniewicz Editors
Uncertainty
Modelling
in Data
Science
Advances in Intelligent Systems and Computing
Volume 832
Series editor
Janusz Kacprzyk, Polish Academy of Sciences, Warsaw, Poland
e-mail: [email protected]
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 such as: computational intelligence, soft computing
including neural networks, fuzzy systems, evolutionary computing and the fusion of these paradigms,
social intelligence, ambient intelligence, computational neuroscience, artificial life, virtual worlds and
society, cognitive science and systems, Perception and Vision, DNA and immune based systems,
self-organizing and adaptive systems, e-Learning and teaching, human-centered and human-centric
computing, recommender systems, intelligent control, robotics and mechatronics including human-machine
teaming, knowledge-based paradigms, learning paradigms, machine ethics, intelligent data analysis,
knowledge management, intelligent agents, intelligent decision making and support, intelligent network
security, trust management, interactive entertainment, Web intelligence and multimedia.
The publications within “Advances in Intelligent Systems and Computing” are primarily 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: [email protected]
Members
Rafael Bello Perez, Universidad Central “Marta Abreu” de Las Villas, Santa Clara, Cuba
e-mail: [email protected]
Emilio S. Corchado, University of Salamanca, Salamanca, Spain
e-mail: [email protected]
Hani Hagras, University of Essex, Colchester, UK
e-mail: [email protected]
László T. Kóczy, Széchenyi István University, Győr, Hungary
e-mail: [email protected]
Vladik Kreinovich, University of Texas at El Paso, El Paso, USA
e-mail: [email protected]
Chin-Teng Lin, National Chiao Tung University, Hsinchu, Taiwan
e-mail: [email protected]
Jie Lu, University of Technology, Sydney, Australia
e-mail: [email protected]
Patricia Melin, Tijuana Institute of Technology, Tijuana, Mexico
e-mail: [email protected]
Nadia Nedjah, State University of Rio de Janeiro, Rio de Janeiro, Brazil
e-mail: [email protected]
Ngoc Thanh Nguyen, Wroclaw University of Technology, Wroclaw, Poland
e-mail: [email protected]
Jun Wang, The Chinese University of Hong Kong, Shatin, Hong Kong
e-mail: [email protected]
Olgierd Hryniewicz
Editors
Uncertainty Modelling
in Data Science
123
Editors
Sébastien Destercke Przemyslaw Grzegorzewski
CNRS, Heudiasyc Faculty of Mathematics and Information
Sorbonne universités, Université Science
de technologie de Compiègne Warsaw University of Technology
Compiegne, France Warsaw, Poland
This Springer imprint is published by the registered company Springer Nature Switzerland AG
The registered company address is: Gewerbestrasse 11, 6330 Cham, Switzerland
Preface
This volume contains the peer-reviewed papers presented at the 9th International
Conference on Soft Methods in Probability and Statistics (SMPS 2018), which
was held in conjunction with the 5th International Conference on Belief Functions
(BELIEF 2018) on 17–21 September 2018 in Compiègne, France. The series of
biannual International Conference on Soft Methods in Probability and Statistics
started in Warsaw in 2002. It then successfully took place in Oviedo (2004), Bristol
(2006), Toulouse (2008), Oviedo/Mieres (2010), Konstanz (2012), Warsaw (2014)
and Rome (2016). SMPS and BELIEF 2018 were organized by the Heudiasyc
laboratory at the Université de Technologie de Compiègne.
Over the last decades, the interest for extensions and alternatives to probability
and statistics has significantly grown in areas as diverse as reliability,
decision-making, data mining and machine learning, optimization, etc. This interest
comes from the need to enrich existing models, in order to include different facets
of uncertainty such as ignorance, vagueness, randomness, conflict or imprecision.
Frameworks such as rough sets, fuzzy sets, fuzzy random variables, random sets,
belief functions, possibility theory, imprecise probabilities, lower previsions,
desirable gambles all share this goal, but have emerged from different needs. By
putting together the BELIEF and SMPS conferences, we hope to increase the
interactions and discussions between the two communities and to converge towards
a more unified view of uncertainty theories.
We also think that the advances, results and tools presented in this volume are
important in the ubiquitous and fast-growing fields of data science, machine
learning and artificial intelligence. Indeed, an important aspect of some of the
learned predictive models is the trust one places in them. Modelling carefully and
with principled methods, the uncertainty associated to the data and the models is
one of the means to increase this trust, as the model will then be able to distinguish
reliable predictions from less reliable ones. In addition, extensions such as fuzzy
sets can be explicitly designed to provide interpretable predictive models,
facilitating user interaction and increasing their trust.
v
vi Preface
Programme Committee
vii
viii Organization
ix
x Contents
1 Introduction
n+1−j
S Xn+1 (t) = , for t ∈ (xj , xj+1 ), j = 0, ..., n. (3)
n+1
The difference between the upper and lower survival functions, called impreci-
sion, is non-zero because of the limited inferential assumptions made, and reflects
the amount of information in the data.
This paper is organized as follows. Section 2 introduces the main idea of
imprecise predictive inference based on ALT and log-rank test. The main novelty
of our approach is that the imprecision results from a classical nonparametric
test, which is the log-rank test, integrated with the Arrhenius function to link
different stress levels. In Sect. 3 we explain why we do not use a single log-rank
test on all stress levels. In Sect. 4 our method is illustrated in two examples.
Section 5 presents some concluding remarks.
Imprecise Statistical Inference for Accelerated Life Testing Data 3
In this section we present new predictive inference based on ALT data and
the log-rank test. The proposed new method consists of two steps. First, the
pairwise log-rank test is used between the stress level Ki and K0 , to get the
intervals [γ i , γ i ] of values γi for which we do not reject the null hypothesis that
the data transformed from level i to level 0, and the original data from level 0,
come from the same underlying distribution, where i = 1, ..., m. With these m
pairs (γ i , γ i ), we define γ = min γ i and γ = max γ i .
Second, we apply the data transformation using γ (γ) for all levels to get
transformed data at level 0 which leads to NPI lower (upper) survival function
S (S). Note that each observation at an increased stress level is transformed
to an interval at level 0, where the interval tends to be larger if a data point
was originally from a higher stress level. If the model fits really well, we expect
most γ i to be quite similar, and also most γ i . The NPI lower survival function is
attained when all data observations at increased stress levels are transformed to
the normal stress level using γ, and the NPI upper survival function results from
the use of γ. If the model fits poorly, γ i are likely to differ a lot, or γ i differ a lot,
or both. Hence, in case of poor model fit, the resulting interval [γ, γ] tends to be
wider than in the case of good model fit. A main novelty of our method is that
imprecision results from pairwise comparisons via a classical test, we comment
further on this in the next section.
In our novel method discussed in Sect. 2, we use pairwise log-rank tests between
stress level Ki and K0 . An alternative would be to use one log-rank test for
the data at all stress levels combined. We now explain why this would not lead
to a sensible method of imprecise statistical inference. Suppose we would test
the null hypothesis that data from all stress levels, transformed using parameter
value γa , originate from the same underlying distribution. Let [γ a , γ a ] be the
interval of such values γa for which this hypothesis is not rejected. If the model
fits very well, we would expect γ a to be close to the γ from Sect. 2 and also γ a
to be close to γ. If however, the model fits poorly, the [γ a , γ a ] interval may be
very small or even empty. Therefore, this leads to less imprecision if the model
fits poorly, and that is the reason why we do the pairwise levels and take the
minimum and the maximum of γ i and γ i , respectively. Then, we are interested in
prediction of one future observation at the normal stress level K0 . So, using the
observations transformed from the increased stress levels K1 , ..., Km as well as
the original data obtained at the normal stress level K0 , we apply NPI to derive
lower and upper survival functions for as described in Sect. 1. The examples in
Sect. 4 illustrate the proposed method of Sect. 2 as well as the problem if we
would use the combined approach for all levels.
4 A. A. H. Ahmadini and F. P. A. Coolen
4 Examples
In this section we present two examples. In example 1 we simulated data at
all levels that correspond to the model for the link function we assume for the
analysis. In example 2 we change these data such that the assumed link function
will not provide a good fit anymore. Together, these examples illustrate our novel
imprecise method, from Sect. 2, as well as the problem that could occur if we
used the log-rank test on all stress levels combined, as discussed in Sect. 3.
4.1 Example 1
The method proposed in Sect. 2 is illustrated in an example, which presents
the temperature-accelerated lifespan test. Data are simulated at three temper-
atures. The normal temperature condition was K0 = 283 and the increased
temperatures stress levels were K1 = 313 and K2 = 353 Kelvin. Ten obser-
vations were simulated from a fully specified model, using the Arrhenius link
function in combination with a Weibull distribution at each temperature. The
Arrhenius parameter γ was set at 5200, and the Weibull distribution at K0 had
shape parameter 3 and scale parameter 7000. This model keeps the same shape
parameter at each temperature, but the scale parameter are linked by the Arrhe-
nius relation, which led to scale parameter 1202.942 at K1 and 183.0914 at K2 .
Imprecise Statistical Inference for Accelerated Life Testing Data 5
Ten units were tested at each temperature, for a total of 30 units used in the
study. The failure times, in hours, are given in Table 1.
To illustrate the log-rank test method using these data, we assume the Arrhe-
nius link function for the data. Note that our method does not assume a paramet-
ric distribution at each stress level. The pairwise log-rank test is used between
K1 and K0 and between K2 and K0 to derive the intervals [γ i , γ i ] of values γi for
which we do not reject the null hypothesis with regard to the well-mixed data
transformation. The resulting intervals [γ i , γ i ] are giving in the first two rows of
Table 2, for three test significance levels. Of course, for larger significance level
the intervals become wider.
According to the accepted intervals in Table 2, we can obtain the NPI lower
and upper survival functions by taking from the pairwise stress level K1 to K0 or
K2 to K0 always the minimum of the γ i and the maximum of the γ i with levels
of significance 0.99, 0.95 and 0.90 values. So, we take γ = min γ i = 3901.267 and
the γ = max γ i = 6563.545 of the pairwise K1 , K0 with 0.99 significance level,
γ = min γ i = 4254.053 and the γ = max γ i = 6251.168 of the pairwise K1 , K0
with 0.95 significance level, and γ = min γ i = 4486.491 and the γ = max γ i =
6017.435 of the pairwise K1 , K0 with 0.90 significance level then transformed
the data to the normal stress level, see Fig. 1(a). In this figure, the lower survival
function S is labeled as S (γ i ) and the upper survival function S is labeled as S
(γ i ). This figure shows that higher significance levels leads to more imprecision
for the NPI lower and upper survival functions.
To illustrate the effect of using the single log-rank test for all stress levels
simultaneously as discussed in Sect. 3, the final row in Table 2 provides the inter-
val [γ a , γ a ] of values γa for all the stress levels together. From this interval we
can again obtain the lower and upper survival functions using NPI, these are
presented in Fig. 1(b). In this example, the data were simulated precisely with
the link function as assumed in our method, so there is not much difference
between the lower and upper survival functions for corresponding significance
levels in Figs. 1(a) and (b). Example 2 will illustrate what happens if the model
does not fit well.
4.2 Example 2
To illustrate our method in case the model does not fit the data well, and also
to show what would have happened if we had used the joint log-rank test in our
method instead of the pairwise tests, we use the same data as in Example 1,
but we change some of these. In Scenario 1 (indicated as Ex 2-1 in Fig. 1), we
multiple the data at level K1 by 1.4. In Scenario 2 (Ex 2-2), we do the same and
in addition we multiply the data at level K2 by 0.8. The resulting data values
are given in the last two columns in Table 1.
For these two scenarios, we have repeated the analysis as described in Exam-
ple 1. The resulting intervals of γ values are given in Tables 3 and 4. Note that
for significance level 0.90 in Scenario 2 the null hypothesis of the joint log-rank
test would be rejected for all values γa , hence we report an empty interval, so
6 A. A. H. Ahmadini and F. P. A. Coolen
clearly our method would not work if we had used this joint test instead of the
pairwise tests.
The NPI lower and upper survival functions in Figs. 1(c) and (e), using our
method as discussed in Sect. 2, have more imprecision. Note that the lower sur-
vival function is identical in both scenarios as the same γ is used, this is because
the increased values at K1 have resulted in smaller values for γ 1 and γ 1 and the
γ in our method is equal to the γ 1 in these cases. In Scenario 2, the observations
at level K2 have decreased, leading to larger γ 2 and γ 2 values, and this leads to
the upper survival functions increasing in comparison to Scenario 1.
If we would have used the joint long-rank test instead of the pairwise tests,
as discussed in Sect. 3, then imprecision would have decreased in these two sce-
narios, as can be seen from Figs. 1(d) and (f). Note that in Fig. 1(f) there are no
lower and upper survival functions corresponding to the use of the joint log-rank
test for significance level 0.90, as this leads to an empty interval of γa values. As
mentioned in Sect. 3, if the model does not fit well, then we are going to sooner
reject the null hypothesis for all the three levels together, see Tables 3 and 4. So
we have a smaller range of values for which we do not reject the null hypoth-
esis. But if the model fits poorly, we actually want a larger range of values, so
increased imprecision. It is obvious that this is achieved by taking the minimum
of the γ i and the maximum of the γ i of the pairwise tests, hence this is our
proposed method in Sect. 2. This is illustrated by Figs. 1(a), (c) and (e).
5 Concluding Remarks
This paper has presented an exploration of the use of a novel statistical method
providing imprecise semi-parametric inference for ALT data, where the impreci-
sion is related to the log-rank test statistics. The proposed method applies the
use of the log-rank test to compare the survival distribution of pairwise stress
levels, in combination with the Arrhenius model finding the accepted interval of
γ values according to the null hypothesis. We explored imprecision through the
use of nonparametric test for the parameter of the link function between different
stress levels, which enabled us to transform the observations at increased stress
levels to interval-valued observations at the normal stress level and achieve fur-
ther robustness. We consider nonparametric predictive inference at the normal
stress level combined with the Arrhenius model linking observations at different
stress levels. We showed why, in our method, we use the imprecision from com-
bined pairwise log-rank tests, and not from a single log-rank test on all stress
levels together. The latter would lead to less imprecision if the model fits poorly,
while our proposed method then leads to more imprecision. In this paper, to illus-
trate basic idea of our novel method, we assumed that failure data are available
at all stress levels including the normal stress level. This may not be realistic. If
there are no failure data at the normal stress level, or only right-censored obser-
vations, then we can apply our method using a higher stress level as the basis
for the combinations, so transform data to that stress level. Then the combined
data at that level could be transformed all together to the normal stress level.
The log-rank test in this approach could be replaced by other comparison tests,
where even the use of tests based on imprecise probability theory [7] could be
explored. This is left as an interesting topic for future research.
References
1. Augustin, T., Coolen, F., de Cooman, G., Troffaes, M.: Introduction to Imprecise
Probabilities. Wiley, Chichester (2014)
2. Coolen, F.: Nonparametric predictive inference. In: International Encyclopedia of
Statistical Science, pp. 968–970. Springer, Berlin (2011)
3. Gehan, E.: A generalized Wilcoxon test for comparing arbitrarily singly-censored
samples. Biometrika 52, 203–224 (1965)
4. Mantel, N.: Evaluation of survival data and two new rank order statistics arising in
its consideration. Cancer Chemother. Rep. 50, 163–170 (1966)
5. Nelson, W.: Accelerated Testing: Statistical Models, Test Plans, and Data Analysis.
Wiley, New Jersey (1990)
6. Peto, R., Peto, J.: Asymptotically efficient rank invariant test procedures. J. R.
Stat. Soc. Ser. A 135, 185–207 (1972)
7. Benavoli, A., Mangili, F., Corani, G., Zaffalon, M., Ruggeri, F.: A Bayesian Wilcoxon
signed-rank test based on the Dirichlet process. In: Proceedings of the 30th Inter-
national Conference on Machine Learning (ICML 2014), pp. 1–9 (2014)
Descriptive Comparison of the Rating
Scales Through Different Scale Estimates:
Simulation-Based Analysis
1 Introduction
The Likert-type scales are frequently used in designing questionnaires to rate
characteristics or attributes that cannot be numerically measured (like satis-
faction, perceived quality, perception...). Although they are easy to answer and
they do not require a special training to use them, respondents often do not
find accurate answers to items and the available statistical methodology to ana-
lyze the data from these questionnaires is rather limited. This is mainly due to
the fact that Likert scales are discrete with a very small number of responses
to choose for each item (often 4 to 7). To overcome this concern, Hesketh et
al. [5] proposed the so-called fuzzy rating scale to allow a complete freedom and
expressiveness in responding, without respondents being constrained to choose
among a few pre-specified responses.
c Springer Nature Switzerland AG 2019
S. Destercke et al. (Eds.): SMPS 2018, AISC 832, pp. 9–16, 2019.
https://doi.org/10.1007/978-3-319-97547-4_2
10 I. Arellano et al.
By drawing the fuzzy number that best represents the respondent’s valua-
tion, the fuzzy rating scale captures the logical imprecision associated with such
variables. Moreover, this fuzzy rating scale allows us to have a rich continuous
scale of measurement, unlike the case of a posterior numerical or fuzzy encoding
(the latter encoding Likert points with fuzzy numbers from a linguistic scale,
and usually made by trained experts).
In previous studies (see Gil et al. [3], Lubiano et al. [6–8]) we have con-
firmed that the results when fuzzy rating scales are considered sometimes differ
importantly from the conclusions drawn from numerically or fuzzy linguistically
encoded Likert values.
As differences can often be even clearer from the dispersion than for the
location perspective, this paper aims to examine, by means of simulation devel-
opments, how location-based ‘scale’ estimates are affected by the considered scale
of measurement.
2 Preliminaries
A (bounded) fuzzy number is a mapping U : R → [0, 1] such that for all
α ∈ [0, 1], the α-level set Uα = {x ∈ R : U (x) ≥ α} if α ∈ (0, 1], and U 0 =
cl{x ∈ R : U (x) > 0} (with ‘cl’ denoting the closure of the set) is a nonempty
compact interval.
In dealing with fuzzy number-valued data, distances will be computed by
considering two different metrics introduced by Diamond and Kloeden [1]: the
2-norm metric ρ2 and the 1-norm metric ρ1 , which for fuzzy numbers U and V
are given by
1
ρ2 (U , V ) = (inf Uα − inf Vα )2 + (sup Uα − sup Vα )2 dα,
2 [0,1]
1 α − inf Vα | + | sup U
α − sup Vα | dα.
ρ1 (U , V ) = | inf U
2 [0,1]
and the sample 1-norm median is the fuzzy number such that for each α
In De la Rosa de Sáa et al. [2] one can find together the most commonly
used location-based scale estimates, namely: the sample Fréchet-type ρ2 -
xn , Me(
Standard Deviation and, for D ∈ {ρ1 , ρ2 } and M ∈ { xn )}, the sam-
ple D-Average Distance Deviation and the sample D-Median Distance
Deviation, which are respectively given by
1 2
n
xn ) =
ρ2 -SD( n ) ,
xi , x
ρ2 (
n i=1
1
n
xn , M ) =
D-ADD( D( xn , M ) = Mei D(
xi , M ), D-MDD( xi , M ) .
n i=1
In fact, fuzzy data will be generated by simulating the four real-valued ran-
dom variables X1 , X2 , X3 and X4 , so that the R×[0, ∞)×[0, ∞)×[0, ∞)-valued
random vector (X1 , X2 , X3 , X4 ) will provide us with the 4-tuples (x1 , x2 , x3 , x4 )
with x1 = center and x2 = radius of the core, and x3 = lower and x4 = upper
spread of the fuzzy number. To each generated 4-tuple (x1 , x2 , x3 , x4 ) we asso-
ciate the fuzzy number Trax1 , x2 , x3 , x4 .
According to the simulation procedure, data have been generated from ran-
dom fuzzy numbers with a bounded reference set and abstracting and mimicking
what we have observed in real-life examples employing the fuzzy rating scale
(FRS). More concretely, fuzzy data have been generated such that
12 I. Arellano et al.
– 100·ω1 % of the data have been obtained by first considering a simulation from
a simple random sample of size 4 from a beta β(p, q) distribution, ordering the
corresponding 4-tuple, and finally computing the values xi . The values of p
and q vary in most cases to cover different distributions (namely, symmetrical
weighting central values, symmetrical weighting extreme values, and asym-
metric ones). In most of the comparative studies involving simulations, the
values from the beta distribution are re-scaled and translated to an interval
[l0 , u0 ] different from [0, 1].
– 100 · ω2 % of the data have been obtained considering a simulation of four
random variables Xi = (u0 − l0 ) · Yi + l0 as follows:
Y1 ∼ β(p, q),
Y2 ∼ Uniform0, min{1/10, Y1 , 1 − Y 1 } ,
Y3 ∼ Uniform0, min{1/5, Y1 − Y2 } ,
Y4 ∼ Uniform 0, min{1/5, 1 − Y1 − Y2 } .
– 100 · ω3 % of the data have been obtained considering a simulation of four
random variables Xi = (u0 − l0 ) · Yi + l0 as follows:
Y1 ∼ β(p, q),
⎧
⎨ Exp(200) if Y1 ∈ [0.25, 0.75]
Y2 ∼ Exp(100 + 4 Y1 ) if Y1 < 0.25
⎩
Exp(500 − 4 Y1 ) otherwise
γ(4, 100) if Y1 − Y2 ≥ 0.25
Y3 ∼
γ(4, 100 + 4 Y1 ) otherwise
γ(4, 100) if Y1 + Y2 ≥ 0.25
Y4 ∼
γ(4, 500 − 4 Y1 ) otherwise.
5 Results
First, FRS data will be simulated in accordance with the above described realistic
simulation procedure. Later, fuzzy data based on a fuzzy rating scale can fairly be
associated/classified in accordance with labels in a Likert scale (more concretely,
with their numerical encoding). This process is to be called “Likertization”.
Furthermore, the associated Likert values could also be later encoded by means
of values from a fuzzy linguistic scale.
Descriptive Comparison of the Rating Scales 13
For carrying out the Likertization, the “minimum distance Likertization cri-
terion” will be employed (see Fig. 2):
Fig. 2. Minimum distance criterion scheme when the reference interval equals [1, k]
In this way, if the considered Likert scale is a k-point one, given a metric
D between fuzzy data and U the free fuzzy response to be classified, then U
is
associated with the integer κ(U ) such that
) = arg
κ(U min , 1{j} ).
D(U
j∈{1,...,k}
1 1
0 0
0 20 40 60 80 100 0 20 40 60 80 100
Table 1. % of simulated samples of size n for which the Euclidean distance between
the sample scale estimate D associated with the FRS and the one associated with
either the NEL (numerically encoded Likert) or the FLS (fuzzy linguistic scale) with
k = 4 different values is greater than ε ∈ {1, 5, 10, 15} and (from top to bottom)
β(p, q) ≡ β(1, 1), β(.75, .75), β(4, 2), and β(6, 1)
% D(FRS) − D(S) > ε (k = 4, β(p, q) ≡ β(1, 1))
D ε = 1 ε = 5 ε = 10 ε = 15
n S = NEL S = FLS S = NEL S = FLS S = NEL S = FLS S = NEL S = FLS
ρ
2 SD(
xn )
ρ2 n )
xn , x
ADD(
xn ))
ρ1
ADD(
xn , Me(
ρ2
MDD( n )
xn , x
xn ))
ρ1
MDD(
xn , Me(
% D(FRS) − D(S) > ε (k = 4, β(p, q) ≡ β(0.75, 0.75))
D ε = 1 ε = 5 ε = 10 ε = 15
n S = NEL S = FLS S = NEL S = FLS S = NEL S = FLS S = NEL S = FLS
ρ
2 SD(
xn )
ρ2
ADD( n )
xn , x
xn ))
ρ1
ADD(
xn , Me(
ρ2
MDD( n )
xn , x
xn ))
ρ1 xn , Me(
MDD(
% D(FRS) − D(S) > ε (k = 4, β(p, q) ≡ β(4, 2))
D ε = 1 ε = 5 ε = 10 ε = 15
n S = NEL S = FLS S = NEL S = FLS S = NEL S = FLS S = NEL S = FLS
ρ
2 SD(
xn )
ρ2
ADD( n )
xn , x
xn ))
ρ1
ADD(
xn , Me(
ρ2
MDD( n )
xn , x
xn ))
ρ1
MDD(
xn , Me(
% D(FRS) − D(S) > ε (k = 4, β(p, q) ≡ β(6, 1))
D ε = 1 ε = 5 ε = 10 ε = 15
n S = NEL S = FLS S = NEL S = FLS S = NEL S = FLS S = NEL S = FLS
ρ
2 SD(
xn )
ρ2
ADD( n )
xn , x
xn ))
ρ1 xn , Me(
ADD(
ρ2
MDD( n )
xn , x
xn ))
ρ1 xn , Me(
MDD(
Descriptive Comparison of the Rating Scales 15
Table 2. % of simulated samples of size n for which the Euclidean distance between
the sample scale estimate D associated with the FRS and the one associated with
either the NEL (numerically encoded Likert) or the FLS (fuzzy linguistic scale) with
k = 5 different values is greater than ε ∈ {1, 5, 10, 15} and (from top to bottom)
β(p, q) ≡ β(1, 1), β(.75, .75), β(4, 2), and β(6, 1)
% D(FRS) − D(S) > ε (k = 5, β(p, q) ≡ β(1, 1))
D ε = 1 ε = 5 ε = 10 ε = 15
n S = NEL S = FLS S = NEL S = FLS S = NEL S = FLS S = NEL S = FLS
ρ
2 SD(
xn )
ρ2 n )
xn , x
ADD(
xn ))
ρ1
ADD(
xn , Me(
ρ2
MDD( n )
xn , x
xn ))
ρ1
MDD(
xn , Me(
% D(FRS) − D(S) > ε (k = 5, β(p, q) ≡ β(0.75, 0.75))
D ε = 1 ε = 5 ε = 10 ε = 15
n S = NEL S = FLS S = NEL S = FLS S = NEL S = FLS S = NEL S = FLS
ρ
2 SD(
xn )
ρ2
ADD( n )
xn , x
xn ))
ρ1
ADD(
xn , Me(
ρ2
MDD( n )
xn , x
xn ))
ρ1 xn , Me(
MDD(
% D(FRS) − D(S) > ε (k = 5, β(p, q) ≡ β(4, 2))
D ε = 1 ε = 5 ε = 10 ε = 15
n S = NEL S = FLS S = NEL S = FLS S = NEL S = FLS S = NEL S = FLS
ρ
2 SD(
xn )
ρ2
ADD( n )
xn , x
xn ))
ρ1
ADD(
xn , Me(
ρ2
MDD( n )
xn , x
xn ))
ρ1
MDD(
xn , Me(
% D(FRS) − D(S) > ε (k = 5, β(p, q) ≡ β(6, 1))
D ε = 1 ε = 5 ε = 10 ε = 15
n S = NEL S = FLS S = NEL S = FLS S = NEL S = FLS S = NEL S = FLS
ρ
2 SD(
xn )
ρ2
ADD( n )
xn , x
xn ))
ρ1
ADD(
xn , Me(
ρ2
MDD( n )
xn , x
xn ))
ρ1 xn , Me(
MDD(
16 I. Arellano et al.
The percentages have been quantified over 1000 samples of n ∈ {10, 30, 100}
FRS simulated (with different betas) data with reference interval [0, 100] (this
last fact being irrelevant for the study). On the basis of Tables 1 and 2 we cannot
get very general conclusions, but we can definitely assert that scale measures
mostly vary more from the FRS-based data to the encoded Likert ones.
Furthermore, one can state some approximate behaviour patterns, such as
– for almost all situations, the robust scale estimate (the last one) provides
us with much higher percentages than non-robust ones; more concretely,
ρ1
-MDD( xn , Me(
xn )) is almost generally more sensitive to the change in the
rating scale type; this is especially clear for small samples;
– distances are uniformly lower for k = 5 than for k = 4 when the midpoint
of the 1-level is beta distributed with (p, q) ∈ {(1, 1), (0.75, 0.75), (4, 2)};
when (p, q) = (6, 1) such a conclusion is appropriate for robust estimates
and ε ∈ {1, 5}, but there is no clear conclusion for non-robust estimates or
greater values of ε.
Acknowledgements. The research is this paper has been partially supported by the
Spanish Ministry of Economy, Industry and Competitiveness Grant MTM2015-63971-
P. Its support is gratefully acknowledged.
References
1. Diamond, P., Kloeden, P.: Metric spaces of fuzzy sets. Fuzzy Sets Syst. 35, 241–249
(1990)
2. De la Rosa de Sáa, S., Lubiano, S., Sinova, S., Filzmoser, P.: Robust scale estimators
for fuzzy data. Adv. Data Anal. Classif. 11, 731–758 (2017)
3. Gil, M.A., Lubiano, M.A., De la Rosa de Sáa, S., Sinova, B.: Analyzing data from
a fuzzy rating scale-based questionnaire: a case study. Psicothema 27, 182–191
(2015)
4. Herrera, F., Herrera-Viedma, E., Martı́nez, L.: A fuzzy linguistic methodology to
deal with unbalanced linguistic term sets. IEEE Trans. Fuzzy Syst. 16(2), 354–370
(2008)
5. Hesketh, T., Pryor, R., Hesketh, B.: An application of a computerized fuzzy graphic
rating scale to the psychological measurement of individual differences. Int. J. Man-
Mach. Stud. 29, 21–35 (1988)
6. Lubiano, M.A., De la Rosa de Sáa, S., Montenegro, M., Sinova, B., Gil, M.A.:
Descriptive analysis of responses to items in questionnaires. Why not using a fuzzy
rating scale? Inf. Sci. 360, 131–148 (2016)
7. Lubiano, M.A., Montenegro, M., Sinova, B., De la Rosa de Sáa, S., Gil, M.A.:
Hypothesis testing for means in connection with fuzzy rating scale-based data:
algorithms and applications. Eur. J. Oper. Res. 251, 918–929 (2016)
8. Lubiano, M.A., Salas, A., Gil, M.A.: A hypothesis testing-based discussion on the
sensitivity of means of fuzzy data with respect to data shape. Fuzzy Sets Syst.
328, 54–69 (2017)
9. Puri, M.L., Ralescu, D.A.: Fuzzy random variables. J. Math. Anal. Appl. 114,
409–422 (1986)
10. Sinova, B., Gil, M.A., Colubi, A., Van Aelst, S.: The median of a random fuzzy
number. The 1-norm distance approach. Fuzzy Sets Syst. 200, 99–115 (2012)
Central Moments of a Fuzzy Random
Variable Using the Signed Distance:
A Look Towards the Variance
I. ABORIGINALS OF BOTH.
Europe. America.
lb. lb.
Mammoth
Buffalo. Bison *1800
White Bear. Ours blanc
Carribou. Renne
Bear. Ours 153.7 *410
Elk. Elan. Original palmated
Red deer. Cerf 288.8 *273
Fallow Deer. Daim 167.8
Wolf. Loup 69.8
Roe. Chevreuil 56.7
Glutton. Glouton. Carcajou
Wild cat. Chat sauvage †30
Lynx. Loup cervier 25.
Beaver. Castor 18.5 *45
Badger. Blaireau 13.6
Red fox. Renard 13.5
Gray fox. Isatis
Otter. Loutre 8.9 †12
Monax. Marmotte 6.5
Vison. Fouine 2.8
Hedgehog. Herisson 2.2
Marten. Marte 1.9 †6
oz.
Water rat. Rat d'eau 7.5
Weasel. Belette 2.2 oz.
Flying squirrel. Polatouche 2.2 †4
Shrew mouse. Musaraigne 1.
EUROPE.
lb.
Sanglier. Wild boar 280.
Mouflon. Wild sheep 56.
Bouquetin. Wild goat
Lievre. Hare 7.6
Lapin. Rabbit 3.4
Putois. Polecat 3.3
Genette 3.1
Desman. Muskrat oz.
Ecureuil. Squirrel 12.
Hermine. Ermin 8.2
Rat. Rat 7.5
Loirs 3.1
Lerot. Dormouse 1.8
Taupe. Mole 1.2
Hampster .6
Zisel
Leming
Souris. Mouse .6
AMERICA
lb.
Tapir 534.
Elk, round horned †450.
Puma
Jaguar 218.
Cabiai 109.
Tamanoir 109.
Tammandua 65.4
Cougar of North-America 75.
Cougar of South-America 59.4
Ocelot
Pecari 46.3
Jaguaret 43.6
Alco
Lama
Paco
Paca 32.7
Serval
Sloth. Unau 27.25
Saricovienne
Kincajou
Tatou Kabassou 21.8
Urson. Urchin
Raccoon. Raton 16.5
Coati
Coendou 16.3
Sloth. Aï 13.
Sapajou Ouarini
Sapajou Coaita 9.8
Tatou Encubert
Tatou Apar
Tatou Cachiea 7.
Little Coendou 6.5
Opossum. Sarigu
Tapeti
Margay
Crabier
Agouti 4.2
Sapajou Saï 3.5
Tatou Cirquinçon
Tatou Tatouate 3.3
Mouffette Squash
Mouffette Chinche
Mouffette Conepate
Scunk
Mouffette. Zorilla
Whabus. Hare. Rabbit
Aperea
Akouchi
Ondatra. Muskrat
Pilori
Great gray squirrel †2.7
Fox squirrel of Virginia †2.625
Surikate 2.
Mink †2.
Sapajou. Sajou 1.8
Indian pig. Cochon d'Inde 1.6
Sapajou Saïmiri 1.5
Phalanger
Coqualain
Lesser gray squirrel †1.5
Black squirrel †1.5
oz.
Red squirrel 10.
Sagoin Saki
Sagoin Pinche
Sagoin Tamarin
Sagoin Ouistiti 4.4
Sagoin Marakine
Sagoin Mico
Cayopollin
Fourmillier
Marmose
Sarigue of Cayenne
Tucan
Red mole
Ground squirrel 4.
I have not inserted in the first table the Phoca,[9] nor leather-winged
bat, because the one living half the year in the water, and the other
being a winged animal, the individuals of each species may visit both
continents.
Of the animals in the first table, Monsieur de Buffon himself informs
us, [XXVII. 130, XXX. 213,] that the beaver, the otter, and shrew
mouse, though of the same species, are larger in America than in
Europe. This should therefore have corrected the generality of his
expressions, XVIII. 145, and elsewhere, that the animals common to
the two countries, are considerably less in America than in Europe,
"et cela sans aucune exception." He tells us too, [Quadrup. VIII.
334, edit. Paris, 1777,] that on examining a bear from America, he
remarked no difference, "dans la forme de cet ours d'Amerique
comparé a celui d'Europe," but adds from Bartram's journal, that an
American bear weighed four hundred pounds, English, equal to three
hundred and sixty-seven pounds French; whereas we find the
European bear examined by Mons. D'Aubenton, [XVII. 82,] weighed
but one hundred and forty-one pounds French. That the palmated
elk is larger in America than in Europe, we are informed by Kalm,[10]
a naturalist, who visited the former by public appointment, for the
express purpose of examining the subjects of natural history. In this
fact Pennant concurs with him. [Barrington's Miscellanies.] The same
Kalm tells us[11] that the black moose, or renne of America, is as
high as a tall horse; and Catesby,[12] that it is about the bigness of a
middle-sized ox. The same account of their size has been given me
by many who have seen them. But Monsieur D'Aubenton says[13]
that the renne of Europe is about the size of a red deer. The weasel
is larger in America than in Europe, as may be seen by comparing its
dimensions as reported by Monsieur D'Aubenton[14] and Kalm. The
latter tells us,[15] that the lynx, badger, red fox, and flying squirrel,
are the same in America as in Europe; by which expression I
understand, they are the same in all material circumstances, in size
as well as others; for if they were smaller, they would differ from the
European. Our gray fox is, by Catesby's account,[16] little different in
size and shape from the European fox. I presume he means the red
fox of Europe, as does Kalm, where he says,[17] that in size "they do
not quite come up to our foxes." For proceeding next to the red fox
of America, he says, "they are entirely the same with the European
sort;" which shows he had in view one European sort only, which
was the red. So that the result of their testimony is, that the
American gray fox is somewhat less than the European red; which is
equally true of the gray fox of Europe, as may be seen by comparing
the measures of the Count de Buffon and Monsieur D'Aubenton.[18]
The white bear of America is as large as that of Europe. The bones
of the mammoth which has been found in America, are as large as
those found in the old world. It may be asked, why I insert the
mammoth, as if it still existed? I ask in return, why I should omit it,
as if it did not exist? Such is the economy of nature, that no instance
can be produced, of her having permitted any one race of her
animals to become extinct; of her having formed any link in her
great work so weak as to be broken. To add to this, the traditionary
testimony of the Indians, that this animal still exists in the northern
and western parts of America, would be adding the light of a taper
to that of the meridian sun. Those parts still remain in their
aboriginal state, unexplored and undisturbed by us, or by others for
us. He may as well exist there now, as he did formerly where we
find his bones. If he be a carnivorous animal, as some anatomists
have conjectured, and the Indians affirm, his early retirement may
be accounted for from the general destruction of the wild game by
the Indians, which commences in the first instant of their connection
with us, for the purpose of purchasing match-coats, hatchets, and
firelocks, with their skins. There remain then the buffalo, red deer,
fallow deer, wolf, roe, glutton, wild cat, monax, bison, hedgehog,
marten, and water-rat, of the comparative sizes of which we have
not sufficient testimony. It does not appear that Messieurs de Buffon
and D'Aubenton have measured, weighed, or seen those of America.
It is said of some of them, by some travellers, that they are smaller
than the European. But who were these travellers? Have they not
been men of a very different description from those who have laid
open to us the other three quarters of the world? Was natural
history the object of their travels? Did they measure or weigh the
animals they speak of? or did they not judge of them by sight, or
perhaps even from report only? Were they acquainted with the
animals of their own country, with which they undertake to compare
them? Have they not been so ignorant as often to mistake the
species? A true answer to these questions would probably lighten
their authority, so as to render it insufficient for the foundation of an
hypothesis. How unripe we yet are, for an accurate comparison of
the animals of the two countries, will appear from the work of
Monsieur de Buffon. The ideas we should have formed of the sizes
of some animals, from the information he had received at his first
publications concerning them, are very different from what his
subsequent communications give us. And indeed his candor in this
can never be too much praised. One sentence of his book must do
him immortal honor. "J'aime autant une personne qui me releve
d'une erreur, qu'une autre qui m'apprend une verité, parce qu'en
effet une erreur corrigée est une verité."[19] He seems to have
thought the cabiai he first examined wanted little of its full growth.
"Il n'etoit pas encore tout-a-fait adulte."[20] Yet he weighed but
forty-six and a half pounds, and he found afterwards,[21] that these
animals, when full grown, weigh one hundred pounds. He had
supposed, from the examination of a jaguar,[22] said to be two years
old, which weighed but sixteen pounds twelve ounces, that when he
should have acquired his full growth, he would not be larger than a
middle-sized dog. But a subsequent account[23] raises his weight to
two hundred pounds. Further information will, doubtless, produce
further corrections. The wonder is, not that there is yet something in
this great work to correct, but that there is so little. The result of this
view then is, that of twenty-six quadrupeds common to both
countries, seven are said to be larger in America, seven of equal
size, and twelve not sufficiently examined. So that the first table
impeaches the first member of the assertion, that of the animals
common to both countries, the American are smallest, "et cela sans
aucune exception." It shows it is not just, in all the latitude in which
its author has advanced it, and probably not to such a degree as to
found a distinction between the two countries.
Proceeding to the second table, which arranges the animals found in
one of the two countries only, Monsieur de Buffon observes, that the
tapir, the elephant of America, is but of the size of a small cow. To
preserve our comparison, I will add, that the wild boar, the elephant
of Europe, is little more than half that size. I have made an elk with
round or cylindrical horns an animal of America, and peculiar to it;
because I have seen many of them myself, and more of their horns;
and because I can say, from the best information, that, in Virginia,
this kind of elk has abounded much, and still exists in smaller
numbers; and I could never learn that the palmated kind had been
seen here at all. I suppose this confined to the more northern
latitudes.[24] I have made our hare or rabbit peculiar, believing it to
be different from both the European animals of those
denominations, and calling it therefore by its Algonquin name,
Whabus, to keep it distinct from these. Kalm is of the same opinion.
[25] I have enumerated the squirrels according to our own
knowledge, derived from daily sight of them, because I am not able
to reconcile with that the European appellations and descriptions. I
have heard of other species, but they have never come within my
own notice. These, I think, are the only instances in which I have
departed from the authority of Monsieur de Buffon in the
construction of this table. I take him for my ground work, because I
think him the best informed of any naturalist who has ever written.
The result is, that there are eighteen quadrupeds peculiar to Europe;
more than four times as many, to wit, seventy four, peculiar to
America; that the[26] first of these seventy-four weighs more than
the whole column of Europeans; and consequently this second table
disproves the second member of the assertion, that the animals
peculiar to the new world are on a smaller scale, so far as that
assertion relied on European animals for support; and it is in full
opposition to the theory which makes the animal volume to depend
on the circumstances of heat and moisture.
The third table comprehends those quadrupeds only which are
domestic in both countries. That some of these, in some parts of
America, have become less than their original stock, is doubtless
true; and the reason is very obvious. In a thinly-peopled country, the
spontaneous productions of the forests, and waste fields, are
sufficient to support indifferently the domestic animals of the farmer,
with a very little aid from him, in the severest and scarcest season.
He therefore finds it more convenient to receive them from the hand
of nature in that indifferent state, than to keep up their size by a
care and nourishment which would cost him much labor. If, on this
low fare, these animals dwindle, it is no more than they do in those
parts of Europe where the poverty of the soil, or the poverty of the
owner, reduces them to the same scanty subsistence. It is the
uniform effect of one and the same cause, whether acting on this or
that side of the globe. It would be erring, therefore, against this rule
of philosophy, which teaches us to ascribe like effects to like causes,
should we impute this diminution of size in America to any imbecility
or want of uniformity in the operations of nature. It may be affirmed
with truth, that, in those countries, and with those individuals in
America, where necessity or curiosity has produced equal attention,
as in Europe, to the nourishment of animals, the horses, cattle,
sheep, and hogs, of the one continent are as large as those of the
other. There are particular instances, well attested, where individuals
of this country have imported good breeders from England, and
have improved their size by care in the course of some years. To
make a fair comparison between the two countries, it will not answer
to bring together animals of what might be deemed the middle or
ordinary size of then species; because an error in judging of that
middle or ordinary size, would vary the result of the comparison.
Thus Mons. D'Aubenton[27] considers a horse of 4 feet five inches
high and 400 lb. weight French, equal to 4 feet 8.6 inches and 436
lb. English, as a middle-sized horse. Such a one is deemed a small
horse in America. The extremes must therefore be resorted to. The
same anatomist[28] dissected a horse of 5 feet 9 inches height,
French measure, equal to 6 feet 1.7 English. This is near 6 inches
higher than any horse I have seen; and could it be supposed that I
had seen the largest horses in America, the conclusion would be,
that ours have diminished, or that we have bred from a smaller
stock. In Connecticut and Rhode Island, where the climate is
favorable to the production of grass, bullocks have been slaughtered
which weighed 2,500, 2,200, and 2,100 lbs. nett; and those of 1,800
lbs. have been frequent. I have seen a hog[29] weigh 1,050 lbs. after
the blood, bowels, and hair had been taken from him. Before he was
killed, an attempt was made to weigh him with a pair of steel yards,
graduated to 1,200 lbs., but he weighed more. Yet this hog was
probably not within fifty generations of the European stock. I am
well informed of another which weighed 1,100 lbs. gross. Asses have
been still more neglected than any other domestic animal in
America. They are neither fed or housed in the most rigorous season
of the year. Yet they are larger than those measured by Mons.
D'Aubenton,[30] of 3 feet 7¼ inches, 3 feet 4 inches, and 3 feet 2½
inches, the latter weighing only 215.8 lbs. These sizes, I suppose,
have been produced by the same negligence in Europe, which has
produced a like diminution here. Where care has been taken of them
on that side of the water, they have been raised to a size bordering
on that of the horse; not by the heat and dryness of the climate, but
by good food and shelter. Goats have been also much neglected in
America. Yet they are very prolific here, bearing twice or three times
a year, and from one to five kids at a birth. Mons. de Buffon has
been sensible of a difference in this circumstance in favor of
America.[31] But what are their greatest weights, I cannot say. A
large sheep here weighs 100 lbs. I observe Mons. D'Aubenton calls a
ram of 62 lbs. one of the middle size.[32] But to say what are the
extremes of growth in these and the other domestic animals of
America, would require information of which no one individual is
possessed. The weights actually known and stated in the third table
preceding will suffice to show, that we may conclude on probable
grounds, that, with equal food and care, the climate of America will
preserve the races of domestic animals as large as the European
stock from which they are derived; and, consequently, that the third
member of Mons. de Buffon's assertion that the domestic animals
are subject to degeneration from the climate of America, is as
probably wrong as the first and second were certainly so.
That the last part of it is erroneous, which affirms that the species of
American quadrupeds are comparatively few, is evident from the
tables taken together. By these it appears that there are an hundred
species aboriginal in America. Mons. de Buffon supposes about
double that number existing on the whole earth.[33] Of these
Europe, Asia, and Africa, furnish suppose one hundred and twenty-
six; that is, the twenty-six common to Europe and America, and
about one hundred which are not in America at all. The American
species, then, are to those of the rest of the earth, as one hundred
to one hundred and twenty-six, or four to five. But the residue of the
earth being double the extent of America, the exact proportion
would have been but as four to eight.
Hitherto I have considered this hypothesis as applied to brute
animals only, and not in its extension to the man of America,
whether aboriginal or transplanted. It is the opinion of Mons. de
Buffon that the former furnishes no exception to it.[34]
"Quoique le sauvage du nouveau monde soit à peu près de même stature que l'homme de notre
monde, cela ne suffit pas pour qu'il puisse faire une exception au fait général du rapetissement de
la nature vivante dans tout ce continent; le sauvage est foible et petit par les organes de la
génération; il n'a ni poil, ni barbe, and nulle ardeur pour sa femelle. Quoique plus léger que
l'Européen, parce qu'il a plus d'habitude à courir, il est cependant beaucoup moins fort de corps; il
est aussi bien moins sensible, et cependant plus craintif et plus lâche; il n'a nulle vivacité, nulle
activité dans l'ame; celle du corps est moins un exercise, un mouvement volontaire qu'une
nécessité d'action causée par le besoin; ôtez lui la faim et la soif, vous détruirez en même tems le
principe actif de tous ses mouvemens; il demeurera stupidement en repos sur ses jambes ou
couché pendant des jours entiers. Il ne faut pas aller chercher plus loin à cause de la vie dispersée
des sauvages et de leur éloignement pour la société; la plus précieuse étincelle du feu de la nature
leur a été refusée; ils manquent d'ardeur pour leur femelle, et par consequent d'amour pour leur
semblables; ne connoissant pas l'attachment le plus vif, le plus tendre de tous, leurs autres
sentimens de ce genre, sont froids et languissans; ils aiment foiblement leurs pères et leurs enfans;
la société la plus intime de toutes, celle de la même famille, n'a donc chez eux que de foibles liens;
la société d'une famille à l'autre n'en a point de tout; dès lors nulle réunion, nulle république, nulle
état social. La physique de l'amour fait chez eux le moral des mœurs; leur cœur est glacé, leur
societé et leur empire dur. Ils ne regardent leurs femmes que comme des servantes de peine ou
des bêtes de somme qu'ils chargent, sans ménagement, du fardeau de leur chasse, et qu'ils
forcent, sans pitié, sans reconnoissance, à des ouvrages qui souvent sont au dessus de leurs
forces; ils n'ont que peu d'enfans; ils en out peu de soin; tout se ressent de leur premier defaut; ils
sont indifférents parce qu'ils sont peu puissants, et cette indifference pour le sexe est la tache
originelle qui flétrit la nature, qui l'empeche de s'épanouir, et qui detruisant les germes de la vie,
coupe en même temps la racine de société. L'homme ne fait donc point d'exception ici. La nature
en lui refusant les puissances de l'amour l'a plus maltraité et plus rapetissé qu'aucun des animaux."
An afflicting picture, indeed, which for the honor of human nature, I am glad to believe has no original.
Of the Indian of South America I know nothing; for I would not honor with the appellation of
knowledge, what I derive from the fables published of them. These I believe to be just as true as the
fables of Æsop. This belief is founded on what I have seen of man, white, red, and black, and what has
been written of him by authors, enlightened themselves, and writing among an enlightened people. The
Indian of North America being more within our reach, I can speak of him somewhat from my own
knowledge, but more from the information of others better acquainted with him, and on whose truth
and judgment I can rely. From these sources I am able to say, in contradiction to this representation,
that he is neither more defective in ardor, nor more impotent with his female, than the white reduced to
the same diet and exercise; that he is brave, when an enterprise depends on bravery; education with
him making the point of honor consist in the destruction of an enemy by stratagem, and in the
preservation of his own person free from injury; or, perhaps, this is nature, while it is education which
teaches us to[35] honor force more than finesse; that he will defend himself against a host of enemies,
always choosing to be killed, rather than to surrender,[36] though it be to the whites, who he knows will
treat him well; that in other situations, also, he meets death with more deliberation, and endures
tortures with a firmness unknown almost to religious enthusiasm with us; that he is affectionate to his
children, careful of them, and indulgent in the extreme; that his affections comprehend his other
connections, weakening, as with us, from circle to circle, as they recede from the centre; that his
friendships are strong and faithful to the uttermost[37] extremity; that his sensibility is keen, even the
warriors weeping most bitterly on the loss of their children, though in general they endeavor to appear
superior to human events; that his vivacity and activity of mind is equal to ours in the same situation;
hence his eagerness for hunting, and for games of chance. The women are submitted to unjust
drudgery. This I believe is the case with every barbarous people. With such, force is law. The stronger
sex imposes on the weaker. It is civilization alone which replaces women in the enjoyment of their
natural equality. That first teaches us to subdue the selfish passions, and to respect those rights in
others which we value in ourselves. Were we in equal barbarism, our females would be equal drudges.
The man with them is less strong than with us, but their women stronger than ours; and both for the
Welcome to our website – the ideal destination for book lovers and
knowledge seekers. With a mission to inspire endlessly, we offer a
vast collection of books, ranging from classic literary works to
specialized publications, self-development books, and children's
literature. Each book is a new journey of discovery, expanding
knowledge and enriching the soul of the reade
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