Next Article in Journal
Pharmaceutical Applications of Molecular Tweezers, Clefts and Clips
Next Article in Special Issue
Investigation of Direct Model Transferability Using Miniature Near-Infrared Spectrometers
Previous Article in Journal
Carotenoids: How Effective Are They to Prevent Age-Related Diseases?
Previous Article in Special Issue
Detecting Zn(II) Ions in Live Cells with Near-Infrared Fluorescent Probes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Calibration Transfer Based on Affine Invariance for NIR without Transfer Standards

1
School of Computer Science and Engineering, Northeastern University, Shenyang 110819, China
2
College of Information Science and Engineering, Northeastern University, Shenyang 110819, China
3
Institute of Automation, Chinese Academy of Sciences, Beijing 100190, China
*
Author to whom correspondence should be addressed.
Molecules 2019, 24(9), 1802; https://doi.org/10.3390/molecules24091802
Submission received: 26 March 2019 / Revised: 4 May 2019 / Accepted: 6 May 2019 / Published: 9 May 2019

Abstract

:
Calibration transfer is an important field for near-infrared (NIR) spectroscopy in practical applications. However, most transfer methods are constructed with standard samples, which are expensive and difficult to obtain. Taking this problem into account, this paper proposes a calibration transfer method based on affine invariance without transfer standards (CTAI). Our method can be utilized to adjust the difference between two instruments by affine transformation. CTAI firstly establishes a partial least squares (PLS) model of the master instrument to obtain score matrices and predicted values of the two instruments, and then the regression coefficients between each of the score vectors and predicted values are computed for the master instrument and the slave instrument, respectively. Next, angles and biases are calculated between the regression coefficients of the master instrument and the corresponding regression coefficients of the slave instrument, respectively. Finally, by introducing affine transformation, new samples are predicted based on the obtained angles and biases. A comparative study between CTAI and the other five methods was conducted, and the performances of these algorithms were tested with two NIR spectral datasets. The obtained experimental results show clearly that, in general CTAI is more robust and can also achieve the best Root Mean Square Error of test sets (RMSEPs). In addition, the results of statistical difference with the Wilcoxon signed rank test show that CTAI is generally better than the others, and at least statistically the same.

1. Introduction

With the characteristics of high efficiency, low cost and non-destructivity, near-infrared (NIR) spectroscopy has been widely used in control of food and pharmaceutical quality [1,2,3,4]. Multivariate calibration methods are commonly used to obtain quantitative or qualitative information from near-infrared spectra, such as principal component regression (PCR) [5,6] and partial least squares (PLS) [7,8,9,10]. Since changes of the instruments and measurement conditions may result in poor applicability of the model. Recalibration can be utilized to solve this problem, but recalibration is time consuming and takes an immense amount of work. In order to reduce consumption of the recalibration, calibration transfer has been widely studied and applied [11]. There are two main situations about calibration transfer: (1) The uniform calibration model is used to predict spectra being measured on multiple instruments; (2) the new spectra are measured on the same instrument after a period of time.
A number of related methods for calibration model transfer have been proposed, which are divided into two categories. Ones require transfer standards and ones not require transfer standards. The first category of methods has the characteristic that a set of samples are separately measured on the master and slave instrument. A great variety of transfer methods with standard samples have been proposed. For examples, SBC [12,13] assumes a linear relationship between predicted values of different instruments. First, the regression coefficient between the spectra and the response values on the master instrument is calculated. Then the predicted values of the master and slave setting are computed based on the regression coefficient. Finally, a linear equation is fitted between the predicted values. PDS proposed by Wang et al. is employed to correct the spectral differences [14]. In PDS [15,16,17,18], each wavelength of the master instrument is related to the wavelength window of the slave instrument, and a band transfer matrix is finally formed based on the regression coefficients of each window. The observation is consistent with this assumption that in various transfer methods the spectral correlation between master and slave is limited to smaller regions. The keys to PDS are the selection of window size and the number of standard samples. Due to the construction of multiple regression models, a huge amount of calculations are desired. The calibration model transfer for near-infrared spectra based on canonical correlation analysis [19] is proposed by Liang et al. The PLS model is built using the master instrument calibration set, and a part of the calibration set of master and slave instrument is taken as standard samples. Then, the features extracted respectively by canonical correlation analysis (CCA) [20,21]. The relationship between master and slave data is established with ordinary least squares (OLS) [22,23], and the test set is finally corrected. For CCA, SBC and PDS, a good result can be achieved with standard samples, but standard samples are difficult to obtain in some cases. For the transfer methods such as calibration transfer via extreme learning machine auto-encoder (TEAM) [24] method, calibration transfer by generalized least squares (GLSW) [25] method and spectral space transform (SST) [26,27] and so on, standard samples are also required, although the principles of these methods are different.
The second category is the methods without transfer standards. For examples, multiplicative scatter correction (MSC) [28,29,30] proposed by Bouveresse et al. first calculates the mean spectra of the calibration set as the reference spectra, then the linear relationship is found between every spectra and the reference spectra, and the slope and bias are obtained; finally, the slope and bias are utilized to correct slave spectra. While the standard samples are not required in MSC, it is difficult to handle complex situations. MSC is a transfer method using pre-processing techniques, and more pre-processing approaches include finite impulse response (FIR) [31] filtering and multivariate filtering via orthogonal signal correction (OSC) [32,33], etc. TCR [34] is also a standard-free method which combines transfer component analysis (TCA) [35] and ordinary least squares (OLS). The basic idea of TCA is to project the data of two instruments in a Reproducing Kernel Hilbert Space, where the data are distributed as close as possible at the same time preserving the key attributes of the original data. TCR is a robust model with good generalization abilities, but does not achieve more accurate predictions. Other techniques belonging to this category include kernel principal component analysis (KPCA) [36,37], domain generalization via invariant feature representation (DICA) [38] and so on.
Different from the above methods, this paper studies the relationship of regression coefficients between the feature vector and predicted values on two spectrometers. Samples of the calibration transfer method based on affine invariance without transfer standards (CTAI) are shown in Figure 1A. The response values of the slave spectrometer are not required, and the map is not necessary between master and slave samples. The samples are further processed under the PLS model. The spectral features and prediction values are respectively obtained, and the processed samples are shown in Figure 1B. We obtain the linear models between the feature vector and the predicted values respectively. According to the linear models of two instruments, the relationship between the predicted values is further obtained. Firstly, the PLS model is built on the master instrument; secondly, the score matrices and predicted values are extracted according to the PLS model, respectively; further, the angles and biases are calculated between two regression coefficients; finally, the prediction values are corrected by affine transformation. If the concentration information of the master spectra and the slave spectra are in the same range, CTAI can achieve more accurate predicted results and more robust model even without standard samples compared with other methods. The predictive performance of CTAI is verified by two near-infrared (NIR) datasets.

2. Results and Discussion

2.1. Analysis of the Corn Dataset

