Food evaluation value prediction method, food evaluation value prediction model generation method, and program

Abstract

This technology improves the accuracy of predicting food quality values. [Solution] The prediction method includes the steps of: acquiring a plurality of analytical parameters that indicate the analysis results of the food to be predicted obtained by a taste sensor; and inputting the acquired plurality of analytical parameters to a prediction model to obtain a food evaluation value that quantifies one of the senses of taste, texture, and smell of the food. The analytical parameters obtained by the taste sensor are obtained based on the sensor intensity at multiple time points measured between the start of measurement and the end of measurement, and include at least one of the following: intensity before the end of measurement, maximum intensity, intensity at a predetermined ratio to the maximum intensity, time to reach the maximum intensity, time to reach the intensity at a predetermined ratio to the maximum intensity, duration for which the intensity is equal to or greater than a predetermined ratio to the maximum intensity, area determined by the intensity over a predetermined period, and the slope of the line connecting the intensity at the first time point and the intensity at the second time point.

Description

[Technical Field]

[0001] This disclosure relates to a technology for predicting food evaluation values ​​using machine learning. [Background technology]

[0002] Food evaluation values, which quantify the senses of taste, texture, and smell of food, are defined by sensory evaluation based on human sensory characteristics. Sensory evaluation has the advantage of directly quantifying human senses, but in order to achieve highly reproducible evaluations, it is necessary to train a large number of trained evaluators (panelists), which is costly and time-consuming.

[0003] In contrast, Patent Document 1 mentions that research and development are being conducted to evaluate texture by physical measurement using instruments, and discloses that it is considered in principle difficult to evaluate human texture by measuring the uniaxial reaction force when food is pressed. [Prior art documents] [Patent Documents]

[0004] [Patent Document 1] Japanese Patent Publication No. 2019-207123 [Overview of the Initiative] [Problems that the invention aims to solve]

[0005] This disclosure provides a technology to improve the prediction accuracy of food evaluation values. [Means for solving the problem]

[0006] The food evaluation value prediction method of this disclosure includes the steps of: obtaining a plurality of analytical parameters that indicate the analysis results of a food to be predicted by an analytical instrument including a taste sensor; and obtaining the food evaluation value by inputting the obtained plurality of analytical parameters into a prediction model that takes the plurality of analytical parameters indicating the analysis results of a food by an analytical instrument including a taste sensor as input and outputs a food evaluation value that quantifies one of the senses of taste, texture and smell of the food. The analytical parameters obtained by the taste sensor include at least one of the following, obtained based on the sensor intensity at multiple time points measured between the start of measurement and the end of measurement: intensity before the end of measurement, maximum intensity, intensity at a predetermined ratio to the maximum intensity, time to reach the maximum intensity, time to reach an intensity at a predetermined ratio to the maximum intensity, maintenance time when the intensity is equal to or greater than a predetermined ratio to the maximum intensity, and area determined by the intensity over a predetermined period. [Brief explanation of the drawing]

[0007] [Figure 1] This is a block diagram showing a food evaluation value prediction model generation system and a food evaluation value prediction system. [Figure 2] This figure shows the parameters (food evaluation values) of the average TI curve plotted using the Overbosch method. [Figure 3] This is a schematic diagram illustrating the measurement state of a taste sensor and the response value (mV) obtained from the sensor. [Figure 4] This is an explanatory diagram relating to multiple analytical parameters obtained by the taste sensor of the first embodiment. [Figure 5] This is a flowchart showing the method for generating a predictive model. [Figure 6] This figure shows the food evaluation values ​​of the first embodiment: Comparative Example 1 for RMinc of the TI method, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation. [Figure 7] This figure shows the food evaluation values ​​of the first embodiment: Comparative Example 1 for Tend of the TI method, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation. [Figure 8] This figure shows the food evaluation values ​​of the first embodiment: Comparative Example 1 for RMdec of the TI method, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation. [Figure 9] This figure shows the food evaluation values ​​of the first embodiment: Comparative Example 1 for Rinc of the TI method, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation. [Figure 10] This figure shows the food evaluation values ​​of the first embodiment: Comparative Example 1 for Rdec of the TI method, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation. [Figure 11] This figure shows the food evaluation values ​​of the first embodiment: Comparative Example 1 for DUR of the TI method, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation. [Figure 12] This figure shows the food evaluation values ​​of the first embodiment: Comparative Example 1 for Dpla of the TI method, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation. [Figure 13] This figure shows the food evaluation values ​​of the first embodiment: Comparative Example 1 for Dinc of the TI method, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation. [Figure 14] This figure shows the food evaluation values ​​of the first embodiment: Comparative Example 1 for Ddec of the TI method, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation. [Figure 15] This figure shows the food evaluation values ​​of the first embodiment: Comparative Example 1 for AUC of the TI method, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation. [Figure 16] This figure shows the food evaluation values ​​of the first embodiment: Comparative Example 1 for Apla of the TI method, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation. [Figure 17] This figure shows the food evaluation values ​​of the first embodiment: Comparative Example 1 for Ainc of the TI method, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation. [Figure 18]Food evaluation value of the first embodiment: A figure showing the accuracy of the regression equations and the analysis parameters used as inputs for the regression equations in Comparative Example 1 and Example 1 for Adec of the TI method. [Figure 19] Food evaluation value of the first embodiment: A figure showing the accuracy of the regression equations and the analysis parameters used as inputs for the regression equations in Comparative Example 1 and Example 1 for astringency of the QDA method. [Figure 20] Food evaluation value of the first embodiment: A figure showing the accuracy of the regression equations and the analysis parameters used as inputs for the regression equations in Comparative Example 1 and Example 1 for citrus of the QDA method. [Figure 21] Food evaluation value of the first embodiment: A figure showing the accuracy of the regression equations and the analysis parameters used as inputs for the regression equations in Comparative Example 1 and Example 1 for herbal medicine of the QDA method. [Figure 22] Food evaluation value of the first embodiment: A figure showing the accuracy of the regression equations and the analysis parameters used as inputs for the regression equations in Comparative Example 1 and Example 1 for kire of the QDA method. [Figure 23] Food evaluation value of the first embodiment: A figure showing the accuracy of the regression equations and the analysis parameters used as inputs for the regression equations in Comparative Example 1 and Example 1 for lastingness of the QDA method. [Figure 24] Food evaluation value of the first embodiment: A figure showing the accuracy of the regression equations and the analysis parameters used as inputs for the regression equations in Comparative Example 1 and Example 1 for rapidness of the QDA method. [Figure 25] Food evaluation value of the first embodiment: A figure showing the accuracy of the regression equations and the analysis parameters used as inputs for the regression equations in Comparative Example 1 and Example 1 for spread of the QDA method. [Figure 26] Food evaluation value of the first embodiment: A figure showing the accuracy of the regression equations and the analysis parameters used as inputs for the regression equations in Comparative Example 1 and Example 1 for thickness of the QDA method. [Figure 27]This figure shows the food evaluation value of the first embodiment: Comparative Example 1 for woody using the QDA method, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation. [Figure 28] This figure shows the food evaluation values ​​of the second embodiment: Comparative Example 1, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation for Rinc of the TI method. [Figure 29] This figure shows the food evaluation values ​​of the second embodiment: Comparative Example 1, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation for RMinc of the TI method. [Figure 30] This figure shows the food evaluation values ​​of the second embodiment: Comparative Example 1, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation for the Rdec of the TI method. [Figure 31] This figure shows the food evaluation values ​​of the second embodiment: Comparative Example 1, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation for RMdec of the TI method. [Figure 32] This figure shows the food evaluation values ​​of the second embodiment: Comparative Example 1, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation for Tstart of the TI method. [Figure 33] This figure shows the food evaluation values ​​of the second embodiment: Comparative Example 1 for Dinc of the TI method, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation. [Figure 34] This figure shows the food evaluation values ​​of the second embodiment: Comparative Example 1 for Tmax of the TI method, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation. [Figure 35] This figure shows the food evaluation values ​​of the second embodiment: Comparative Example 1, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation for Dpla of the TI method. [Figure 36] This figure shows the food evaluation values ​​of the second embodiment: Comparative Example 1 for Ddec of the TI method, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation. [Figure 37]This figure shows the food evaluation values ​​of the second embodiment: Comparative Example 1 for Tend of the TI method, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation. [Figure 38] This figure shows the food evaluation values ​​of the second embodiment: Comparative Example 1, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation for the DUR of the TI method. [Figure 39] This figure shows the food evaluation values ​​of the second embodiment: Comparative Example 1 for Ainc of the TI method, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation. [Figure 40] This figure shows the food evaluation values ​​of the second embodiment: Comparative Example 1, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation for Apla of the TI method. [Figure 41] This figure shows the food evaluation values ​​of the second embodiment: Comparative Example 1, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation for Adec of the TI method. [Figure 42] This figure shows the food evaluation values ​​of the second embodiment: Comparative Example 1 for AUC of the TI method, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation. [Figure 43] This figure shows the food evaluation value of the second embodiment: the initial intensity of the QDA method, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation. [Figure 44] This figure shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1 and Example 1 regarding the food evaluation value of the second embodiment: lastingness of the QDA method. [Figure 45] This figure shows the food evaluation values ​​of the second embodiment: Comparative Example 1, the accuracy of the regression equation for Example 1, and the analytical parameters used as input to the regression equation for bitterness using the QDA method. [Figure 46] This figure shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1 and Example 1 regarding the astringency of the food evaluation value of the second embodiment using the QDA method. [Figure 47]This figure shows the food evaluation values ​​of the second embodiment: Comparative Example 1, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation for the monotonic (Flat) QDA method. [Figure 48] This figure shows the food evaluation values ​​of the second embodiment: Comparative Example 1, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation for stickiness (coating) using the QDA method. [Figure 49] This figure shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1 and Example 1 regarding the food evaluation value of the second embodiment: aroma (caramel) using the QDA method. [Figure 50] This diagram illustrates the proportion of each parameter selected as an important variable. [Modes for carrying out the invention]

