# On the Classification of Colored Textures From a Texture-Ranking Experiment: Observers Ability of Discrimination Quantification

Amadou Sawadogo, Dominique Lafon-Pham, Simplice Dossou-Gbété

<figure><img src="image_1.png" /></figure>

## To cite this version:

Amadou Sawadogo, Dominique Lafon-Pham, Simplice Dossou-Gbété. On the Classification of Colored Textures From a Texture-Ranking Experiment: Observers Ability of Discrimination Quantification. International Journal of Statistics and Probability, 2020, 10 (1), pp.1. 10.5539/ijsp.v10n1p1. hal-04662050

HAL Id: hal-04662050

https://imt-mines-ales.hal.science/hal-04662050v1

Submitted on 25 Jul 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# On the Classification of Colored Textures From a Texture-Ranking Experiment: Observers Ability of Discrimination Quantification

Amadou Sawadogo<fnref n="1" />, Dominique Lafon<fnref n="2" />, Simplice Dossou-Gbété<fnref n="3" />

Correspondence: Amadou Sawadogo, Department of Mathematics and Computer, University of Félix Houphouët Boigny, Postal address 22 BP 582 Abidjan 22, Ivory Coast. E-mail: amadou.sawadogo@gmail.com

Received: January 14, 2020

Accepted: November 9, 2020

Online Published: November 30, 2020

doi:10.5539/ijsp.v10n1p1

URL: https://doi.org/10.5539/ijsp.v10n1p1

## Abstract

This paper is devoted to the construction of Fechnerian scales on a physical dimension of investigated colored textures. For this purpose, we considered the extension of the Mallows-Bradley-Terry model for the analysis of the data collected from a contrast-sorting experiment. A likelihood ratio test procedure was proposed in order to choose between the two following hypotheses: discrimination and non-discrimination between the investigated stimuli. In addition, post-hoc analyzes allowed us to find out which of the stimuli differ from the others. Our findings indicate that the subjective attribute of visual contrast appears to be a psycho physical scale that maps to the physical scale corresponding to the Michelson contrast. Mainly, the estimates of the model index of discrimination parameter of the stimuli show that the ability of the observers to discriminate between the textures according to the visual contrast varies with respect to the color ranges and the textures types. According to the luminance contrasts ability of discrimination, the Isotropic texture type is the best, followed by the Random-dots texture type, then by the Horizontal grating type and the Vertical grating type is the least. The Fechnerian scales on the physical dimension of the Michelson contrast of the colored textures depend on the chromaticness of the colored textures phases and the texture types. The psycho physical method of identification would be the best when determining the related thresholds.

**Keywords:** Fechnerian scale, physical scale, colored texture, rank data, Michelson contrast, visual contrast, bridge sampling, likelihood ratio test, post-hoc test

## 1. Introduction

In contrast with the past, sensory properties of materials are subject to growing attention both in an hedonic and utilitarian point of view. In industrial sectors such as luxury, cosmetics or transports, the visual aspect of objects is a criterion of evaluation and decision in customer's choices. Effect colors (metallic and pearlescent) are now commonly used. Depending on the angle of view, they produce different visual feelings.

The classical color descriptors (Fairchild, 2013; Hunt & Pointer, 2011), e.g., CIE Lab, CIE Luv, CIECAM02 etc, do not provide a complete representation of these effect colors and more generally on surfaces exhibiting a colored texture. It is now clearly established that our visual system is concerned not only with individual colors but also with contrasts. However, it remains unclear how to characterize visual similarities between textures belonging to the same type.

The purpose of this study was to construct Fechnerian scales on a physical dimension of stimuli capable of describing the behavior of the average observer with respect to the perceived contrast of colored textures types. The similarity relations between isochromatic colored textures poorly different regarding to their luminance level were investigated. To this end, the Michelson contrast is considered to be the physical continuum in the physical domain of the colored textures. The colored textures consisted of four types of texture (Random-dots textures, Isotropic textures, Horizontal gratings and Vertical grating) and three color ranges (red, green and yellow). Then, we conducted a softcopy experiment consisting in ranking the colored textures according to the attribute of visual contrast. Ties were permitted in the ranking experiment. The data thus collected are rank data with ties.

In the literature, many models of rank data have been proposed over the years. They could be categorized into four classes (Marden, 1995; Fligner & Verducci, 1986; Fligner & Verducci, 1988): (i) paired-comparison models, (ii) distance based-models, (iii) multistage models and (iv) order-statistics models. Applications of the models to real life situations can be found in the literature; for example, in (Marden, 1995; Critchlow & Fligner, 1991). In the statistical modeling of ranking

<footnote n="1">¹ Department of Mathematics and Computer, University of Félix Houphouët Boigny, Abidjan, Ivory Coast</footnote>

<footnote n="2">² Ecole des Mines d'Alès, Site Hélioparc, Pau, France</footnote>

<footnote n="3">³ Université de Pau et des Pays de l'Adour, E2S UPPA, CNRS, LMAP, PAU, France. International Chair in Mathematical Physics and Applications (ICMPA-UNESCO Chair), University of Abomey-Calavi, Benin</footnote>
data, the paired comparison approach is widely used and relies upon a clear formulation of probabilities distribution models. The basic idea behind this approach, in particular the model considered for the analysis of the rank data is that, the observer reports their preferences only after having an unambiguous ranking, by starting with paired comparisons, with possibility of ties. In other words, the multiple pairwise judgments are consistent. Consistency means transitivity. Hence, each pair of stimuli will have a little triple of probability attached. Rao and Kupper (1967) and Davidson (1970) suggested particular parametrizations of the paired comparisons when extending the Bradley-Terry model to account for ties. The parametrization of the paired-comparisons assumed were Davidson's model (Davidson, 1970). The latter is defined as follows: choose either stimulus *i* or stimulus *j* respectively with probabilities $\theta_{ij} = \pi_i/(\pi_i + \pi_j + v\sqrt{\pi_i\pi_j})$, $\theta_{ji} = \pi_j/(\pi_i + \pi_j + v\sqrt{\pi_i\pi_j})$ and do not choose with probability $\theta_{ij2} = v\sqrt{\pi_i\pi_j}/(\pi_i + \pi_j + v\sqrt{\pi_i\pi_j})$ with $v \ge 0$, where $\pi_i$ denotes the "worth", an index of relative preference, of the *i*th stimulus and $v \ge 0$ an index of discrimination which is, in our context, interpreted as a measure of difficulty of the population of observers to discriminate between the set of stimuli. From this point of view, we made use of the extension (Sawadogo et al., 2017) of the Mallows-Bradley-Terry ranking model for one block comparison consisting of all the stimuli of interest to analyze the collected ranking data.

This paper is organized as follows. Section 2 describes the apparatus, the investigated stimuli, the raw data set and indicates how they had been collected. Section 3 recalls the generalization of the Mallows-Bradley-Terry ranking model. In section 4, a likelihood ratio test has been first performed to answer the question of discernability of the stimuli. Secondly, post-hoc analyzes based on post-hoc tests have been applied to find out which stimuli differ from the others. Thirdly, we have addressed the construction of Fechnerian scales on the physical dimension of stimuli using different psycho physical methods and indicate how such scales differ. A discussion has been carried out where our findings are compared with recent studies in Section 5.

The data used for this study were collected by a softcopy experiment.

## 2. Apparatus, Stimuli and Procedure

### 2.1 Apparatus