The training errors, prediction errors, cross-validation errors, biases and the correlation coefficients for the predicted vs. actual results about the PLS model of the corn dataset are shown in Table 1. Large correlation coefficients and small biases can be seen in all results. The results reflect a good linear relationship between the spectra and measured values of the corn dataset. There are no significant differences between Root Mean Square Error of calibration set (RMSEC), Minimum Root Mean Square Error of Cross-Validation (RMSECV) and Root Mean Square Error of test set (RMSEP), indicating that there is no over-fitting and under-fitting phenomenon, which can explain the reasonable selection of the number of latent variables. Moreover, we can see that RMSEPm of the PLS on the instrument m5spec are smaller than the RMSEPm of the instrument mp6spec. For most calibration transfer methods, it is important that the master instrument has more accurate prediction results. Thus, m5spec as the master instrument and mp6spec as the slave instrument is a more reasonable choice.
In order to more fully assess the predicted performance of CTAI, the methods MSC, TCR, CCA, SBC and PDS are tested. In this work, when PDS was performed, PLS was utilized to compute the transformation function. For the PLS model, the optimal number of latent variables is shown in Table 1. The optimal dimensionality of the subspace in TCR is 4, 6, 10 and 10. In addition, optimal window sizes of PDS are all 3. We set the standard samples in range [5,30]. When the model is stable, the number of standard samples is selected for modeling based on the smallest RMSEC criteria.
As shown in Table 2, we can see the correlation coefficients rpre and corresponding ppre values, which indicate the prediction values between the master instrument and the slave instrument are linearly correlated. We can also see that the tpre is greater than the t critical value. We then know the bias adjustment in predicted results should be implemented. Furthermore, the RMSE of prediction without any correction for the slave instrument shows more error of prediction than the master instrument. The corrected results of CTAI result in a significant reduction in RMSE of prediction. The same situation can be found between y m and y ˜ n in Table 2. The absolute value of t in each component is 15.437, 19.657, 19.408 and 8.762, respectively. The critical value of t is 2.131, and all results are greater than it. It is further proved that the adjustment of bias is very important. For the corn dataset, the effect of correction in CTAI is vividly described in Figure 2. It can be seen that the corrected predicted values of CTAI more close to the straight line, and RMSEP is greatly reduced.
Moreover, the results listed in Table 3 and Table 4 show the difference between the 16 predictive corn samples by different methods. In general, the results of CTAI exhibit the best performance for prediction compared to other five methods. When moisture is used as the property, CTAI achieves the lowest RMSEP (0.21095). More specifically, the RMSEP improvements provided by CTAI with respect to MSC, TCR, CCA, SBC and PDS are as high as 87.35%, 46%, 9.48%, 50.45% and 12.96%, respectively. Though there are no statistically significant differences, CTAI is greatly improved in predictive accuracy compared with CCA and TCR. There is a significant difference at the 95% confidence level between CTAI and MSC, SBC and PDS. When oil is used as the property, it can be seen that there is no significant difference between RMSEC and RMSEP in different transfer methods, so the over-fitting phenomenon does not appear. CTAI also produces the lowest RMSECV (0.08141) and RMSEC (0.08233). The results by Wilcoxon signed rank test reveal that CTAI is significantly different from MSC and TCR and has similar performance compared with CCA, SBC and PDS. It is noticeable that the RMSEP improvement rates of CTAI compared with CCA, SBC and PDS are 27.98%, 1.52% and 13.28%, respectively. Other properties are similar with the property of oil; CTAI achieves better predictive performance.
In order to compare the predictive stability of various methods, Figure 3, Figure 4, Figure 5 and Figure 6 show the plots of measured vs. predicted values for the calibration set and the test set. If the model predicts better, the point will be closer to the straight line. When moisture is used as the property, it is observed from Figure 3 that CTAI is in general closer to the straight line than the other models. It confirms that the CTAI achieves the best overall performance. When oil is used as the property, it is clear that CTAI provides satisfactory results not only in the calibration set but also in the test set. It reconfirmed that CTAI achieves more accurate prediction results. In addition, the standard error has also achieves good results in CTAI compared with others. From the discussion above, one can easily conclude that CTAI can achieve the best performance in all models and has better generalization ability.

2.2. Analysis of the Wheat Dataset

The RMSEP of the PLS model is listed in Table 1. We can see that the predicted performance of the instrument B1 is better than B3 and the instrument B3 is better than B2. Thus, three combinations (B1-B2; B1-B3; B3-B2) of the instruments B1, B2 and B3 are used to analyze the wheat dataset. The first instrument of every combination stands for master instrument and the second instrument stands for slave instrument. For PLS model, the optimal number of latent variables is 14, 15 and 15, respectively, and the corresponding optimal dimensionality of the subspace in TCR is 17, 12 and 17, respectively. Moreover, the optimal number of window sizes for B1-B2, B1-B3 and B3-B2 is 3, 9 and 13, respectively.
For the three combinations of instruments (B1-B2; B1-B3; B3-B2), we can see between y m and y ˜ n the correlation coefficients rpre are large and ppre are close to zero in Table 2. Hence, there is a linear relationship between the predicted values of the two instruments for wheat dataset. For all combinations, the absolute value of t is greater than tcritical_value. So there is a significant bias between uncorrected predicted values of the slave instrument and predicted values of the master instrument. So we can correct the predicted values of the slave instrument by affine transformation. The experimental results show that the prediction performance of CTAI is significantly enhanced. We found the same phenomenon for the uncorrected prediction values of the slave instrument relative to the master instrument actual values. Furthermore, for the predicted performance of CTAI, Figure 7 shows the difference between uncorrected and corrected predicted values for B1-B2, B1-B3 and B3-B2. It can be seen that CTAI plays an important role in the correction of predicted values.
In addition, Table 3 lists the results of different methods for calibration set and test set. For the B1-B2, CTAI produces the lowest RMSEP (0.41419) and the second lowest RMSEC (0.55682). For PDS and CCA, it is worth noting that RMSEP is significantly larger than RMSEC. Therefore, the predictive performance of PDS and CCA are poor under this setting. Further, a statistical testing is utilized to evaluate the RMSEP difference between the CTAI and other methods for the wheat dataset. The Wilcoxon signed rank sum test was performed and at the significance level alpha = 0.05. It can be seen from Table 4 that there is a statistically significant difference compared with CCA, SBC and PDS. In addition, the improvement rates of prediction provided by CTAI for MSC and TCR are up to 55.07% and 52.32%, respectively. For the combination (B1-B3), CTAI displays the lowest RMSEP (0.68215), followed by TCR (0.72996) and SBC (0.79294). For PDS, we can see that under-fitting still existed under this setting, and for CCA, this phenomenon also exists, but it is not particularly serious. The results by Wilcoxon signed rank test show that CTAI is significantly different from MSC, CCA, SBC and PDS (shown in Table 4). Compared with TCR, RMSEP improvement rates of CTAI can reach 6.55%. For the last combination, both RMSEP and RMSEC achieve the best predicted results. Further, except for PDS, the differences between CTAI and other models are statistically significant at the 95% confidence level. Compared with PDS, the RMSEP improvements of CTAI are as high as 79.05%. It is also worth noting that there is no under-fitting phenomenon in PDS under the current setting, but the predicted results are still poor. Therefore, the predictive performance of PDS is worse for wheat datasets under the current model.
To further display the predictive abilities of different models, the correlation between measured and predicted values obtained in Figure 8, Figure 9 and Figure 10. Zero differences between measured and predicted values result in points over the straight line of the plot. It can be seen that good correlations are found between expected and predicted concentrations, which confirm the good performance of CTAI. CTAI achieved the lowest standard error for three combinations. Moreover, the predictive abilities of PDS and CCA are poor for wheat dataset. For SBC, PDS and CCA, they require standard samples and TCR requires reference values of the slave instrument samples, both of which are expensive and difficult to obtain. Obviously, this means that CTAI shows much more outstanding performance.