[0008] Embodiments of this disclosure will be described below with reference to the drawings. Figure 1 is a block diagram showing the food evaluation value prediction model generation system 1 and the food evaluation value prediction system 2.

[0009] The food evaluation value prediction model generation method described herein can be executed by the food evaluation value prediction model generation system 1. Furthermore, the food evaluation value prediction method described herein can be executed by the food evaluation value prediction system 2. Both the food evaluation value prediction model generation system 1 and the food evaluation value prediction system 2 can be implemented using a computer.

[0010] As shown in Figure 1, the food evaluation value prediction model generation system 1 constructs a prediction model 20 using machine learning with training data D1. The food evaluation value prediction system 2 uses the prediction model 20 to predict a food evaluation value (estimated value) that quantifies one of the senses of taste, texture, and smell of the food, based on multiple analytical parameters that show the analysis results of the food to be predicted by analytical instruments 3 such as taste sensors. The food evaluation value prediction system 2 uses the prediction model 20 that was previously generated by machine learning by the food evaluation value prediction model generation system 1. Details of the food evaluation value prediction model generation system 1 and the food evaluation value prediction system 2 will be described later.

[0011] <First Embodiment> [Food evaluation values ​​(sensory evaluation values, bitterness evaluation values)] The food evaluation values ​​in this disclosure are numerical values ​​that quantify one of the senses of taste, texture, and smell of food. In the first embodiment, bitterness among the senses of taste is used as an example. Nine types of bittering agents and bittering materials were used as food samples to evaluate only bitterness. Specifically, these were Caffeine, Jamaica quassia, Naringin, Wormwood, Kuding tea, Iso alpha acids, Gentiana lutea (root), Quina, and Amur cork tree. A selection panel of 50 people compared each sample with a control (0.04% Caffeine aqueous solution) and serially diluted samples, and selected the one they felt was more bitter (two-point discrimination method). The inflection point of the sigmoid curve plotted based on the selectivity of each concentration was defined as the subjective equivalent concentration of bitterness, and in subsequent evaluations and analyses, samples were prepared and provided at that concentration. Each sample was evaluated using Quantitative Descriptive Analysis (QDA) and Time-Intensity Analysis (TI), and the following 25 food evaluation values ​​(bitterness evaluation values) were obtained.

[0012] [Food evaluation values ​​showing the qualitative characteristics of bitterness based on Quantitative Descriptive Analysis (QDA method)] A panel of 14 selected members determined 10 characteristic parameters and their definitions suitable for the sample. Through trial evaluations, appropriate evaluation frequencies and sample quantities were determined. After sufficient evaluation training, the intensity of each characteristic parameter was evaluated on a 15cm line scale. The evaluation was repeated twice. Each sample was assigned a three-digit random number, and the liquid color was blinded during the test. The following shows the 10 food evaluation values ​​(characteristic items) and their definitions. For example, the definition of the characteristic item "astringency" is "irritation to the oral cavity." • Astringency: Irritation in the oral cavity (the higher the intensity, the stronger the irritation) • Rapidity of expression: The speed of the timing (higher intensity means faster). • Sharp: The gradient of the rising angle (the higher the strength, the greater the gradient) • Kire: The gradient of the angle of the drop (the greater the strength, the steeper the gradient) • Lastingness: The length of the lingering effect (longer for stronger effects) • Thickness: The flavor profile from the middle to the last stage (the higher the thickness, the stronger the flavor profile). • Spread: Spreads throughout the entire mouth (the stronger the intensity, the wider the spread). • Herbal medicine: Herbal, mossy bitterness (the stronger the concentration, the more herbal or mossy the taste). • Woody: A woody, pencil lead-like bitterness (the stronger the intensity, the more woody or pencil lead-like the taste). • Citrus: Grapefruit-like bitterness (the higher the intensity, the stronger the grapefruit-like taste)

