Hiroharu Natsume,Shogo Okamoto,Hikaru Nagano
TDS Similarity: Outlier Analysis Using a Similarity Index to Compare Time-Series Responses of Temporal Dominance of Sensations Tasks
temporal dominance of sensations,individual analysis,similarity index,discrete time,strawberry
Temporal dominance of sensations (TDS) methods are used to record temporally developing sensations while eating food samples. Results of TDS tasks are typically discussed using averages across multiple trials and panels, and few methods have been developed to analyze differences between individual trials. We defined a similarity index between two time-series responses of TDS tasks. This index adopts a dynamic level to determine the importance of the timing of attribute selection. With a small dynamic level, the index focuses on the duration for attributes to be selected rather than on the timing of the attribute selection. With a large dynamic level, the index focuses on the temporal similarity between two TDS tasks. We performed an outlier analysis based on the developed similarity index using the results of TDS tasks performed in an earlier study. Certain samples were categorized as outliers irrespective of the dynamic level, whereas the categorization of a few samples depended on the level. The similarity index developed in this study achieved individual analyses of TDS tasks, including outlier detection, and adds new analysis techniques to TDS methods.
TDS Similarity: Outlier Analysis Using a Similarity Index to Compare Time-Series Responses of Temporal Dominance of Sensations Tasks

Temporal dominance of sensations (TDS) methods are used to record temporally developing sensations while eating food samples. Results of TDS tasks are typically discussed using averages across multiple trials and panels, and few methods have been developed to analyze differences between individual trials. We defined a similarity index between two time-series responses of TDS tasks. This index adopts a dynamic level to determine the importance of the timing of attribute selection. With a small dynamic level, the index focuses on the duration for attributes to be selected rather than on the timing of the attribute selection. With a large dynamic level, the index focuses on the temporal similarity between two TDS tasks. We performed an outlier analysis based on the developed similarity index using the results of TDS tasks performed in an earlier study. Certain samples were categorized as outliers irrespective of the dynamic level, whereas the categorization of a few samples depended on the level. The similarity index developed in this study achieved individual analyses of TDS tasks, including outlier detection, and adds new analysis techniques to TDS methods.

The temporal dominance of sensations (TDS) method [1,2] enables us to record the temporal evolution of multiple types of sensations, including gustatory, olfactory, and textural feelings while eating food samples. Such time evolution is caused by the physical, chemical, and thermal changes of foods in the mouth, part of which is brought about by saliva [3] and mouth and tongue motion [4], temporal changes in sensory responses including adaptation at peripheral and central nervous levels [5], and cross-modal interactions [6,7,8,9]. For example, the habituation of olfactory sensations is prominent in humans [5,10] and may affect the time-series responses while eating foods. TDS methods have garnered attention in food science since the 2010s because they are more cost-effective than previous time-series sensory evaluation methods such as the time-intensity method [11,12]. The properties of the data acquired through TDS methods and analysis techniques for TDS methods continue to be intensively studied by many researchers [13,14,15]. TDS methods are largely used for food and beverages; however, recent studies have applied the method to other types of stimuli and modalities. Examples include audio [16], visual [17], audiovisual [18], and haptic stimuli [19,20]. Furthermore, there exist modified TDS methods, including the temporal check-all-that-apply method [21,22]. In addition, the TDS method can be used to evaluate temporal and overall food liking or preferences [23,24] and to investigate their temporal drivers. The TDS method is primarily used to determine the tendency of entire panels, and the average of multiple trials from all panels is often discussed. For example, the TDS curve [2], a typical analysis and visualization method for the TDS method, plots the proportion of the panels that select each sensory attribute as dominant over time. Such proportions are called dominance proportions. Inevitably, other analysis methods based on TDS curves, including sensory trajectories [25] and casual analyses [26,27], deal with the average response among the panels as well. These average responses are influenced by outlier samples; however, methods to detect outliers have not been developed for TDS methods. Meyners et al. proposed a method to test whether a panel is statistically different from the other panels using a randomization test [28]. Lepage et al. developed a method to test whether panels differ from each other in their sensory responses and abilities to distinguish food products based on dominance durations that are time periods for descriptive attributes to be selected as dominant [29]. These methods assume that a panel repeatedly evaluates the same foods. Therefore, the methods developed by Meyners [28] and Lepage et al. [29] cannot be applied to individual TDS trials. It is possible to simulate time-series responses of TDS tasks with semi-Markov chains [30,31,32]. Cardot et al. [30] classified the whole panel into several semi-Markov processes according to the characteristics of their responses. However, this corresponded to a classification analysis of panels, and was not intended to discuss individual TDS trials. Individual differences and effects of mental states exist in food preferences and perception [33,34,35,36,37]; hence, a method for analyzing differences among multiple TDS trials is required. However, to date no method has been proposed to investigate the similarity or dissimilarity between different TDS trials. Owing to the quantitative similarities or dissimilarities between multiple TDS trials, several analytical methods that require distance scales between samples can be developed. Representative examples of such methods include outlier analysis and hierarchical clustering. In this study, we developed an index to evaluate the similarity between two TDS trials. This index compares the time series responses of two trials discretized into R intervals using the Euclidean distance. The value of R is the level of dynamics, and indicates the extent to which the temporal aspect should be considered necessary. For small R values, the focus is on the similarity of attribute selection during the trial periods between two trials rather than on the temporal similarity between the trials. We applied this index to the result of TDS tasks on strawberries [38] and examined how the similarity index depends on R. As an essential application of the similarity index, we performed an outlier analysis of TDS tasks for strawberries. To the best of our knowledge, outlier analysis methods have not yet been developed for TDS methods. Even well-controlled experiments cannot avoid outliers in samples. Our study can provide the TDS method with a common data-preprocessing method to detect outliers. The present study is an extension of an earlier study [39] in which the concept of dissimilarity between TDS trials was introduced; however, the variable R was not used. In the present study, the effects of the R values are intensively investigated.