3. Materials and Methods

3.1. Dataset Description

3.1.1. Corn Dataset

The corn dataset, which contains 80 samples, was measured on three NIR spectrometers (m5, mp5 and mp6). Each sample consists of four components: Moisture, oil, protein, and starch. The wavelength range is 1100–2400 nm with interval 2 nm (700 channels). The spectra measured in m5spec were used as the master spectra, and the spectra measured by mp6spec were used as the secondary spectra. The data can be obtained from http://www.eigenvector.com/data/Corn/. The dataset was divided into a calibration set of 64 samples and a test set of 16 samples based on Kennard-Stone (KS) algorithm. The NIR spectra are shown in Figure 11A, which represents the difference between m5 and mp6.

3.1.2. Wheat Dataset

The wheat dataset was used as the shootout data for the International Diffuse Conference 2016, and the protein content was chosen as the property. Related information about the wheat dataset at http://www.idrc-chambersburg.org/content.aspx?page_id=22&club_id=409746&module_id=191116 can be easily accessed. 248 samples of the wheat dataset from three different NIR instrument manufacturers (B1, B2 and B3) were analyzed. According to KS algorithm, 198 samples were chosen as the calibration set and the remainder of samples formed the test set. The wavelength range is 570–1100 nm with an interval of 0.5 nm. The spectral difference between B1 and B2 is shown in Figure 11B. The spectral difference between B1 and B3 is shown in Figure 11C. The spectral difference between B2 and B3 is shown in Figure 11D.

3.2. Determination of the Optimal Parameters

Latent variables of PLS in CTAI are allowed to take values in the set [1,15], and it is determined by the 10-fold cross-validation. The optimal number of latent variables is selected only when the lowest RMSECV.
Five methods were used for comparison, where the latent variable range and parameter optimization all of SBC, CCA, PDS and MSC in PLS are consistent with CTAI. In particular, the window size in PDS is searched for from 3 to 16 in increments of 2, and is selected by 5-fold cross-validation. In addition, the dimensionality of the TCA space in TCR is estimated in the range [1,24] and the optimization criteria are consistent as described in [24].

3.3. Model Performance Evaluation