[0013] [Food evaluation values ​​showing the time-dependent characteristics of bitterness using the Time-Intensity (TI) method] A selection panel of 13 participants conducted trial evaluations to determine appropriate sample quantities, swallowing times, etc. After sufficient evaluation training, the temporal intensity of "bitterness" was evaluated using a slider bar on a 10cm line scale. The evaluation was repeated three times. Figure 2 shows the parameters (food evaluation values) of the average TI curve plotted using the Overbosch method. In Figure 2, the vertical axis represents intensity, and the horizontal axis represents time. As shown in Figure 2, in the first embodiment, there are 15 food evaluation values ​​that indicate the temporal characteristics of bitterness using the TI method. The 15 food evaluation values ​​and their definitions are as follows. • Rinc: Slope of increase • RMinc: Maximum slope of increase • Rdec: Slope of decrease • RMdec: Maximum slope of decrease • Tstart: The first time bitterness was detected. • Dinc: Increased time • Tmax: Time at which the intensity first reaches its maximum. • Dpla: Duration of the plateau of maximum intensity ·Ddec: Decrease time • Tend: The time when the bitterness is no longer noticeable. • DUR: Total time spent feeling bitterness Ainc: Area under the increasing curve • Apla: Area under the curve of the maximum intensity plateau • Adec: Area under the decreasing curve • AUC: Area under the curve of the entire curve

[0014] [Example of analysis parameter: Taste sensor] An example of analytical instrument 3 is a taste sensor. An example of analytical parameters obtained by a taste sensor is described below. A taste sensor is a sensor that obtains response information to a lipid membrane. Each sample was analyzed using eight types of taste sensors. The specific items of the eight types of taste sensors and their identifiers are as follows. ·Acidic bitterness:C00 • Hydrochloride bitter taste: BT0 • Basic salt bitterness: AN0 • Astringency: AE1 • Acidity: CA0 • Saltiness: CT0 • Umami: AAE • Sweetness: GL1 For each taste sensor, a first-stage initial taste measurement, a second-stage first aftertaste measurement, and a third-stage second aftertaste measurement were performed, and the response value over a 30-second measurement period was measured in 1-second increments at each stage. Figure 3 is a schematic explanatory diagram showing the measurement state of the taste sensor and the response value (mV) obtained from the sensor. The measurement method is as schematically shown in Figure 3 for membrane potential (sensor output), first by immersing the taste sensor in a reference solution to obtain a baseline membrane potential (sensor output). Next, when the taste sensor is immersed in the sample, the membrane potential changes due to interaction with the taste substance, and the amount or portion of the change corresponds to the initial taste (sometimes referred to as test_1). Next, the taste sensor is washed with the reference solution. Next, the taste sensor is immersed in the reference solution and the change in membrane potential is measured. The amount or portion of the change corresponds to the first aftertaste (sometimes referred to as CPA1). In the first embodiment, the taste sensor is further washed with the reference solution, and then the taste sensor is immersed in the reference solution and the change in membrane potential is measured. The amount or part of the change corresponds to the second aftertaste (sometimes referred to as CPA2).

[0015] Figure 4 is an explanatory diagram of the multiple analytical parameters obtained by the taste sensor of the first embodiment. As shown in Figure 4, the analytical parameters are obtained based on the sensor intensity at multiple time points (0 seconds, 1 second, 2 seconds, ..., 29 seconds, 30 seconds) measured between the start and end of the measurement (30 seconds). The 10 basic analytical parameters obtained by the taste sensor and their definitions are as follows. ·Imax: Maximum intensity • Tmax: Time to reach maximum intensity • Plateau Intensity: 90% Imax (Maximum intensity plateau) • Plateau_Reach_Time: Time to reach 90% Imax (maximum intensity plateau) • Plateau_Duration: Duration of the plateau of maximum intensity • Intensity_at_0_Sec: Intensity at 0 seconds (start of measurement) • Intensity_at_30_Sec: Intensity at 30 seconds (end of measurement) • AUC_total: Area under the curve of the entire curve ·AUC_pre_plateau: Area under the curve before Plateau_Reach_Time AUC_plateau: Area under the curve of Plateau_Duration The above 10 parameters are basic parameters. Since the initial taste, the first aftertaste, and the second aftertaste each have the above 10 reference parameters, there are a total of 30 parameters. For example, the maximum intensity of the initial taste is denoted as Imax_test_1, the maximum intensity of the first aftertaste is denoted as Imax_CPA_1, and the maximum intensity of the second aftertaste is denoted as Imax_CPA_2. Furthermore, for each basic parameter, combination parameters can be generated that represent the difference (def) and ratio (ratio) of two combinations of the initial taste, first aftertaste, and second aftertaste. For example, the ratio of the initial taste to the second aftertaste in the AUC_total of the taste sensor (C00) can be expressed as TasteSensor_SB2C00_AUC_total_ratio_test_1_vs_CPA_2. The difference between the initial taste and the second aftertaste in the AUC_total of the taste sensor (C00) can be expressed as TasteSensor_SB2C00_AUC_total_diff_test_1_vs_CPA_2. Based on one basic parameter, there are three combination parameters: the initial taste and first aftertaste, the initial taste and second aftertaste, and the first aftertaste and second aftertaste. Each combination has two types: a difference and a ratio. Therefore, based on one basic parameter, 3 × 2 = 6 combination parameters can be generated.

[0016] In summary, the analytical parameters obtained by the taste sensor may include intensities prior to the end of the measurement (30 seconds). An example of this is Intensity_at_0_Sec [intensity at 0 seconds (start of measurement)]. In some cases, Plateau_Intensity [90% Imax (plateau of maximum intensity)] can also be used as an example. The analytical parameters obtained by the taste sensor may include an intensity that is a predetermined percentage of the maximum intensity. In this example, if the predetermined percentage is set to 90%, Plateau_Intensity[90%Imax (plateau of maximum intensity)] can be exemplified. In the first embodiment, the predetermined percentage is 90%, but any value greater than 0% and less than 100%, such as 80%, can be appropriately set. The analytical parameters obtained from the taste sensor can exemplify the time it takes to reach a predetermined percentage of the maximum intensity. An example of this is Plateau_Reach_Time [time to reach 90%Imax (plateau of maximum intensity)]. The analytical parameters obtained by the taste sensor may include the duration at which the intensity is equal to or greater than a predetermined percentage of the maximum intensity. An example of this is Plateau_Duration [Duration of the maximum intensity plateau]. The analytical parameters obtained by the taste sensor may include an area determined by the intensity over a predetermined period. Examples include at least one of AUC_total [area under the curve of the entire curve], AUC_pre_plateau [area under the curve before Plateau_Reach_Time], and AUC_plateau [area under the curve of Plateau_Duration]. Conventionally, a single taste sensor only recorded the intensity at the end of the measurement (30 seconds), resulting in only the first and first aftertaste measurements. Specifically, only Intensity_at_30_Sec_test_1 and Intensity_at_30_Sec_CPA_1 were known. However, by continuously measuring the intensity at multiple time points during the measurement period, as described above, a large number of analytical parameters could be obtained. This allows for improved prediction accuracy in machine learning. Examples and comparative examples demonstrating the effect will be described later.