Stimuli were presented on a color CRT monitor (Philips) with a refresh rate of 85 hertz and using the resolution 1600×1200 pixels. Prior to the experiments we had characterized and calibrated accurately the CRT monitor in order to control the colorimetric attributes of the displayed colors. The white point of the monitor has been set to the standard illuminant $D_{65}$. A colorimetric calibration with a PR650-spectroradiometer (Photo Research Inc. USA) of the CRT monitor allowed to display colored textures corresponding to predetermined colors physically expressed in the standard CIE1931 (Commission internationale de l'éclairage), notation (X, Y, Z). The observers sat 60 centimeters away from the screen. The experimental environment has been strictly controlled. It was important to minimize spatial color interferences on the perception of the stimuli displayed on the monitor screen. The experiment were carried out under diffuse lighting conditions. The light source in the room was a courtesy light with a fluorescent tube (Just Normlicht Color Control Daylight 5000 kelvins) and a diffuser of light. The experiment took place in a blind room with the same neutral grey color for walls, ceiling and floor in order to avoid possible local variations of the ambient light. The monitor background was hidden by grey surfaces in order to reduce visual distractions (only the screen and the mouse were distinguishable from the background).

### 2.2 Stimuli and Contrast-Sorting Experiment

The experiment consisted in ranking isochromatic textures by a total of forty observers. The observers had a normal or corrected-to-normal visual acuity and a normal color vision. They were informed about the final purpose of the experiment. Prior to the start of the experiment, each observer had been trained in order to have him or her understand what we mean by visual contrast: we used series of colored textures of the same natural scene that exhibit different levels of visual contrast. After making sure that each observer understood the task, the experiment began. Each of the forty observers completed 12 trials. An isochromatic texture is characterized by its color range and the type of texture it belongs to. A color belonging to a given color range is described by its chromaticity coordinates and luminance (x, y, Y) in the well-known three-dimensional CIE XYZ color space (Bonton et al., 2000): (x, y) denotes its chromaticity coordinates and Y its luminance. A colored texture belonging to a given type of texture consists of two isochromatic colors i.e., described by the two tristimulus coordinates ($x_0, y_0, Y_0$) and ($x_0, y_0, Y_0 + \Delta Y$): a luminance value $Y_0$ for each color range has been chosen as the reference for each color. Then we applied a luminance discrepancy $\Delta Y = Y - Y_0$. A Michelson contrast level M of a given colored texture is defined by $M = |\Delta Y|/(Y_0 + 2 * |\Delta Y|)$. A sample of the textures types is presented in Figure 1. A trial consisted in displaying simultaneously on the CRT monitor screen, 20 colored textures belonging to the same type of texture pattern and the same color range; then ranking them according to the appreciation of the visual contrast perceived between the two components of each colored texture (from the lowest visual contrast corresponding to a very little difference between the luminances of the two components to the highest visual contrast). The observers were asked to use integer from 1 to 20; rank 1 is assigned to the lowest visual contrast texture. Ties were permitted and expressed by assigning the same rank to the stimuli exhibiting the same visual contrast according to the
observers. The colored textures were randomly displayed on a uniform gray field for each trial. After an observer had completed a trial, he or she keyed his responses and finally clicked to ask for a new trial. The observers were free to use as much time as they found necessary for the set of 12 trials. For the 12 trials, the textures were presented to in a pseudo random order in order to avoid a familiarization: two successive trials present to the observers are necessarily composed of stimuli with different color ranges and different textures types. During the ranking experiment, it appeared that each observer took an average of 90 minutes for a complete round. The textures consisted of square uncompressed Bitmap images (200×200 pixels) with a color depth of 24 bits per pixel. They were displayed as a square of 5 centimeters × 5 centimeters without any resizing (5 degrees). Three color ranges (red, green and yellow) and four types of texture pattern (Random-dots texture, Isotropic texture, Horizontal grating and Vertical grating texture) had been investigated. The chromaticity coordinates and luminance of the considered color ranges in each color range are given in Table 1.

<table><caption>Table 1. The chromaticity coordinates and luminance of the three color ranges considered in the CIE 1931 standard color space</caption><thead><tr><th>Tristimulus values</th><th>Red</th><th>Green</th><th>Yellow</th></tr></thead><tbody><tr><td>x<sub>0</sub></td><td>0.51</td><td>0.29</td><td>0.39</td></tr><tr><td>y<sub>0</sub></td><td>0.33</td><td>0.54</td><td>0.49</td></tr><tr><td>Y<sub>0</sub> (candela per square meter)</td><td>15.56</td><td>20.28</td><td>84.83</td></tr><tr><td>λ (nanometers)</td><td>675</td><td>550</td><td>570</td></tr></tbody></table>

<figure><img src="image_3.png" /><figcaption>Figure 1. The investigated types of colored textures consisting of isochromatic colors: equal chromaticity coordinates (x; y) and different luminance coordinates Y</figcaption></figure>

The textures are identified according to their color ranges and the type of textures they belong to by two letters: the first letter denotes the color range and the latter the type of the textures. The Red, Green and Yellow color ranges are respectively identified by the following letters R, G and Y. The types of textures, namely Random-dot textures, Isotropic textures, Horizontal gratings and Vertical gratings are respectively identified by the letters R, I, H and V. Therefore, the indicator GR denotes the set of colored textures in the Green color range belonging to the Random-dots texture type; any colored texture of the set GR is identified by the prefix GR.

The collected data had been organized as follows: the data set consisted in all the elements of the set $\{R, G, Y\} \times \{R, I, H, V\}$, each with 40 rows (the observers) and 20 columns (the colored textures).

Finally, each data is described by the observers identifiers, the color range, the type of texture, the luminance discrepancies
ΔY in colored textures, the Michelson contrast M levels associated to the luminance discrepancies and the ranks assigned by the observers to the colored textures.

The statistical model used for the analysis of the ranking data is described and defined in the section below.

### 3. The Extension of the Mallows-Bradley-Terry Ranking Model

As mentioned in Introduction, the paired comparison approach is the one considered. The basic idea behind this approach is that the observer reports their preferences only after having an unambiguous ranking, by starting with paired comparisons. We mean by unambiguous ranking, transitivity in the set of paired comparison.

Suppose that the paired comparison experiment involves *q* stimuli, where independent comparisons are made for the pair of stimuli {i, j}. It is assumed that ties are not permitted, and that the order of item presentation is unimportant. The basic parameters are then the $q(q-1)/2$ quantities $\theta_{ij}$, where for $i < j$, $\theta_{ij}$ denotes the probability that *i* is preferred to *j* in a paired comparison of these two stimuli. The Bradley-Terry model is given by $\theta_{ij} = \pi_i/(\pi_i + \pi_j)$ where $\pi_i > 0$ is the probability that the stimulus *i* is ranked lowest when the entire set of *q* stimuli are submitted for ranking.

In the most well studied case, there is simply one block consisting of all the stimuli of interest, and the resulting comparison is just a complete ranking of these stimuli, provided by each of *n* observers. To model the probability of such a ranking *r* of *q* stimuli, Mallows (1957) proposes a natural extension of the Bradley-Terry paired comparison model. According to the Mallows-Bradley-Terry ranking model, the probability of the rank vector *r*

$$ \mathbb{P}(\mathbf{r}; \boldsymbol{\pi}) = c(\boldsymbol{\pi}) \prod_{j=1}^{q} \pi_j^{q-r(j)}, \qquad (1) $$

where $c(\boldsymbol{\pi})$ stands for the normalizing constant. This result is well-known and called the Mallows-Bradley-Terry ranking model (Critchlow & Fligner, 1991).

Sawadogo et al. (2017) adapted the Mallows-Bradley-Terry ranking model to ties by considering the Davidson (1970) extension of the Bradley-Terry model. To account for ties in a paired comparisons experiment, Davidson proposed an extension of the Bradley-Terry model as follows: choose either stimulus *i* or stimulus *j* respectively with probabilities $\theta_{ij} = \pi_i/(\pi_i + \pi_j + v \sqrt{\pi_i\pi_j})$, $\theta_{ji} = \pi_j/(\pi_i + \pi_j + v \sqrt{\pi_i\pi_j})$ and do not choose with probability $\theta_{ij2} = v \sqrt{\pi_i\pi_j}/(\pi_i + \pi_j + v \sqrt{\pi_i\pi_j})$ with $v \ge 0$, where $\pi_i$ denotes the "worth", an index of relative preference, of the *i*th stimulus and $v \ge 0$ an index of discrimination which is, in our context, interpreted as a measure of difficulty of the population of observers to discriminate between the set of stimuli. Notice that the Bradley-Terry model is obtained from the Davidson's when $v = 0$.

Let us consider a set of *q* items labeled by distinct numbers from 1 to *q*. The *q* items are submitted to an observer for ranking with ties. A ranking is obtained by making independently all the pairwise comparisons until having an unambiguous ranking. The resulting ranks vector *r* is a finite sequence of integers of length *q* where each value belongs to the set of consecutive integers from 1 to *q*. Then the *j*-th component of *r*, say *r*(*j*), is the rank assigned to the item labeled *j*. In what follows, we assume the following definitions which will be used to define the extension of the Mallows-Bradley-Terry ranking model that account for ties.

**Definition 1** Given a ranking $\mathbf{r} = \{\mathbf{r}(j)\}_{j=1,2,\dots,q}$ of the set $\mathcal{A}(q)$ of all rankings of $q$ stimuli in which ties are possible and $\mathbf{r}_+ = \max_{j=1,2,\dots,q} \mathbf{r}(j)$. The pattern of ties associated with the ranking with ties $\mathbf{r}$, denoted $\lambda_r$, is a vector of length $\mathbf{r}_+$ defined by

$$ \lambda_r = (\lambda_r(1), \dots, \lambda_r(k), \dots, \lambda_r(\mathbf{r}_+)), $$

where $\lambda_r(k) = \text{card}\{\mathbf{r} = k\}$ denotes the cardinality of the set $\{\mathbf{r} = k\} = \{j \in \{1, 2, \dots, q\}; \mathbf{r}(j) = k\}, \forall k = 1, 2, \dots, \mathbf{r}_+.$

**Definition 2** Given a ranking $\mathbf{r} = \{\mathbf{r}(j)\}_{j=1,2,\dots,q}$ of the set $\mathcal{A}(q)$. The mid-rank of stimulus *i*, denoted $\mathbf{r}^*(i)$, is defined by

$$ \mathbf{r}^*(i) = \frac{\lambda_r\{\mathbf{r}(i)\} + 1}{2} + \sum_{k=1}^{r(i)-1} \lambda_r(k), \quad \forall i \in \{1, \dots, q\}. $$

The extension of the Mallows-Bradley-Terry model for ranking with ties is defined as follows.

**Definition 3** The generalization of the Mallows-Bradley-Terry model for ranking with ties is defined by the following probability mass function,

$$ \mathbb{P}(\mathbf{r}; \boldsymbol{\pi}, v) = c(\boldsymbol{\pi}, v)v^{\frac{1}{2}\left\{\sum_{k=1}^{r_+} \lambda_r^2(k) - q\right\}} \prod_{j=1}^{q} \pi_j^{q-r^*(j)}, \forall \mathbf{r} \in \mathcal{A}(q) \qquad (2) $$
with

$$ \log\{c(\pi, v)\} = -\log\left[\sum_{r} v^{\frac{1}{2}} \left\{\sum_{k=1}^{r_{+}} \lambda_{r}^{2}(k) - q\right\} \prod_{j=1}^{q} \pi_{j}^{q-r_{*}(j)}\right], $$

where $\mathcal{A}(q)$ denotes the sample space.

The parameter $\pi_j, j \in \{1, 2, \dots, q\}$ denote the "worth", an index of relative preference, of the j-th stimulus, $\pi_j > 0, \sum_{j=1}^q \pi_j = 1$. Whereas, the parameter $v$ is an index of discrimination and interpreted as a measure of difficulty of the complete population of the observers to discriminate between the stimuli being compared. One can remark that if $r$ is a ranking without ties, the model given by Equation (2) is exactly the classical Mallows-Bradley-Terry ranking model defined above by Equation (1) because one will have $r_+ = q$ and $\lambda_r(k) = 1, \forall k = 1, \dots, q$ leading to $r^*(i) = r(i), \forall i = 1, \dots, q$. The Fechnerian scale provides a continuous unidimensional dimensionless variable for preference. It expresses the degree of preference that an experimental observer shows.

The structural property $\pi_j > 0, \sum_{j=1}^q \pi_j = 1$ on the parameters $\pi_j, j = 1, 2, \dots, q$ allows re-parameterization of the model as follows: $\theta_j = \log(\pi_j/\pi_q)$ for $j = 1, 2, \dots, q$ and $\gamma = \log(v)$ with $v > 0$. Following the re-parameterization stated above, the Mallows-Bradley-Terry model for ranking with ties given by Equation (2) can be rewritten in a manner given by the definition, viz.

**Definition 4** The Mallows-Bradley-Terry model for ranking with ties is defined as follows

$$ \mathbb{P}(\mathbf{r}; \theta, \gamma) = c(\theta, \gamma) \exp \left[ \frac{1}{2} \left\{ \sum_{\alpha=1}^{r_+} \lambda_r^2(\alpha) - q \right\} \gamma + \sum_{j=1}^{q} \left\{ q - \mathbf{r}^*(j) \right\} \theta_j \right], \forall \mathbf{r} \in \mathcal{A}(q) \quad (3) $$

where $c(\theta, \gamma) = \left( \sum_s \exp \left[ \frac{1}{2} \left\{ \sum_{\alpha=1}^{s_+} \lambda_s^2(\alpha) - q \right\} \gamma + \sum_{j=1}^{q} \left\{ q - s^*(j) \right\} \theta_j \right] \right)^{-1}$ denotes the normalizing constant.

The advantage of this re-parameterization is that one obtains a regular curved exponential family model (David & Mark, 2005) with $q$ dimensions and the statistic $[(q - \mathbf{r}^*(j))_{j \in \{1, 2, \dots, q-1\}}$, $1/2\{\sum_{\alpha=1}^{r_+} \lambda_r^2(\alpha) - q\}$] is a sufficient statistic for the model parameter $(\theta, \gamma)$ vector.

Throughout the remainder of the paper, we will use the extension of the Mallows-Bradley-Terry ranking model in the form given by Equation (3).

**Proposition 1** The log-likelihood function $\mathcal{L}_n(\theta, \gamma)$ of the canonical parameters $(\theta, \gamma)$, given an iid sample of rankings with ties $\mathbf{r}_i, i \in \{1, 2, \dots, n\}$ is

$$ \mathcal{L}_n(\theta, \gamma) = \sum_{j=1}^{q} \theta_j \sum_{i=1}^{n} \left\{ q - \mathbf{r}_i^*(j) \right\} + \frac{1}{2} \gamma \left[ -nq + \sum_{i=1}^{n} \left\{ \sum_{\alpha=1}^{r_{i+}} \lambda_{\mathbf{r}_i}^2(\alpha) \right\} \right] + n \log\{c(\theta, \gamma)\}. \quad (4) $$

# 4. Estimation of the Stimuli Worth Parameters and Population Index of Discrimination

## 4.1 Likelihood Ratio Testing: Discrimination Versus Non-discrimination

From the Davidson's model stated above, the probability of do not choose $\theta_{i/2} = v \sqrt{\pi_i \pi_j}/(\pi_i + \pi_j + v \sqrt{\pi_i \pi_j})$ with $v \ge 0$ tends to 1 when $v$ tends to infinity. The setting $v = 0$ corresponds to a certain discrimination. The non-discrimination means that there is an equal probability of choosing either one or the other of both the two stimuli, whatever the concerned stimuli; which leads to $v > 0, \forall j \in \{1, 2, \dots, q\}, \pi_j = 1/q$. Thus, the question of discrimination of the stimuli is answered by making an hypothesis testing: the hypothesis of discrimination of the textures versus the hypothesis of non-discrimination of the textures. We were interested in testing the null hypothesis $H_0: \pi_i = 1/q, \forall i \in \{1, 2, \dots, q\}, v > 0$ which corresponds to non-discrimination between the stimuli, versus the alternative hypothesis $H_1: \exists i; \pi_i \ne 1/q, v > 0$ which means discrimination between them. A likelihood ratio test based on the ratio

$$ \Lambda_n = \frac{\max_{\gamma} \mathcal{L}_n(\theta_0, \gamma)}{\mathcal{L}_n(\tilde{\theta}_0, \tilde{\gamma})}, \quad (5) $$

where $\mathcal{L}_n(\theta, \gamma) = \exp\{l_n(\theta, \gamma)\}$ denotes the likelihood of the model parameter vector $(\theta, \gamma)$, is obtained by using $-2 \log \Lambda_n$ as the test statistic. Under $H_0$, the statistic $-2 \log \Lambda_n$ has an approximate $\chi^2(q-1)$ distribution with $q-1$ degrees of freedom when $n$ tends to infinity. The null hypothesis $H_0$ is rejected if $-2 \log \Lambda_n > \chi_{1-\alpha}^2(q-1)$ where the quantities $\Lambda_n$ and
$\chi_{1-\alpha}^2(q-1)$ are the statistic $\Lambda_n$ value and the $(1-\alpha)$ quantile of the chi-square distribution with $q-1$ degrees of freedom respectively. Equivalently, $H_0$ is rejected if $\alpha > p$-value.

**Proposition 2** Given $(\hat{\theta}, \hat{\gamma})$ (resp. $(\theta_0, \hat{\gamma}_0)$) the Maximum Likelihood Estimates of the canonical parameter vector $(\theta, \gamma)$ (resp. $(\theta_0, \gamma)$ under $H_0$). Then

$$ \log \Lambda_n = n \log \left\{ \frac{c(\hat{\theta}, \hat{\gamma})}{c(\theta_0, \hat{\gamma}_0)} \right\} + \sum_{i=1}^{n} \left[ \frac{1}{2} \sum_{\alpha=1}^{r_{i,+}} \lambda_{r_i}^2(\alpha) - q \right] \left\{ \hat{\gamma}_0 - \hat{\gamma} \right\} + \sum_{j=1}^{q} \left\{ q - \mathbf{r}_i^*(j) \right\} \left\{ \hat{\theta}_{0,j} - \hat{\theta}_j \right\}. $$

It results from $\theta_{0,j} = 0, \forall j \in \{1, 2, \dots, q\}$ that

$$ \log \Lambda_n = n \log \left\{ \frac{c(\hat{\theta}, \hat{\gamma})}{c(\theta_0, \hat{\gamma}_0)} \right\} + \sum_{i=1}^{n} \left[ \frac{1}{2} \sum_{\alpha=1}^{r_{i,+}} \lambda_{r_i}^2(\alpha) - q \right] \left\{ \hat{\gamma}_0 - \hat{\gamma} \right\} - \sum_{j=1}^{q} \left\{ q - \mathbf{r}_i^*(j) \right\} \hat{\theta}_j, \quad (6) $$

given $\mathbf{r}_1, \mathbf{r}_2, \dots, \mathbf{r}_n$.

The computation of the ratio $c(\hat{\theta}, \hat{\gamma})/c(\theta_0, \hat{\gamma}_0)$ of the normalizing constants involved in Equation (6) is required and involved sums with a huge number of terms. There is no closed form analytical expression for the calculation of this ratio. To overcome this difficulty, we use the path sampling technique; see (e.g., Hunter and Handcock, 2012. This technique is a particular case of a Monte Carlo method called bridge sampling (e.g., Meng and Wong, 1996). Hereafter, the following results describe how it works in our case.

**Lemma 1** Let $\varphi: [0, 1] \rightarrow \mathbb{R}^q$ be a smooth mapping defined as $\varphi(u) = (\theta(u), \gamma(u))$, $\forall u \in [0, 1]$ such that $\varphi(0) = (\theta_0, \hat{\gamma}_0)$ and $\varphi(1) = (\hat{\theta}, \hat{\gamma})$. Then,

$$ \mathbb{E}_{\varphi(u)} \left\{ \frac{d}{du} \log \mathbb{P}(\mathbf{r}; \varphi(u)) \right\} = \frac{d}{du} \sum_{\mathbf{r}} \mathbb{P}(\mathbf{r}; \varphi(u)) = 0. $$

**Proposition 3** One has,

$$ \log \left\{ \frac{c(\varphi(1))}{c(\varphi(0))} \right\} = \mathbb{E} \left\{ \mathrm{\textit{Z}} \cdot \left( \frac{d\varphi}{du} \circ U \right) \right\}, \quad (7) $$

where

$$ \mathrm{\textit{Z}}(\mathbf{r}) = \left[ \left\{ q - \mathbf{r}^*(j) \right\}_{j \in \{1, 2, \dots, q-1\}}, \frac{1}{2} \left\{ \sum_{k=1}^{r_+} \lambda_r^2(k) - q \right\} \right] $$

is a sufficient statistic for the parameter vector $\varphi(u) = (\theta(u), \gamma(u))$. The expectation $\mathbb{E}$ is taken with respect to the joint distribution of the random vector $(U, \mathbf{r})$, where $U$ is uniform $[0, 1]$ and $\mathbf{r} \mid U = u$ is distributed according to $\mathbb{P}(\mathbf{r}; \varphi)$ with parameter $\varphi(U)$.

**Corollary 1** Equation (7) suggests that $\log\left\{\frac{c(\varphi(1))}{c(\varphi(0))}\right\}$ could be estimated by drawing a sample $(u_1, \mathbf{r}_1), \dots, (u_N, \mathbf{r}_N)$ from the joint distribution $(U, \mathbf{r})$, then calculating the sample average

$$ \frac{1}{N} \sum_{k=1}^{N} \mathrm{\textit{Z}}(\mathbf{r}_k) \frac{d\varphi}{du}(u_k). $$

In our implementation of the path sampling method, we considered $\varphi(u) = (\theta_0, \hat{\gamma}_0) + u \cdot [(\hat{\theta}, \hat{\gamma}) - (\theta_0, \hat{\gamma}_0)]$ which corresponds to a linear path from $(\theta_0, \hat{\gamma}_0)$ to $(\hat{\theta}, \hat{\gamma})$.

## 4.2 Rank Data Analysis and Fechnerian Scales Constructing

The data analysis by means of a fixed effects additive model showed that the visual contrast is a psychophysical scale corresponding to the physical scale determined by the Michelson contrast of the colored textures (Sawadogo et al., 2014). Furthermore, one has graphically noticed that the ability of the observers to discriminate between the colored textures according to the visual contrast varies according to the color ranges and the textures types. But, we have quantified neither the ability of discrimination of the set of the observers nor established a classification of the texture types and/or the color ranges. In the present paper, our objective is to provide some answers to these questions and make assumptions about the related Fechnerian sensory scales.

The model parameter vector $(\pi, \nu)$ has been estimated by means of Monte Carlo MM-algorithm method; see, Sawadogo et al. (2017). A bootstrap-t confidence interval has been implemented for the calculation of the confidence intervals
$CI_{1-\alpha}(\pi, \nu)$ of the model parameter vector $(\pi, \nu)$ at the significance level of $\alpha = 5\%$. For this purpose, we wrote some own codes using the R software (R Core Team, 2020). The results of the model parameter vector $(\pi, \nu)$ estimates are available in Appendix 2. In multiple regression, the multiple correlation coefficient $R^2$ is often used to measure the percentage of variation in the data that is explained by the regression model. An analogous quantity, in our context, is the percentage of nonuniformity $R^2$ in the data that is explained by the model, see, e.g., Marden, 1995. See Appendice 3 for the definition of the $R^2$ and its values in our data. The results show that the model fits well all the data.

Figure 2, Figure 3 and Figure 4 depict the variations of the Bradley-Terry-Luce parameters $(\pi_j)_{j=1,2,\dots,q}$ of the textures with respect to its Michelson contrast values (M). The curves in the mentioned Figures are obtained by making use of the loess() function of R Software. The estimates of the colored textures worth parameters $\pi_j, j \in \{1, 2, \dots, q\}$ increase when its corresponding Michelson contrast levels decrease whatever the signs of their luminance discrepancies $\Delta Y$. Thus, an ordering relation based on the Michelson contrasts of the colored textures always corresponds to the ordering obtained during the psychophysical experiments whatever the type of colored texture and the color range. This suggests to consider the visual contrast as a psycho physical scale corresponding to the physical scale determined by the textures' Michelson contrasts M. This result is coherent with one found in (Sawadogo et al, 2014).

The likelihood ratio tests of $H_0$ versus $H_1$ result in the reject of $H_0$ at a significance level of $\alpha = 0.05$ for all the observed data since the $p$-values (see, Appendice 4) are all less than $\alpha = 0.05$. It comes that, all the observers were capable to discriminate between all the colored textures according to the visual contrast at a significance level of 5%.

From the likelihood ratio tests, significant differences in luminance contrasts are detected between the colored textures whatever the color range and the texture type. Given a color range and a texture type, it is necessary to find out which colored textures differ from the others (Dzhafarov & Perry, 2014). Therefore, the hypotheses to test are $H_0' : m_k = m_j, k \neq j$ vs $H_1' : m_k \neq m_j$, where $m_k$ represents the medians of the stimuli $k$. A Conover post-hoc test (Conover, 1999) with $p$-values adjustment according to the Holm procedure (Holm, 1979) has been used. To this end, we make use of the PMCMRplus (Pohlert, 2019) package available in R software. The results of the multiple paired comparisons tests are available in Appendix 5: box-whisker plot for each of the colored textures belonging to a given color range and texture type is plotted. The range of the whiskers indicate the extremes. Letter symbols are depicted on top of each box. Different letters indicate significant differences between groups determined using $p$-values at the selected level of $\alpha = 5\%$.

The estimates of the model index of discrimination parameter $v$ differ from a colored texture to another whatever the type of the color range(see, Figure 5). This explains that the observers ability to discriminate between the colored textures according to the visual contrast varies according to the color ranges and the texture types. The probability of do not choose tends to 1 when $v$ tends to infinity; this probability tends to 0 when $v$ tends to 0: the smaller the values of $v$, the better the ability of discrimination between the textures. The estimates of the index of discrimination $v$ are generally small in the green color range, indicating that discrimination is best in this color range. Moreover, one can classify either the textures with respect to the pattern they belong to for a given color range or the color ranges given the texture type basing on the index of discrimination $v$. For instance, in the green color range, the Isotropic textures are the best discriminated with a value of $\tilde{v} = 0.6027$, secondly the Random-dot textures ($\tilde{v} = 0.7273$), thirdly the Horizontal gratings ($\tilde{v} = 1.2114$) and finally Vertical gratings ($\tilde{v} = 1.3243$). This holds true for both of the red color range and the yellow color ones. Whatever the texture type, the sense of perception of the set of observers which determined their ability of discrimination is sharper in the green color range than the other color range. Then, follows the red color range and finally the yellow color range, except in the type of Isotropic textures. In this texture type, the yellow color range comes before the red one.

The construction of Fechnerian scales on physical dimensions of stimuli relies on the concept of threshold, namely absolute threshold and discrimination threshold. Fundamentally, there are two kinds of task that are used to obtain thresholds (Dzhafarov & Colonius, 1999; Pelli & Farell, 1995): classification and adjustment. In the former, the experimenter controls the stimulus and the observer makes a judgment based on the resulting percept whereas in the latter the observer adjusts the stimulus to satisfy a perceptual criterion specified by the experimenter. Since the task of adjustment is subjective, classification tasks would be suitable for the determination of those thresholds. In the literature, three kinds of classification are widely used: yes/no, two-alternative forced choice (2afc), and identification. All three call for the observer to classify stimuli (or their subjective responses). Applying these three kinds of methods to determine the thresholds, the identification one would give best results; see, e.g., Pelli & Farell, 1995. The drawback of the 2afc is that because the observer has only two alternatives and thus will be right half the time even if the luminance contrast is invisible, a relatively large number of trials is required to obtain a good threshold estimate. This method would give good results after the identification task. The observer in a yes/no experiment can't avoid introducing an internal subjective criterion in deciding whether each faint ambiguous percept deserves a 'yes' or a 'no'. This task would be used as a last resort.
<figure><img src="image_4.png" /><figcaption>Figure 2. Variation of the Bradley-Terry-Luce parameters ($\pi_j$)$_{j=1,2,...,q}$ with respect to its Michelson contrast levels M in the Red color range</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>Figure 3. Variation of the Bradley-Terry-Luce parameters ($\pi_j$)$_{j=1,2,...,q}$ with respect to its Michelson contrast levels M in the Green color range</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>Figure 4. Variation of the Bradley-Terry-Luce parameters ($\pi_j$)$_{j=1,2,...,q}$ with respect to its Michelson contrast levels M in the Yellow color range</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>Figure 5. Luminance contrasts ability of discrimination versus investigated colored textures patterns</figcaption></figure>

**5. Discussion**

The findings of this study clearly show that the ability of discrimination of the luminance contrast in the chromatic textures varies according to the color ranges and the textures patterns. The observers ability of discrimination is sharper in the green color range, followed by the red color range and finally the yellow color range, except in the type of Isotropic textures. In this type of texture, the visual perception is highest in the green color range, followed by yellow color range and finally the red one. Thus, it is sharper to detect luminance contrast in green color textures than the others ranges. This result supports the hypothesis that color information is processed independently from many luminance-based tasks (Livingston & Hubel, 1987).

Whatever the color range, the ability of discrimination is best in the type of Isotropic textures; followed successively by the Random-dots textures type and the Horizontal gratings type, and finally the Vertical gratings ones. It results that luminance
contrast detection in chromatic images does not depend only on the luminance and chromatic attributes. This result is in agreement with the ones obtained by (Gegenfurtner & Kiper, 1990) who studied *the contrast detection in luminance and chromatic noise*. They show that it is incorrect to assume that psychophysical performance is always determined by the activity of just a luminance and a chromatic mechanism. Therein, they demonstrated the existence of additional mechanisms that determine performance in detection of contrast. They concluded that the possible involvement of these mechanisms must be considered when one interprets the results of other visual tasks. Our data support the hypothesis that the importance of color to visual performance can be revealed if chromatic and luminance variations are both present in visual targets (Gur & Akri, 1992). Moreover, Gur & Akri have shown that when color and luminance contrasts are combined, human contrast sensitivity is enhanced at all spatial frequencies tested and that there is a substantial facilitatory interaction between the chromatic and achromatic systems.

A large-sample test based on the likelihood ratio shows that the sense of perception of the set of observers were sharp enough to detect the existing luminance contrast difference. Furthermore, one has been interested in testing whether $m_i = m_j$ across pairs of stimuli (see, e.g., Dzhafarov & Perry, 2014) via post-hoc tests. Those tests allowed us to find out which colored textures differ from the others. The data analysis gives rise to a sensory scale that shows a monotonic functional relationship with the Michelson contrast on which the ranking experiment is based: the shapes of the curves in Figs. 2-4 of the Bradley-Terry-Luce parameters ($\pi$) versus the Michelson contrast (M) of the colored textures indicate higher preference for lower contrast. We recall that the observers have ranked the stimuli according to their appreciation of the perceived visual contrast, from the lowest contrast to the highest contrast. They were asked to use integer from 1 to 20; rank 1 is assigned to the lowest visual contrast texture. Therefore, in the context of our experimentation, preference parameter $\pi_i$ is referred to luminance contrast sensitivity for stimuli $i$. Thus, higher preference means lower luminance contrast sensitivity. This is what the shapes of the curves in Figs. 2-4 render. In short, the Fechnerian scales on the physical dimension of the colored textures depend on the chromaticness of the colored phases of textures and the texture types. The psycho physical method of identification would be the best when determining the related thresholds.

During the ranking experiment, it appeared that each observer took an average of 90 minutes for a complete round, which is a bit long, even with the remarkable ability of chromatic adaptation that has the human visual system. Indeed, fatigue impairs observer's perception. The color space used for the physical description of the colored textures is the CIE 1931 standard observer using 2 degrees and the formula of contrast used is the Michelson one.

## Acknowledgements

The authors are grateful to the referees.

## References

Bonton, P., Fernandez-Maloigne, C., & Trémau, A. (2000). *Image numérique couleur (De l'acquisition au traitement)*. Dunod, Paris.

Chaudhuri, A. (2011). *Fundamentals of sensory perception*. Don mills, Ont: Oxford University Press.

Conover, W. J. (1999). *Practical Nonparametric Statistics* (3rd ed.). Wiley.

Critchlow, E. D., & Fligner, A. M. (1991). Paired comparison, triple comparison, and ranking experiments as generalized linear models, and their implementation on GLIM. *Psychometrika*, *56*, 517–533.

CIE (1926). *Commission internationale de l'Eclairage Proceedings, 1924*. Cambridge: Cambridge University Press.

Davidson, R. R. (1970). On extending the Bradley-Terry Model to accommodate ties in paired comparison experiments. *Journal of the American Statistical Association*, *65*, 317–328. https://doi.org/10.1080/01621459.1970.10481082

Dzhafarov, E. N., & Colonius, H. (1999). Fechnerian metrics in unidimensional and multidimensional stimulus spaces. *Psychonomic Bulletin Review*, *6*, 239-268. https://doi.org/10.3758/BF03212329

Dzhafarov, E. N., & Colonius, H. (2006). Reconstructing distances among objects from their dis-criminability. *Psychometrika*, *71*, 365-386. https://doi.org/10.1007/s11336-003-1126-9

Dzhafarov, E. N., & Perry, L. (2014). Perceptual matching and sorites: Experimental study of an ancient Greek paradox. *Attention, Perception, & Psychophysics*, *76*, 2441-2464. https://doi.org/10.3758/s13414-014-0711-6

Fairchild, M. D. (2013). *Colour appearance models*. Hoboken, Nj: Wiley.

Falmagne, J-C. (1985). *Elements of Psychophysical Theory*. Oxford University Press, Inc. https://doi.org/10.1002/9781118653128

Fligner, M. A., & Verducci, J. S. (1986). Distance Based Ranking Models. *Journal of the Royal Statistical Society*, *48*,
359–369. https://doi.org/10.1111/j.2517-6161.1986.tb01420.x

Fligner, M. A., & Verducci, J. S. (1988). Multistage Ranking Models. *Journal of the Royal Statistical Society, 83*, 892–901. https://doi.org/10.1080/01621459.1988.10478679

Gur, M., & Akri, V. (1992). Isoluminant stimuli may not expose the full contribution of color to visual functioning: Spatial contrast sensitivity measurements indicate interaction between color and luminance processing. *Vision Res.*, 32, 1253-1262. https://doi.org/10.1016/0042-6989(92)90220-D

Gegenfurtner, K. R., & Kiper, D. C. (1992). Contrast detection in luminance and chromatic noise. *J. Opt. Soc. Am. A*, *9*, 1880-1888. https://doi.org/10.1364/JOSAA.9.001880

Holm, S. (1979). A Simple Sequentially Rejective Multiple Test Procedure. *Scandinavian Journal of Statistics*, **6**, 65–70.

Hunt, R. W. G., & Pointer, M. R. (2011). *Measuring Colour* (4 ed.). Chichester, UK: John Wiley & Sons. Hunter, D. R., & Handcock, M. S. (2012). Inference in curved exponential family models for networks. *Journal of Computational and Graphical Statistics*, *15*, 565–583. https://doi.org/10.1002/9781119975595

Iversion, G., & Luce, R. D. (1998). *The Representational Measurement Approach to Psychophysical and Judgmental Problems, Measurement, Judgment, and decision Making*, by Academic Press. https://doi.org/10.1016/B978-012099975-0/50003-8

Livingtone, M., & Hubel, D. H. (1987). Psychophysical evidence for separate channels for the perception of form, color, movement and depth. *Neuroscience*, **7**, 3416-3468. https://doi.org/10.1523/JNEUROSCI.07-11-03416.1987

Marden, J. I. (1995). *Analysing and Modeling Rank Data*. London: Chapman and Hall.

Meng, X.-L., & Wong, W. H. (1996). Simulating ratios of normalizing constants via a simple identity: A theoretical exploration. *Statistica Sinica*, **6**, 831-860.

Palmer, S. E. (2002). *Vision Science: photons to phenomenology (third printing)*. The MIT Press. ISBN 0-262-16183-4.

Pelli, D. G., & Farell, B. (1995). *Psychophysical methods*. In Handbook of optics (2nd ed.), Vol. 1 ted. M. Bass, E.W. Van Stryland, D. R. Williams, & W. L. Wolfe, pp. 29.1-29.13. McGraw-Hili, New York.

Pohlert, T. (2019). *PMCMRplus: Calculate Pairwise Multiple Comparisons of Mean Rank Sums Extended*. R package version 1.4.2.

R Core Team. (2020). R: A Language and Environment for Statistical Computing, R Foundation for Statistical Computing, Vienna.

Robert, C. (1996). *Méthodes de Monte Carlo par chaînes de Markov*, Economica.

Reinhart, E., Khan, E. A., Akyez, A. O., & Johnson, G. (2008). *Imaging: Fundamentals and applications*. A K Peeters/CRC Press.

Sawadogo, A., Lafon, D., & Dossou-Gbété, S. (2014). Statistical analysis of rank data from a visual matching of colored textures. *Journal of Applied Statistics*, *41*(11), 2462-2482. https://doi.org/10.1080/02664763.2014.920775

Sawadogo, A., Dossou-Gbété, S., & Lafon, D. (2017). Ties in one block comparison experiments: a generalization of the Mallows-Bradley-Terry ranking model. *Journal of Applied Statistics*, *44*(14), 2621-2644. https://doi.org/10.1080/02664763.2016.1259400

Wandell, B. A. (1995). *Foundation of vision*. Sinauer Associates, Incorporated.

# Appendix

## 1. Proofs of Propositions

### 1.1 Proofs of Proposition 1

Let $r_i$, $i \in \{1, 2, \dots, n\}$ be a sample of rankings with ties independant and identically distributed as the Mallows-Bradley-Terry distribution of ranking with ties $P(r; \theta, \gamma)$. The log-likehood of the canonical parameter vector $(\theta, \gamma)$ is given
by

$$
\begin{align*}
\mathcal{L}_n(\boldsymbol{\theta}, \gamma) &= \sum_{i=1}^{n} \mathbb{P}(\boldsymbol{r}_i; \boldsymbol{\theta}, \gamma) \\
&= \sum_{j=1}^{q} \theta_j \sum_{i=1}^{n} \{q - \boldsymbol{r}_i^*(j)\} + \frac{1}{2}\gamma [-nq + \sum_{i=1}^{n} \{\sum_{k=1}^{r_i^*} \lambda_{\boldsymbol{r}_i}^2(k)\}] + n \log\{c(\boldsymbol{\theta}, \gamma)\},
\end{align*}
$$

$$
\text{with } \log\{c(\boldsymbol{\theta}, \gamma)\} = -\log\left(\sum_s \exp\left[\frac{1}{2}\left\{\sum_{k=1}^{s_+} \lambda_s^2(k) - q\right\}\gamma + \sum_{j=1}^q \left\{q - s^*(j)\right\}\theta_j\right]\right).
$$

1.2 Proofs of Proposition 3

Indeed, the model defined by Equation (3) can be rewritten

$$
\mathbb{P}(\boldsymbol{r}; \varphi(u)) = \exp\left[ \boldsymbol{z}(\boldsymbol{r}) \cdot \varphi(u) - \log\{c(\varphi(u))\} \right].
$$

We have,

$$
\frac{d}{du} \mathbb{P}(\boldsymbol{r}; \varphi(u)) = {}^{t}\nabla_{\varphi} \mathbb{P}(\boldsymbol{r}, \varphi(u)). \frac{d\varphi(u)}{du},
$$

and,

$$
{}^{t}\nabla_{\varphi} \mathbb{P}(\boldsymbol{r}; \varphi(u)) = \left[ \boldsymbol{z}(\boldsymbol{r}) - \frac{1}{c(\varphi(u))} {}^{t}\nabla_{\varphi} c(\varphi(u)) \right] \mathbb{P}(\boldsymbol{r}; \varphi(u)).
$$

Then,

$$
\frac{d}{du} \mathbb{P}(\boldsymbol{r}, \varphi(u)) = \left[ \boldsymbol{z}(\boldsymbol{r}) - \frac{1}{c(\varphi(u))} {}^{t}\nabla_{\varphi} c(\varphi(u)) \right] \frac{d\varphi(u)}{du} \cdot \mathbb{P}(\boldsymbol{r}; \varphi(u)).
$$

It results that, using Lemma 1,

$$
\sum_{\boldsymbol{r}} \frac{d}{du} \mathbb{P}(\boldsymbol{r}; \varphi(u)) = 0
$$

is equivalent to,

$$
\frac{1}{c(\varphi(u))} {}^{t}\nabla_{\varphi} c(\varphi(u)) \frac{d\varphi(u)}{du} = \mathbb{E}_{\varphi(u)} \left\{ {}^{t}\boldsymbol{z} \frac{d\varphi(u)}{du} \right\}.
$$

By integrating this last equality on the closed interval $[0, 1]$, one obtains

$$
\log \left\{ \frac{c(\varphi(1))}{c(\varphi(0))} \right\} = \int_{0}^{1} \mathbb{E}_{\varphi(u)} \left\{ {}^{t}\mathbf{Z} \frac{d\varphi(u)}{du} \right\} du. \quad (8)
$$

One has,

$$
\mathbb{E}_{\varphi(u)} \left\{ {}^{t}\mathbf{Z} \frac{d\varphi(u)}{du} \right\} = \mathbb{E} \left\{ {}^{t}\mathbf{Z} \left( \frac{d\varphi}{du} \circ U \right) \mid U = u \right\}
$$

thus,

$$
\begin{align*}
\int_0^1 \mathbb{E}_{\varphi(u)} \{ {}^t \mathbf{Z} \frac{d\varphi(u)}{du} \} du &= \mathbb{E} \left[ \mathbb{E} \left\{ {}^t \mathbf{Z} \left( \frac{d\varphi}{du} \circ U \right) \mid U = u \right\} \right] \\
&= \mathbb{E} \left\{ {}^t \mathbf{Z} \left( \frac{d\varphi}{du} \circ U \right) \right\}.
\end{align*}
$$

In fine

$$
\log \left\{ \frac{c(\varphi(1))}{c(\varphi(0))} \right\} = \mathbb{E} \left\{ {}^{t}\mathbf{Z} \left( \frac{d\varphi}{du} \circ U \right) \right\}.
$$

The mathematical expectation $\mathbb{E}$ is taken with respect to the joint distribution of the random vector $(U, \boldsymbol{r})$, where $U$ is distributed as the uniform distribution on the interval $[0, 1]$ and $\boldsymbol{r} \mid U = u$ distributed as $\mathbb{P}(\boldsymbol{r}; \varphi)$ of parameter $\varphi(U)$.
## 2. Estimates of the Bradley-Terry-Luce Parameters and Indexes of Discrimination for the Colored Textures Belonging to the Red Color Range and the Green One

<table><caption>Table 2. Estimates of the model parameter vector $(\pi, \nu)$ for each type of texture in the Red color range</caption><thead><tr><th colspan="2">RR</th><th colspan="2">RI</th><th colspan="2">RH</th><th colspan="2">RV</th></tr><tr><th>π̂</th><th>CI<sub>1-α</sub>(π)</th><th>π̂</th><th>CI<sub>1-α</sub>(π)</th><th>π̂</th><th>CI<sub>1-α</sub>(π)</th><th>π̂</th><th>CI<sub>1-α</sub>(π)</th></tr></thead><tbody><tr><td>0.0021</td><td>(0.0015,0.0056)</td><td>0.0194</td><td>(0.0174,0.0213)</td><td>0.0246</td><td>(0.0210,0.0281)</td><td>0.0023</td><td>(0.0031,0.0077)</td></tr><tr><td>0.0825</td><td>(0.0661,0.0989)</td><td>0.0715</td><td>(0.0643,0.0788)</td><td>0.0565</td><td>(0.0499,0.0631)</td><td>0.0799</td><td>(0.0693,0.0905)</td></tr><tr><td>0.0306</td><td>(0.0161,0.0451)</td><td>0.0412</td><td>(0.0380,0.0445)</td><td>0.0422</td><td>(0.0373,0.0470)</td><td>0.0429</td><td>(0.0338,0.0521)</td></tr><tr><td>0.1348</td><td>(0.0807,0.1889)</td><td>0.1015</td><td>(0.0817,0.1213)</td><td>0.0845</td><td>(0.0756,0.0935)</td><td>0.0828</td><td>(0.0697,0.0958)</td></tr><tr><td>0.0610</td><td>(0.0446,0.0775)</td><td>0.0574</td><td>(0.0555,0.0593)</td><td>0.0583</td><td>(0.0561,0.0604)</td><td>0.0665</td><td>(0.0618,0.0713)</td></tr><tr><td>0.1000</td><td>(0.0755,0.1245)</td><td>0.0754</td><td>(0.0709,0.0799)</td><td>0.0779</td><td>(0.0650,0.0908)</td><td>0.0869</td><td>(0.0751,0.0986)</td></tr><tr><td>0.0007</td><td>(0.0017,0.0030)</td><td>0.0153</td><td>(0.0135,0.0171)</td><td>0.0181</td><td>(0.0109,0.0254)</td><td>0.0015</td><td>(0.0014,0.0044)</td></tr><tr><td>0.0531</td><td>(0.0423,0.0639)</td><td>0.0534</td><td>(0.0464,0.0605)</td><td>0.0543</td><td>(0.0515,0.0571)</td><td>0.0597</td><td>(0.0542,0.0652)</td></tr><tr><td>0.0214</td><td>(0.0081,0.0346)</td><td>0.0333</td><td>(0.0295,0.0371)</td><td>0.0388</td><td>(0.0322,0.0455)</td><td>0.0387</td><td>(0.0314,0.0460)</td></tr><tr><td>0.0272</td><td>(0.0134,0.0410)</td><td>0.0364</td><td>(0.0343,0.0386)</td><td>0.0460</td><td>(0.043,0.049)</td><td>0.0422</td><td>(0.0333,0.0512)</td></tr><tr><td>0.0985</td><td>(0.0787,0.1184)</td><td>0.0799</td><td>(0.0752,0.0846)</td><td>0.0735</td><td>(0.0653,0.0817)</td><td>0.0806</td><td>(0.0658,0.0954)</td></tr><tr><td>0.0683</td><td>(0.0577,0.0789)</td><td>0.0578</td><td>(0.0520,0.0636)</td><td>0.0637</td><td>(0.0537,0.0737)</td><td>0.0741</td><td>(0.0645,0.0836)</td></tr><tr><td>0.0057</td><td>(0.0015,0.0128)</td><td>0.0220</td><td>(0.0191,0.0249)</td><td>0.0283</td><td>(0.0210,0.0357)</td><td>0.0044</td><td>(0.0033,0.0122)</td></tr><tr><td>0.0543</td><td>(0.0387,0.0700)</td><td>0.0540</td><td>(0.0470,0.0609)</td><td>0.0507</td><td>(0.0434,0.0581)</td><td>0.0689</td><td>(0.0573,0.0805)</td></tr><tr><td>0.0181</td><td>(0.0069,0.0294)</td><td>0.0309</td><td>(0.0264,0.0354)</td><td>0.0399</td><td>(0.0382,0.0417)</td><td>0.0361</td><td>(0.0239,0.0483)</td></tr><tr><td>0.0178</td><td>(0.0051,0.0305)</td><td>0.0325</td><td>(0.0278,0.0372)</td><td>0.0333</td><td>(0.0273,0.0392)</td><td>0.0314</td><td>(0.0193,0.0434)</td></tr><tr><td>0.0256</td><td>(0.0116,0.0395)</td><td>0.0352</td><td>(0.0318,0.0386)</td><td>0.0408</td><td>(0.0377,0.0439)</td><td>0.0284</td><td>(0.0179,0.0389)</td></tr><tr><td>0.0471</td><td>(0.0368,0.0573)</td><td>0.0486</td><td>(0.0473,0.0499)</td><td>0.0536</td><td>(0.0481,0.0591)</td><td>0.0579</td><td>(0.0508,0.0649)</td></tr><tr><td>0.1286</td><td>(0.0615,0.1957)</td><td>0.1003</td><td>(0.0869,0.1137)</td><td>0.0773</td><td>(0.0639,0.0907)</td><td>0.0829</td><td>(0.0747,0.0911)</td></tr><tr><td>0.0227</td><td>(0.0152,0.0302)</td><td>0.0339</td><td>(0.0280,0.0398)</td><td>0.0376</td><td>(0.0364,0.0388)</td><td>0.0319</td><td>(0.0243,0.0395)</td></tr><tr><th>ν̂</th><th>CI<sub>1-α</sub>(ν)</th><th>ν̂</th><th>CI<sub>1-α</sub>(ν)</th><th>ν̂</th><th>CI<sub>1-α</sub>(ν)</th><th>ν̂</th><th>CI<sub>1-α</sub>(ν)</th></tr><tr><td>1.2364</td><td>(0.9177,1.5550)</td><td>1.0242</td><td>(0.9862,1.0621)</td><td>1.3215</td><td>(1.210,1.433)</td><td>1.5291</td><td>(1.3469,1.7113)</td></tr></tbody></table>

<table><caption>Table 3. Estimates of the model parameter vector $(\pi, \nu)$ for each type of texture in the Green color range</caption><thead><tr><th colspan="2">GR</th><th colspan="2">GI</th><th colspan="2">GH</th><th colspan="2">GV</th></tr><tr><th>π̂</th><th>CI<sub>1-α</sub>(π)</th><th>π̂</th><th>CI<sub>1-α</sub>(π)</th><th>π̂</th><th>CI<sub>1-α</sub>(π)</th><th>π̂</th><th>CI<sub>1-α</sub>(π)</th></tr></thead><tbody><tr><td>0.0006</td><td>(0.0005,0.0018)</td><td>0.0010</td><td>(0.000,0.002)</td><td>0.0120</td><td>(0.0040,0.0201)</td><td>0.0090</td><td>(0.0005,0.0175)</td></tr><tr><td>0.0446</td><td>(0.0288,0.0603)</td><td>0.0282</td><td>(0.0211,0.0353)</td><td>0.0607</td><td>(0.0530,0.0684)</td><td>0.0676</td><td>(0.0573,0.0780)</td></tr><tr><td>0.0114</td><td>(0.0080,0.0149)</td><td>0.0093</td><td>(0.0062,0.0123)</td><td>0.0339</td><td>(0.0279,0.0400)</td><td>0.0303</td><td>(0.0203,0.0404)</td></tr><tr><td>0.3574</td><td>(0.2072,0.5075)</td><td>0.6137</td><td>(0.5606,0.6667)</td><td>0.1245</td><td>(0.0931,0.1559)</td><td>0.1022</td><td>(0.0742,0.1303)</td></tr><tr><td>0.0318</td><td>(0.0211,0.0425)</td><td>0.0191</td><td>(0.0139,0.0242)</td><td>0.0519</td><td>(0.0460,0.0577)</td><td>0.0538</td><td>(0.0495,0.0581)</td></tr><tr><td>0.0954</td><td>(0.0441,0.1466)</td><td>0.0466</td><td>(0.0328,0.0604)</td><td>0.0812</td><td>(0.0761,0.0863)</td><td>0.0963</td><td>(0.0801,0.1124)</td></tr><tr><td>0.0002</td><td>(1e-04,4e-04)</td><td>0.0003</td><td>(1e-04,8e-04)</td><td>0.0063</td><td>(0.0009,0.0117)</td><td>0.0082</td><td>(0.0015,0.0150)</td></tr><tr><td>0.0213</td><td>(0.0153,0.0273)</td><td>0.0141</td><td>(0.0090,0.0193)</td><td>0.0454</td><td>(0.0424,0.0485)</td><td>0.0456</td><td>(0.0382,0.0530)</td></tr><tr><td>0.0124</td><td>(0.0077,0.0171)</td><td>0.0095</td><td>(0.0062,0.0127)</td><td>0.0433</td><td>(0.0368,0.0498)</td><td>0.0370</td><td>(0.0302,0.0437)</td></tr><tr><td>0.0141</td><td>(0.0103,0.0178)</td><td>0.0109</td><td>(0.0074,0.0145)</td><td>0.0445</td><td>(0.0412,0.0479)</td><td>0.0412</td><td>(0.0322,0.0503)</td></tr><tr><td>0.0618</td><td>(0.0399,0.0836)</td><td>0.0348</td><td>(0.0253,0.0443)</td><td>0.0710</td><td>(0.0642,0.0778)</td><td>0.0822</td><td>(0.0672,0.0973)</td></tr><tr><td>0.0661</td><td>(0.0384,0.0937)</td><td>0.0378</td><td>(0.0235,0.0520)</td><td>0.0703</td><td>(0.0610,0.0796)</td><td>0.0782</td><td>(0.0652,0.0911)</td></tr><tr><td>0.0013</td><td>(0.0006,0.0032)</td><td>0.0021</td><td>(0.0005,0.0036)</td><td>0.0155</td><td>(0.0073,0.0236)</td><td>0.0178</td><td>(0.0072,0.0285)</td></tr><tr><td>0.0314</td><td>(0.0199,0.0429)</td><td>0.0193</td><td>(0.0134,0.0251)</td><td>0.0559</td><td>(0.0480,0.0639)</td><td>0.0549</td><td>(0.0464,0.0633)</td></tr><tr><td>0.0083</td><td>(0.0048,0.0118)</td><td>0.0067</td><td>(0.0031,0.0103)</td><td>0.0391</td><td>(0.0327,0.0454)</td><td>0.0314</td><td>(0.0252,0.0377)</td></tr><tr><td>0.0046</td><td>(0.0013,0.0079)</td><td>0.0044</td><td>(0.0016,0.0072)</td><td>0.0195</td><td>(0.0106,0.0284)</td><td>0.0230</td><td>(0.0141,0.0320)</td></tr><tr><td>0.0089</td><td>(0.0049,0.0128)</td><td>0.0060</td><td>(0.0031,0.0089)</td><td>0.0273</td><td>(0.0195,0.0350)</td><td>0.0282</td><td>(0.0176,0.0387)</td></tr><tr><td>0.0250</td><td>(0.0174,0.0326)</td><td>0.0167</td><td>(0.0097,0.0237)</td><td>0.0550</td><td>(0.0493,0.0607)</td><td>0.0522</td><td>(0.0447,0.0598)</td></tr><tr><td>0.1863</td><td>(0.0932,0.2795)</td><td>0.1040</td><td>(0.0274,0.1805)</td><td>0.1070</td><td>(0.0917,0.1222)</td><td>0.1082</td><td>(0.0630,0.1533)</td></tr><tr><td>0.0174</td><td>(0.0125,0.0222)</td><td>0.0157</td><td>(0.0121,0.0193)</td><td>0.0357</td><td>(0.0319,0.0394)</td><td>0.0325</td><td>(0.0261,0.0390)</td></tr><tr><th>ν̂</th><th>CI<sub>1-α</sub>(ν)</th><th>ν̂</th><th>CI<sub>1-α</sub>(ν)</th><th>ν̂</th><th>CI<sub>1-α</sub>(ν)</th><th>ν̂</th><th>CI<sub>1-α</sub>(ν)</th></tr><tr><td>0.7273</td><td>(0.5892,0.8653)</td><td>0.6027</td><td>(0.4043,0.8011)</td><td>1.2114</td><td>(1.0712,1.3516)</td><td>1.3243</td><td>(1.1424,1.5063)</td></tr></tbody></table>
<table><caption>Table 4. Estimates of the model parameter vector ($\pi$, $\nu$) for each type of texture in the Yellow color range</caption><thead><tr><th colspan="2">YR</th><th colspan="2">YI</th><th colspan="2">YH</th><th colspan="2">YV</th></tr><tr><th>$\hat{\pi}$</th><th>$CI_{1-\alpha}(\pi)$</th><th>$\hat{\pi}$</th><th>$CI_{1-\alpha}(\pi)$</th><th>$\hat{\pi}$</th><th>$CI_{1-\alpha}(\pi)$</th><th>$\hat{\pi}$</th><th>$CI_{1-\alpha}(\pi)$</th></tr></thead><tbody><tr><td>0.0069</td><td>(0.0060, 0.0198)</td><td>0.0006</td><td>(0.0001, 0.0012)</td><td>0.0222</td><td>(0.0187, 0.0257)</td><td>0.0242</td><td>(0.0195, 0.0290)</td></tr><tr><td>0.1001</td><td>(0.0650, 0.1352)</td><td>0.0574</td><td>(0.0373, 0.0774)</td><td>0.0697</td><td>(0.0648, 0.0746)</td><td>0.0634</td><td>(0.0593, 0.0676)</td></tr><tr><td>0.0373</td><td>(0.0177, 0.0570)</td><td>0.0174</td><td>(0.0112, 0.0236)</td><td>0.0516</td><td>(0.0465, 0.0568)</td><td>0.0592</td><td>(0.0555, 0.0629)</td></tr><tr><td>0.1435</td><td>(0.063, 0.224)</td><td>0.3663</td><td>(0.2572, 0.4754)</td><td>0.0852</td><td>(0.0758, 0.0946)</td><td>0.0632</td><td>(0.0605, 0.0659)</td></tr><tr><td>0.0681</td><td>(0.0565, 0.0797)</td><td>0.0385</td><td>(0.0281, 0.0489)</td><td>0.0663</td><td>(0.0593, 0.0732)</td><td>0.0619</td><td>(0.0581, 0.0658)</td></tr><tr><td>0.0864</td><td>(0.0754, 0.0974)</td><td>0.0468</td><td>(0.0328, 0.0607)</td><td>0.0710</td><td>(0.064, 0.078)</td><td>0.0752</td><td>(0.0696, 0.0808)</td></tr><tr><td>0.0043</td><td>(0.0050, 0.0136)</td><td>0.0004</td><td>(2e-04, 1e-03)</td><td>0.0191</td><td>(0.0144, 0.0237)</td><td>0.0145</td><td>(0.0090, 0.0199)</td></tr><tr><td>0.0572</td><td>(0.0436, 0.0707)</td><td>0.0293</td><td>(0.0208, 0.0377)</td><td>0.0606</td><td>(0.0542, 0.0671)</td><td>0.0605</td><td>(0.0532, 0.0678)</td></tr><tr><td>0.0108</td><td>(0.0015, 0.0231)</td><td>0.0027</td><td>(0.0016, 0.0037)</td><td>0.0328</td><td>(0.0264, 0.0392)</td><td>0.0341</td><td>(0.0315, 0.0366)</td></tr><tr><td>0.0179</td><td>(0.0035, 0.0324)</td><td>0.0049</td><td>(0.0027, 0.0071)</td><td>0.0367</td><td>(0.0329, 0.0405)</td><td>0.0427</td><td>(0.0400, 0.0455)</td></tr><tr><td>0.1277</td><td>(0.0777, 0.1776)</td><td>0.2119</td><td>(0.1642, 0.2595)</td><td>0.0803</td><td>(0.0730, 0.0877)</td><td>0.0625</td><td>(0.0593, 0.0656)</td></tr><tr><td>0.0537</td><td>(0.0361, 0.0712)</td><td>0.0298</td><td>(0.0230, 0.0365)</td><td>0.0601</td><td>(0.0559, 0.0643)</td><td>0.0662</td><td>(0.0625, 0.0699)</td></tr><tr><td>0.0083</td><td>(0.0037, 0.0204)</td><td>0.0021</td><td>(0.0008, 0.0034)</td><td>0.0271</td><td>(0.0222, 0.0320)</td><td>0.0292</td><td>(0.0253, 0.0331)</td></tr><tr><td>0.0477</td><td>(0.0302, 0.0652)</td><td>0.0202</td><td>(0.0150, 0.0253)</td><td>0.0484</td><td>(0.0427, 0.0541)</td><td>0.0660</td><td>(0.0609, 0.0712)</td></tr><tr><td>0.0090</td><td>(0.0046, 0.0225)</td><td>0.0021</td><td>(0.0011, 0.0032)</td><td>0.0303</td><td>(0.0278, 0.0328)</td><td>0.0310</td><td>(0.0248, 0.0373)</td></tr><tr><td>0.0193</td><td>(0.0040, 0.0346)</td><td>0.0072</td><td>(0.0053, 0.0092)</td><td>0.0365</td><td>(0.0319, 0.0410)</td><td>0.0457</td><td>(0.0418, 0.0497)</td></tr><tr><td>0.0260</td><td>(0.0076, 0.0445)</td><td>0.0105</td><td>(0.0084, 0.0126)</td><td>0.0401</td><td>(0.0345, 0.0457)</td><td>0.0477</td><td>(0.0447, 0.0507)</td></tr><tr><td>0.0356</td><td>(0.0204, 0.0509)</td><td>0.0146</td><td>(0.0103, 0.0189)</td><td>0.0447</td><td>(0.0426, 0.0468)</td><td>0.0547</td><td>(0.0508, 0.0586)</td></tr><tr><td>0.1200</td><td>(0.0817, 0.1584)</td><td>0.1220</td><td>(0.0169, 0.2270)</td><td>0.0857</td><td>(0.0679, 0.1035)</td><td>0.0643</td><td>(0.0560, 0.0727)</td></tr><tr><td>0.0202</td><td>(0.0149, 0.0255)</td><td>0.0154</td><td>(0.0117, 0.0191)</td><td>0.0316</td><td>(0.0290, 0.0343)</td><td>0.0336</td><td>(0.0292, 0.0379)</td></tr><tr><td>$\hat{\bar{v}}$</td><td>$CI_{1-\alpha}(\bar{v})$</td><td>$\hat{\bar{v}}$</td><td>$CI_{1-\alpha}(\bar{v})$</td><td>$\hat{\bar{v}}$</td><td>$CI_{1-\alpha}(\bar{v})$</td><td>$\hat{\bar{v}}$</td><td>$CI_{1-\alpha}(\bar{v})$</td></tr><tr><td>1.4452</td><td>(1.1983, 1.6922)</td><td>0.8474</td><td>(0.6464, 1.0483)</td><td>1.5043</td><td>(1.4669, 1.5416)</td><td>1.5678</td><td>(1.4259, 1.7097)</td></tr></tbody></table>

### 3. Percentage of Nonuniformity in the Data Explained by the Model

We define,

$$R^2 = 1 - \frac{GOF}{Denom}, \quad GOF = 2 \sum_r N_r \log\left(\frac{N_r}{n\mathbb{P}(r; \hat{\theta})}\right) \quad \text{and} \quad Denom = 2 \sum_r N_r \log\left(\frac{N_r}{n/K}\right)$$

with, $n$ the number of observers, $N_r$ the observed frequency of the ranking with ties $r$ and $K = \text{card}(\mathcal{A}(q))$. If the model exactly fits the data, $R^2 = 1$, and if it performs no better than the uniform, $R^2 = 0$.

<table><caption>Table 5. Percentage of nonuniformity in the data</caption><thead><tr><th>Data</th><th>RR</th><th>RI</th><th>RH</th><th>RV</th></tr></thead><tbody><tr><td>$R^2(\%)$</td><td>83.23</td><td>83.76</td><td>82.20</td><td>80.73</td></tr></tbody></table>

<table><caption>Table 6. Percentage of nonuniformity in the data.</caption><thead><tr><th>Data</th><th>GR</th><th>GI</th><th>GH</th><th>GV</th></tr></thead><tbody><tr><td>$R^2(\%)$</td><td>87.81</td><td>89.16</td><td>84.38</td><td>81.86</td></tr></tbody></table>

<table><caption>Table 7. Percentage of nonuniformity in the data.</caption><thead><tr><th>Data</th><th>YR</th><th>YI</th><th>YH</th><th>YV</th></tr></thead><tbody><tr><td>$R^2(\%)$</td><td>77.11</td><td>77.58</td><td>76.95</td><td>76.48</td></tr></tbody></table>
## 4. Estimates of the Parameter $\gamma_0$, Statistic Values and p-values, Associated With the Observed Data for the Colored Textures

<table><caption>Table 8. Estimates of $\gamma_0$ involved in the large-sample test based on the likelihood ratio for the colored textures belonging to the red color ranges</caption><thead><tr><th>Textures</th><th>RR</th><th>RI</th><th>RH</th><th>RV</th></tr></thead><tbody><tr><td>$\hat{\gamma}_0$</td><td>0.1380</td><td>-0.0212</td><td>0.2461</td><td>0.2601</td></tr><tr><td>$CI_{1-\alpha}(\gamma)$</td><td>(0.1020,0.1523)</td><td>(-0.0360,0.0231)</td><td>(0.2160,0.2510)</td><td>(0.1600,0.3500)</td></tr></tbody></table>

<table><caption>Table 9. Estimates of $\gamma_0$ involved in the large-sample test based on the likelihood ratio for the colored textures belonging to the Green color ranges</caption><thead><tr><th>Textures</th><th>GR</th><th>GI</th><th>GH</th><th>GV</th></tr></thead><tbody><tr><td>$\hat{\gamma}_0$</td><td>-0.1040</td><td>-0.2001</td><td>0.2366</td><td>0.2324</td></tr><tr><td>$CI_{1-\alpha}(\gamma)$</td><td>(-1.0023,0.2000)</td><td>(-0.2812,0.0575)</td><td>(0.2022,0.2552)</td><td>(0.2012,0.2500)</td></tr></tbody></table>

<table><caption>Table 10. Estimates of $\gamma_0$ involved in the large-sample test based on the likelihood ratio for the colored textures belonging to the Yellow color ranges</caption><thead><tr><th>Textures</th><th>YR</th><th>YI</th><th>YH</th><th>YV</th></tr></thead><tbody><tr><td>$\hat{\gamma}_0$</td><td>0.2573</td><td>0.0784</td><td>0.2512</td><td>0.2490</td></tr><tr><td>$CI_{1-\alpha}(\gamma)$</td><td>(0.1600,0.3100)</td><td>(0.0000,0.0910)</td><td>(0.2400,0.2723)</td><td>(0.2000,0.2632)</td></tr></tbody></table>

<table><caption>Table 11. Statistic values and p-values associated with the observed data for the colored textures belonging to Red color range</caption><thead><tr><th>Textures</th><th>-2 log λ<sub>n</sub></th><th>p-value</th></tr></thead><tbody><tr><td>RR</td><td>1105.1861</td><td>0.0000</td></tr><tr><td>RI</td><td>936.9275</td><td>0.0000</td></tr><tr><td>RH</td><td>738.4551</td><td>0.0000</td></tr><tr><td>RV</td><td>1004.1607</td><td>0.0000</td></tr></tbody></table>

<table><caption>Table 12. Statistic values and p-values associated with the observed data for the colored textures belonging to Green color range</caption><thead><tr><th>Textures</th><th>-2 log λ<sub>n</sub></th><th>p-value</th></tr></thead><tbody><tr><td>GR</td><td>1555.7132</td><td>0.0000</td></tr><tr><td>GI</td><td>1723.2944</td><td>0.0000</td></tr><tr><td>GH</td><td>987.3761</td><td>0.0000</td></tr><tr><td>GV</td><td>1012.5445</td><td>0.0000</td></tr></tbody></table>

<table><caption>Table 13. Statistic values and p-values associated with the observed data for the colored textures belonging to Yellow color range</caption><thead><tr><th>Textures</th><th>-2 log λ<sub>n</sub></th><th>p-value</th></tr></thead><tbody><tr><td>YR</td><td>1108.5782</td><td>0.0000</td></tr><tr><td>YI</td><td>2001.4651</td><td>0.0000</td></tr><tr><td>YH</td><td>762.6826</td><td>0.0000</td></tr><tr><td>YV</td><td>640.1447</td><td>0.0000</td></tr></tbody></table>
# 5. Results of the Conover Post-Hoc Tests

## 5.1 Results of the Conover Post-Hoc Tests for the Investigated Textures Type in Red Color Range

<table><caption>Table 14. Summary statistics for each colored texture and a letter symbol: different letters indicate significant differences between colored textures in the set RR based on the selected level of α = 0.05</caption><thead><tr><th>RR</th><th>Median</th><th>1st.Quartile</th><th>3rd.Quartile</th><th>Sig. group</th></tr></thead><tbody><tr><td>RRLCA</td><td>19.00</td><td>18.375</td><td>19.000</td><td>ab</td></tr><tr><td>RRLCO</td><td>5.00</td><td>4.000</td><td>5.625</td><td>cdef</td></tr><tr><td>RRLCR</td><td>11.50</td><td>11.000</td><td>14.000</td><td>ghi</td></tr><tr><td>RRLCU</td><td>2.00</td><td>1.375</td><td>3.000</td><td>c</td></tr><tr><td>RRLFE</td><td>6.75</td><td>6.000</td><td>7.625</td><td>cdefj</td></tr><tr><td>RRLGA</td><td>3.75</td><td>3.000</td><td>4.000</td><td>cd</td></tr><tr><td>RRLK</td><td>19.50</td><td>19.000</td><td>20.000</td><td>a</td></tr><tr><td>RRLMN</td><td>8.00</td><td>7.500</td><td>9.000</td><td>efgj</td></tr><tr><td>RRLMO</td><td>16.00</td><td>14.000</td><td>16.500</td><td>bhk</td></tr><tr><td>RRLNB</td><td>13.75</td><td>12.000</td><td>14.500</td><td>ghi '</td></tr><tr><td>RRLNV</td><td>3.00</td><td>3.000</td><td>4.250</td><td>cde</td></tr><tr><td>RRLRB</td><td>6.50</td><td>6.000</td><td>7.000</td><td>cdef</td></tr><tr><td>RRLSC</td><td>18.00</td><td>17.000</td><td>19.000</td><td>abk</td></tr><tr><td>RRLSR</td><td>9.00</td><td>8.000</td><td>9.000</td><td>defgj</td></tr><tr><td>RRLTC</td><td>16.25</td><td>14.000</td><td>17.000</td><td>abhk</td></tr><tr><td>RRLTR</td><td>15.00</td><td>14.500</td><td>17.000</td><td>abhk</td></tr><tr><td>RRLV</td><td>14.00</td><td>12.875</td><td>15.000</td><td>hik</td></tr><tr><td>RRLY</td><td>10.00</td><td>9.000</td><td>10.125</td><td>fgij</td></tr><tr><td>RRLZN</td><td>2.75</td><td>1.500</td><td>3.625</td><td>c</td></tr><tr><td>RRLZR</td><td>12.00</td><td>10.500</td><td>12.000</td><td>ghij</td></tr></tbody></table>

<figure><img src="image_8.png" /><figcaption>Figure 6. Box-whisker plot for each of the 20 random-dots textures in the Red color range</figcaption></figure>
<table><caption>Table 15. Summary statistics for each colored texture and a letter symbol: different letters indicate significant differences between colored textures in the set RI based on the selected level of α = 0.05</caption><thead><tr><th>RI</th><th>Median</th><th>1st.Quartile</th><th>3rd.Quartile</th><th>Sig. group</th></tr></thead><tbody><tr><td>RILCA</td><td>18.50</td><td>18.000</td><td>19.000</td><td>ab</td></tr><tr><td>RILCO</td><td>5.00</td><td>5.000</td><td>5.000</td><td>cdef</td></tr><tr><td>RILCR</td><td>12.00</td><td>11.000</td><td>13.000</td><td>ghij</td></tr><tr><td>RILCU</td><td>1.50</td><td>1.000</td><td>2.000</td><td>c</td></tr><tr><td>RILFE</td><td>6.75</td><td>6.000</td><td>7.500</td><td>cdefg</td></tr><tr><td>RILGA</td><td>4.00</td><td>3.500</td><td>4.000</td><td>cde</td></tr><tr><td>RILK</td><td>20.00</td><td>19.000</td><td>20.000</td><td>a</td></tr><tr><td>RILMN</td><td>8.00</td><td>7.875</td><td>9.000</td><td>dfgh</td></tr><tr><td>RILMO</td><td>16.00</td><td>14.000</td><td>17.000</td><td>abijk</td></tr><tr><td>RILNB</td><td>13.25</td><td>12.750</td><td>14.500</td><td>hijk</td></tr><tr><td>RILNV</td><td>3.00</td><td>3.000</td><td>3.500</td><td>cde</td></tr><tr><td>RILRB</td><td>7.00</td><td>6.000</td><td>8.000</td><td>defg</td></tr><tr><td>RILSC</td><td>18.00</td><td>16.875</td><td>18.625</td><td>abk</td></tr><tr><td>RILSR</td><td>9.00</td><td>7.500</td><td>9.000</td><td>dfgh</td></tr><tr><td>RILTC</td><td>16.50</td><td>15.000</td><td>17.250</td><td>abik</td></tr><tr><td>RILTR</td><td>15.00</td><td>13.750</td><td>16.250</td><td>abik</td></tr><tr><td>RILV</td><td>14.00</td><td>12.500</td><td>15.125</td><td>bijk</td></tr><tr><td>RILY</td><td>10.00</td><td>9.000</td><td>11.000</td><td>fghj</td></tr><tr><td>RILZN</td><td>2.00</td><td>1.500</td><td>2.000</td><td>ce</td></tr><tr><td>RILZR</td><td>12.00</td><td>11.000</td><td>13.000</td><td>ghij</td></tr></tbody></table>

<figure><img src="image_9.png" /><figcaption>Figure 7. Box-whisker plot for each of the 20 isotropic textures in the Red color range</figcaption></figure>
<table><caption>Table 16. Summary statistics for each colored texture and a letter symbol: different letters indicate significant differences between colored textures belonging to the set RH based on the selected level of α = 0.05</caption><thead><tr><th>RH</th><th>Median</th><th>1st.Quartile</th><th>3rd.Quartile</th><th>Sig. group</th></tr></thead><tbody><tr><td>RHLCA</td><td>18.75</td><td>17.000</td><td>19.000</td><td>ab</td></tr><tr><td>RHLCO</td><td>7.00</td><td>5.000</td><td>9.000</td><td>cde</td></tr><tr><td>RHLCR</td><td>12.25</td><td>10.000</td><td>15.625</td><td>fghi</td></tr><tr><td>RHLCU</td><td>2.00</td><td>1.375</td><td>3.625</td><td>c</td></tr><tr><td>RHLFE</td><td>6.50</td><td>5.000</td><td>10.000</td><td>cde</td></tr><tr><td>RHLGA</td><td>3.75</td><td>3.000</td><td>5.000</td><td>c</td></tr><tr><td>RHLK</td><td>19.00</td><td>18.000</td><td>20.000</td><td>a</td></tr><tr><td>RHLMN</td><td>9.00</td><td>7.000</td><td>11.000</td><td>def</td></tr><tr><td>RHLMO</td><td>15.50</td><td>11.500</td><td>18.000</td><td>abghi</td></tr><tr><td>RHLNB</td><td>12.75</td><td>9.875</td><td>14.000</td><td>dfgh</td></tr><tr><td>RHLNV</td><td>3.50</td><td>2.500</td><td>5.375</td><td>ce</td></tr><tr><td>RHLRB</td><td>6.25</td><td>4.500</td><td>7.250</td><td>ce</td></tr><tr><td>RHLSC</td><td>17.00</td><td>15.875</td><td>18.000</td><td>abi</td></tr><tr><td>RHLSR</td><td>9.75</td><td>7.250</td><td>11.625</td><td>defg</td></tr><tr><td>RHLTC</td><td>15.00</td><td>11.375</td><td>18.000</td><td>abghi</td></tr><tr><td>RHLTR</td><td>16.00</td><td>14.500</td><td>17.000</td><td>abhi</td></tr><tr><td>RHLV</td><td>14.00</td><td>11.500</td><td>16.000</td><td>bfghi</td></tr><tr><td>RHLY</td><td>9.00</td><td>7.000</td><td>11.625</td><td>defg</td></tr><tr><td>RHLZN</td><td>2.75</td><td>2.000</td><td>4.125</td><td>c</td></tr><tr><td>RHLZR</td><td>11.25</td><td>9.500</td><td>14.000</td><td>dfgh</td></tr></tbody></table>

<figure><img src="image_10.png" /><figcaption>Figure 8. Box-whisker plot for each of the 20 horizontal gratings in the Red color range</figcaption></figure>
<table><caption>Table 17. Summary statistics for each colored texture and a letter symbol: different letters indicate significant differences between colored textures belonging to the set RV based on the selected level of α = 0.05</caption><thead><tr><th>RV</th><th>Median</th><th>1st.Quartile</th><th>3rd.Quartile</th><th>Sig. group</th></tr></thead><tbody><tr><td>RVLCA</td><td>19.00</td><td>18.000</td><td>19.000</td><td>a</td></tr><tr><td>RVLCO</td><td>4.00</td><td>3.000</td><td>7.625</td><td>bc</td></tr><tr><td>RVLCR</td><td>11.75</td><td>11.000</td><td>14.250</td><td>def</td></tr><tr><td>RVLCU</td><td>4.00</td><td>3.000</td><td>7.125</td><td>bc</td></tr><tr><td>RVLFE</td><td>6.25</td><td>4.750</td><td>8.000</td><td>bcg</td></tr><tr><td>RVLGA</td><td>4.00</td><td>2.750</td><td>6.000</td><td>b</td></tr><tr><td>RVLK</td><td>19.00</td><td>18.000</td><td>20.000</td><td>a</td></tr><tr><td>RVLMN</td><td>8.00</td><td>6.875</td><td>9.000</td><td>bcdg</td></tr><tr><td>RVLMO</td><td>15.50</td><td>12.750</td><td>17.000</td><td>ef</td></tr><tr><td>RVLNB</td><td>13.50</td><td>12.000</td><td>15.000</td><td>def</td></tr><tr><td>RVLNV</td><td>4.00</td><td>3.000</td><td>6.250</td><td>bc</td></tr><tr><td>RVLRB</td><td>6.00</td><td>4.000</td><td>7.125</td><td>bc</td></tr><tr><td>RVLSC</td><td>18.00</td><td>17.375</td><td>19.000</td><td>ah</td></tr><tr><td>RVLSR</td><td>8.00</td><td>6.000</td><td>9.000</td><td>bcg</td></tr><tr><td>RVLTC</td><td>16.00</td><td>12.750</td><td>18.000</td><td>efh</td></tr><tr><td>RVLTR</td><td>15.00</td><td>13.000</td><td>16.000</td><td>aeh</td></tr><tr><td>RVLV</td><td>14.50</td><td>12.000</td><td>16.625</td><td>aeh</td></tr><tr><td>RVLY</td><td>10.00</td><td>9.375</td><td>10.250</td><td>cdfg</td></tr><tr><td>RVLZN</td><td>4.00</td><td>3.000</td><td>6.500</td><td>bc</td></tr><tr><td>RVLZR</td><td>12.00</td><td>10.000</td><td>13.000</td><td>defg</td></tr></tbody></table>

<figure><img src="image_11.png" /><figcaption>Figure 9. Box-whisker plot for each of the 20 vertical gratings in the Red color range</figcaption></figure>
## 5.2 Results of the Conover Post-Hoc Tests for the Investigated Textures Type in the Green Color Range

<table><caption>Table 18. Summary statistics for each colored texture and a letter symbol: different letters indicate significant differences between colored textures belonging to the set GR based on the selected level of α = 0.05</caption><thead><tr><th>GR</th><th>Median</th><th>1st.Quartile</th><th>3rd.Quartile</th><th>Sig. group</th></tr></thead><tbody><tr><td>GRLCA</td><td>19.00</td><td>18.375</td><td>19.000</td><td>ab</td></tr><tr><td>GRLCO</td><td>6.00</td><td>5.500</td><td>6.000</td><td>cde</td></tr><tr><td>GRLCR</td><td>13.00</td><td>12.000</td><td>15.000</td><td>fgh</td></tr><tr><td>GRLCU</td><td>1.50</td><td>1.000</td><td>2.000</td><td>c</td></tr><tr><td>GRLFE</td><td>7.00</td><td>7.000</td><td>7.625</td><td>de</td></tr><tr><td>GRLGA</td><td>3.00</td><td>3.000</td><td>4.000</td><td>cd</td></tr><tr><td>GRLK</td><td>20.00</td><td>19.000</td><td>20.000</td><td>a</td></tr><tr><td>GRLMN</td><td>9.25</td><td>9.000</td><td>11.000</td><td>efgi</td></tr><tr><td>GRLMO</td><td>13.75</td><td>12.000</td><td>15.000</td><td>fgh</td></tr><tr><td>GRLNB</td><td>13.00</td><td>11.875</td><td>14.000</td><td>fghi</td></tr><tr><td>GRLNV</td><td>4.00</td><td>3.875</td><td>5.000</td><td>cd</td></tr><tr><td>GRLRB</td><td>5.00</td><td>4.000</td><td>5.000</td><td>cd</td></tr><tr><td>GRLSC</td><td>18.00</td><td>17.500</td><td>18.000</td><td>ab</td></tr><tr><td>GRLSR</td><td>8.00</td><td>8.000</td><td>8.125</td><td>dei</td></tr><tr><td>GRLTC</td><td>15.00</td><td>13.750</td><td>17.000</td><td>bfh</td></tr><tr><td>GRLTR</td><td>17.00</td><td>16.000</td><td>17.000</td><td>abh</td></tr><tr><td>GRLV</td><td>15.00</td><td>14.000</td><td>16.000</td><td>bfh</td></tr><tr><td>GRLY</td><td>9.25</td><td>8.875</td><td>10.000</td><td>egi</td></tr><tr><td>GRLZN</td><td>2.00</td><td>1.500</td><td>2.000</td><td>c</td></tr><tr><td>GRLZR</td><td>11.00</td><td>10.000</td><td>11.625</td><td>efgi</td></tr></tbody></table>

<figure><img src="image_12.png" /><figcaption>Figure 10. Box-whisker plot for each of the 20 random-dots textures in the Red color range</figcaption></figure>
<table><caption>Table 19. Summary statistics for each colored texture and a letter symbol: different letters indicate significant differences between colored textures belonging to the set GI based on the selected level of α = 0.05</caption><thead><tr><th>GI</th><th>Median</th><th>1st.Quartile</th><th>3rd.Quartile</th><th>Sig. group</th></tr></thead><tbody><tr><td>GILCA</td><td>19.00</td><td>18.000</td><td>19.000</td><td>ab</td></tr><tr><td>GILCO</td><td>6.00</td><td>5.875</td><td>6.000</td><td>cdef</td></tr><tr><td>GILCR</td><td>13.50</td><td>12.000</td><td>15.000</td><td>ghij</td></tr><tr><td>GILCU</td><td>1.00</td><td>1.000</td><td>1.000</td><td>k</td></tr><tr><td>GILFE</td><td>7.50</td><td>7.000</td><td>8.625</td><td>cdlm</td></tr><tr><td>GILGA</td><td>3.00</td><td>3.000</td><td>3.625</td><td>efk</td></tr><tr><td>GILK</td><td>20.00</td><td>19.000</td><td>20.000</td><td>a</td></tr><tr><td>GILMN</td><td>10.00</td><td>8.000</td><td>13.000</td><td>cgilm</td></tr><tr><td>GILMO</td><td>13.00</td><td>12.000</td><td>15.000</td><td>ghil</td></tr><tr><td>GILNB</td><td>13.00</td><td>11.000</td><td>14.000</td><td>ghilm</td></tr><tr><td>GILNV</td><td>4.75</td><td>4.000</td><td>5.000</td><td>defk</td></tr><tr><td>GILRB</td><td>4.50</td><td>4.000</td><td>5.000</td><td>defk</td></tr><tr><td>GILSC</td><td>18.00</td><td>17.000</td><td>18.000</td><td>abj</td></tr><tr><td>GILSR</td><td>8.00</td><td>7.000</td><td>9.000</td><td>cdem</td></tr><tr><td>GILTC</td><td>15.25</td><td>13.125</td><td>17.250</td><td>bghj</td></tr><tr><td>GILTR</td><td>17.00</td><td>15.000</td><td>17.000</td><td>abhj</td></tr><tr><td>GILV</td><td>15.00</td><td>13.875</td><td>16.000</td><td>abghj</td></tr><tr><td>GILY</td><td>9.00</td><td>8.375</td><td>10.000</td><td>cdilm</td></tr><tr><td>GILZN</td><td>2.00</td><td>2.000</td><td>2.000</td><td>fk</td></tr><tr><td>GILZR</td><td>11.00</td><td>10.000</td><td>11.250</td><td>cgilm</td></tr></tbody></table>

<figure><img src="image_13.png" /><figcaption>Figure 11. Box-whisker plot for each of the 20 isotropic textures in the Green color range</figcaption></figure>
<table><caption>Table 20. Summary statistics for each colored texture and a letter symbol: different letters indicate significant differences between colored textures belonging to the set GH based on the selected level of α = 0.05</caption><thead><tr><th>GH</th><th>Median</th><th>1st.Quartile</th><th>3rd.Quartile</th><th>Sig. group</th></tr></thead><tbody><tr><td>GHLCA</td><td>19.00</td><td>18.375</td><td>19.000</td><td>ab</td></tr><tr><td>GHLCO</td><td>6.00</td><td>6.000</td><td>7.125</td><td>cdef</td></tr><tr><td>GHLCR</td><td>15.00</td><td>12.000</td><td>15.625</td><td>aghi</td></tr><tr><td>GHLCU</td><td>1.50</td><td>1.000</td><td>2.000</td><td>c</td></tr><tr><td>GHLFE</td><td>8.00</td><td>7.000</td><td>12.125</td><td>degj</td></tr><tr><td>GHLGA</td><td>3.50</td><td>3.000</td><td>4.000</td><td>cf</td></tr><tr><td>GHLK</td><td>19.75</td><td>18.500</td><td>20.000</td><td>b</td></tr><tr><td>GHLMN</td><td>10.75</td><td>8.000</td><td>13.250</td><td>dghj</td></tr><tr><td>GHLMO</td><td>13.25</td><td>9.875</td><td>15.000</td><td>dghj</td></tr><tr><td>GHLNB</td><td>13.00</td><td>8.750</td><td>14.000</td><td>dghj</td></tr><tr><td>GHLNV</td><td>5.00</td><td>4.000</td><td>6.750</td><td>cef</td></tr><tr><td>GHLRB</td><td>5.00</td><td>4.000</td><td>6.250</td><td>cef</td></tr><tr><td>GHLSC</td><td>18.00</td><td>17.000</td><td>18.000</td><td>ab</td></tr><tr><td>GHLSR</td><td>8.00</td><td>7.375</td><td>9.625</td><td>defj</td></tr><tr><td>GHLTC</td><td>13.75</td><td>11.000</td><td>17.000</td><td>ghij</td></tr><tr><td>GHLTR</td><td>17.00</td><td>16.000</td><td>17.000</td><td>abi</td></tr><tr><td>GHLV</td><td>15.00</td><td>14.000</td><td>16.000</td><td>abhi</td></tr><tr><td>GHLY</td><td>9.00</td><td>7.000</td><td>10.625</td><td>defj</td></tr><tr><td>GHLZN</td><td>2.00</td><td>1.500</td><td>3.250</td><td>c</td></tr><tr><td>GHLZR</td><td>10.00</td><td>8.750</td><td>11.125</td><td>degj</td></tr></tbody></table>

<figure><img src="image_14.png" /><figcaption>Figure 12. Box-whisker plot for each of the 20 horizontal gratings in the Green color range</figcaption></figure>
<table><caption>Table 21. Summary statistics for each colored texture and a letter symbol: different letters indicate significant differences between colored textures belonging to the set GV based on the selected level of α = 0.05</caption><thead><tr><th>GV</th><th>Median</th><th>1st.Quartile</th><th>3rd.Quartile</th><th>Sig. group</th></tr></thead><tbody><tr><td>GVLCA</td><td>19.00</td><td>18.000</td><td>19.000</td><td>ab</td></tr><tr><td>GVLCO</td><td>6.00</td><td>5.000</td><td>6.000</td><td>cde</td></tr><tr><td>GVLCR</td><td>14.00</td><td>12.000</td><td>15.625</td><td>afgh</td></tr><tr><td>GVLCU</td><td>2.50</td><td>2.000</td><td>3.500</td><td>c</td></tr><tr><td>GVLFE</td><td>7.00</td><td>7.000</td><td>9.875</td><td>dei</td></tr><tr><td>GVLGA</td><td>3.00</td><td>2.000</td><td>5.000</td><td>c</td></tr><tr><td>GVLK</td><td>19.75</td><td>18.500</td><td>20.000</td><td>b</td></tr><tr><td>GVLMN</td><td>10.00</td><td>8.500</td><td>11.625</td><td>dfi</td></tr><tr><td>GVLMO</td><td>14.00</td><td>11.375</td><td>15.625</td><td>fgi</td></tr><tr><td>GVLNB</td><td>12.00</td><td>9.500</td><td>14.000</td><td>fi</td></tr><tr><td>GVLNV</td><td>4.00</td><td>3.000</td><td>5.000</td><td>ce</td></tr><tr><td>GVLRB</td><td>5.00</td><td>4.000</td><td>6.000</td><td>ce</td></tr><tr><td>GVLSC</td><td>18.00</td><td>16.500</td><td>18.000</td><td>abh</td></tr><tr><td>GVLSR</td><td>8.00</td><td>7.500</td><td>9.625</td><td>dei</td></tr><tr><td>GVLTC</td><td>15.00</td><td>12.750</td><td>17.000</td><td>fgh</td></tr><tr><td>GVLTR</td><td>16.50</td><td>15.375</td><td>17.000</td><td>abgh</td></tr><tr><td>GVLV</td><td>15.00</td><td>14.000</td><td>16.000</td><td>abfgh</td></tr><tr><td>GVLY</td><td>9.00</td><td>8.500</td><td>10.000</td><td>dei</td></tr><tr><td>GVLZN</td><td>3.00</td><td>1.875</td><td>4.000</td><td>c</td></tr><tr><td>GVLZR</td><td>11.00</td><td>9.875</td><td>11.000</td><td>dfi</td></tr></tbody></table>

<figure><img src="image_15.png" /><figcaption>Figure 13. Box-whisker plot for each of the 20 vertical gratings in the Green color range</figcaption></figure>
## 5.3 Results of the Conover Post-Hoc Tests for the Investigated Textures Type in the Yellow Color Range

<table><caption>Table 22. Summary statistics for each colored texture and a letter symbol: different letters indicate significant differences between colored textures belonging to the set YR based on the selected level of $\alpha = 0.05$</caption><thead><tr><th>YR</th><th>Median</th><th>1st.Quartile</th><th>3rd.Quartile</th><th>Sig. group</th></tr></thead><tbody><tr><td>YRLCA</td><td>19.00</td><td>17.875</td><td>19.000</td><td>ab</td></tr><tr><td>YRLCO</td><td>4.00</td><td>3.000</td><td>5.000</td><td>cd</td></tr><tr><td>YRLCR</td><td>10.00</td><td>9.000</td><td>11.125</td><td>efgh</td></tr><tr><td>YRLCU</td><td>2.75</td><td>2.000</td><td>3.125</td><td>c</td></tr><tr><td>YRLFE</td><td>6.00</td><td>4.375</td><td>7.000</td><td>cde</td></tr><tr><td>YRLGA</td><td>5.00</td><td>3.000</td><td>6.000</td><td>cd</td></tr><tr><td>YRLK</td><td>20.00</td><td>19.000</td><td>20.000</td><td>a</td></tr><tr><td>YRLMN</td><td>7.00</td><td>6.500</td><td>8.000</td><td>cdef</td></tr><tr><td>YRLMO</td><td>17.00</td><td>16.000</td><td>17.500</td><td>abi</td></tr><tr><td>YRLNB</td><td>14.00</td><td>13.000</td><td>15.000</td><td>bgi</td></tr><tr><td>YRLNV</td><td>3.00</td><td>2.000</td><td>4.000</td><td>c</td></tr><tr><td>YRLRB</td><td>7.25</td><td>7.000</td><td>8.000</td><td>cdef</td></tr><tr><td>YRLSC</td><td>16.00</td><td>16.000</td><td>18.000</td><td>abi</td></tr><tr><td>YRLSR</td><td>9.25</td><td>9.000</td><td>10.000</td><td>defh</td></tr><tr><td>YRLTC</td><td>17.25</td><td>16.000</td><td>18.000</td><td>abi</td></tr><tr><td>YRLTR</td><td>14.00</td><td>13.000</td><td>15.000</td><td>bgi</td></tr><tr><td>YRLV</td><td>12.00</td><td>11.375</td><td>13.125</td><td>ghi</td></tr><tr><td>YRLY</td><td>11.00</td><td>9.875</td><td>11.625</td><td>fgh</td></tr><tr><td>YRLZN</td><td>3.00</td><td>2.000</td><td>4.000</td><td>c</td></tr><tr><td>YRLZR</td><td>12.25</td><td>11.000</td><td>14.000</td><td>ghi</td></tr></tbody></table>

<figure><img src="image_16.png" /><figcaption>Figure 14. Box-whisker plot for each of the 20 random-dots textures in the Yellow color range</figcaption></figure>
<table><caption>Table 23. Summary statistics for each colored texture and a letter symbol: different letters indicate significant differences between colored textures belonging to the set YI based on the selected level of α = 0.05</caption><thead><tr><th>YI</th><th>Median</th><th>1st.Quartile</th><th>3rd.Quartile</th><th>Sig. group</th></tr></thead><tbody><tr><td>YILCA</td><td>19.00</td><td>18.000</td><td>19.000</td><td>a</td></tr><tr><td>YILCO</td><td>4.00</td><td>4.000</td><td>5.000</td><td>bcd</td></tr><tr><td>YILCR</td><td>9.75</td><td>9.000</td><td>10.000</td><td>efg</td></tr><tr><td>YILCU</td><td>1.50</td><td>1.000</td><td>2.000</td><td>b</td></tr><tr><td>YILFE</td><td>6.00</td><td>5.000</td><td>6.000</td><td>bcde</td></tr><tr><td>YILGA</td><td>5.00</td><td>5.000</td><td>6.000</td><td>bcde</td></tr><tr><td>YILK</td><td>20.00</td><td>19.000</td><td>20.000</td><td>a</td></tr><tr><td>YILMN</td><td>8.00</td><td>7.000</td><td>8.000</td><td>cef</td></tr><tr><td>YILMO</td><td>17.00</td><td>16.000</td><td>17.125</td><td>ah</td></tr><tr><td>YILNB</td><td>15.00</td><td>14.000</td><td>16.000</td><td>ahi</td></tr><tr><td>YILNV</td><td>2.00</td><td>2.000</td><td>2.500</td><td>bd</td></tr><tr><td>YILRB</td><td>7.00</td><td>7.000</td><td>8.000</td><td>cdef</td></tr><tr><td>YILSC</td><td>17.00</td><td>15.875</td><td>18.000</td><td>ah</td></tr><tr><td>YILSR</td><td>9.00</td><td>9.000</td><td>10.000</td><td>efg</td></tr><tr><td>YILTC</td><td>17.00</td><td>15.500</td><td>18.000</td><td>ah</td></tr><tr><td>YILTR</td><td>14.00</td><td>12.500</td><td>15.000</td><td>ghi</td></tr><tr><td>YILV</td><td>12.00</td><td>12.000</td><td>13.250</td><td>ghi</td></tr><tr><td>YILY</td><td>11.00</td><td>10.500</td><td>11.000</td><td>fgi</td></tr><tr><td>YILZN</td><td>3.00</td><td>2.375</td><td>3.000</td><td>bcd</td></tr><tr><td>YILZR</td><td>13.00</td><td>12.000</td><td>14.125</td><td>ghi</td></tr></tbody></table>

<figure><img src="image_17.png" /><figcaption>Figure 15. Box-whisker plot for each of the 20 isotropic textures in the Yellow color range</figcaption></figure>
<table><caption>Table 24. Summary statistics for each colored texture and a letter symbol: different letters indicate significant differences between colored textures belonging to the set YH based on the selected level of α = 0.05</caption><thead><tr><th>YH</th><th>Median</th><th>1st.Quartile</th><th>3rd.Quartile</th><th>Sig. group</th></tr></thead><tbody><tr><td>YHLCA</td><td>18.25</td><td>17.000</td><td>19.000</td><td>a</td></tr><tr><td>YHLCO</td><td>5.00</td><td>4.000</td><td>6.250</td><td>bc</td></tr><tr><td>YHLCR</td><td>9.00</td><td>8.000</td><td>11.000</td><td>bdef</td></tr><tr><td>YHLCU</td><td>3.00</td><td>2.000</td><td>4.000</td><td>c</td></tr><tr><td>YHLFE</td><td>5.00</td><td>4.000</td><td>6.125</td><td>bc</td></tr><tr><td>YHLGA</td><td>5.00</td><td>3.500</td><td>7.000</td><td>bc</td></tr><tr><td>YHLK</td><td>19.00</td><td>17.875</td><td>20.000</td><td>a</td></tr><tr><td>YHLMN</td><td>7.00</td><td>5.875</td><td>8.000</td><td>bcde</td></tr><tr><td>YHLMO</td><td>16.00</td><td>14.375</td><td>17.000</td><td>aghi</td></tr><tr><td>YHLNB</td><td>15.00</td><td>12.875</td><td>17.000</td><td>afghi</td></tr><tr><td>YHLNV</td><td>3.00</td><td>2.000</td><td>5.000</td><td>c</td></tr><tr><td>YHLRB</td><td>7.75</td><td>5.625</td><td>8.000</td><td>bcd</td></tr><tr><td>YHLSC</td><td>17.50</td><td>15.000</td><td>18.625</td><td>ag</td></tr><tr><td>YHLSR</td><td>10.00</td><td>9.000</td><td>11.125</td><td>defh</td></tr><tr><td>YHLTC</td><td>17.00</td><td>15.875</td><td>17.625</td><td>agi</td></tr><tr><td>YHLTR</td><td>14.00</td><td>11.500</td><td>16.250</td><td>afghi</td></tr><tr><td>YHLV</td><td>13.00</td><td>11.000</td><td>14.000</td><td>fghi</td></tr><tr><td>YHLY</td><td>11.00</td><td>10.000</td><td>13.000</td><td>efhi</td></tr><tr><td>YHLZN</td><td>3.00</td><td>2.000</td><td>4.125</td><td>c</td></tr><tr><td>YHLZR</td><td>12.25</td><td>11.750</td><td>14.000</td><td>fghi</td></tr></tbody></table>

<figure><img src="image_18.png" /><figcaption>Figure 16. Box-whisker plot for each of the 20 horizontal gratings in the Yellow color range</figcaption></figure>
<table><caption>Table 25. Summary statistics for each colored texture and a letter symbol: different letters indicate significant differences between colored textures belonging to the set YV based on the selected level of α = 0.05</caption><thead><tr><th>YV</th><th>Median</th><th>1st.Quartile</th><th>3rd.Quartile</th><th>Sig. group</th></tr></thead><tbody><tr><td>YVLCA</td><td>19.0</td><td>17.00</td><td>19.00</td><td>ab</td></tr><tr><td>YVLCO</td><td>6.0</td><td>4.00</td><td>7.50</td><td>cd</td></tr><tr><td>YVLCR</td><td>9.0</td><td>6.00</td><td>10.00</td><td>cdef</td></tr><tr><td>YVLCU</td><td>6.0</td><td>3.25</td><td>8.25</td><td>cd</td></tr><tr><td>YVLFE</td><td>7.5</td><td>5.50</td><td>9.00</td><td>cde</td></tr><tr><td>YVLGA</td><td>4.0</td><td>2.25</td><td>6.00</td><td>c</td></tr><tr><td>YVLK</td><td>19.0</td><td>17.00</td><td>20.00</td><td>a</td></tr><tr><td>YVLMN</td><td>7.5</td><td>6.00</td><td>9.75</td><td>cde</td></tr><tr><td>YVLMO</td><td>16.5</td><td>15.50</td><td>18.00</td><td>abg</td></tr><tr><td>YVLNB</td><td>14.5</td><td>13.25</td><td>15.25</td><td>abgh</td></tr><tr><td>YVLNV</td><td>6.0</td><td>3.00</td><td>8.25</td><td>cd</td></tr><tr><td>YVLRB</td><td>4.5</td><td>2.75</td><td>8.00</td><td>cd</td></tr><tr><td>YVLSC</td><td>17.0</td><td>15.25</td><td>18.00</td><td>abg</td></tr><tr><td>YVLSR</td><td>5.0</td><td>2.00</td><td>9.00</td><td>cd</td></tr><tr><td>YVLTC</td><td>17.0</td><td>16.50</td><td>19.00</td><td>abg</td></tr><tr><td>YVLTR</td><td>13.0</td><td>12.00</td><td>14.00</td><td>bfgh</td></tr><tr><td>YVLV</td><td>12.0</td><td>11.00</td><td>13.00</td><td>efgh</td></tr><tr><td>YVLY</td><td>11.0</td><td>6.00</td><td>11.00</td><td>defh</td></tr><tr><td>YVLZN</td><td>5.0</td><td>3.00</td><td>8.50</td><td>cd</td></tr><tr><td>YVLZR</td><td>14.0</td><td>12.00</td><td>15.00</td><td>bfgh</td></tr></tbody></table>

<figure><img src="image_19.png" /><figcaption>Figure 17. Box-whisker plot for each of the 20 vertical gratings in the Yellow color range</figcaption></figure>

## Copyrights

Copyright for this article is retained by the author(s), with first publication rights granted to the journal.

This is an open-access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/4.0/).