In this experiment, root mean squared error RMSE is employed as indicators for parameter selection and model evaluation. Furthermore, RMSEC is the training error, RMSECV denotes the cross-validation error and RMSEP indicates the prediction error of the test set. The RMSE calculation method is written as:
RMSE = ( y y ^ ) T ( y y ^ ) / n
where y ^ is the predict value, y is the measured value and n represents the number of samples.
Bias and standard error (SE) are also utilized as reference indicators for model evaluation. The bias and SE are as follows:
{ bias = i n ( y i y ^ i ) / n SE   =   ( ( y y ^ ) T ( y y ^ ) bias ) / n
Moreover, the Pearson correlation coefficient and corresponding test is used to determine if there is a linear relationship between the master instrument and the slave instrument. One-Sample t-Test is also utilized to determine whether a bias adjustment in predicted results should be implemented [11].
In order to compare CTAI and other methods further, another important parameter (h) is cited in order to compare the rate of improvement, defined as follows:
h = ( 1 RMSEP RMSEP other ) × 100 %
where RMSEP represents the prediction error of CTAI and RMSEP other represents the others.
In addition, the Wilcoxon signed rank sum test at the 95% confidence level is used to determine whether there is a significant difference between CTAI and the others.

3.4. Computational Environment

All experimental procedures were implemented on a personal computer by python language, software version python 2.7, and run on an acer notebook with a 2.60 GHz Intel (R) Core (TM) i5-3230M CPU, 8 GB RAM and a Microsoft Windows 7 operating system (Acer Incorporated, Taiwan, China). Normalization and cross-validation are performed using the sklearn package. The Wilcoxon signed rank test is implemented using the scipy package and other programs are implemented by the individual.

3.5. Calibration Transfer

3.5.1. Notation

In the following text, matrices are represented by bold capital letters (e.g., X), column vectors by bold lower case letters (e.g., y) and scalars by italic letters (e.g., a). The transposition operation is indicated by superscript T.

3.5.2. Overview of PLS

PLS is used to establish the linear relationship between the input space and the response space. The purpose of the PLS model is to ensure the optimal number of latent variables. The latent variables are linear combinations of the primitive variables. The latent variables are calculated in this way so that they contain a maximum of relevant information concerning the relation between X and y . Mathematically, this is shown by the following objective function.
H = argmax w           cov X w ,   y subject   to       | | w | | 2 = 1
where w represents the weight vector. This objective is a maximization problem under one constraint, which can be settled in virtue of the Lagrange multiplier method.
Assuming a PLS model is built between spectral matrix X n × p and concentration vector y n × 1 , the model is named PLS1 (n denotes the number of samples and p represents the optimal numbers of latent variables). In the algorithm, the first weighting vector must be the primary eigenvector of the matrix X T y y T X . From the second latent variable on, it requires the following latent variables to be orthogonal (uncorrelated) to the former ones. Hence, the following weighting vectors will be the dominant eigenvectors of the matrix X T y y T X ; also, repeat a sequence of the steps until convergence. The PLS1 is built using the following model:
{ X n × p = T n × A ( P p × A ) T + E n × p y n × 1 = T n × A ( Q 1 × A ) T + F n × 1
where T is the score matrix and P and Q represent the X-loading matrix and y-loading vector, respectively; E and F denote the matrix of residuals; A is the optimal number of principal components over the master instrument PLS model.
Finally, the regression coefficient β of the model can be written as follows:
β = W ( P T W ) - 1 Q T
where W = [ w 1 , w 2 , , w A ] represents the weight matrix.

3.5.3. Affine Transformation

This paper focuses on the rotation and translation properties of two-dimensional affine transformation [39]. After transformation, the original line is still a straight line and the original parallel line is still parallel. Affine transformation is a transformation of coordinates. Based on Figure 12, the derivation is written as follows:
Point P in the original coordinate system (black) is (x, y). A counterclockwise rotation of the point P is equivalent to clockwise rotation of the coordinate system. Thus, the point P in the black coordinate system is equivalent with the point P in the red coordinate system after the rotation. Based on this conclusion, we can determine the coordinates of the point P by simple stereo geometry, and then add the offset of the X-axis and the Y-axis based on this position; the formula is as follows:
{ x = x cos θ y sin θ + Δ x y = y cos θ + x sin θ + Δ y
where θ is the angle of rotation, Δ x is the offset on the X axis and Δ y is the offset on the Y axis; x and y are coordinate in the new coordinate system.

3.5.4. Calibration Transfer Method based on Affine Transformation

Based on the inputs and outputs { X m , y m } from the master instrument, and the inputs { X s } from the slave instrument, our task is to predict the unknown outputs { y ^ s } in the slave instrument. We assume that X m and X s are the spectra of two similar substances, and y m and y ^ s are in the same range. Due to the difference between two instruments, the observed spectral data are different. The observations from the perspective of the master instrument model are as follows:
{ y ^ m = F ( X m , β m ) = i = 1 A t i m q i m y ˜ s = F ( X s , β m ) = i = 1 A t ˜ i s q i m
where F is the linear prediction function, which is obtained by partial least squares in this paper; β m is the coefficient of the master model and y ^ m , t i m and q i m are the predicted values, the i-th column score vector and the loading vector, respectively. Accordingly, y ˜ s and t ˜ i s are the biased predicted values and the i-th biased column score vector for the slave instrument, respectively.
Therefore, the score vectors and predicted values both of the two instruments are different. As a result, there is a certain bias that needs to be corrected in the coefficient between the score vector and predicted values.
When correcting the bias, direct calculation will produce large errors. In order to solve this problem, we need to transform the score vectors and predicted values of the master and slave instrument into the range [0, 1] and thus keep the same scale between different values. The corresponding equations are given as follows:
{ t m - norm = ( t i m min ( t i m ) ) / ( max ( t i m ) min ( t i m ) ) y ^ m - norm = ( y ^ m min ( y ^ m ) ) / ( max ( y ^ m ) - min ( y ^ m ) ) t ˜ s - norm = ( t ˜ s min ( t ˜ i s ) ) / ( max ( t ˜ i s ) min ( t ˜ i s ) ) y ˜ s - norm = ( y ˜ s min ( y ˜ s ) ) / ( max ( y ˜ s ) - min ( y ˜ s ) )
where t i m - norm and y ^ m - norm are the normalized score vector and the predicted values of the master instrument, respectively; t ˜ i s - norm and y ˜ s - norm are the normalized and biased score vector and predicted values, respectively.
Two linear regression equations between score vector and predicted values are as follows:
{ y ^ m - norm = t i m - norm tan θ i m + b i m y ˜ s - norm = t ˜ i s - norm tan θ ˜ i s + b ˜ i s
where tan θ i m and tan θ ˜ i s are the regression coefficients (slopes) computed on the two instrument; b i m and b ˜ i s are the intercepts.
In order to more intuitively reflect the difference between two instruments, it can be better understood from Figure 13. The blue line is the regression coefficient between the score vector and predicted values. The black and red coordinate systems are the observations of the master and slave instrument, and there is a difference from different observations.
The unknown angles and biases between two instruments are solved as follows:
Firstly, the regression coefficient β m , the weight W m and loading P m matrix of PLS are obtained.
Secondly, a linear regression both of master and slave instrument is performed and slopes and intercepts are determined, respectively.
On the grounds of the PLS model, the score matrices and predicted values are calculated as shown below:
{ T m = X m W m ( P m W m ) - 1 ,       y ^ m = X m β m T ˜ s = X s W m ( P m W m ) - 1     ,         y ˜ s = X s β m
where T m and T ˜ s represent the score matrices of two instruments.
The score matrix T m , predicted values y ^ m , the score matrix T ˜ s and predicted values y ˜ s are pre-processed using Equation (9).
According to score vector of each column and predicted values, the least square is used to compute the corresponding slopes and intercepts, respectively. The equations are as follows:
{ min θ i m , b i m y ^ m - norm T aug m * [ tan θ i m b i m ]       min θ ˜ i s , b ˜ i s y ˜ s - norm T aug s * [ tan θ ˜ i s b ˜ i s ] 2 2
where T aug m is an augmented matrix [ t i m - norm , 1 ] ; T aug s is an augmented matrix [ t ˜ i s - norm , 1 ] ; 1 is the column vector with all ones.
Finally, the angle and biases between the two instruments are obtained. The equations for calculating the angles and biases are as follows:
{ Δ θ i = θ i m θ ˜ i s Δ b i = b i m b ˜ i s
where Δ θ i is the angle of the two coefficients; Δ b i is the corresponding bias.
The score matrix and predicted values of the test set are extracted by Equation (11).
The angles and biases obtained by Equation (13) are brought into the affine transformation to correct the predicted values. Since the rotation angle is relative to the origin of the coordinate, each sample needs to be adjusted before rotation. The equation is shown as follows:
U ^ i = U ˜ i M i
where the matrix M i = [ λ t cos Δ θ i λ t sin Δ θ i 0 λ y sin Δ θ i λ y cos Δ θ i 0 0 b i m 1 ] , U ˜ i = [ t ˜ i s - test , y ˜ s - test , 1 ] and U ^ i = [ t ^ i s - test , y ^ s - test , 1 ] . In addition, λ t = [ ( t ˜ i s - test min ( t ˜ i s ) ) / ( max ( t ˜ i s ) min ( t ˜ i s ) ) + min ( t ˜ i s ) ] × ( max ( t ˜ i s ) min ( t ˜ i s ) ) and λ y = [ ( y ˜ i s - test min ( y ˜ s ) ) / ( max ( y ˜ s ) min ( y ˜ s ) ) + min ( y ˜ s ) ] × ( max ( y ˜ s ) min ( y ˜ s ) ) represent the corresponding scaling factors for feature vector and predicted values, respectively; t ˜ i s - test and y ˜ s - test are biased score vector and predicted values of the test set, respectively; y ^ s - test is corrected predicted values; t ^ i s - test is corrected score vector.
Each column score vector and predicted values are solved separately, and a prediction matrix is obtained. The mean of the prediction matrix is the final predicted values.
Therefore, according to the expansion of the predicted values, β s is as follows:
β s = ( ( X s ) T X s ) 1 ( X s ) T ( i A ( t ˜ i s - test * λ t * sin Δ θ i + ( y ˜ s - test b ˜ i s * 1 ) * λ y * cos Δ θ i + b i m ) / A )

3.5.5. Summary of CTAI

Given calibration set of the master ( X cal m , y cal m ) , calibration set of the slave X cal s and test set ( X test s , y test s ) .
  • The PLS model is built on the calibration set ( X cal m , y cal m ) and the coefficient β m ; the weight matrix W m and the loading matrix P m can be obtained.
  • Modeling of affine transformation; it consists of the two datasets ( X cal m , y cal m ) and X cal s .
    (a)
    Computing ( T cal m , y ^ cal m ) and ( T ˜ cal s , y ˜ cal s ) of master and slave instrument by Equation (11).
    (b)
    ( T cal m , y ^ cal m ) and ( T ˜ cal s , y ˜ cal s ) are normalized separately by Equation (9).
    (c)
    ( tan θ i m , b i m ) and ( tan θ ˜ i s , b ˜ i s ) are calculated by Equation (12).
    (d)
    Computing Δ θ i angle and Δ b i bias between master and slave instrument by Equation (13).
  • Prediction.
    (a)
    ( T ˜ test s , y ˜ test s ) is obtained by Equation (11).
    (a)
    The matrix M i is introduced to correct predicted values by Equation (14).
    (c)
    The corrected prediction values are accumulated. The mean values are the last result.

4. Conclusions

In this study, the relationship of regression coefficients between feature vector and predicted values on different instruments was investigated and CTAI was proposed for calibration transfer based on affine invariance without transfer standards (CTAI). Based on the PLS model of the master instrument, the score matrix and the predicted values of the master spectra, the pseudo score matrix and the pseudo predicted values of the slave spectra are obtained. Then, angles and biases between the coefficients of the master instrument and the corresponding coefficients of the slave instrument are computed. Finally, new samples are corrected by affine transformation. Different transfer methods are tested with two NIR datasets, CTAI achieves the lowest RMSEP and standard error, and the results of statistical difference indicate that CTAI is generally better than other methods, which proves that CTAI is successfully used to correct the difference on different instruments. Hence, the proposed method may provide an efficient way for calibration transfer when standard samples are unavailable in practical applications.

Author Contributions

Conceptualization, Y.Z. and P.S.; methodology, Y.Z. and P.S.; software, Z.Z.; validation, Z.Z.; formal analysis, P.S., Y.Z., Z.Z. and J.Y.; data curation S.G.; writing—original draft preparation, Z.Z.; writing—review and editing, J.Y., P.S. and S.P.; visualization, S.P.; supervision, P.S. and Y.Z.; project administration, Z.Z.; funding acquisition, Y.Z.

Funding

This research was funded by National Natural Science Foundation of China (Grant no. 61601104), Natural Science Foundation of Hebei Province (Grant no. F2017501052) and the Basic Science Research Fund of Northeast University at Qin Huang Dao (Grant no. XNB201611).

Conflicts of Interest

No conflict of interest exits in the submission of this manuscript, and the manuscript is approved by all authors for publication. I would like to declare on behalf of my co-authors that the work described was original research that has not been published previously, and not under consideration for publication elsewhere, in whole or in part. All the authors listed have approved the manuscript that is enclosed.

References

  1. Huang, H.; Yu, H.; Xu, H.; Ying, Y. Near infrared spectroscopy for on/in-line monitoring of quality in foods and beverages: A review. J. Food Eng. 2008, 87, 303–313. [Google Scholar] [CrossRef]
  2. Roggo, Y.; Chalus, P.; Maurer, L.; Lema-Martinez, C.; Edmond, A.; Jent, N. A review of near infrared spectroscopy and chemometrics in pharmaceutical technologies. J. Pharm. Biomed. Anal. 2007, 44, 683–700. [Google Scholar] [CrossRef]
  3. Martinez, J.C.; Guzmán-Sepúlveda, J.R.; Bolañoz Evia, G.R.; Córdova, T.; Guzmán-Cabrera, R. Enhanced Quality Control in Pharmaceutical Applications by Combining Raman Spectroscopy and Machine Learning Techniques. Int. J. Thermophys. 2018, 39, 79. [Google Scholar] [CrossRef]
  4. Porep, J.U.; Kammerer, D.R.; Carle, R. On-line application of near infrared (NIR) spectroscopy in food production. Trends Food Sci. Tech. 2015, 46, 211–230. [Google Scholar] [CrossRef]
  5. Geladi, P.; Esbensen, K. Regression on multivariate images: Principal component regression for modeling, prediction and visual diagnostic tools. J. Chemom. 1991, 5, 97–111. [Google Scholar] [CrossRef]
  6. Næs, T.; Martens, H. Principal component regression in NIR analysis: View-points, background details and selection of components. J. Chemom. 1988, 2, 155–167. [Google Scholar] [CrossRef]
  7. Wold, S.; Sjöström, M.; Eriksson, L. PLS-regression: A basic tool of chemometrics. Chemom. Intell. Lab. Syst. 2001, 58, 109–130. [Google Scholar] [CrossRef]
  8. Sijmen, D.J. SIMPLS: An alternative approach to partial least squares regression. Chemom. Intell. Lab. Syst. 1993, 18, 251–263. [Google Scholar]
  9. Geladi, P.; Kowalski, B.R. Partial least-squares regression: A tutorial. Anal. Chim. Acta 1986, 185, 1–17. [Google Scholar] [CrossRef]
  10. Matthew, B.; Rayens, W. Partial least squares for discrimination. J. Chemometrics. 2012, 30, 446–452. [Google Scholar]
  11. Workman, J.J. A Review of Calibration Transfer Practices and Instrument Differences in Spectroscopy. Appl. Spectrosc. 2018, 72, 340–365. [Google Scholar] [CrossRef]
  12. Bouveresse, E.; Hartmann, C.; Massart, D.L.; Last, I.R.; Prebble, K.A. Standardization of near-infrared spectrometric instruments. Anal. Chem. 1996, 68, 982–990. [Google Scholar] [CrossRef]
  13. Feudale, R.N.; Woody, N.A.; Tan, H.; Myles, A.J.; Brown, S.D.; Ferré, J. Transfer of multivariate calibration models: A review. Chemom. Intell. Lab. Syst. 2002, 64, 181–192. [Google Scholar] [CrossRef]
  14. Wang, Y.; Veltkamp, D.J.; Kowalski, B.R. Multivariate instrument standardization. Anal. Chem. 1991, 63, 2750–2756. [Google Scholar] [CrossRef]
  15. Wang, Y.; Michael, J.L.; Kowalski, B.R. Improvement of multivariate calibration through instrument standardization. Anal. Chem. 1992, 64, 562–564. [Google Scholar] [CrossRef]
  16. Bouveresse, E.; Massart, D. Improvement of the piecewise direct standardisation procedure for the transfer of NIR spectra for multivariate calibration. Chemom. Intell. Lab. Syst. 1996, 32, 201–213. [Google Scholar] [CrossRef]
  17. Wang, Z.; Thomas, D.; Kowalski, B.R. Additive background correction in multivariate instrument standardization. Anal. Chem. 1995, 67, 2379–2385. [Google Scholar] [CrossRef]
  18. Tan, H.-W.; Brown, S.D. Wavelet hybrid direct standardization of near-infrared multivariate calibrations. J. Chemometrics. 2001, 15, 647–663. [Google Scholar] [CrossRef]
  19. Fan, W.; Liang, Y.; Yuan, D.; Wang, J. Calibration model transfer for near-infrared spectra based on canonical correlation analysis. Anal. Chim. Acta 2008, 623, 22–29. [Google Scholar] [CrossRef] [PubMed]
  20. Zheng, K.; Zhang, X.; Iqbal, J.; Fan, W.; Wu, Ti.; Du, Y.; Liang, Y. Calibration transfer of near-infrared spectra for extraction of informative components from spectra with canonical correlation analysis. J. Chemometrics. 2014, 28, 773–784. [Google Scholar] [CrossRef]
  21. Melzer, T.; Reiter, M.; Bischof, H. Appearance models based on kernel canonical correlation analysis. Pattern Recognit. 2003, 36, 1961–1971. [Google Scholar] [CrossRef]
  22. Leng, L.; Zhang, T.; Kleinman, L.; Zhu, W. Ordinary least square regression, orthogonal regression, geometric mean regression and their applications in aerosol science. J. Phys. Conf. Ser. 2007, 78, 012084. [Google Scholar] [CrossRef] [Green Version]
  23. Donald, C.; Orcutt, G.H. Application of least squares regression to relationships containing auto-correlated error terms. J. Amer. Stat. Assoc. 1949, 44, 32–61. [Google Scholar]
  24. Chen, W.-R.; Bin, J.; Lu, H.-M.; Zhang, Z.-M.; Liang, Y.-Z. Calibration transfer via an extreme learning machine auto-encoder. Analyst 2016, 141, 1973–1980. [Google Scholar] [CrossRef] [PubMed]
  25. Wise, B.M.; Martens, H.; Høy, M.; Bro, R.; Brockhoff, P.B. Calibration Transfer by Generalized Least Squares. In Proceedings of the Seventh Scandinavian Symposium on Chemometrics (SSC7), Copenhagen, Denmark, 19–23 August 2001. [Google Scholar]
  26. Du, W.; Chen, Z.-P.; Zhong, L.-J.; Wang, S.-X.; Yu, R.-Q.; Nordon, A.; Littlejohn, D.; Holden, M. Maintaining the predictive abilities of multivariate calibration models by spectral space transformation. Anal. Chim. Acta. 2011, 690, 64–70. [Google Scholar] [CrossRef] [PubMed]
  27. Chen, Z.P.; Li, L.M.; Yu, R.Q.; Littlejohn, D.; Nordon, A.; Morris, J.; Dann, A.S.; Jeffkins, P.A.; Richardson, M.D.; Stimpson, S.L. Systematic prediction error correction: A novel strategy for maintaining the predictive abilities of multivariate calibration models. Analyst. 2010, 136, 98–106. [Google Scholar] [CrossRef] [PubMed]
  28. Kramer, K.E.; Morris, R.E.; Rose-Pehrsson, S.L. Comparison of two multiplicative signal correction strategies for calibration transfer without standards. Chemom. Intell. Lab. Syst. 2008, 92, 33–43. [Google Scholar] [CrossRef]
  29. Preisner, O.; Lopes, J.A.; Guiomar, R.; Machado, J.; José, C.M. Fourier transform infrared (FT-IR) spectroscopy in bacteriology: towards a reference method for bacteria discrimination. Anal. Bioanal. Chem. 2007, 387, 1739–1748. [Google Scholar] [CrossRef] [PubMed]
  30. Isaksson, T.; Næs, T. The effect of multiplicative scatter correction (MSC) and linearity improvement in NIR spectroscopy. Appl. Spectrosc. 1988, 42, 1273–1284. [Google Scholar] [CrossRef]
  31. Blank, T.B.; Sum, S.T.; Brown, S.D.; Monfre, S.L. Transfer of near-infrared multivariate calibrations without standards. Anal. Chem. 1996, 68, 2987–2995. [Google Scholar] [CrossRef]
  32. Wold, S.; Antti, H.; Lindgren, F.; Öhman, J. Orthogonal signal correction of near infrared spectra. Chemom. Intell. Lab. Syst. 1998, 44, 175–185. [Google Scholar] [CrossRef]
  33. Sjöblom, J.; Svensson, O.; Josefson, M.; Kullberg, H.; Wold, S. An evaluation of orthogonal signal correction applied to calibration transfer of near infrared spectra. Chemom. Intell. Lab. Syst. 1998, 44, 229–244. [Google Scholar] [CrossRef]
  34. Malli, B.; Birlutiu, A.; Natschläger, T. Standard-free calibration transfer-An evaluation of different techniques. Chemom. Intell. Lab. Syst. 2017, 161, 49–60. [Google Scholar] [CrossRef]
  35. Pan, S.J.; Tsang, I.; Kwok, J.; Yang, Q. Domain adaptation via transfer component analysis. IEEE Trans. Neural Netw. 2011, 22, 199–210. [Google Scholar] [CrossRef] [PubMed]
  36. Schölkopf, B.; Smola, A.; Müller, K.R. Kernel principal component analysis. In Artificial Neural Networks — ICANN’97, Proceeding of 7th International Conference Lausanne, Lausanne, Switzerland, 8–10 October 1997; Springer: Berlin/Heidelberg, Germany, 1997; pp. 583–588. [Google Scholar]
  37. Schölkopf, B.; Smola, A.; Müller, K.-R. Nonlinear component analysis as a kernel eigenvalue problem. Neural Comput. 1998, 10, 1299–1319. [Google Scholar] [CrossRef]
  38. Muandet, K.; Balduzzi, D.; Schölkopf, B. Domain generalization via invariant feature representation. In Proceedings of the 30th International Conference on Machine Learning (ICML-13), Atlanta, GA, USA, 16–21 June 2013; pp. 10–18. [Google Scholar]
  39. Bloomenthal, J.; Jon, R. Homogeneous coordinates. Visual Computer. 1994, 11, 15–26. [Google Scholar] [CrossRef] [Green Version]
Sample Availability: Samples are not available from the authors.
Figure 1. Data setting of the calibration transfer method based on affine invariance without transfer standards (CTAI). We assume the data to be available in (A), and the data after being processed based on PLS model of the master instrument is shown in (B).
Figure 1. Data setting of the calibration transfer method based on affine invariance without transfer standards (CTAI). We assume the data to be available in (A), and the data after being processed based on PLS model of the master instrument is shown in (B).
Molecules 24 01802 g001
Figure 2. The relationship between the uncorrected and the corrected predict values for corn dataset by (A) moisture, (B) oil, (C) protein and (D) starch. The blue and red dots represent the uncorrected and the corrected predicted results for each sample, respectively.
Figure 2. The relationship between the uncorrected and the corrected predict values for corn dataset by (A) moisture, (B) oil, (C) protein and (D) starch. The blue and red dots represent the uncorrected and the corrected predicted results for each sample, respectively.
Molecules 24 01802 g002
Figure 3. Moisture content predicted for corn dataset as determined by (A) CTAI, (B) MSC, (C) TCR, (D) CCA, (E) SBC and (F) PDS. The blue and red dots represent the results for each sample in the train set and test set, respectively.
Figure 3. Moisture content predicted for corn dataset as determined by (A) CTAI, (B) MSC, (C) TCR, (D) CCA, (E) SBC and (F) PDS. The blue and red dots represent the results for each sample in the train set and test set, respectively.
Molecules 24 01802 g003
Figure 4. Oil content predicted for corn dataset as determined by (A) CTAI, (B) MSC, (C) TCR, (D) CCA, (E) SBC and (F) PDS. The blue and red dots represent the results for each sample in the train set and test set, respectively.
Figure 4. Oil content predicted for corn dataset as determined by (A) CTAI, (B) MSC, (C) TCR, (D) CCA, (E) SBC and (F) PDS. The blue and red dots represent the results for each sample in the train set and test set, respectively.
Molecules 24 01802 g004
Figure 5. Protein content predicted for corn dataset as determined by (A) CTAI, (B) MSC, (C) TCR, (D) CCA, (E) SBC and (F) PDS. The blue and red dots represent the results for each sample in the train set and test set, respectively.
Figure 5. Protein content predicted for corn dataset as determined by (A) CTAI, (B) MSC, (C) TCR, (D) CCA, (E) SBC and (F) PDS. The blue and red dots represent the results for each sample in the train set and test set, respectively.
Molecules 24 01802 g005
Figure 6. Starch content predicted for corn dataset as determined by (A) CTAI, (B) MSC, (C) TCR, (D) CCA, (E) SBC and (F) PDS. The blue and red dots represent the results for each sample in the train set and test set, respectively.
Figure 6. Starch content predicted for corn dataset as determined by (A) CTAI, (B) MSC, (C) TCR, (D) CCA, (E) SBC and (F) PDS. The blue and red dots represent the results for each sample in the train set and test set, respectively.
Molecules 24 01802 g006
Figure 7. The relationship between the uncorrected and the corrected predict values for wheat dataset by (A) B1-B2, (B) B1-B3 and (C) B3-B2. The blue and red dots represent the uncorrected and the corrected predicted results for each sample, respectively.
Figure 7. The relationship between the uncorrected and the corrected predict values for wheat dataset by (A) B1-B2, (B) B1-B3 and (C) B3-B2. The blue and red dots represent the uncorrected and the corrected predicted results for each sample, respectively.
Molecules 24 01802 g007
Figure 8. Protein content predicted between instruments B1 and B2 for wheat dataset as determined by (A) CTAI, (B) MSC, (C) TCR, (D) CCA, (E) SBC and (F) PDS. The blue and red dots represent the results for each sample in the train set and test set, respectively.
Figure 8. Protein content predicted between instruments B1 and B2 for wheat dataset as determined by (A) CTAI, (B) MSC, (C) TCR, (D) CCA, (E) SBC and (F) PDS. The blue and red dots represent the results for each sample in the train set and test set, respectively.
Molecules 24 01802 g008
Figure 9. Protein content predicted between instruments B1 and B3 for wheat dataset as determined by (A) CTAI, (B) MSC, (C) TCR, (D) CCA, (E) SBC and (F) PDS. The blue and red dots represent the results for each sample in the train set and test set, respectively.
Figure 9. Protein content predicted between instruments B1 and B3 for wheat dataset as determined by (A) CTAI, (B) MSC, (C) TCR, (D) CCA, (E) SBC and (F) PDS. The blue and red dots represent the results for each sample in the train set and test set, respectively.
Molecules 24 01802 g009
Figure 10. Protein content predicted between instruments B3 and B2 for wheat dataset as determined by (A) CTAI, (B) MSC, (C) TCR, (D) CCA, (E) SBC and (F) PDS. The blue and red dots represent the results for each sample in the train set and test set, respectively.
Figure 10. Protein content predicted between instruments B3 and B2 for wheat dataset as determined by (A) CTAI, (B) MSC, (C) TCR, (D) CCA, (E) SBC and (F) PDS. The blue and red dots represent the results for each sample in the train set and test set, respectively.
Molecules 24 01802 g010
Figure 11. (A) Spectral differences between m5 and mp6 of corn samples; (B) spectral differences between B1 and B2 of wheat samples; (C) spectral differences between B1 and B3 of wheat samples; (D) spectral differences between B2 and B3 of wheat samples.
Figure 11. (A) Spectral differences between m5 and mp6 of corn samples; (B) spectral differences between B1 and B2 of wheat samples; (C) spectral differences between B1 and B3 of wheat samples; (D) spectral differences between B2 and B3 of wheat samples.
Molecules 24 01802 g011
Figure 12. Derivation of affine transformation. In the coordinate system, the counterclockwise rotation of P is equivalent to the clockwise rotation of the coordinate system.
Figure 12. Derivation of affine transformation. In the coordinate system, the counterclockwise rotation of P is equivalent to the clockwise rotation of the coordinate system.
Molecules 24 01802 g012
Figure 13. Theory of CTAI. tan θ is the coefficient between the feature vector and the predicted values. The angles and deviations observed under different instruments are different. We correct the predicted value of the slave instrument with the rotation and translation of affine transformation.
Figure 13. Theory of CTAI. tan θ is the coefficient between the feature vector and the predicted values. The angles and deviations observed under different instruments are different. We correct the predicted value of the slave instrument with the rotation and translation of affine transformation.
Molecules 24 01802 g013
Table 1. Summary of the partial least squares (PLS) models and properties.
Table 1. Summary of the partial least squares (PLS) models and properties.
InstrumentReference ValuesRMSECmRMSEPmRMSECVmin (LV)Biasmrmpm
m5specmoisture0.005990.007640.01066(14)0.00080.999732.6 × 10−24
m5specoil0.026860.056640.05049(15)−0.013270.93321.3 × 10−7
m5specprotein0.05070.100660.11012(15)0.028140.976321 × 10−10
m5specstarch0.095390.189930.19227(15)0.017890.974641.6 × 10−10
mp6specmoisture0.099910.156370.14775(10)−0.026780.920834.2 × 10−7
mp6specoil0.060520.090980.09872(12)0.018680.876978.2 × 10−6
mp6specprotein0.101010.133380.15043(12)0.021280.966591.1 × 10−9
mp6specstarch0.276360.267230.35978(9)0.021240.931361.6 × 10−7
B1protein0.32880.332540.50337(15)0.009060.985082.3 × 10−38
B2protein0.216360.837550.32441(15)−0.131240.84857.2 × 10−15
B3protein0.302880.515670.43896(15)−0.0340.960093.2 × 10−28
RMSECm: Root Mean Square Error of calibration set; RMSEPm: Root Mean Square Error of test set; RMSECVmin: Minimum Root Mean Square Error of Cross-Validation; LV: The optimal number of latent variables is selected only with the lowest RMSECV; rm: Pearson correlation coefficient for predicted vs. actual values; pm: p values corresponding to the Pearson correlation coefficient is obtained by test.
Table 2. Summary of the relevant results between uncorrected and CTAI corrected.
Table 2. Summary of the relevant results between uncorrected and CTAI corrected.
Instrument Reference Valuesm5spec*-mp6specB1*-B2B1*-B3B3*-B2
MoistureOilProteinStarchProtein
y ^ m vs   y ˜ s RMSEPupre1.607050.79892.067972.117430.698942.925411.23368
RMSEPpre0.212550.069220.131950.333580.315370.626320.65398
kpre0.64980.771290.945530.825270.888090.762900.86909
rpre0.816440.895980.962860.921970.975940.876950.93715
ppre1.1 × 10−42.6 × 10−62.3 × 10−93.8 × 10−72 × 10−336.8 × 10−171.3 × 10−23
tpre−15.42919.335−19.1478.8382.29210.684-3.826
y m vs   y ˜ s RMSEPu1.607620.815322.096652.102910.719772.900111.08008
RMSEP0.210950.082330.166140.347140.414190.682150.38446
k0.651910.532970.987360.793290.968980.856930.93896
r0.819220.788580.958440.914870.967700.895170.97796
p1.0 × 10−42.8 × 10−45.1 × 10−96.9 × 10−72.2 × 10−301.8 × 10−182.5 × 10−34
t−15.43719.657−19.4088.7622.25610.649−3.701
tcritical_value2.1312.1312.1312.1312.012.012.01
*: The master instrument; RMSEPupre: RMSEP of uncorrected slave instrument relative to primary instrument prediction; RMSEPpre: RMSEP of CTAI corrected slave instrument relative to primary instrument prediction; kpre: The slope between predicted values of uncorrected slave instrument and primary prediction; rpre: Correlation coefficient of uncorrected slave prediction relative to master prediction; ppre: p values corresponding to the Pearson correlation coefficient are obtained by test; tpre: The result of One-Sample t-Test between uncorrected slave prediction and master prediction; RMSEPu: RMSEP of uncorrected slave instrument relative to primary actual values; RMSEP: RMSEP of CTAI corrected slave instrument relative to primary actual values; k: The slope between predicted values of uncorrected slave instrument and primary actual values; r: Pearson correlation coefficient of uncorrected slave prediction relative to primary actual values; p: p values corresponding to the Pearson correlation coefficient are obtained by test; t: The result of One-Sample t-Test between uncorrected slave prediction and master actual values; tcritical_value: The t critical value for n–1 degrees of freedom at the significance level alpha = 0.05.
Table 3. Summary of Root Mean Square Error of test set (RMSEP) and Root Mean Square Error of calibration set (RMSEC) of different methods. The m5spec was used as the master spectra, and the mp6spec was used as the secondary spectra for corn dataset. The protein content was chosen as the property for wheat dataset.
Table 3. Summary of Root Mean Square Error of test set (RMSEP) and Root Mean Square Error of calibration set (RMSEC) of different methods. The m5spec was used as the master spectra, and the mp6spec was used as the secondary spectra for corn dataset. The protein content was chosen as the property for wheat dataset.
MethodCTAIMSCTCRCCASBCPDS
moistureRMSEC0.226461.928390.618730.15996(14a)0.18506(5a)0.14742(17a)
RMSEP0.210951.66890.390660.23304(14a)0.42574(5a)0.24238(17a)
oilRMSEC0.081411.216470.145430.15764(6a)0.08423(23a)0.10794(28a)
RMSEP0.082331.232090.142250.11432(6a)0.08361(23a)0.09495(28a)
proteinRMSEC0.172471.772940.282970.27860(14a)0.17422(6a)0.24662(23a)
RMSEP0.166141.800870.352230.39535(14a)0.19101(6a)0.28193(23a)
starchRMSEC0.395171.891651.210930.33937(10a)0.38426(23a)0.62099(23a)
RMSEP0.347141.931290.798520.85704(10a)0.36969(23a)0.78977(23a)
B1*-B2RMSEC0.556821.311530.992461.11889(5a)0.48509(6a)1.3676(7a)
RMSEP0.414190.921940.868812.68469(5a)0.4677(6a)4.09019(7a)
B1*-B3RMSEC0.818952.916950.846820.68529(15a)1.00007(8a)0.57858(5a)
RMSEP0.682152.405870.729961.10564(15a)0.79294(8a)1.33547(5a)
B3*-B2RMSEC0.547531.250960.769721.57073(14a)0.56236(5a)2.1039(8a)
RMSEP0.384461.384680.636892.29856(14a)0.53534(5a)1.83564(8a)
a: Number of standard samples; the number of samples for slave instrument with labels is 20 in TCR.
Table 4. RMSEP comparison of CTAI and other methods, RMSEP improvements and p values by the Wilcoxon signed rank test (α = 0.05). The m5spec was used as the master spectra, and the mp6spec was used as the secondary spectra for corn dataset. The protein content was chosen as the property for wheat dataset.
Table 4. RMSEP comparison of CTAI and other methods, RMSEP improvements and p values by the Wilcoxon signed rank test (α = 0.05). The m5spec was used as the master spectra, and the mp6spec was used as the secondary spectra for corn dataset. The protein content was chosen as the property for wheat dataset.
MSCTCRCCASBCPDS
h(%)ph(%)ph(%)ph(%)ph(%)p
moisture87.354.3 × 10 −4460.539.480.4350.450.0112.960.04
oil93.314.3 × 10 −442.120.0127.980.321.520.2313.280.46
protein90.774.3 × 10 −452.830.0957.970.0313.020.2341.060.01
starch82.024.3 × 10 −456.520.2359.490.836.090.0256.040.75
B1*-B255.070.1152.320.7984.575.3 × 10 −911.442.6 × 10 −989.879.2 × 10 −3
B1*-B371.647.5 × 10 −106.550.1138.31.8 × 10 −513.971 × 10 −548.929.8 × 10 −5
B3*-B272.233.1 × 10−939.634.6 × 10−383.270.0228.187.5 × 10 −1079.050.06

Share and Cite

MDPI and ACS Style

Zhao, Y.; Zhao, Z.; Shan, P.; Peng, S.; Yu, J.; Gao, S. Calibration Transfer Based on Affine Invariance for NIR without Transfer Standards. Molecules 2019, 24, 1802. https://doi.org/10.3390/molecules24091802

AMA Style

Zhao Y, Zhao Z, Shan P, Peng S, Yu J, Gao S. Calibration Transfer Based on Affine Invariance for NIR without Transfer Standards. Molecules. 2019; 24(9):1802. https://doi.org/10.3390/molecules24091802

Chicago/Turabian Style

Zhao, Yuhui, Ziheng Zhao, Peng Shan, Silong Peng, Jinlong Yu, and Shuli Gao. 2019. "Calibration Transfer Based on Affine Invariance for NIR without Transfer Standards" Molecules 24, no. 9: 1802. https://doi.org/10.3390/molecules24091802

APA Style

Zhao, Y., Zhao, Z., Shan, P., Peng, S., Yu, J., & Gao, S. (2019). Calibration Transfer Based on Affine Invariance for NIR without Transfer Standards. Molecules, 24(9), 1802. https://doi.org/10.3390/molecules24091802

Article Metrics

Back to TopTop