Tactile Evaluation Feedback System for Multi-Layered Structure Inspired by Human Tactile Perception Mechanism

Tactile sensation is one type of valuable feedback in evaluating a product. Conventionally, sensory evaluation is used to get direct subjective responses from the consumers, in order to improve the product’s quality. However, this method is a time-consuming and costly process. Therefore, this paper proposes a novel tactile evaluation system that can give tactile feedback from a sensor’s output. The main concept of this system is hierarchically layering the tactile sensation, which is inspired by the flow of human perception. The tactile sensation is classified from low-order of tactile sensation (LTS) to high-order of tactile sensation (HTS), and also to preference. Here, LTS will be correlated with physical measures. Furthermore, the physical measures that are used to correlate with LTS are selected based on four main aspects of haptic information (roughness, compliance, coldness, and slipperiness), which are perceived through human tactile sensors. By using statistical analysis, the correlation between each hierarchy was obtained, and the preference was derived in terms of physical measures. A verification test was conducted by using unknown samples to determine the reliability of the system. The results showed that the system developed was capable of estimating preference with an accuracy of approximately 80%.


Introduction
A product's tactile sensation is one of the imperative factors in making the decision to purchase a product, besides its functionality and usability [1]. Moreover, McCabe and Nowlis also prove that tactile cues significantly affect impulse purchasing [2]. Apart from that, with the presence of tactile cue, the confidence level in product evaluations increases [3]. Today product developers increasingly recognize the importance of tactile sensation as one of the added values against rivals in a competitive market [2]. For example, in development of the interiors of automobiles especially, for the parts that are directly touched by the consumer, such as seats, armrests and internal panels, the product designers choose materials from fabric to leather and plastic, and then finish them with unique kinds of texture that influence people's subjective responses and the value of the product itself [4,5]. Generally, sensory evaluation as the mainstream method is performed to evaluate subjective responses, because this method provides direct responses from the consumer. However, this method may not be reliable