[0017] <Food Evaluation Value Prediction Model Generation System 1, Food Evaluation Value Prediction System 2> As shown in Figure 1, the food evaluation value prediction system 2 includes a prediction target analysis parameter acquisition unit 21 and a prediction unit 22. The food evaluation value prediction model generation system 1 includes a training data acquisition unit 10, training data D1, and a learning unit 11. The training data D1 is stored in memory 1a. The training data acquisition unit 10 and the learning unit 11 are implemented by processor 1b. The prediction target analysis parameter acquisition unit 21 and the prediction unit 22 are implemented by processor 2b. In the first embodiment, processors 1b and 2b in one device implement each unit, but the system is not limited to this. For example, each process may be distributed using a network, and multiple processors 1b and 2b may be configured to execute the processing of each unit. That is, one or more processors 1b and 2b execute the processing.

[0018] [Predictive model generation] The training data acquisition unit 10 takes multiple analytical parameters representing the analysis results of food by the analytical instrument 3 as input (explanatory variables) and acquires training data D1 with food evaluation values ​​as output (dependent variable). In the first embodiment, the multiple branching parameters in the training data D1 include analytical parameters obtained by the taste sensor (indicated as TS in Figure 1). The number of analytical parameters obtained by the taste sensor includes 8 types of taste sensors × (10 basic parameters + 6 × 10 combination parameters) = approximately 560 parameters. In one model, there is one food evaluation value per food as the output (dependent variable) for the input (group of explanatory variables). In Figure 1, if the food evaluation value X1 is, for example, Tstart, then the food evaluation value (Tstart) of food 1 is X 1_1 This is written as follows, and the food evaluation value (Tstart) of food 2 is X 1_2 Similarly, the food evaluation value (Tstart) of food 9 is X 1_9This is how it is written. In other words, in the first embodiment, in order to generate a predictive model 20 that predicts food evaluation values ​​(Tstart), there are nine rows of data, each row taking the above-mentioned multiple analysis parameters as input and outputting one food evaluation value. Although not shown in Figure 1, there are 25 types of food evaluation values, so training data D1 is obtained similarly for each food evaluation value, and one predictive model 20 is generated for each type of food evaluation value. There are 25 types of predictive models 20 generated.

[0019] The learning unit 11 shown in Figure 1 generates a predictive model 20 that outputs food evaluation values ​​by performing machine learning using the training data D1. In the first embodiment, the predictive model 20 uses a regression equation based on PLS (partial least squares) regression analysis, but it is not limited to this as long as it is a supervised machine learning model. For example, various models such as regression trees, random forests, support vector machines, neural networks, and ensembles can be used. In the first embodiment, since there are only nine samples for a large number of explanatory variables, PLS regression analysis, which has a high predictive value even with a small number of samples for a large number of explanatory variables, is used. However, if a large number of samples can be prepared, it is not limited to PLS regression analysis, and other regression methods can also be used.

[0020] Furthermore, in the first embodiment, PLS regression analysis, which is robust against multicollinearity, is used. However, in order to reduce multicollinearity as much as possible, a first removal processing unit 12 that performs a first removal process and a second removal processing unit 13 that performs a second removal process are provided. The first removal process is a process of determining which analytical parameters to be removed are selected from among a plurality of parameters obtained by the analytical instrument 3, and whose correlation coefficient with the food evaluation value is lower than a predetermined threshold. The second removal process is a process of determining which of the analytical parameters obtained by the analytical instrument 3, which have a relatively high correlation coefficient with the food evaluation value, will be excluded. The generation flow of the prediction model 20, including the first removal process, the second removal process, and PLS regression analysis, will be explained with reference to Figure 5.

[0021] Figure 5 is a flowchart illustrating the method for generating a predictive model. As shown in Figure 5, the first elimination process is initiated. Specifically, in step ST1, the correlation coefficient between each of the multiple analytical parameters (approximately 500) obtained by multiple analytical instruments 3 (multiple taste sensors) and the food evaluation value is calculated. In the next step ST2, analytical parameters whose absolute value of the correlation coefficient with the food evaluation value is lower than a predetermined threshold (0.5) are determined to be excluded and deleted. This completes the first elimination process. In one example, approximately 50 explanatory variables remained at the end of the first elimination process.

[0022] Next, the second elimination process is initiated. Specifically, in step ST3, all possible combinations (pairs) of multiple analytical parameters are generated. In the next step, ST4, the pairs of analytical parameters generated in step ST3 are sorted in descending order of their correlation coefficient with the food evaluation value. In the next step, ST5, one of the first or second analytical parameters constituting a pair is selected as the analytical parameter to be excluded and deleted. More specifically, the determination of the analytical parameter to be deleted is performed randomly. If there are many analytical parameters, the analytical parameter to be deleted may be determined using the RV coefficient with the food evaluation value obtained by MFA analysis. The process in step ST5 is repeated until the number of remaining analytical parameters (explanatory variables) is 30 or less (step ST6: NO). The second elimination process makes it possible to reduce the number of explanatory variables to 30 or less.

[0023] In step ST6, if the number of remaining analysis parameters (explanatory variables) is 30 or less (step ST6: YES), then in the next step ST7, PLS regression analysis is performed. A detailed explanation of PLS ​​regression analysis is omitted. Briefly, it involves repeatedly performing PLS mapping based on the remaining explanatory variables to select the optimal factors and generate a regression equation, and then generating a regression equation using only explanatory variables with high variable weights (VIP) (e.g., thresholds of 0.9 or higher). The termination conditions for the iteration include either the number of explanatory variables being less than or equal to a predetermined value (5 in the first embodiment) or the number of explanatory variables being greater than or equal to the number of selected factors. This generates a predictive model 20 (regression equation) with 5 or fewer analysis parameters as input. In some cases, the number of analysis parameters may be greater than 5.

[0024] [Food Evaluation Value Prediction] The prediction target analysis parameter acquisition unit 21 acquires multiple analysis parameters obtained by the analysis instrument 3. The analysis parameters to be acquired are five or fewer analysis parameters selected during the generation stage of the prediction model 20, and these differ for each prediction model 20. The analysis parameters are measured by the analysis instrument 3 in the same manner as the analysis parameters included in the training data D1 acquired by the training data acquisition unit 10.

[0025] As shown in Figure 1, the prediction model 20 constituting the prediction unit 22 is a model pre-trained to take multiple analysis parameters representing the analysis results of food by the analytical instrument 3 as input and output food evaluation values, as described above. The prediction unit 22 inputs the multiple analysis parameters acquired by the prediction target analysis parameter acquisition unit 21 into the prediction model 20 and outputs food evaluation values. The multiple analysis parameters input into the prediction model 20 do not include analysis parameters that were determined to be excluded in the first exclusion process and the second exclusion process. Furthermore, the number of multiple analysis parameters input into the prediction model 20 in the prediction stage is smaller than the number of multiple analysis parameters obtained by the multiple analytical instruments 3 in the prediction model 20 generation stage.

[0026] [Effects of this disclosure] Comparative Example 1 To demonstrate the effects of the above disclosure, Comparative Example 1 and Example 1 are shown below. In Comparative Example 1, the multiple analytical parameters obtained by the analytical instrument 3 (taste sensor) used in the generation stage of the prediction model 20 are only analytical parameters obtained from the taste sensor, and there are only two types for each taste sensor: Intensity_at_30_Sec_test_1 and Intensity_at_30_Sec_CPA_1. Therefore, the analytical parameters of Comparative Example 1 are 2 types of parameters × 8 types of taste sensors = 16 analytical parameters.