In this section, the typical task of the TDS method and the calculation of dominance proportions are introduced. More information can be acquired from [1,13]. As shown in Figure 1, the tasks of the TDS method are performed using a graphical user interface on a computer. The interface has a maximum of 10–12 buttons with labels of attribute words representing tastes, odors, and mouth feels. In addition, there are two buttons to start and stop the recording. When a panel performs a task, the participant places the food sample in their mouth and simultaneously presses the start button. The panel then selects buttons with attribute words. These buttons correspond to the most dominant sensation that the panel experiences at each moment. Multiple buttons cannot be selected at the same time. After selecting one attribute button, that button is kept selected until another button is selected. The same attribute words can be selected more than one time or never. When the food is swallowed, the panel presses the stop button, and the task is terminated. The computer application records the time when the panel selects the buttons. Individual panels repeat the tasks for the same food sample several times. The result of an example TDS task is shown in Figure 2. The response to attribute i in the jth TDS trial is expressed by a binary function of continuous time t. When attribute i is selected at t in jth trial, ; otherwise, . The continuous time t is the normalized time ranging from 0 to 1. The commencement and end of the individual task correspond to and , respectively. We use an additional function , which equals 1 only before the first attribute button is selected. In other words, after a certain button has been selected, it remains 0 until the end. Hence, at each moment t, the following equation holds where q is the number of attributes used in the TDS task. Multiple trials are averaged using the following equation: where n is the number of trials (the number of panels × repetitions) and represents the proportion with which attribute i is selected at t among all the trials. After being smoothed over time, this proportion is referred to as a dominance proportion [1].

A method for calculating the similarity of the results of two arbitrary TDS tasks is described in this section. The results of the TDS tasks, which are continuous functions, are converted to discrete functions before computing the similarity. This process follows where R is the number of discretized time intervals and is the average value of on the kth discretization interval (), is equivalent to the dominance duration in the kth interval. The interval used to calculate the dominance durations is . When , has the value only for , which represents the dominance duration for the entire period of the trial. This conversion from continuous to discrete forms can be applied to the dominance proportion as well, which is expressed as . We defined the similarity index S between two TDS trials a and b with a certain R value:The second term of (4) corresponds to the average Euclidean distance in a -dimensional space between TDS trials a and b, which is normalized by such that its maximum value is 1. In the -dimensional space, the coordinates of the ith axis are and . The distance corresponds to the dissimilarity of the two trials. The similarity is the complement of the dissimilarity for the 1. In addition, can be applied to the dominance proportion P. For example, indicates the similarity of trial to the centroid of entire trials. We defined the similarity between the two trials; however, the dissimilarity or distance in the second term of (4) may be preferable for different analysis purposes. Certain aspects of the above computational process are similar to the method proposed by Pineau et al. for hypothesis testing of two sets of TDS curves [13].