Concept
In order to develop this evaluation method, human subjective responses (dependent variables) need to be collected and correlated with physical measures (independent variables) from designated sensors. This method is a basic method of quantifying or translating subjective responses to physical measures, and many previous works used this approach [5,11,13,[18][19][20][21][22][23][24][25]. One of the extensions of this method is to construct a hierarchy of stages of subjective responses, for example, by classifying subjective responses into two categories, which are sensorial and affective responses [21]. All of these approaches have a similar reasoning behind them. It is based on the biological structure of human perception and cognition. When a surface is being touched, it stimulates mechanoreceptors and thermoreceptors beneath the skin [26], and the signals are mapped onto outer areas of the brain. Here, psychophysical responses of the stimuli are processed, for example, softness, roughness, warmness and so on of the object. Then, these responses trigger other areas of the brain where they are synthesized, and compared to experiences, to form affective responses [27]. In addition, Chen et al. [21] and Nagano et al. [28] suggested dividing into another two layers of affective responses which were labelled affective and preferential layers, for better understanding of complex human affective responses.
Based on the above discussion, this paper employs the concept of hierarchically layering the subjective responses and correlates the psychophysical responses with physical measures as shown in Figure 1. Note that this paper addresses subjective responses and physical measures as below. For subjective responses, psychophysical response is labelled as low-order of tactile sensation (LTS).
Based on the above discussion, this paper employs the concept of hierarchically layering the subjective responses and correlates the psychophysical responses with physical measures as shown in Figure 1. Note that this paper addresses subjective responses and physical measures as below. For subjective responses, psychophysical response is labelled as low-order of tactile sensation (LTS).
Affective response as high-order of tactile sensation (HTS) and preferential response as preference. LTS is a layer that includes terms which directly express physical properties in the context of tactile sense, such as cool, smooth, soft, and so on. On the other hand, HTS is a layer that is comprised of emotional and affective expressions, such as fit, refreshing, luxury, and so on. Lastly, preference is a layer that describes personal preferences which are strongly related to decision-making process while purchasing a product, such as prefer, pleasure, comfort, and so on. Conventionally, haptic information can be generalized into four main aspects which are roughness, compliance, coldness, and slipperiness [21,29]. Roughness is the most studied aspect compared to others because it is the most vital aspect in differentiating between surface textures. Previous works [11,30] showed that the vibration which is evoked during dynamic interaction between skin and the object's surface is considered to have a major role in roughness perception. Compliance has a number of ways to physically express it, and the main streams are object's stiffness and surface deformation/Young's modulus [29]. In a detailed investigation [31], the researchers discovered that 90 percent of the information in perceiving compliance was associated with the perception of surface deformation. Coldness is perceived due to the rate of heat loss from the skin to the object when touching it and it is mainly influenced by the material's thermal properties, such as thermal conductivity, specific heat, and density [32]. Lastly, slipperiness is perceived through two channels: kinesthetic (when an opposite force works in the opposite direction of the motion) and cutaneous (when skin stretches tangentially to the surface) [33]. These are related to the friction that occurs during the interaction between two surfaces: skin and surface of material. In conclusion, this paper measures four physical measures which correspond to four main aspects of haptic information; vibration for roughness, bulk displacement for surface deformation for compliance, thermal properties for coldness, and friction for slipperiness. The specific measurement methods that are used in this work are reported here.
This hierarchical structure of tactile perception process mentioned above may not perfectly correspond to the process in our brain; however, it may help product development in understanding the overall structure of tactile perception mechanism in human. The correlation between LTS-HTS and HTS-Preference will be achieved from sensory evaluation data (refer Section 5). Then, the relationship between LTS and physical measures will be obtained by correlating sensory evaluation data of LTS and physical measures data. Once each relation between hierarchies is obtained, we may estimate the human preference for unknown samples by measuring the physical quantities. Conventionally, haptic information can be generalized into four main aspects which are roughness, compliance, coldness, and slipperiness [21,29]. Roughness is the most studied aspect compared to others because it is the most vital aspect in differentiating between surface textures. Previous works [11,30] showed that the vibration which is evoked during dynamic interaction between skin and the object's surface is considered to have a major role in roughness perception. Compliance has a number of ways to physically express it, and the main streams are object's stiffness and surface deformation/Young's modulus [29]. In a detailed investigation [31], the researchers discovered that 90 percent of the information in perceiving compliance was associated with the perception of surface deformation. Coldness is perceived due to the rate of heat loss from the skin to the object when touching it and it is mainly influenced by the material's thermal properties, such as thermal conductivity, specific heat, and density [32]. Lastly, slipperiness is perceived through two channels: kinesthetic (when an opposite force works in the opposite direction of the motion) and cutaneous (when skin stretches tangentially to the surface) [33]. These are related to the friction that occurs during the interaction between two surfaces: skin and surface of material. In conclusion, this paper measures four physical measures which correspond to four main aspects of haptic information; vibration for roughness, bulk displacement for surface deformation for compliance, thermal properties for coldness, and friction for slipperiness. The specific measurement methods that are used in this work are reported here.

Method
This hierarchical structure of tactile perception process mentioned above may not perfectly correspond to the process in our brain; however, it may help product development in understanding the overall structure of tactile perception mechanism in human. The correlation between LTS-HTS and HTS-Preference will be achieved from sensory evaluation data (refer Section 5). Then, the relationship between LTS and physical measures will be obtained by correlating sensory evaluation data of LTS and physical measures data. Once each relation between hierarchies is obtained, we may estimate the human preference for unknown samples by measuring the physical quantities.

Method
A sensory evaluation with semantic differential method was carried out. 15 adults with age between their twenties and forties were asked to touch freely with their hands and evaluate 26 samples of door armrest on a seven-point unipolar scale. All participants had given their informed consent for inclusion before they participated in the study, and it was approved on 7 July 2016 by the School  . The unipolar scale is used to avoid translation problems between opposite adjectives [34]. Moreover, the scale is only defined at the endpoints to prevent varying interpretations of verbal anchors and unevenness between anchors [35]. The sensory evaluation included 22 items of adjectives listed and classified into LTS, HTS and preference as shown in Table 1. The adjectives were selected by referring to previous works [20,22,28,[36][37][38] and discussion with door-armrest developers. In order to exclude visual effects from the sensory evaluation, the experiment was carried out as a blind test with samples' details undisclosed. Table 2 shows the outer layer samples which are made of synthetic leather, genuine leather, fabric, polyvinyl chloride and resin, and the cross-sectional structure types of each sample; there are five types as shown in Figure 2. A sensory evaluation with semantic differential method was carried out. 15 adults with age between their twenties and forties were asked to touch freely with their hands and evaluate 26 samples of door armrest on a seven-point unipolar scale. All participants had given their informed consent for inclusion before they participated in the study, and it was approved on 7 July 2016 by the School of Science and Technology, Keio University Research Ethics Committee . The unipolar scale is used to avoid translation problems between opposite adjectives [34]. Moreover, the scale is only defined at the endpoints to prevent varying interpretations of verbal anchors and unevenness between anchors [35]. The sensory evaluation included 22 items of adjectives listed and classified into LTS, HTS and preference as shown in Table 1. The adjectives were selected by referring to previous works [20,22,28,[36][37][38] and discussion with door-armrest developers. In order to exclude visual effects from the sensory evaluation, the experiment was carried out as a blind test with samples' details undisclosed. Table 2 shows the outer layer samples which are made of synthetic leather, genuine leather, fabric, polyvinyl chloride and resin, and the cross-sectional structure types of each sample; there are five types as shown in Figure 2. A sensory evaluation with semantic differential method was carried out. 15 adults with age between their twenties and forties were asked to touch freely with their hands and evaluate 26 samples of door armrest on a seven-point unipolar scale. All participants had given their informed consent for inclusion before they participated in the study, and it was approved on 7 July 2016 by the School of Science and Technology, Keio University Research Ethics Committee . The unipolar scale is used to avoid translation problems between opposite adjectives [34]. Moreover, the scale is only defined at the endpoints to prevent varying interpretations of verbal anchors and unevenness between anchors [35]. The sensory evaluation included 22 items of adjectives listed and classified into LTS, HTS and preference as shown in Table 1. The adjectives were selected by referring to previous works [20,22,28,[36][37][38] and discussion with door-armrest developers. In order to exclude visual effects from the sensory evaluation, the experiment was carried out as a blind test with samples' details undisclosed. Table 2 shows the outer layer samples which are made of synthetic leather, genuine leather, fabric, polyvinyl chloride and resin, and the cross-sectional structure types of each sample; there are five types as shown in Figure 2. A sensory evaluation with semantic differential method was carried out. 15 adults with age between their twenties and forties were asked to touch freely with their hands and evaluate 26 samples of door armrest on a seven-point unipolar scale. All participants had given their informed consent for inclusion before they participated in the study, and it was approved on 7 July 2016 by the School of Science and Technology, Keio University Research Ethics Committee . The unipolar scale is used to avoid translation problems between opposite adjectives [34]. Moreover, the scale is only defined at the endpoints to prevent varying interpretations of verbal anchors and unevenness between anchors [35]. The sensory evaluation included 22 items of adjectives listed and classified into LTS, HTS and preference as shown in Table 1. The adjectives were selected by referring to previous works [20,22,28,[36][37][38] and discussion with door-armrest developers. In order to exclude visual effects from the sensory evaluation, the experiment was carried out as a blind test with samples' details undisclosed. Table 2 shows the outer layer samples which are made of synthetic leather, genuine leather, fabric, polyvinyl chloride and resin, and the cross-sectional structure types of each sample; there are five types as shown in Figure 2. A sensory evaluation with semantic differential method was carried out. 15 adults with age between their twenties and forties were asked to touch freely with their hands and evaluate 26 samples of door armrest on a seven-point unipolar scale. All participants had given their informed consent for inclusion before they participated in the study, and it was approved on 7 July 2016 by the School of Science and Technology, Keio University Research Ethics Committee . The unipolar scale is used to avoid translation problems between opposite adjectives [34]. Moreover, the scale is only defined at the endpoints to prevent varying interpretations of verbal anchors and unevenness between anchors [35]. The sensory evaluation included 22 items of adjectives listed and classified into LTS, HTS and preference as shown in Table 1. The adjectives were selected by referring to previous works [20,22,28,[36][37][38] and discussion with door-armrest developers. In order to exclude visual effects from the sensory evaluation, the experiment was carried out as a blind test with samples' details undisclosed. Table 2 shows the outer layer samples which are made of synthetic leather, genuine leather, fabric, polyvinyl chloride and resin, and the cross-sectional structure types of each sample; there are five types as shown in Figure 2. A sensory evaluation with semantic differential method was carried out. 15 adults with age between their twenties and forties were asked to touch freely with their hands and evaluate 26 samples of door armrest on a seven-point unipolar scale. All participants had given their informed consent for inclusion before they participated in the study, and it was approved on 7 July 2016 by the School of Science and Technology, Keio University Research Ethics Committee . The unipolar scale is used to avoid translation problems between opposite adjectives [34]. Moreover, the scale is only defined at the endpoints to prevent varying interpretations of verbal anchors and unevenness between anchors [35]. The sensory evaluation included 22 items of adjectives listed and classified into LTS, HTS and preference as shown in Table 1. The adjectives were selected by referring to previous works [20,22,28,[36][37][38] and discussion with door-armrest developers. In order to exclude visual effects from the sensory evaluation, the experiment was carried out as a blind test with samples' details undisclosed. Table 2 shows the outer layer samples which are made of synthetic leather, genuine leather, fabric, polyvinyl chloride and resin, and the cross-sectional structure types of each sample; there are five types as shown in Figure 2. A sensory evaluation with semantic differential method was carried out. 15 adults with age between their twenties and forties were asked to touch freely with their hands and evaluate 26 samples of door armrest on a seven-point unipolar scale. All participants had given their informed consent for inclusion before they participated in the study, and it was approved on 7 July 2016 by the School of Science and Technology, Keio University Research Ethics Committee . The unipolar scale is used to avoid translation problems between opposite adjectives [34]. Moreover, the scale is only defined at the endpoints to prevent varying interpretations of verbal anchors and unevenness between anchors [35]. The sensory evaluation included 22 items of adjectives listed and classified into LTS, HTS and preference as shown in Table 1. The adjectives were selected by referring to previous works [20,22,28,[36][37][38] and discussion with door-armrest developers. In order to exclude visual effects from the sensory evaluation, the experiment was carried out as a blind test with samples' details undisclosed. Table 2 shows the outer layer samples which are made of synthetic leather, genuine leather, fabric, polyvinyl chloride and resin, and the cross-sectional structure types of each sample; there are five types as shown in Figure 2. A sensory evaluation with semantic differential method was carried out. 15 adults with age between their twenties and forties were asked to touch freely with their hands and evaluate 26 samples of door armrest on a seven-point unipolar scale. All participants had given their informed consent for inclusion before they participated in the study, and it was approved on 7 July 2016 by the School of Science and Technology, Keio University Research Ethics Committee . The unipolar scale is used to avoid translation problems between opposite adjectives [34]. Moreover, the scale is only defined at the endpoints to prevent varying interpretations of verbal anchors and unevenness between anchors [35]. The sensory evaluation included 22 items of adjectives listed and classified into LTS, HTS and preference as shown in Table 1. The adjectives were selected by referring to previous works [20,22,28,[36][37][38] and discussion with door-armrest developers. In order to exclude visual effects from the sensory evaluation, the experiment was carried out as a blind test with samples' details undisclosed. Table 2 shows the outer layer samples which are made of synthetic leather, genuine leather, fabric, polyvinyl chloride and resin, and the cross-sectional structure types of each sample; there are five types as shown in Figure 2. A sensory evaluation with semantic differential method was carried out. 15 adults with age between their twenties and forties were asked to touch freely with their hands and evaluate 26 samples of door armrest on a seven-point unipolar scale. All participants had given their informed consent for inclusion before they participated in the study, and it was approved on 7 July 2016 by the School of Science and Technology, Keio University Research Ethics Committee . The unipolar scale is used to avoid translation problems between opposite adjectives [34]. Moreover, the scale is only defined at the endpoints to prevent varying interpretations of verbal anchors and unevenness between anchors [35]. The sensory evaluation included 22 items of adjectives listed and classified into LTS, HTS and preference as shown in Table 1. The adjectives were selected by referring to previous works [20,22,28,[36][37][38] and discussion with door-armrest developers. In order to exclude visual effects from the sensory evaluation, the experiment was carried out as a blind test with samples' details undisclosed. Table 2 shows the outer layer samples which are made of synthetic leather, genuine leather, fabric, polyvinyl chloride and resin, and the cross-sectional structure types of each sample; there are five types as shown in Figure 2. A sensory evaluation with semantic differential method was carried out. 15 adults with age between their twenties and forties were asked to touch freely with their hands and evaluate 26 samples of door armrest on a seven-point unipolar scale. All participants had given their informed consent for inclusion before they participated in the study, and it was approved on 7 July 2016 by the School of Science and Technology, Keio University Research Ethics Committee . The unipolar scale is used to avoid translation problems between opposite adjectives [34]. Moreover, the scale is only defined at the endpoints to prevent varying interpretations of verbal anchors and unevenness between anchors [35]. The sensory evaluation included 22 items of adjectives listed and classified into LTS, HTS and preference as shown in Table 1. The adjectives were selected by referring to previous works [20,22,28,[36][37][38] and discussion with door-armrest developers. In order to exclude visual effects from the sensory evaluation, the experiment was carried out as a blind test with samples' details undisclosed. Table 2 shows the outer layer samples which are made of synthetic leather, genuine leather, fabric, polyvinyl chloride and resin, and the cross-sectional structure types of each sample; there are five types as shown in Figure 2. A sensory evaluation with semantic differential method was carried out. 15 adults with age between their twenties and forties were asked to touch freely with their hands and evaluate 26 samples of door armrest on a seven-point unipolar scale. All participants had given their informed consent for inclusion before they participated in the study, and it was approved on 7 July 2016 by the School of Science and Technology, Keio University Research Ethics Committee . The unipolar scale is used to avoid translation problems between opposite adjectives [34]. Moreover, the scale is only defined at the endpoints to prevent varying interpretations of verbal anchors and unevenness between anchors [35]. The sensory evaluation included 22 items of adjectives listed and classified into LTS, HTS and preference as shown in Table 1. The adjectives were selected by referring to previous works [20,22,28,[36][37][38] and discussion with door-armrest developers. In order to exclude visual effects from the sensory evaluation, the experiment was carried out as a blind test with samples' details undisclosed. Table 2 shows the outer layer samples which are made of synthetic leather, genuine leather, fabric, polyvinyl chloride and resin, and the cross-sectional structure types of each sample; there are five types as shown in Figure 2. A sensory evaluation with semantic differential method was carried out. 15 adults with age between their twenties and forties were asked to touch freely with their hands and evaluate 26 samples of door armrest on a seven-point unipolar scale. All participants had given their informed consent for inclusion before they participated in the study, and it was approved on 7 July 2016 by the School of Science and Technology, Keio University Research Ethics Committee . The unipolar scale is used to avoid translation problems between opposite adjectives [34]. Moreover, the scale is only defined at the endpoints to prevent varying interpretations of verbal anchors and unevenness between anchors [35]. The sensory evaluation included 22 items of adjectives listed and classified into LTS, HTS and preference as shown in Table 1. The adjectives were selected by referring to previous works [20,22,28,[36][37][38] and discussion with door-armrest developers. In order to exclude visual effects from the sensory evaluation, the experiment was carried out as a blind test with samples' details undisclosed. Table 2 shows the outer layer samples which are made of synthetic leather, genuine leather, fabric, polyvinyl chloride and resin, and the cross-sectional structure types of each sample; there are five types as shown in Figure 2. A sensory evaluation with semantic differential method was carried out. 15 adults with age between their twenties and forties were asked to touch freely with their hands and evaluate 26 samples of door armrest on a seven-point unipolar scale. All participants had given their informed consent for inclusion before they participated in the study, and it was approved on 7 July 2016 by the School of Science and Technology, Keio University Research Ethics Committee . The unipolar scale is used to avoid translation problems between opposite adjectives [34]. Moreover, the scale is only defined at the endpoints to prevent varying interpretations of verbal anchors and unevenness between anchors [35]. The sensory evaluation included 22 items of adjectives listed and classified into LTS, HTS and preference as shown in Table 1. The adjectives were selected by referring to previous works [20,22,28,[36][37][38] and discussion with door-armrest developers. In order to exclude visual effects from the sensory evaluation, the experiment was carried out as a blind test with samples' details undisclosed. Table 2 shows the outer layer samples which are made of synthetic leather, genuine leather, fabric, polyvinyl chloride and resin, and the cross-sectional structure types of each sample; there are five types as shown in Figure 2. A sensory evaluation with semantic differential method was carried out. 15 adults with age between their twenties and forties were asked to touch freely with their hands and evaluate 26 samples of door armrest on a seven-point unipolar scale. All participants had given their informed consent for inclusion before they participated in the study, and it was approved on 7 July 2016 by the School of Science and Technology, Keio University Research Ethics Committee . The unipolar scale is used to avoid translation problems between opposite adjectives [34]. Moreover, the scale is only defined at the endpoints to prevent varying interpretations of verbal anchors and unevenness between anchors [35]. The sensory evaluation included 22 items of adjectives listed and classified into LTS, HTS and preference as shown in Table 1. The adjectives were selected by referring to previous works [20,22,28,[36][37][38] and discussion with door-armrest developers. In order to exclude visual effects from the sensory evaluation, the experiment was carried out as a blind test with samples' details undisclosed. Table 2 shows the outer layer samples which are made of synthetic leather, genuine leather, fabric, polyvinyl chloride and resin, and the cross-sectional structure types of each sample; there are five types as shown in Figure 2. A sensory evaluation with semantic differential method was carried out. 15 adults with age between their twenties and forties were asked to touch freely with their hands and evaluate 26 samples of door armrest on a seven-point unipolar scale. All participants had given their informed consent for inclusion before they participated in the study, and it was approved on 7 July 2016 by the School of Science and Technology, Keio University Research Ethics Committee . The unipolar scale is used to avoid translation problems between opposite adjectives [34]. Moreover, the scale is only defined at the endpoints to prevent varying interpretations of verbal anchors and unevenness between anchors [35]. The sensory evaluation included 22 items of adjectives listed and classified into LTS, HTS and preference as shown in Table 1. The adjectives were selected by referring to previous works [20,22,28,[36][37][38] and discussion with door-armrest developers. In order to exclude visual effects from the sensory evaluation, the experiment was carried out as a blind test with samples' details undisclosed. Table 2 shows the outer layer samples which are made of synthetic leather, genuine leather, fabric, polyvinyl chloride and resin, and the cross-sectional structure types of each sample; there are five types as shown in Figure 2. A sensory evaluation with semantic differential method was carried out. 15 adults with age between their twenties and forties were asked to touch freely with their hands and evaluate 26 samples of door armrest on a seven-point unipolar scale. All participants had given their informed consent for inclusion before they participated in the study, and it was approved on 7 July 2016 by the School of Science and Technology, Keio University Research Ethics Committee . The unipolar scale is used to avoid translation problems between opposite adjectives [34]. Moreover, the scale is only defined at the endpoints to prevent varying interpretations of verbal anchors and unevenness between anchors [35]. The sensory evaluation included 22 items of adjectives listed and classified into LTS, HTS and preference as shown in Table 1. The adjectives were selected by referring to previous works [20,22,28,[36][37][38] and discussion with door-armrest developers. In order to exclude visual effects from the sensory evaluation, the experiment was carried out as a blind test with samples' details undisclosed. Table 2 shows the outer layer samples which are made of synthetic leather, genuine leather, fabric, polyvinyl chloride and resin, and the cross-sectional structure types of each sample; there are five types as shown in Figure 2. A sensory evaluation with semantic differential method was carried out. 15 adults with age between their twenties and forties were asked to touch freely with their hands and evaluate 26 samples of door armrest on a seven-point unipolar scale. All participants had given their informed consent for inclusion before they participated in the study, and it was approved on 7 July 2016 by the School of Science and Technology, Keio University Research Ethics Committee . The unipolar scale is used to avoid translation problems between opposite adjectives [34]. Moreover, the scale is only defined at the endpoints to prevent varying interpretations of verbal anchors and unevenness between anchors [35]. The sensory evaluation included 22 items of adjectives listed and classified into LTS, HTS and preference as shown in Table 1. The adjectives were selected by referring to previous works [20,22,28,[36][37][38] and discussion with door-armrest developers. In order to exclude visual effects from the sensory evaluation, the experiment was carried out as a blind test with samples' details undisclosed. Table 2 shows the outer layer samples which are made of synthetic leather, genuine leather, fabric, polyvinyl chloride and resin, and the cross-sectional structure types of each sample; there are five types as shown in Figure 2.

Principal Component Analysis
From the survey data of sensory evaluation, principal component analysis with varimax rotation was carried out using statistical analysis software (SPSS Ver. 22, IBM, Armonk, NY, USA) to reduce the number of variables by grouping the adjectives that have strong correlation to each other into an independent semantic variable or principal component (PC). Previous research [21,39,40] has shown that the adjectives can be grouped, and principal component analysis is one of the methods used. By reducing the number of variables, they become practical and easy for interpreting the data. Tables 3 and 4 show the results of principal component analysis for LTS and HTS, respectively. All components that have loadings higher than 0.50 are in bold (this is an arbitrary limit).

Principal Component Analysis
From the survey data of sensory evaluation, principal component analysis with varimax rotation was carried out using statistical analysis software (SPSS Ver. 22, IBM, Armonk, NY, USA) to reduce the number of variables by grouping the adjectives that have strong correlation to each other into an independent semantic variable or principal component (PC). Previous research [21,39,40] has shown that the adjectives can be grouped, and principal component analysis is one of the methods used. By reducing the number of variables, they become practical and easy for interpreting the data. Tables 3 and 4 show the results of principal component analysis for LTS and HTS, respectively. All components that have loadings higher than 0.50 are in bold (this is an arbitrary limit).

Principal Component Analysis
From the survey data of sensory evaluation, principal component analysis with varimax rotation was carried out using statistical analysis software (SPSS Ver. 22, IBM, Armonk, NY, USA) to reduce the number of variables by grouping the adjectives that have strong correlation to each other into an independent semantic variable or principal component (PC). Previous research [21,39,40] has shown that the adjectives can be grouped, and principal component analysis is one of the methods used. By reducing the number of variables, they become practical and easy for interpreting the data. Tables 3 and 4 show the results of principal component analysis for LTS and HTS, respectively. All components that have loadings higher than 0.50 are in bold (this is an arbitrary limit).

Principal Component Analysis
From the survey data of sensory evaluation, principal component analysis with varimax rotation was carried out using statistical analysis software (SPSS Ver. 22, IBM, Armonk, NY, USA) to reduce the number of variables by grouping the adjectives that have strong correlation to each other into an independent semantic variable or principal component (PC). Previous research [21,39,40] has shown that the adjectives can be grouped, and principal component analysis is one of the methods used. By reducing the number of variables, they become practical and easy for interpreting the data. Tables 3 and 4 show the results of principal component analysis for LTS and HTS, respectively. All components that have loadings higher than 0.50 are in bold (this is an arbitrary limit).

Principal Component Analysis
From the survey data of sensory evaluation, principal component analysis with varimax rotation was carried out using statistical analysis software (SPSS Ver. 22, IBM, Armonk, NY, USA) to reduce the number of variables by grouping the adjectives that have strong correlation to each other into an independent semantic variable or principal component (PC). Previous research [21,39,40] has shown that the adjectives can be grouped, and principal component analysis is one of the methods used. By reducing the number of variables, they become practical and easy for interpreting the data. Tables 3 and 4 show the results of principal component analysis for LTS and HTS, respectively. All components that have loadings higher than 0.50 are in bold (this is an arbitrary limit).

Principal Component Analysis
From the survey data of sensory evaluation, principal component analysis with varimax rotation was carried out using statistical analysis software (SPSS Ver. 22, IBM, Armonk, NY, USA) to reduce the number of variables by grouping the adjectives that have strong correlation to each other into an independent semantic variable or principal component (PC). Previous research [21,39,40] has shown that the adjectives can be grouped, and principal component analysis is one of the methods used. By reducing the number of variables, they become practical and easy for interpreting the data. Tables 3 and 4 show the results of principal component analysis for LTS and HTS, respectively. All components that have loadings higher than 0.50 are in bold (this is an arbitrary limit).

Principal Component Analysis
From the survey data of sensory evaluation, principal component analysis with varimax rotation was carried out using statistical analysis software (SPSS Ver. 22, IBM, Armonk, NY, USA) to reduce the number of variables by grouping the adjectives that have strong correlation to each other into an independent semantic variable or principal component (PC). Previous research [21,39,40] has shown that the adjectives can be grouped, and principal component analysis is one of the methods used. By reducing the number of variables, they become practical and easy for interpreting the data. Tables 3 and 4 show the results of principal component analysis for LTS and HTS, respectively. All components that have loadings higher than 0.50 are in bold (this is an arbitrary limit).

Principal Component Analysis
From the survey data of sensory evaluation, principal component analysis with varimax rotation was carried out using statistical analysis software (SPSS Ver. 22, IBM, Armonk, NY, USA) to reduce the number of variables by grouping the adjectives that have strong correlation to each other into an independent semantic variable or principal component (PC). Previous research [21,39,40] has shown that the adjectives can be grouped, and principal component analysis is one of the methods used. By reducing the number of variables, they become practical and easy for interpreting the data. Tables 3 and 4 show the results of principal component analysis for LTS and HTS, respectively. All components that have loadings higher than 0.50 are in bold (this is an arbitrary limit). There were six principal components for LTS extracted with 91.7% of the total variance. PC1 was interpreted as "dampness" dimension, with high loadings on damp and wet, PC2 was "coldness" Sensors 2017, 17, 2601 6 of 14 dimension, with high loadings on chilly and cold. PC3 was "micro-roughness" dimension, with high loadings on smooth and silky. PC4 was "macro-roughness" dimension, with high loadings on rough and bumpy. PC5 was "hardness" dimension, with high loadings on tough and hard. Lastly, PC6 was "hollowness" dimension, with high loadings on brittle and hollow. On the other hand, there were four principal components for HTS extracted with 92.3% of the total variance. PC1 was identified as "embracingness" dimension, with high loadings on fit and embraceable. PC2 was "refreshingness" dimension, with high loadings on reviving and refreshing. PC3 was "excitingness" dimension, with high loadings on exciting and exhilarating. Lastly, PC4 was "expensiveness" dimension, with high loadings on cheap (negative sign) and luxury. Here, luxury was seen to have slightly loading on PC1, in other words, luxury had a combination of not cheap and fit.

Data Collection of Physical Measures
As mentioned in Section 2, four physical measures (vibration, bulk displacement for surface deformation, thermal property, and friction) for each sample will be acquired in this section by a proposed tactile sensor for vibration and commercialized tactile sensors for others.

Vibration
During interaction between skin and object, vibration is one of the physical effects that are evoked. There are four kinds of mechanoreceptors in human glabrous skin that perceived vibration or mechanical stimuli: fast adapting, FA I (Meissner corpuscle), slow adapting, SA I (Merkel's disc), FA II (Pacinian corpuscle) and SA II (Ruffini ending) [41]. Furthermore, each mechanoreceptor has its peculiar frequency band of vibrating stimuli, and the perceptible frequency range of human is up to 1000 Hz [42].
This paper referred to previous research on the method of collecting and indexing vibrational data. Asaga et al. had collected vibrational data by tracing on the surface of samples with a piezoelectric element. Then, the vibrational data was compared with mechanoreceptors properties. As a result, two vibratory stimuli values, I FA I and I FA II , which correspond to the firing status of FA I and FA II, were determined and used to quantify roughness [11]. Here, a 15 mm × 22 mm × 3 mm acrylic resin plate with a piezoelectric element attached to was fabricated as shown in Figure 3. Piezoelectric element is mostly used for actuating or sensing vibration in numerous research as it has simple mechanism so that it is easy to implement in any design [7,[43][44][45][46]. The sensor was placed 45 • to the door armrest sample and traced on with a velocity of 50 mm/s under a load of 0.49 N. Then, two values of vibratory stimuli, I FA I and I FA II [V 2 ·Hz] which corresponded to mechanoreceptor FA I and FA II were estimated and will be used in correlating with LTS in the next section. placed 45° to the door armrest sample and traced on with a velocity of 50 mm/s under a load of 0.49 N. Then, two values of vibratory stimuli, IFA I and IFA II [V 2 ·Hz] which corresponded to mechanoreceptor FA I and FA II were estimated and will be used in correlating with LTS in the next section.

Bulk Displacement
Compliance or softness is perceived when skin is pressed on an object, and both the finger and the object deform as well as change their profile/pressure distributions [7]. However, the deformation of skin will not be measured because previous research [31] proved that 90 percent of the information in perceiving compliance is associated with the perception of surface deformation. Consequently, this paper proposes to measure bulk displacement when a fixed force is applied. By using an indentation hardness tester (TK-HS100, Tokushu-Keisoku. Co., Ltd, Yokohama, Japan) as shown in Figure 4, a sample was pressed with loads from 5 N to 30 N, with an interval of 5 N and the corresponding bulk displacements, d [mm] were measured.

Thermal Property
Coldness or warmness is perceived when heat is transferred from or to our skin when we touch them [47]. Perception of temperature is attributed to the thermal property between skin and an object [39]. By using Thermo Labo II B (FR-07, Kato Tech. Co., Ltd., Kyoto, Japan) which is a heat flux sensor, the silicone rubber surfaced sensor with dimension of 35  122 mm was preheated to 33 °C which was the average finger skin temperature, and then it was placed on the sample as shown in Figure 5. The peak heat transfer speed, qmax [-] was determined. Note that the sample was left in a room with temperature of 23 °C.

Bulk Displacement
Compliance or softness is perceived when skin is pressed on an object, and both the finger and the object deform as well as change their profile/pressure distributions [7]. However, the deformation of skin will not be measured because previous research [31] proved that 90 percent of the information in perceiving compliance is associated with the perception of surface deformation. Consequently, this paper proposes to measure bulk displacement when a fixed force is applied. By using an indentation hardness tester (TK-HS100, Tokushu-Keisoku. Co., Ltd, Yokohama, Japan) as shown in Figure 4, a sample was pressed with loads from 5 N to 30 N, with an interval of 5 N and the corresponding bulk displacements, d [mm] were measured. placed 45° to the door armrest sample and traced on with a velocity of 50 mm/s under a load of 0.49 N. Then, two values of vibratory stimuli, IFA I and IFA II [V 2 ·Hz] which corresponded to mechanoreceptor FA I and FA II were estimated and will be used in correlating with LTS in the next section.

Bulk Displacement
Compliance or softness is perceived when skin is pressed on an object, and both the finger and the object deform as well as change their profile/pressure distributions [7]. However, the deformation of skin will not be measured because previous research [31] proved that 90 percent of the information in perceiving compliance is associated with the perception of surface deformation. Consequently, this paper proposes to measure bulk displacement when a fixed force is applied. By using an indentation hardness tester (TK-HS100, Tokushu-Keisoku. Co., Ltd, Yokohama, Japan) as shown in Figure 4, a sample was pressed with loads from 5 N to 30 N, with an interval of 5 N and the corresponding bulk displacements, d [mm] were measured.

Thermal Property
Coldness or warmness is perceived when heat is transferred from or to our skin when we touch them [47]. Perception of temperature is attributed to the thermal property between skin and an object [39]. By using Thermo Labo II B (FR-07, Kato Tech. Co., Ltd., Kyoto, Japan) which is a heat flux sensor, the silicone rubber surfaced sensor with dimension of 35  122 mm was preheated to 33 °C which was the average finger skin temperature, and then it was placed on the sample as shown in Figure 5. The peak heat transfer speed, qmax [-] was determined. Note that the sample was left in a room with temperature of 23 °C.

Thermal Property
Coldness or warmness is perceived when heat is transferred from or to our skin when we touch them [47]. Perception of temperature is attributed to the thermal property between skin and an object [39]. By using Thermo Labo II B (FR-07, Kato Tech. Co., Ltd., Kyoto, Japan) which is a heat flux sensor, the silicone rubber surfaced sensor with dimension of φ35 × 122 mm was preheated to 33 • C which was the average finger skin temperature, and then it was placed on the sample as shown in Figure 5. The peak heat transfer speed, q max [-] was determined. Note that the sample was left in a room with temperature of 23 • C.

Friction Force
Slipperiness or stickiness is a perception when skin slides over on an object's surface, and the skin stretches and adheres to the surface [39]. Furthermore, according to previous research, this perception is mainly attributed to friction forces or friction coefficients [48][49][50]. Hence, by using Built-up Static-Dynamic Friction Measuring Device (TL201Ts, Trinity-Lab. Co., Ltd., Tokyo, Japan) with a skin-like urethane pad as shown in Figure 6, frictional force was measured and its variance, Afric was computed to represent the magnitude of fluctuation of the frictional force. The sensor was placed on a sample with a preload of 1.47 N and traced with a velocity of 5 mm/s.

Correlation between Subjective Responses and Physical Measures
As mentioned in Section 2, this paper suggests hierarchically layering the subjective responses and correlating between physical measures-LTS, LTS-HTS, and HTS-Preference. For this purpose a commonly used technique, multiple linear regression analysis, is implemented [18,[51][52][53][54]. This method is very simple and easy to interpret since there are multiple correlations between upper and lower layers [21,22,55].
Before conducting multiple linear regression analysis between LTS (dependent variables) and the physical measures (independent variables), principal component analysis with varimax rotation was performed to group physical measures that had strong correlation, and ensured no multicollinearity between independent variables. The result is as shown in Table 5. There were four principal components extracted with 94.9% of the total variance; PC1 was associated with bulk displacements for all load conditions, PC2 with vibratory stimuli values of IFA I and IFA II, PC3 with peak heat transfer speed, qmax, and PC4 with variance of dynamic frictional force, Afric. This result supports the concept of four main aspects of haptic information which are vibration, bulk displacement for surface deformation, thermal property, and friction.

Friction Force
Slipperiness or stickiness is a perception when skin slides over on an object's surface, and the skin stretches and adheres to the surface [39]. Furthermore, according to previous research, this perception is mainly attributed to friction forces or friction coefficients [48][49][50]. Hence, by using Built-up Static-Dynamic Friction Measuring Device (TL201Ts, Trinity-Lab. Co., Ltd., Tokyo, Japan) with a skin-like urethane pad as shown in Figure 6, frictional force was measured and its variance, A fric was computed to represent the magnitude of fluctuation of the frictional force. The sensor was placed on a sample with a preload of 1.47 N and traced with a velocity of 5 mm/s.

Friction Force
Slipperiness or stickiness is a perception when skin slides over on an object's surface, and the skin stretches and adheres to the surface [39]. Furthermore, according to previous research, this perception is mainly attributed to friction forces or friction coefficients [48][49][50]. Hence, by using Built-up Static-Dynamic Friction Measuring Device (TL201Ts, Trinity-Lab. Co., Ltd., Tokyo, Japan) with a skin-like urethane pad as shown in Figure 6, frictional force was measured and its variance, Afric was computed to represent the magnitude of fluctuation of the frictional force. The sensor was placed on a sample with a preload of 1.47 N and traced with a velocity of 5 mm/s.

Correlation between Subjective Responses and Physical Measures
As mentioned in Section 2, this paper suggests hierarchically layering the subjective responses and correlating between physical measures-LTS, LTS-HTS, and HTS-Preference. For this purpose a commonly used technique, multiple linear regression analysis, is implemented [18,[51][52][53][54]. This method is very simple and easy to interpret since there are multiple correlations between upper and lower layers [21,22,55].
Before conducting multiple linear regression analysis between LTS (dependent variables) and the physical measures (independent variables), principal component analysis with varimax rotation was performed to group physical measures that had strong correlation, and ensured no multicollinearity between independent variables. The result is as shown in Table 5. There were four principal components extracted with 94.9% of the total variance; PC1 was associated with bulk displacements for all load conditions, PC2 with vibratory stimuli values of IFA I and IFA II, PC3 with peak heat transfer speed, qmax, and PC4 with variance of dynamic frictional force, Afric. This result supports the concept of four main aspects of haptic information which are vibration, bulk displacement for surface deformation, thermal property, and friction.

Correlation between Subjective Responses and Physical Measures
As mentioned in Section 2, this paper suggests hierarchically layering the subjective responses and correlating between physical measures-LTS, LTS-HTS, and HTS-Preference. For this purpose a commonly used technique, multiple linear regression analysis, is implemented [18,[51][52][53][54]. This method is very simple and easy to interpret since there are multiple correlations between upper and lower layers [21,22,55].
Before conducting multiple linear regression analysis between LTS (dependent variables) and the physical measures (independent variables), principal component analysis with varimax rotation was performed to group physical measures that had strong correlation, and ensured no multicollinearity between independent variables. The result is as shown in Table 5. There were four principal components extracted with 94.9% of the total variance; PC1 was associated with bulk displacements for all load conditions, PC2 with vibratory stimuli values of I FA I and I FA II , PC3 with peak heat transfer speed, q max , and PC4 with variance of dynamic frictional force, A fric . This result supports the concept of four main aspects of haptic information which are vibration, bulk displacement for surface deformation, thermal property, and friction. By using principal components scores, multiple linear regression with a stepwise method was conducted three times using statistical analysis software (SPSS Ver. 22, IBM, Armonk, NY, USA), and Figure 7 shows the result: (1) LTS (as dependent variables) and physical measures (as independent variables); (2) HTS (as dependent variables) and LTS (as independent variables); and (3) Preference (as dependent variables) and HTS (as independent variables). The grayscale gradient of the lines corresponds to the standardized coefficients. Moreover, the solid and dash lines indicate the positive and negative coefficients, respectively. R 2 is a multiple coefficient of determination that the amount of variance in the regression model, and it is a better indicator in representing strength of the relationship between independent variables and the dependent variable [56]. However, there was one component in HTS (the "excitingness" dimension) that could not successfully construct its multiple regression with components in LTS because the correlations between the "excitingness" dimension with components in LTS were not significantly strong. In order to build estimating equations for preference from physical measures, the "excitingness" dimension was not included during the construction of multiple regressions between Preference and HTS.  By using principal components scores, multiple linear regression with a stepwise method was conducted three times using statistical analysis software (SPSS Ver. 22, IBM, Armonk, NY, USA), and Figure 7 shows the result: (1) LTS (as dependent variables) and physical measures (as independent variables); (2) HTS (as dependent variables) and LTS (as independent variables); and (3) Preference (as dependent variables) and HTS (as independent variables). The grayscale gradient of the lines corresponds to the standardized coefficients. Moreover, the solid and dash lines indicate the positive and negative coefficients, respectively. R 2 is a multiple coefficient of determination that the amount of variance in the regression model, and it is a better indicator in representing strength of the relationship between independent variables and the dependent variable [56]. However, there was one component in HTS (the "excitingness" dimension) that could not successfully construct its multiple regression with components in LTS because the correlations between the "excitingness" dimension with components in LTS were not significantly strong. In order to build estimating equations for preference from physical measures, the "excitingness" dimension was not included during the construction of multiple regressions between Preference and HTS. Furthermore, estimating equations for preference (prefer and pleasure) from principal com ponents of physical measures were derived by substituting multiple regression equations obt ained from the analysis. Hence, the feedback of subjective responses on tactile aspect for do or armrest can be acquired from physical measures by using the equations as follows: Furthermore, estimating equations for preference (prefer and pleasure) from principal components of physical measures were derived by substituting multiple regression equations obtained from the analysis. Hence, the feedback of subjective responses on tactile aspect for door armrest can be acquired from physical measures by using the equations as follows:

Verification Experiment
In this section, the equations obtained were verified by using unknown samples of a, b, and c (refer Table 6). First, a sensory evaluation test for each sample with the same condition and participants in Section 3.1 was carried out; however, only adjectives in the preference layer were asked. Next, similar physical measures as in Section 4 for unknown samples were measured. Consequently, principal component scores for each principal component in the physical measures layer were computed by using principal component loadings. Then, the evaluation scores of P prefer and P pleasure were calculated by using estimating Equations (1) and (2), respectively. Lastly, the actual and estimated scores were compared in Figure 8. Gray plots are the 26 samples that are used in the process to derive the estimating equations, and black plots are the unknown samples. Dashed lines indicate one-to-one relationships. (2)

Verification Experiment
In this section, the equations obtained were verified by using unknown samples of a, b, and c (refer Table 6). First, a sensory evaluation test for each sample with the same condition and participants in Section 3.1 was carried out; however, only adjectives in the preference layer were asked. Next, similar physical measures as in Section 4 for unknown samples were measured. Consequently, principal component scores for each principal component in the physical measures layer were computed by using principal component loadings. Then, the evaluation scores of Pprefer and Ppleasure were calculated by using estimating Equations (1) and (2), respectively. Lastly, the actual and estimated scores were compared in Figure 8. Gray plots are the 26 samples that are used in the process to derive the estimating equations, and black plots are the unknown samples. Dashed lines indicate one-to-one relationships.  The accuracy of this tactile evaluation feedback system was determined by calculating the percent error for each sample using the following equation.
Percent error = |Estimated score − Actual score| Actual score ×100%, As a result, the maximum percent error when estimating prefer and pleasure were 20.8% and 16.2%, respectively. In other words, the developed system has the accuracy of 79% for prefer and 84% for pleasure in giving feedback on tactile evaluation.
Similarly, the percent errors when estimating prefer and pleasure for unknown samples a, b and c were computed, and the result was concluded in Table 7. Hence, by comparing the percentage error for prefer and pleasure, both have smaller percent errors compared to the developed system. Thus, this system can give estimation on unknown product's evaluation successfully.  (2)

Verification Experiment
In this section, the equations obtained were verified by using unknown samples of a, b, and c (refer Table 6). First, a sensory evaluation test for each sample with the same condition and participants in Section 3.1 was carried out; however, only adjectives in the preference layer were asked. Next, similar physical measures as in Section 4 for unknown samples were measured. Consequently, principal component scores for each principal component in the physical measures layer were computed by using principal component loadings. Then, the evaluation scores of Pprefer and Ppleasure were calculated by using estimating Equations (1) and (2), respectively. Lastly, the actual and estimated scores were compared in Figure 8. Gray plots are the 26 samples that are used in the process to derive the estimating equations, and black plots are the unknown samples. Dashed lines indicate one-to-one relationships.  The accuracy of this tactile evaluation feedback system was determined by calculating the percent error for each sample using the following equation.
Percent error = |Estimated score − Actual score| Actual score ×100%, As a result, the maximum percent error when estimating prefer and pleasure were 20.8% and 16.2%, respectively. In other words, the developed system has the accuracy of 79% for prefer and 84% for pleasure in giving feedback on tactile evaluation.
Similarly, the percent errors when estimating prefer and pleasure for unknown samples a, b and c were computed, and the result was concluded in Table 7. Hence, by comparing the percentage error for prefer and pleasure, both have smaller percent errors compared to the developed system. (2)

Verification Experiment
In this section, the equations obtained were verified by using unknown samples of a, b, and c (refer Table 6). First, a sensory evaluation test for each sample with the same condition and participants in Section 3.1 was carried out; however, only adjectives in the preference layer were asked. Next, similar physical measures as in Section 4 for unknown samples were measured. Consequently, principal component scores for each principal component in the physical measures layer were computed by using principal component loadings. Then, the evaluation scores of Pprefer and Ppleasure were calculated by using estimating Equations (1) and (2), respectively. Lastly, the actual and estimated scores were compared in Figure 8. Gray plots are the 26 samples that are used in the process to derive the estimating equations, and black plots are the unknown samples. Dashed lines indicate one-to-one relationships.  The accuracy of this tactile evaluation feedback system was determined by calculating the percent error for each sample using the following equation.
Percent error = |Estimated score − Actual score| Actual score ×100%, As a result, the maximum percent error when estimating prefer and pleasure were 20.8% and 16.2%, respectively. In other words, the developed system has the accuracy of 79% for prefer and 84% for pleasure in giving feedback on tactile evaluation.
Similarly, the percent errors when estimating prefer and pleasure for unknown samples a, b and c were computed, and the result was concluded in Table 7. Hence, by comparing the percentage error for prefer and pleasure, both have smaller percent errors compared to the developed system. (2)

Verification Experiment
In this section, the equations obtained were verified by using unknown samples of a, b, and c (refer Table 6). First, a sensory evaluation test for each sample with the same condition and participants in Section 3.1 was carried out; however, only adjectives in the preference layer were asked. Next, similar physical measures as in Section 4 for unknown samples were measured. Consequently, principal component scores for each principal component in the physical measures layer were computed by using principal component loadings. Then, the evaluation scores of Pprefer and Ppleasure were calculated by using estimating Equations (1) and (2), respectively. Lastly, the actual and estimated scores were compared in Figure 8. Gray plots are the 26 samples that are used in the process to derive the estimating equations, and black plots are the unknown samples. Dashed lines indicate one-to-one relationships.  The accuracy of this tactile evaluation feedback system was determined by calculating the percent error for each sample using the following equation.
Percent error = |Estimated score − Actual score| Actual score ×100%, As a result, the maximum percent error when estimating prefer and pleasure were 20.8% and 16.2%, respectively. In other words, the developed system has the accuracy of 79% for prefer and 84% for pleasure in giving feedback on tactile evaluation.
Similarly, the percent errors when estimating prefer and pleasure for unknown samples a, b and c were computed, and the result was concluded in Table 7. Hence, by comparing the percentage error for prefer and pleasure, both have smaller percent errors compared to the developed system. The accuracy of this tactile evaluation feedback system was determined by calculating the percent error for each sample using the following equation.
Percent error = |Estimated score − Actual score| Actual score × 100%, As a result, the maximum percent error when estimating prefer and pleasure were 20.8% and 16.2%, respectively. In other words, the developed system has the accuracy of 79% for prefer and 84% for pleasure in giving feedback on tactile evaluation.
Similarly, the percent errors when estimating prefer and pleasure for unknown samples a, b and c were computed, and the result was concluded in Table 7. Hence, by comparing the percentage error for prefer and pleasure, both have smaller percent errors compared to the developed system. Thus, this system can give estimation on unknown product's evaluation successfully.

Discussion and Conclusions
The above physical measuring and statistical analysis results have shown that the proposed concepts of the developed assessment system can be considered adequate with slight errors. First of all, the concept of hierarchy stages of subjective responses has helped us to easily interpret the main aspect of tactile sensation that is related to the preference. From Figure 7, the "embracingness" dimension seems to have high correlation to both prefer and pleasure, compared to other HTS components. Furthermore, the "hardness" dimension shows a strong correlation to the "embracingness" dimension. Hence, the biggest influence on the preference layer can be concluded to be the "hardness" dimension in the case of tactile assessment of door armrests. However, the structure is provisional based on the kind of the object. In another study on tactile assessment of film and board materials for confectionery packaging [21], "roughness" seemed to be the most important factor to the affective layer (equivalent to HTS in this paper).
Besides, this paper suggests using the final-end product in the tactile assessment of a product that has layers of different materials. This argument is supported by the result obtained that shows the "hardness" dimension is the most important aspect in tactile assessment of door armrests. The perception of "hardness" involves both the kinesthetic and cutaneous systems [29]. This perception may not be evaluated accurately by just using only the outer layer of the sample, as the product had a layered structure.
Furthermore, this paper suggests correlating LTS with four main aspects of haptic information (roughness, compliance, coldness, and slipperiness). The physical measures for each aspects of haptic information are selected based on the kind of stimuli that evoked the receptors, and also the physical effects that occurred during the interaction between the skin and the object. Other than the physical measures that are selected in this paper (i.e. vibration for roughness, bulk displacement for surface deformation for compliance, thermal property for coldness, and friction for slipperiness), there are many other possible physical measures that are used in other studies. For example, Chen et al. chose three dimensional pictures or topography of a surface's texture for roughness [21]. However, the result is not quite convincing. In addition, the author mentioned that there is a need for further study on the measurement of roughness and it is probably related to vibration, as mentioned in other research [30,57]. Thus, the concept proposed in this paper may help in selecting suitable physical measures.
More work is required to find the appropriate physical measures to correlate with LTS, because in this paper, several coefficient of determinations obtained from the multiple regression analysis between LTS (dependent variable) and physical measures (independent variable) are less than 0.5 (arbitrary lower limit for strong correlation), especially in the case of "micro-roughness", "macro-roughness" and "hollowness". Moreover, physical measure of vibration, which was expected to have correlation with roughness, was found to be not significantly correlated.
In addition, the other work is to classify people by clustering them according to their preference and then construct each group's hierarchy structure of tactile sensation. This may help product developers in targeting their market. Before that, there is a need to increase the number of participants and vary the cohorts of people, for example, broaden the age groups.