[0027] Example 1 In Example 1, the analytical parameters obtained by the analytical instrument 3 (taste sensor) used in the generation stage of the prediction model 20 are added to the analytical parameters obtained by the taste sensor in Comparative Example 1, and as described above, a total of approximately 560 analytical parameters are input from the eight types of taste sensors.

[0028] Comparative Example 1 and Example 1 for some of the food evaluation values ​​among the 25 types will be explained.

[0029] Figure 6 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1 and Example 1 for the food evaluation value of the first embodiment: RMinc of the TI method. Regarding the accuracy of the regression equation, Q 2 and cumulative R 2 Y and R 2 A value closer to 1 indicates higher prediction accuracy, while a value closer to 0 for RMSE (error) indicates smaller error.

[0030] Figure 7 shows the food evaluation value of the first embodiment: Comparative Example 1 for Tend of the TI method, the accuracy of the regression equation of Example 1, and the analytical parameters that serve as input to the regression equation.

[0031] Figure 8 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1, Example 1, and the food evaluation value of the TI method RMdec in the first embodiment.

[0032] Figure 9 shows Comparative Example 1, the accuracy of the regression equation for Example 1, and the analytical parameters used as input to the regression equation for the food evaluation value: Rinc of the TI method in the first embodiment.

[0033] Figure 10 shows Comparative Example 1, the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1, the Rdec of the TI method for the food evaluation value of the first embodiment.

[0034] Figure 11 shows Comparative Example 1, the accuracy of the regression equation for Example 1, and the analytical parameters used as input to the regression equation for the food evaluation value: DUR of the TI method in the first embodiment.

[0035] Figure 12 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1, Example 1, and the food evaluation value of the TI method Dpla in the first embodiment.

[0036] Figure 13 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1, Example 1, and the food evaluation value of the TI method for Dinc in the first embodiment.

[0037] Figure 14 shows Comparative Example 1, the accuracy of the regression equation for Example 1, and the analytical parameters used as input to the regression equation for the food evaluation value: Ddec of the TI method in the first embodiment.

[0038] Figure 15 shows the food evaluation value of the first embodiment: Comparative Example 1 for the AUC of the TI method, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation.

[0039] Figure 16 shows the food evaluation values ​​of the first embodiment: Comparative Example 1 for Apla of the TI method, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation.

[0040] Figure 17 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1, Example 1, and Ainc of the food evaluation value of the first embodiment:TI method.

[0041] Figure 18 shows Comparative Example 1, the accuracy of the regression equation for Example 1, and the analytical parameters used as input to the regression equation for the food evaluation value of the first embodiment: Adec of the TI method.

[0042] Figure 19 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1 and Example 1 regarding the astringency of the food evaluation value of the first embodiment using the QDA method.

[0043] Figure 20 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1 and Example 1 of the food evaluation value: QDA method for citrus in the first embodiment.

[0044] Figure 21 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1 and Example 1 of the food evaluation value of the first embodiment: herbal medicine using the QDA method.

[0045] Figure 22 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1 and Example 1 regarding the food evaluation value: kire (sharpness) of the QDA method in the first embodiment.

[0046] Figure 23 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1 and Example 1 regarding the lastingness of the food evaluation value of the first embodiment using the QDA method.

[0047] Figure 24 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1 and Example 1 regarding the rapidness of the food evaluation value of the first embodiment using the QDA method.

[0048] Figure 25 shows the spread of the food evaluation value:QDA method in the first embodiment, comparative example 1, the accuracy of the regression equation in example 1, and the analytical parameters used as input to the regression equation.

[0049] FIG. 26 is a diagram showing the accuracy of the regression equations and the analysis parameters serving as the inputs of the regression equations for Comparative Example 1 and Example 1 with respect to the thickness of the food evaluation value: QDA method in the first embodiment.

[0050] FIG. 27 is a diagram showing the accuracy of the regression equations and the analysis parameters serving as the inputs of the regression equations for Comparative Example 1 and Example 1 with respect to the woody of the food evaluation value: QDA method in the first embodiment.

[0051] Regarding Comparative Example 1 and Example 1 shown in FIGS. 6 to 27, Q 2 and the cumulative R 2 Y, and R 2 are all in a relationship where Example 1 is larger than Comparative Example 1, and RMSE is in a relationship where Example 1 is smaller than Comparative Example 1. That is, it can be seen that the prediction accuracy can be improved when using a large number of analysis parameters obtained by continuously measuring the intensities at a plurality of time points during the measurement period by the taste sensor (Example 1) compared to the case of using the analysis parameter indicating the intensity at the end point of the measurement by the taste sensor (30 seconds) (Comparative Example 1).

[0052] <Second Embodiment> [Food Evaluation Value (Sensory Evaluation Value, Sweetness Evaluation Value)] In the second embodiment, sweetness among the tastes is illustrated. Ten types of sweeteners and sweetening materials were used as food samples to evaluate only sweetness. Specifically, these were sucrose, sucralose, acesulfame K, advantame, thaumatin, stevia extract (A), stevia extract (B), stevia extract (C), monk fruit extract (A), and monk fruit extract (B). Twenty selected panelists compared each sample with a control (6% sucrose aqueous solution) and serially diluted samples, selecting the one they perceived as sweeter (two-point discrimination method). The concentration at which the selectivity rate was 50% was defined as the subjective equivalent concentration of sweetness, and samples were prepared to this concentration for subsequent evaluations and analyses. Each sample was evaluated using Quantitative Descriptive Analysis (QDA) and Time-Intensity Analysis (TI), and the following 22 food evaluation values ​​(sweetness evaluation values) were obtained.