We used the data from TDS tasks collected by Shimaoka et al. [38]. They conducted the TDS tasks for strawberries by involving seventeen paid university student panels (fourteen males and three females) in their 20s. They were not experts in the sensory evaluation of foods, and familiarized themselves with the TDS tasks, including how to use graphical user interfaces, before the experiment for several minutes until they felt confident. Individual panels replicated three TDS tasks in which one strawberry was eaten in each task. In total, 51 tasks were conducted (). The participants used eight types of attributes () selected by ten participants through a vote among 218 candidate attributes. The eight attributes were sweet, sour, fruity, green, watery, juicy, aromatic, and light. Figure 3 shows the dominance proportions calculated from all the tasks. Typically, in the earlier phase, sweet, juicy, and fruity were dominantly selected. The attribute aromatic was more prominent in the middle phase than in the other phases. Further, sour was frequently selected in the middle phase, and remained dominant in the last period. Similar results were reported in another study of TDS tasks with strawberries [26].

The similarity index varies with R. In particular, the index is almost different between small R values (such as 1–3) and large R values (such as 100). Here, using strawberries as an example, we investigated how the similarity indices are different between different R values. We calculated the correlation coefficients between the similarity index values between different R values to obtain the influences of R values on similarity scores. We assumed the strawberry data in the previous section to be a set of continuous functions for t: , where 51 and 8 are the numbers of TDS tasks or replications and attributes, as used in the previous experiment [38]. As in (2), is the average of all trials or dominance proportions. For given values of R, is the set of the values of for each trial and : where and P are the symbols of the jth TDS trial and mean dominance proportions, respectively. When , the dominance durations of two trials were compared without considering the temporality. When , the similarity was computed for each of the three intervals, that is, the early, middle, and last phases of the task, and the temporality of the dominance proportions was moderately considered. When , provided that the average eating period of a strawberry was about 30 s [38], two pairs of trials were compared at intervals of approximately 3 s. The correlation coefficient between two sets and is defined as where and are the standard deviations of and , respectively, and and are the mean values of and , respectively. The values for are shown in Table 1. When the values of and were close, the two sets and are correlated more strongly. For example, the indices for and are strongly correlated with ; however, those for and are only weakly correlated with . Figure 4 is the graphical version of Table 1, and shows values with and varying from 1 to 100. Similar to Table 1, the correlation coefficient between and is higher when and are closer to each other. When , the correlation coefficient is greater than for any , indicating that the similarity values for are largely or moderately correlated with those for the other R values.

To test the normality of the value distribution in , we applied the Kolmogorov–Smirnov normality test to . When the hypothesized normality was not denied, we judged potential outliers using the range defined by the mean and standard deviation. Trials, that is, members in , were considered potential outliers if they did not fall within the range of , where is the sample mean and is the sample standard deviation calculated among all the values in . This range contains of all the samples, the remaining with small similarity values being screened out. The range was for considering of the samples. Note that the upper side of the distribution does not need to be considered because trials with high similarity represent the population well. When the samples in are judged to not be subjected to the normal distribution, outlier detection methods should be selected depending on the type of distribution. Typically, a method using box plots and interquartile range can be used for determining outliers among non-normally distributed samples. If the mean is calculated across all samples, it can be biased by potential outliers. To avoid such cases, a method based on the minimum covariance determinant [40] may be used.

We performed an outlier analysis on the TDS data for strawberries for each . Figure 5 shows the distribution of the similarity indices between the individual trials and the average, that is, the TDS curves, for each R value. The p-value of the Kolmogorov–Smirnov test was greater than for all R values; hence, we assumed the normal distribution to detect potential outliers. The detected outliers are shown in Table 2 as “+”. Each trial is identified by the panel ID of an alphabetic letter and the iteration count of the task. Trial C3, the third trial of panel C, was classified as an outlier for any R value, while other trials were classified as such only for small, large, or medium R values. For example, N1, that is, the first trial of panel N, was categorized as a potential outlier only when . The time series for each trial are shown in Figure 6.