[0053] [Food evaluation values ​​showing the qualitative characteristics of sweetness based on Quantitative Descriptive Analysis (QDA method)] A panel of 21 selected members determined seven characteristic parameters and their definitions that were suitable for the sample. Through trial evaluations, appropriate evaluation frequencies and sample quantities were determined. After sufficient evaluation training, the intensity of each characteristic parameter was evaluated on a 15cm line scale. The evaluation was repeated twice. Each sample was assigned a three-digit random number, and the liquid color was blinded during the testing. The following shows the 10 food evaluation values ​​(characteristic items) and their definitions. For example, the definition of the characteristic item "Top Intensity" is shown to be "Top Taste Sensation of Sweetness." • Initial Intensity: The initial taste sensation of sweetness (higher intensity indicates a stronger taste sensation). • Lastingness: A sweet aftertaste lingers, leaving a lingering sensation (the stronger the intensity, the longer the aftertaste). • Bitterness: Maximum bitterness intensity (higher intensity indicates stronger bitterness) • Astringency: Irritation, numbness, or bitterness in the mouth (the higher the intensity, the stronger the irritation). • Flat: Mild, light, bland (the higher the intensity, the milder it is) • Stickiness (Coating): Sticky, clinging, rough, persistent (the stronger the stickiness, the more it feels like it's sticking to the inside of the mouth) • Aroma (Caramel): Darkness of brown sugar, caramel, savory, soy sauce, soybeans, cooked flavor (the higher the intensity, the stronger the brown sugar flavor)

[0054] [Food evaluation values ​​showing the time-dependent characteristics of sweetness using the Time-Intensity (TI) method] A selection panel of 20 participants conducted trial evaluations to determine appropriate sample quantities, swallowing times, etc. Furthermore, after sufficient evaluation training, the temporal intensity of "sweetness" was evaluated using a slider bar on a 10cm line scale. The evaluation was repeated three times. As shown in Figure 2, in the second embodiment, there are 15 food evaluation values ​​that indicate the temporal characteristics of sweetness using the TI method. The 15 food evaluation values ​​and their definitions are as follows. • Rinc: Slope of increase • RMinc: Maximum slope of increase • Rdec: Slope of decrease • RMdec: Maximum slope of decrease • Tstart: The first time sweetness was detected. • Dinc: Increased time • Tmax: Time at which the intensity first reaches its maximum. • Dpla: Duration of the plateau of maximum intensity ·Ddec: Decrease time • Tendency: The time when you no longer taste sweetness. • DUR: Total time spent feeling sweetness Ainc: Area under the increasing curve • Apla: Area under the curve of the maximum intensity plateau • Adec: Area under the decreasing curve • AUC: Area under the curve of the entire curve

[0055] [Example of analysis parameter: Taste sensor] An example of analytical instrument 3 is a taste sensor. An example of analytical parameters obtained by the taste sensor will be described. The taste sensor is a sensor that obtains response information to the lipid membrane. Each sample was analyzed with nine types of taste sensors. The specific items of the nine types of taste sensors and the sensor identifiers are as follows. In the second embodiment, UM2 is added to the eight types of taste sensors in the first embodiment. The measurement method for each taste sensor is the same as in the first embodiment. ·Acidic bitterness:C00 • Hydrochloride bitter taste: BT0 • Basic salt bitterness: AN0 • Astringency: AE1 • Acidity: CA0 • Saltiness: CT0 Umami: AAE, UM2 • Sweetness: GL1 For each taste sensor, a first-stage initial taste measurement, a second-stage first aftertaste measurement, and a third-stage second aftertaste measurement were performed, and the response value over a 30-second measurement period was measured in 1-second increments at each stage. Figure 3 is a schematic explanatory diagram showing the measurement state of the taste sensor and the response value (mV) obtained from the sensor. The measurement method is as schematically shown in Figure 3 for membrane potential (sensor output), first by immersing the taste sensor in a reference solution to obtain a baseline membrane potential (sensor output). Next, when the taste sensor is immersed in the sample, the membrane potential changes due to interaction with the taste substance, and the amount or portion of the change corresponds to the initial taste (sometimes referred to as test_1). Next, the taste sensor is washed with the reference solution. Next, the taste sensor is immersed in the reference solution and the change in membrane potential is measured. The amount or portion of the change corresponds to the first aftertaste (sometimes referred to as CPA1). In the first embodiment, the taste sensor is further washed with the reference solution, and then the taste sensor is immersed in the reference solution and the change in membrane potential is measured. The amount or part of the change corresponds to the second aftertaste (sometimes referred to as CPA2).

[0056] The following describes the multiple analytical parameters obtained by the taste sensor in the second embodiment. Similar to the first embodiment, the analytical parameters are obtained based on the sensor intensity at multiple time points (0 seconds, 1 second, 2 seconds, ..., 29 seconds, 30 seconds) measured between the start and end of the measurement (30 seconds). The 13 basic analytical parameters obtained by the taste sensor and their definitions are as follows. ·Imax: Maximum intensity • Tmax: Time to reach maximum intensity • Plateau Intensity: 90% Imax (Maximum intensity plateau) • Plateau_Reach_Time: Time to reach 90% Imax (maximum intensity plateau) • Plateau_Duration: Duration of the plateau of maximum intensity • Intensity_at_1_Sec: Intensity at 1 second (start of measurement) • Intensity_at_3_Sec: Intensity at 3 seconds (start of measurement) • Intensity_at_30_Sec: Intensity at 30 seconds (end of measurement) • Slope_1_to_3: Slope from 1 second to 3 seconds. • Slope_3_to_Tmax: Slope from 3 seconds to Tmax • Slope_Tmax_to_30: Slope from Tmax to 30 seconds • Slope_3_to_Plateau: Slope from 3 seconds to 90% Imax · Slope_Plateau_to_Tmax: Slope from 90% Imax to Tmax In the first embodiment, a parameter relating to area is used, but in the second embodiment, the parameter relating to area is replaced with a parameter relating to slope. Since the parameter relating to area is the product of intensity and time, it is a parameter that includes intensity. However, a parameter relating to the slope of the line connecting the intensity at the first time point and the intensity at the second time point is a parameter that does not include intensity, thus differentiating it from parameters that include intensity, such as Imax. The 13 parameters listed above are basic parameters. Since the initial taste, the first aftertaste, and the second aftertaste each have the 13 reference parameters listed above, there are a total of 39 parameters. For example, the maximum intensity of the initial taste is denoted as Imax_test_1, the maximum intensity of the first aftertaste is denoted as Imax_CPA_1, and the maximum intensity of the second aftertaste is denoted as Imax_CPA_2. Furthermore, for each basic parameter, it becomes possible to generate combination parameters that represent the difference (def) between two combinations of the initial taste, first aftertaste, and second aftertaste. For example, the difference between the initial taste and second aftertaste in the AUC_total of the taste sensor (C00) can be expressed as TasteSensor_SB2C00_AUC_total_diff_test_1_vs_CPA_2. Based on one basic parameter, there are three combination parameters: the initial taste and first aftertaste, the initial taste and second aftertaste, and the first aftertaste and second aftertaste. Therefore, it becomes possible to generate three combination parameters based on one basic parameter.

[0057] In taste sensor analysis, measurement error correction is performed using interpolation and difference. In this correction, a regression equation is created with the measured value of sample A, designated as a control, as the dependent variable and the measurement order as the independent variable. An interpolated value is calculated that takes into account the state change of the waiting sample, and the difference is calculated from the test data. In situations where a curve is drawn from measured values, as in this method, a smooth curve may not be drawn when the difference is calculated due to the slope and inflection point of the control. Therefore, in this method, the curve drawn using raw data without corrected interpolation and difference processing was used to calculate parameters (Tmax, ...Time, Slope...) based on the curve shape. For the curve, the response values ​​at each time point were arithmetic mean of the data from three trials that were smoothed using Savizchy-Golay, and each parameter was calculated from the resulting average curve.

[0058] Regarding the method for generating the predictive model, the analysis parameters obtained from the taste sensors included in the training data D1 consist of 9 types of taste sensors × (13 basic parameters + 3 × 13 combination parameters) = approximately 468 parameters. Regarding the method for generating the predictive model, before performing the first elimination process in the first embodiment, the coefficient of variation (CV) of each explanatory variable was calculated, and variables whose coefficient of variation (CV) was smaller than a predetermined threshold (e.g., 0.05) were removed from the explanatory variables.

[0059] [Effects of this disclosure] To demonstrate the effects of the present disclosure described above, Comparative Example 1 and Example 1 are shown below. In Comparative Example 1, the analytical parameters obtained by the analytical instrument 3 used in the generation stage of the prediction model 20 are only those obtained from the taste sensor, and for each taste sensor, there are only two types: Intensity_at_30_Sec_test_1 and Intensity_at_30_Sec_CPA_1. Therefore, the analytical parameters in Comparative Example 1 are 2 types of parameters × 9 types of taste sensors = 18 analytical parameters.

[0060] In Example 1, the analytical parameters obtained by the analytical instrument 3 (taste sensor) used in the generation stage of the prediction model 20 are added to the analytical parameters from the taste sensor in Comparative Example 1, and as described above, a total of approximately 468 analytical parameters are input from the nine types of taste sensors.

[0061] Comparative Example 1 and Example 1 for some of the food evaluation values ​​among the 22 types will be explained.

[0062] Figure 28 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1, Example 1, and Rinc of the food evaluation value of the second embodiment:TI method.

[0063] Figure 29 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1, Example 1, and RMinc of the food evaluation value of the second embodiment:TI method.

[0064] Figure 30 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1, Example 1, and Rdec of the food evaluation value of the second embodiment:TI method.

[0065] Figure 31 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1 and Example 1 for the RMdec of the TI method in the second embodiment of food evaluation values.

[0066] Figure 32 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1, Example 1, and Tstart of the TI method for the food evaluation value of the second embodiment.

[0067] Figure 33 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1, Example 1, and Dinc of the food evaluation value of the second embodiment:TI method.

[0068] Figure 34 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1, Example 1, and Tmax of the food evaluation value of the second embodiment:TI method.

[0069] Figure 35 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1 and Example 1 for the food evaluation value of the second embodiment: Dpla of the TI method.

[0070] Figure 36 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1, Example 1, and Ddec of the food evaluation value of the second embodiment:TI method.

[0071] Figure 37 shows the food evaluation value of the second embodiment: Comparative Example 1, the accuracy of the regression equation of Example 1, and the analytical parameters used as input to the regression equation for the Tend of the TI method.

[0072] Figure 38 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1, Example 1, and the food evaluation value of the second embodiment: DUR of the TI method.

[0073] Figure 39 shows the accuracy of the regression equations and the analytical parameters used as inputs to the regression equations for Comparative Example 1, Examples 1 and 2, and for Ainc of the food evaluation value of the second embodiment: TI method.

[0074] Figure 40 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1, Example 1, and Apla of the TI method for the food evaluation value of the second embodiment.

[0075] Figure 41 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1, Example 1, and Adec of the food evaluation value of the second embodiment:TI method.

[0076] Figure 42 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1, Example 1, and the AUC of the TI method for the food evaluation value of the second embodiment.

[0077] Figure 43 shows the accuracy of the regression equations and the analytical parameters used as inputs to the regression equations for Comparative Example 1, Examples 1 and 2, regarding the initial intensity of the top food evaluation value of the QDA method in the second embodiment.

[0078] Figure 44 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1 and Example 1 regarding the food evaluation value of the second embodiment: lastingness of the QDA method.

[0079] Figure 45 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1 and Example 1 regarding the bitterness of the food evaluation value: QDA method in the second embodiment.

[0080] Figure 46 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1 and Example 1 regarding the astringency of the food evaluation value of the second embodiment using the QDA method.

[0081] Figure 47 shows the accuracy of the regression equations and the analytical parameters used as inputs to the regression equations for Comparative Example 1, Examples 1 and 2, and for the monotonic (Flat) food evaluation value of the second embodiment of the QDA method.

[0082] Figure 48 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1 and Example 1 regarding the food evaluation value of the QDA method for coating in the second embodiment.

[0083] Figure 49 shows the accuracy of the regression equation and the analytical parameters used as input to the regression equation for Comparative Example 1 and Example 1 regarding the food evaluation value of the second embodiment: aroma (Caramel) using the QDA method.

[0084] Regarding Comparative Example 1 and Example 1 shown in Figures 28 to 49, Q 2 And, cumulative R 2 Y and R 2 In all cases, Example 1 showed a larger correlation with Comparative Example 1, while the RMSE was smaller for Example 1 than for Comparative Example 1. In other words, it can be seen that using a large number of analytical parameters obtained by continuously measuring the intensity at multiple time points during the measurement period with the taste sensor (Example 1) can improve prediction accuracy compared to using analytical parameters that show the intensity at the end of the measurement with the taste sensor (30 seconds) (Comparative Example 1).

[0085] Therefore, it can be understood that the food evaluation value prediction method of this disclosure is effective whether the food evaluation value for bitterness is as shown in the first embodiment or for sweetness as shown in the second embodiment.

[0086] [Selection rate for each parameter] Figure 50 is an explanatory diagram showing the proportion of each parameter selected as an important variable. As shown in Figure 50, since the contributing components and receptor properties differ for each of the five tastes, the types of important variables selected by regression analysis showed characteristics specific to each taste.

[0087] In sweetness prediction and bitterness prediction, the types of sensors selected as important variables were similar, with selection rates concentrated on sweetness, umami, and bitterness sensors. All three tastes are known to be sensed by G protein-coupled receptors. The concentration of selection rates is thought to be due to the fact that the perception mechanism is significantly different from that of sourness and saltiness, which are ion channel receptors. Furthermore, in both sweetness and bitterness prediction, the bitterness sensor had the highest selection rate. Sweeteners can have not only sweetness but also characteristic bitterness and astringency. Therefore, it is thought that bitterness sensors were selected as an important variable in many cases in sweetness prediction as well.

[0088] In sweetness prediction, the usage rate of CPA_2 (second aftertaste) tended to be lower than in bitterness prediction. Since many bitter substances are lipid-soluble and their adsorption time to receptors and oral lipid membranes contributes to sensory persistence, the value of CPA_2 was selected more often as an important variable. On the other hand, in the sensory evaluation results for sweetness, "initial intensity" and "persistence" showed a high correlation, suggesting that the intensity of the sensation contributes to the aftertaste, and therefore the value of CPA_2, which indicates membrane adsorption capacity, was not effective in prediction.

[0089] [1] As in the embodiment described above, the food evaluation value prediction method includes the steps of: acquiring a plurality of analytical parameters that indicate the analysis results of a food to be predicted by an analytical instrument 3 including a taste sensor; and obtaining a food evaluation value by inputting the acquired plurality of analytical parameters into a prediction model 20 that has been pre-trained to take the plurality of analytical parameters that indicate the analysis results of a food by an analytical instrument 3 including a taste sensor as input and output a food evaluation value that quantifies one of the senses of taste, texture, and smell of the food. The analytical parameters obtained by the taste sensor may include at least one of the following, which are obtained based on the sensor intensity at multiple time points measured between the start of measurement and the end of measurement: intensity before the end of measurement, maximum intensity, intensity at a predetermined ratio to the maximum intensity, time to reach maximum intensity, time to reach intensity at a predetermined ratio to the maximum intensity, duration for which the intensity is equal to or greater than a predetermined ratio to the maximum intensity, area determined by the intensity over a predetermined period, and the slope of the line connecting the intensity at the first time point and the intensity at the second time point. As shown in the comparison between Comparative Example 1 and Example 1, using a large number of analytical parameters obtained by continuously measuring the intensity at multiple time points during the measurement period using a taste sensor (Example 1) improves the prediction accuracy of food evaluation values ​​compared to using analytical parameters that show the intensity at the end of the measurement by the taste sensor (30 seconds) (Comparative Example 1).

[0090] [2] The food evaluation value prediction method described in [1] above may include analytical parameters obtained by a taste sensor, which are obtained by aftertaste measurements from the third stage onward, following the first stage of initial taste measurement and the second stage of aftertaste measurement. By using analytical parameters obtained not only from the first stage of pre-taste measurement and the second stage of aftertaste measurement, but also from aftertaste measurement in the third stage and beyond, the prediction accuracy of machine learning can be improved.

[0091] [3] The food evaluation value prediction method described in [1] or [2] above may include a value corresponding to the difference or ratio between the analytical parameter obtained by the taste sensor and the analytical parameter obtained by the measurement at a predetermined stage among multiple stages of measurement, including pre-taste measurement and after-taste measurement. In this way, the analytical parameters of the taste sensor can be increased, enriching the expression of taste, which in turn improves the prediction accuracy of machine learning.

[0092] [4] As in the embodiment described above, the food evaluation value prediction model generation method includes the step of generating a prediction model 20 that outputs a texture evaluation value, using training data D1 which takes as input multiple analysis parameters indicating the analysis results of food by an analytical instrument 3 including a taste sensor, and outputs a food evaluation value that quantifies one of the senses of taste, texture, and smell of food. The analysis parameters obtained by the taste sensor may include at least one of the following, obtained based on the sensor intensity at multiple time points measured between the start of measurement and the end of measurement: intensity before the end of measurement, maximum intensity, intensity at a predetermined ratio to the maximum intensity, time to reach the maximum intensity, time to reach an intensity at a predetermined ratio to the maximum intensity, duration for which the intensity is equal to or greater than a predetermined ratio to the maximum intensity, area determined by the intensity over a predetermined period, and the slope of the line connecting the intensity at the first time point and the intensity at the second time point.

[0093] [5] As in the above embodiment, the program may cause one or more processors 1b, 2b to execute the method described in any of [1] to [4] above.

[0094] The computer-readable temporary recording medium according to the above embodiment stores the above program.

[0095] Although embodiments of this disclosure have been described above with reference to the drawings, it should be understood that the specific configurations are not limited to these embodiments. The scope of this disclosure is indicated not only by the description of the embodiments above but also by the claims, and further includes all modifications within the meaning and scope equivalent to the claims.

[0096] The structures adopted in each of the above embodiments can be adopted in any other embodiment. The specific configuration of each part is not limited to the embodiments described above, and various modifications are possible without departing from the spirit of this disclosure.

[0097] <Variation> (A) In the above embodiments, bitterness is given as an example of a food evaluation value that quantifies one of the senses of taste, texture, and smell of food, but the invention is not limited to this. For example, this disclosure can also be applied to food evaluation values ​​that quantify tastes such as sweetness, umami, depth of flavor, saltiness, and sourness as tastes of food. Furthermore, since the analytical method using analytical instrument 3 mentioned in the examples can analyze aroma-contributing components, it can also be applied to food evaluation values ​​that quantify smell. Furthermore, it can also be applied to food evaluation values ​​that quantify textures such as chewiness and biteability of food. Since the analytical method using analytical instrument 3 mentioned in the examples can analyze components that contribute to viscosity and hardness, it can also be applied to food evaluation values ​​that quantify texture.

[0098] (B) In the above embodiment, the number of types of food evaluation values ​​is 25, but is not limited to this. If one type of food evaluation value is sufficient, the number of types of food evaluation values ​​may be just one, and the number of types of food evaluation values ​​can also be changed in various ways.

[0099] (C) The sensory evaluation method for food evaluation values ​​is not limited to the above. For example, the number of evaluators can be changed as appropriate.

[0100] (D) In ​​the above embodiment, both the first removal process and the second removal process are performed, but the embodiment is not limited to this. Only the first removal process may be performed, or only the second removal process may be performed, and the prediction accuracy can be improved compared to when no removal process is performed.

[0101] (E) In the above embodiment, the analytical parameters input to the prediction model 20 are parameters obtained from a taste sensor, but are not limited to these. In addition to the parameters obtained from the taste sensor described above, parameters obtained from other analytical instruments 3 may also be added. Examples of other analytical instruments 3 include real-time direct mass spectrometers (DART-MS), near-infrared analyzers (NIRS), nuclear magnetic resonance spectrometers (NMR), and gas chromatography-mass spectrometers (GC-MS). [Explanation of symbols]

[0102] 1b: Processor 2b: Processor 3:Analytical equipment 20: Predictive Models D1: Training data

Claims

1. The steps include obtaining multiple analytical parameters that represent the analysis results of the food to be predicted using analytical instruments including a taste sensor, and The process involves taking multiple analytical parameters representing the results of food analysis performed by an analytical instrument including a taste sensor as input, and inputting the acquired multiple analytical parameters into a pre-trained predictive model that outputs a food evaluation value that quantifies one of the senses of taste, texture, and smell of the food, in order to obtain the food evaluation value. Includes, A method for predicting food evaluation values, wherein the analytical parameters obtained by the taste sensor include at least one of the following, based on the sensor intensity at multiple time points measured between the start and end of the measurement: intensity before the end of the measurement, maximum intensity, intensity at a predetermined ratio to the maximum intensity, time to reach the maximum intensity, time to reach the intensity at a predetermined ratio to the maximum intensity, duration for which the intensity is equal to or greater than a predetermined ratio to the maximum intensity, area determined by the intensity over a predetermined period, and the slope of the line connecting the intensity at the first time point and the intensity at the second time point.

2. The method for predicting food evaluation values ​​according to claim 1, wherein the analytical parameters obtained by the taste sensor include analytical parameters obtained from a third stage and subsequent aftertaste measurements, which are performed after the first stage of initial taste measurement and the second stage of aftertaste measurement.

3. The food evaluation value prediction method according to claim 1 or 2, wherein the analytical parameters obtained by the taste sensor include values ​​corresponding to the difference or ratio between analytical parameters obtained from measurements at a predetermined stage and analytical parameters obtained from measurements at stages other than the predetermined stage, among a plurality of measurements including pre-taste measurement and after-taste measurement.

4. The process includes the step of generating a predictive model that outputs a texture evaluation value, using training data that takes multiple analytical parameters representing the results of food analysis by analytical instruments including a taste sensor as input, and outputs a food evaluation value that quantifies one of the senses of taste, texture, and smell of the food. A method for generating a food evaluation value prediction model, wherein the analytical parameters obtained by the taste sensor include at least one of the following, based on the sensor intensity at multiple time points measured between the start and end of the measurement: intensity before the end of the measurement, maximum intensity, intensity at a predetermined ratio to the maximum intensity, time to reach the maximum intensity, time to reach an intensity at a predetermined ratio to the maximum intensity, duration for which the intensity is equal to or greater than a predetermined ratio to the maximum intensity, area determined by the intensity over a predetermined period, and the slope of the line connecting the intensity at the first time point and the intensity at the second time point.

5. A program that causes one or more processors to execute the method according to claim 1, 2, or 4.

Images