Here, we discuss or interpret the reasons for certain trials being judged as outliers in the analysis in Section 6.2, which corresponds to the semantic validation of the outlier analysis. In trial N1, was selected for almost the complete duration of the task. As shown in Figure 3, on average was the most dominant attribute in the second half of the task. Nevertheless, N1 was classified as an outlier when . This was because the number of attributes selected during the task was too small in trial N1. On average, several types of attributes, such as , , , and , were selected for sustained periods; however, in trial N1 the dominance durations of all attributes except for and were 0, which distinguishes this trial from the average trial. E1 was classified as an outlier at medium R values. In this trial, attributes that exhibited low proportions in Figure 3 were selected over the entire period, such as and . In average trials, was largely selected in the early phase of the task, as shown in Figure 3. However, in trial E1 was selected even in the middle and last parts of the task. This trial was not considered an outlier when focusing on more chronological details with large R values, such as 50 and 100. In attribute selection, C2 and C3 exhibited similar patterns. However, C2 was classified as an outlier only when , whereas C3 was considered an outlier for every R value, as shown in Table 2. For and 100, C2 exhibited small similarity values; however, they were barely within the range of the inliers. In both trials, was selected in the later period, in which was not prominent on average. In the early period, however, typical attributes such as and were selected in both C2 and C3. In C3, and were selected earlier than in the average trial, as shown in Figure 3. Hence, C3 was judged as an outlier with large R values, with which the selected attributes were chronically compared. Further, in C3, , which was selected in the majority of the trials, was not selected at all. This caused C3 to be judged as an outlier with small R, in which dominance durations are more important than timing. As aforementioned, four trials that were detected as outliers with at least one R value seem distinct from the typical panel behavior shown in Figure 3. In the above outlier analysis, different values of R led to different outliers. This means that the developed similarity index changes the level of dynamics considered depending on the R value.

A key feature of the developed similarity index is that it can change the degree of emphasis on timing depending on the value of R. We assume that an appropriate R value depends on the research target, which varies according to the purpose of the analysis and the type of foods to be tested. For example, a low R value may be used for foods with large individual differences in eating time, such as chewy foods. For such foods, closely focusing on time-series is not meaningful. In contrast, a large R value may be preferred if researchers wish to conduct a detailed investigation of the temporal changes of sensations. As a systematic method, we propose setting R based on the correlation coefficients of the similarity values between the different R values. In the example of strawberries, as shown in Figure 4, the correlation coefficients between the index values at and those at the other R values ranging from 1–100 exhibit constant high values. Such R (in this case, 20) is recommendable because it is in accordance with both small R values, where the timings of attribute selection are not deemed important, and large R values, where the timings are considered important. For strawberries, the average task duration was approximately 30 s [38]. Hence, for , one interval corresponds to about 1.5 s. Previous studies have discussed TDS curves by dividing the entire trial period into three or four equal parts; in other words, they used [13,29] or [41]. However, these studies not provide suggestions on the question of how many intervals the entire period should be split into. One notable criterion for determining the R value was suggested in [26], where continuous TDS curves were discretized with an interval of approximately 1 s for causality analysis. The authors referred to the minimum interval between two successive button selections during the TDS tasks. The panels largely switched attributes in more than 1 s. This can be a reasonable method for determining R. It may not be beneficial to discuss and analyze the behaviors of panels when the temporal resolution is smaller than the minimum behavioral response time. Another perspective on how to set an appropriate R value involves dummy outliers. For example, in our case several TDS samples of non-strawberry foods could have been mixed with those for strawberries. The samples of non-strawberry foods could then be perfectly determined as outliers, with the dummy TDS samples used for performance evaluation of the outlier analysis. In this study, only strawberries were used as food samples. Therefore, at present the nature of the index discussed above only holds for strawberries. Furthermore, the data used in this study were collected from panels with a limited background. The generalizability of the index must be confirmed by further testing with a wider variety of foods and panels. To compute the similarity between trials, we used the Euclidean distance. Other popular distance measures include the Manhattan distance and distances based on probability. The effect of the choice of distance scale on the similarity index should be further studied. The similarity index developed in this study treats all attributes as having the same importance. However, it may be practical to use only essential attributes or to weight important attributes in order to compute similarity. Moreover, as demonstrated in this study, outlier analysis based on the developed similarity index can detect potentially outlying samples. However, it does not show how such samples differ from typical samples. Thus, a supplementary method may be desirable to help determine whether samples detected as outliers should be removed.

In the TDS method, the average of all trials is discussed generally, and few analyses of individual trials are conducted. Therefore, in this paper we have proposed a similarity index that can be used to compare individual TDS trials. This index can adjust the importance of the temporal similarity by changing the R parameter, which represents the number of discretizations. Using this index, we performed an outlier analysis for TDS on strawberry panel data. Different trials were classified as outliers with and without an emphasis on time series. These trials were confirmed to be vastly different from the average trial. Outlier analyses can be employed in all experiments using the TDS method to exclude outliers that substantially bias the average consumer responses. In addition, this index can be useful for other individual analyses that use distance measures, such as hierarchical clustering [42]. Challenges to be solved in the future include determining the appropriate R value for different analyses. The method used for outlier analysis can be improved in the future as well. Currently, it does not indicate the reasons behind certain samples being judged as outliers; thus, a method that is better able to explain these reasons is required.