Mathematical Problems in Engineering
Volume 2015 (2015), Article ID 761280, 9 pages
http://dx.doi.org/10.1155/2015/761280

Research Article

Recursive Gaussian Process Regression Model for Adaptive Quality Monitoring in Batch Processes

Le Zhou,¹ Junghui Chen,² and Zhihuan Song¹

¹State Key Laboratory of Industrial Control Technology, Zhejiang University, Hangzhou, Zhejiang 310027, China
²Department of Chemical Engineering, Chung-Yuan Christian University, Chung-Li 320, Taiwan

Received 5 November 2014; Revised 23 December 2014; Accepted 23 December 2014

Academic Editor: Gang Li

Copyright © 2015 Le Zhou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

In chemical batch processes with slow responses and a long duration, it is time-consuming and expensive to obtain sufficient normal data for statistical analysis. With the persistent accumulation of the newly evolving data, the modelling becomes adequate gradually and the subsequent batches will change slightly owing to the slow time-varying behavior. To efficiently make use of the small amount of initial data and the newly evolving data sets, an adaptive monitoring scheme based on the recursive Gaussian process (RGP) model is designed in this paper. Based on the initial data, a Gaussian process model and the corresponding SPE statistic are constructed at first. When the new batches of data are included, a strategy based on the RGP model is used to choose the proper data for model updating. The performance of the proposed method is finally demonstrated by a penicillin fermentation batch process and the result indicates that the proposed monitoring scheme is effective for adaptive modelling and online monitoring.

1. Introduction

In modern industries, batch and semibatch processes are of great importance for the production of high-quality and value-added specialty chemicals, such as semiconductors, pharmaceuticals, and biological products. For product consistency and quality improvement, multivariate statistical process monitoring (MSPM) has been developed and widely used in many batch processes [1–11]. Nomikos and MacGregor [12, 13] are the first to apply multiway principle component analysis (MPCA) and multiway partial least squares (MPLS) to batch process monitoring. MPCA is powerful at analyzing the batch-to-batch variations and monitoring the newly evolving batches with the assumption that the batchwise unfolded data follow a multivariate Gaussian distribution. However, there are still some disadvantages of the conventional MPCA and its extensions. Firstly, MPCA assumes that the length of each batch is equal, which is unlikely to happen in practice. Furthermore, the entire batch data are needed for MPCA for online monitoring. Hence, the unknown future data of a newly evolving batch have to be estimated. Also, it is hard to reveal the timewise variations since the entire batch data are treated as a single object [14]. To solve these problems, variable-wise unfolding approaches were developed [15]. This method provides a straightforward scheme for online monitoring and the equal length of the collected batches is not required.

It is noticed that, however, both batchwise and variable-wise unfolding methods are based on the sufficient batch data. Abundant data are easy to be collected in the rapid manufacturing process, such as wafer etching, spray-drying, and spray-coating. Nevertheless, it is time-consuming and expensive to obtain data from slow chemical processes, such as emulsion polymerization, fermentation, and pharmaceutical and biotechnical products [15]. In this situation, it is not proper to construct monitoring models after all the sufficient data are collected. It is interesting to design a statistical scheme to construct an initial model using limited batch data and then adapt the strategy when newly evolving batch data arrive. Therefore, Zhao et al. [14, 16] proposed an adaptive monitoring model based on ICA [14]. In their article, multiphase ICA was constructed based on limited reference cycles and an adaptive algorithm was formulated with the accumulation of newly available normal batches. However, their adaptive strategy was complex and time-consuming since the adaptive criterion was not straightforward and the updated model needed to be retrained completely.

Process monitoring is crucial not only to process safety but also to quality improvement. When the quality variables are examined, an input-output model is needed to reveal the quality-relevant variations. Since the batch process data are highly nonlinear, dynamic, and seriously interconnected with uncertainties, it is hard to build the relationship between the process variables and the quality variables and detect the abnormal production conditions. Therefore, there has been an increasing interest in quality monitoring of batch processes [13, 16–19]. However, their research was also conducted based on the assumption that the normal batch data were sufficient, when the chemical batch processes are operated in a long batch time. Moreover, the process is slow-responding. Thus, only limited batch runs can be used for modelling. With the restricted batch data, the statistical process analysis would be inaccurate. Also, the time-varying behavior of the newly evolving batches may not be completely described based on the data recorded in the early stage of the operating batches. Hence, an adaptive strategy is important for batch quality monitoring, especially in slow processes.

In this paper, a recursive Gaussian process (RGP) model is designed for adaptive quality monitoring with limited initial batch data. A Gaussian process (GP) model is trained at first using limited batch data to construct the relationship between the process and the quality variables. In the GP model, the measure of the confidence level in the prediction is also taken into account. Hence, not only the mean value but also the variance of the quality prediction is supplied. With the estimation of the variance from GP as a level of confidence, the adaptive strategy is proposed to determine when and how the model should be updated. This updating strategy is straightforward and easy to implement. It is helpful to describe the time-varying dynamics over batches. Using the online updating strategy, an adaptive quality monitoring scheme is constructed.

The rest of the paper is organized as follows. The Gaussian process model is revisited in Section 2. Next, the batch process monitoring based on GP using limited batch data is discussed. The SPE statistic and its confidence bound are constructed then. With the accumulation of newly evolving batches, the updating strategy is developed and the implementation of the recursive GP based adaptive quality monitoring in the batch process is discussed in Section 3. Also, an industrial case is provided to demonstrate the efficiency of the proposed method in Section 4. Finally, conclusions are made.

2. GP Regression Model

In supervised learning, the objective is to infer a distribution over functions given the inputs and outputs . Usually, a parametric representation for the function is assumed. Then the model parameters are inferred instead of . The methods above are called parametric modelling. In nonparametric modelling, without limitations of the model structure, the distribution over functions themselves is inferred. In a GP model, Bayesian inference is used for estimating the distribution of functions.

Given inputs and outputs , a GP model defines a Gaussian prior distribution as [20]

where is the mean value of outputs and is the covariance function in which . In the GP model, the data collected are not assumed to follow the same Gaussian distribution any more. On the contrary, all the samples follow a joint Gaussian distribution. Hence, the joint Gaussian distribution of the training data and the test data can be estimated:

where and . According to the property of conditional Gaussian distribution, the posterior distribution of the test data has the following form:

It is noticed that it is common to use the mean function as since GP is flexible enough to model arbitrary mean value [20]. In the kernel function of GP, any positive definite kernel can be used for GP covariance. Generally, three common kernels, linear kernel, Gaussian kernel, and squared-exponential kernel, are mostly used [21]. In this paper, Gaussian kernel is employed, which is given by [22]

where controls the vertical length scales of the function and reflects the horizontal variations. To estimate the kernel parameters, which are called hyperparameters, the empirical Bayes approach is used. Firstly, the marginal log-likelihood is written as [20]

Then the log-likelihood is maximized with respect to each hyperparameter as [23]

In the GP model, it will take to compute and for hyperparameter gradient estimation. Particularly in the batch process model, the model needs to be updated with the newly evolving batches. With the accumulation of the sample size, the training of the GP model is time-consuming. In this research work, a RGP model needs to be constructed.

3. Adaptive Quality Monitoring Based on RGP in Batch Processes

For a slowly varying batch process, all the batch data cannot sufficiently be collected at the initial operation stage. If the monitoring system is not adapted with the newly evolving batches, it will cause nuisance alarms. An online adaptive strategy is proposed to determine when the model needs to be updated and how to adapt the model without retraining the GP model using all the data set. The update criterion is based on the estimation of the variance from GP and the SPE statistic to judge if the newly evolving data are needed for updating. Then the covariance matrix of GP is recursively updated to provide a better representative of the current state of the system.

3.1. GP Based Batch Process Monitoring with Limited Batches

Using limited batch data at first, the data are collected in the form of three-dimensional array. Assume the measured process variables are summarized as , where variables are collected at times and is the batch number. Simultaneously, the key quality variable is collected as . For simple explanation, only single quality variable is measured, but it is very easy to extend it to multiple quality variables. Also, in this paper, a variable-wise unfolding method is used for data preprocessing [15]. Using the unfolded matrix and , a GP model can be trained as an initial quality monitoring model to identify the early correlations between the process and the quality variables.

The prediction error of the GP model and the corresponding SPE statistic can be calculated as [24]

The confidence limit of SPE is approximated by distribution as , in which and are the parameters of distribution and they are estimated as [25]

The significance level of SPE statistics is 99%. When a new sample of the process variables is collected, the prediction value of quality can be estimated using (3) and the SPE statistic is calculated using (7). Generally, if the SPE statistic exceeds the control limit, it is thought that a fault occurs and vice versa. It is noticed that, however, the accuracy of the statistic is not evaluated because most of the time the collected data are not sufficient in the initial modeling stage. In the GP model, not only the mean value but also the variance of prediction of each sample is provided. With the variance serving as a measure of the confidence for the prediction, the required information can be found and included to enhance the accuracy of the monitoring model. Since the prediction follows a Gaussian distribution, the real value of quality mostly relies on a region covering of the prediction distribution around the mean. Owing to the uncertainty of prediction, the SPE statistic can also be estimated in a corresponding confidence region instead of a single value. According to the criterion, the upper and lower limit of mean can be estimated as . Then the prediction error and the SPE statistic of and are calculated using (7). Similarly, can also be constructed based on the prediction error between mean value and the real value. It is noticed that and are not the real bounds of the SPE statistic because of the quadratic property. However, it is easy to prove that the upper and lower limits can be determined by , , and :

In this way, the confidence bound of the SPE statistic is constructed to evaluate its reliability. If is below the control limit while is above it, the normal state this time becomes questionable. By contrast, in some states there may be false alarms if only is above the control limit while indicates it is normal. Thus, the confidence bound of the SPE statistic is helpful for online updating strategy for quality monitoring in batch processes.

3.2. Updating Strategy and RGP Based on Adaptive Quality Monitoring

After the initial monitoring model is built using the early collected normal batches, with the new batch data evolved, one needs to judge if the model should be adjusted to enhance the predictability of the monitoring model. Based on the discussions in the previous section, we have the following.

(1)If the upper bond of the test data is still below the control limit, there is statistically sufficient evidence to show that it is normal, no more new information is needed, and the original model is good enough without being updated.(2)If the lower bond of the test data is larger than the control limit, chances are that almost all the potential values of predictions are far from the real data. In this situation, it is highly doubted that an abnormal variation has occurred, and it is mainly caused by the unknown disturbances instead of the normal time-varying behavior.(3)If the control limit falls in the region of the statistic, in which it is larger than and smaller than , the cases become indistinct. The result indicates that part of the potential values has exceeded the control limited while others have not. Hence, the part of false alarm results from the failure of initial models and an updating model is needed to adapt the model to the state of current batches.

Figure 1 illustrates the updating strategy of our method and the relationship between and of a test data and the control limit. In Figure 1, the value of -coordinate is the control limit based on the mean prediction error of the normal data. Hence, when the control limit is smaller than , all the test data are defined as faults; when the control limit is larger than , it indicates that almost all the possible prediction values fall in the normal region and all the test data would be normal.

Figure 1: Relationship between the SPE statistics and the control limit.

When the test data are added for updating the model, a recursive method is used. Suppose inputs and outputs are used for training the initial GP model; the new data are added for model updating. Then the new data follow a new joint Gaussian distribution as

where . The updated covariance can be directly calculated using data which are computationally demanded. In this paper, a recursive method is employed using Sherman-Morrison-Woodbury formula to estimate by the covariance of the first samples [26]. The updated covariance is calculated as

where and . Hence, when the new data need to be added, the recursive formula for can be used. After the model has been updated, the corresponding monitoring system can be set up as GP.

To sum up, the procedures of adaptive quality monitoring based on RGP are listed as follows.

Initial Phase

(1)Collect the historical normal batch data even if there are a limited number of operation batches.(2)Unfold the data using variable-wise unfolding and train an initial GP model.(3)Calculate the SPE statistic and the control limit based on the currently available data. Estimate the data regions using (7) and (9).

Adaptive Phase

(4)Calculate the SPE statistic for new batch data one time point by one time point, and judge if each sampling data point needs to be included for updating the model. If yes, go to Step ; otherwise, go to Step .(5)Update the covariance using (11) and calculate the new SPE statistic and the region of the control limit based on new and in the updated model.

Monitoring Phase

(6)Check if the current SPE statistic exceeds the control limit.

For batch-size updating, the region of control limits and is estimated based on all the useful data after one batch run. Then, the active data are chosen in a newly evolving batch to update the model. For sample-size updating, during the batch run, the region of control limit will be changed as soon as the model is updated. Compared with the two updating strategies, the selected active data in the batch-size updating may be similar and redundant. Thus, in this paper, sample-size updating is applied since it is more time-saving.

Case Study. In this section, a fed-batch penicillin benchmark is illustrated to evaluate the performance of the proposed method. In a typical operating process for penicillin fermentation, the microorganisms are generated in certain circumstances during the initial stage. When the substrate has been consumed by the microorganisms, the penicillin is produced as the metabolite. The concentration of microorganisms reaches the stationary phase. During the penicillin production stage, glucose is fed continuously to achieve a high product formation rate due to the catabolized repressor effect [27, 28]. The batch data are carried out using a simulator (PenSim v2.0) developed by the monitoring and control group at the Illinois Institute of Technology (http://www.chee.iit.edu/~cinar). The flow diagram of the penicillin fermentation process is given in Figure 2.

Figure 2: Flow diagram of the penicillin fermentation process.

In this paper, the duration of every batch is 400 hours and the sampling intervals of the process variables are set to be 1 hour. For modelling, a total of 11 key process variables are selected. They are tabulated in Table 1. The concentration of penicillin is chosen as the quality variable. It is noticed that the sampling interval of quality cannot be 1 hour in reality as the concentration needs to be examined offline at the lab generally. Therefore, to mimic the real situation, it is assumed that the quality data are recorded every 20 h. To build the input-output relationship using the unequal size of the process variables and the quality variables, some research has been done [29]. Since it is not the main concern of this paper, the simple downsampling method is used here to equalize the sampling frequency of inputs and outputs. Besides, the measured variables are added by Gaussian noises () to simulate the real process environment. Forty normal batches are generated to demonstrate the proposed method. To simulate the slow time-varying behavior among batches, the culture volume of these normal batches is set to change gradually from 100 L to 120 L. Although the culture volume has changed, it is still considered that these batches are normal and the variations over batches can be tolerated. For constructing initial model, 5 batches data are first used. According to the downsampling method, a GP model can be constructed based on the process variables and the quality variables and . The remaining 35 batches are employed for updating the test. The prediction result of GP training is illustrated in Figure 3. The solid curve represents the prediction of the mean and the red dashed line represents the real data. The shade region denotes of the prediction distribution around the mean value. It can be seen that the shadow is very thin, indicating that the training model is quite accurate and representative.

Table 1: Process variables used in the modeling of the fed-batch penicillin process.

Number
Measured variables

1
Aeration rate (L/h)

2
Agitator power (W)

3
Substrate feed rate (L/h)

4
Substrate feed temperature (K)

5
Dissolved oxygen concentration (g/L)

6
Carbon dioxide concentration (g/L)

7
pH

8
Fermenter temperature (K)

9
Generated heat (kcal)

10
Cooling water flow rate (L/h)

11
Culture volume (L)

Figure 3: Prediction result of the GP model using initial batches.

To test the efficiency of the proposed method, the prediction result using the initial GP model without updating is firstly discussed. With the slow time-varying dynamics, the final quality production has changed slightly from one batch after another. If the initial model is still used to predict all the following bathes, the prediction of mean will drift from the real value gradually (Figure 4). Moreover, even of the prediction distribution cannot cover the real value in the last several batches. In this situation, the SPE statistics based on the prediction error have several false alarms (Figure 5). The false alarm rate is the percentage of regular samples as faulty ones for the normal part (Type I error). By contrast, with the new batches data, the proposed RGP model can judge when the model should be updated and which data should be selected. Thus, the time-varying dynamics are reflected by the updating model and the prediction performance is improved. In Figure 6, it is found that the predicted value is close to the real data and the covariance shadow is narrow, which indicates the result is reliable. Based on the prediction result, the SPE statistic can be calculated for all the 35 bathes. As shown in Figure 7, there are almost no false alarms during testing all the samples. To sum up, the RMSE values and false alarms of these two models listed in Table 2 indicate that the accuracy of predicting the monitoring performance is improved by recursive GP. RMSE (root-mean-square error) is the sample standard deviation of the differences between predicted values and observed values.

Table 2: Mean prediction and the process monitoring result in the penicillin process.

Type
GP model
Recursive GP model

RMSE values
0.0092
0.0016

False alarm
0.2615
0.0240

Figure 4: Prediction result of the test data using the initial GP model.

Figure 5: The SPE statistic of the test data using the initial GP model.

Figure 6: Prediction result of the test data using the recursive GP model.

Figure 7: The SPE statistic of the test data using the recursive GP model.

Finally, for a clearer explanation of the updating model procedures, the 6th batch data are demonstrated particularly. In Figure 8(a), it can be found that the prediction becomes inaccurate staring from the 25th sampling interval if the original model is still used. The quality monitoring results are shown in Figure 9(a). Part of the SPE statistic region of these data has exceeded the control limit. Take the 25th sampling, for instance. is larger than the confidence bound while is still below it. This is caused by the time-varying dynamics instead of a fault. Therefore, this data point should be added for model updating. After adding all the needed data points, the model can be updated, and the prediction result and the quality monitoring performance are given in Figures 8(b) and 9(b). The result reveals the prediction is better and the false alarm has been eliminated. It is also noticed that only 7 data sets from the 6th batch run are treated as the active ones and used for updating the model based on the criterion proposed in Section 3.2, which makes the adaptive strategy time-saving since only the data containing useful information are used for updating.

(a)

(b)

Figure 8: Prediction result of the 6th batch before and after model updating.

(a)

(b)

Figure 9: The SPE statistic of the 6th batch before and after model updating.

4. Conclusion

In this paper, an adaptive quality monitoring scheme is proposed based on the RGP model. It is mainly used for the slow time-varying batch process with limited initial batches. The cases discussed in this paper are commonly seen in fine chemicals. Using the initial normal data, a GP model and the corresponding SPE statistic can be constructed for revealing the variations of early batches. With the accumulation of the newly evolving data, an updating strategy is introduced to judge if the new data should be added so that the model can be enhanced. Then a recursive GP model is used for implementation of updating. The proposed method is efficient since it avoids the computation of the inversion of the covariance matrix in the GP model. It is helpful for online updating. The feasibility and the monitoring performance have been demonstrated by a real penicillin fermentation batch process. Among different batches, the operating state has changed slowly and the proposed method can make an accurate prediction and right monitoring decision over batches. Therefore, the proposed adaptive monitoring scheme is an effective tool for solving the slow time-varying behavior of the process. The proposed method is mainly used for slow-responding chemical processes; the corresponding model changes slightly during the operation. In this situation, only the kernel matrix is easily updated for online modelling. For more general processes, regardless of the fast-responding or the slow-responding processes, in the future work it is desirable to derive a new recursive GP model in which both kernel matrix and kernel parameters are updated simultaneously. Also, in this work, only single quality output is considered. For MIMO (multi-input and multioutput) regression, simply, a few independent GP models can be used if several key quality variables needed to be monitored, but the above method cannot reveal the interactions between different outputs. However, MIMO GP is hard to implement since the covariance matrix of the outputs will increase exponentially when the intracovariance between the outputs is calculated, not to mention a recursive MIMO GP. This part is worthy of being extended.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (61273167) and the National Basic Research Program (973 Program) of China Grant no. 2012CB720505.

References

J. Chen and K.-C. Liu, “On-line batch process monitoring using dynamic PCA and dynamic PLS models,” Chemical Engineering Science, vol. 57, no. 1, pp. 63–75, 2002.

L. H. Chiang, R. D. Braatz, and E. L. Russell, Fault Detection and Diagnosis in Industrial Systems, Springer, Berlin, Germany, 2001.

S. X. Ding, Model-Based Fault Diagnosis Techniques, Springer, New York, NY, USA, 2008.

X.-T. Doan and R. Srinivasan, “Online monitoring of multi-phase batch processes using phase-based multivariate statistical process control,” Computers & Chemical Engineering, vol. 32, no. 1-2, pp. 230–243, 2008.

Z. Ge, F. Gao, and Z. Song, “Batch process monitoring based on support vector data description method,” Journal of Process Control, vol. 21, no. 6, pp. 949–959, 2011.

J. Korbicz, Fault Diagnosis: Models, Artificial Intelligence, Applications, Springer, 2004.

C. Ündey and A. Çinar, “Statistical monitoring of multistage, multiphase batch processes,” IEEE Control Systems Magazine, vol. 22, no. 5, pp. 40–52, 2002.

B. M. Wise, N. B. Gallagher, S. W. Butler, D. D. White Jr., and G. G. Barna, “A comparison of principal component analysis, multiway principal component analysis, trilinear decomposition and parallel factor analysis for fault detection in a semiconductor etch process,” Journal of Chemometrics, vol. 13, no. 3-4, pp. 379–396, 1999.

J. Yu and S. J. Qin, “Multiway gaussian mixture model based multiphase batch process monitoring,” Industrial and Engineering Chemistry Research, vol. 48, no. 18, pp. 8585–8594, 2009.

H. H. Yue, S. J. Qin, R. J. Markle, C. Nauert, and M. Gatto, “Fault detection of plasma etchers using optical emission spectra,” IEEE Transactions on Semiconductor Manufacturing, vol. 13, no. 3, pp. 374–385, 2000.

L. Zhou, J. Chen, Z. Song, Z. Ge, and A. Miao, “Probabilistic latent variable regression model for process-quality monitoring,” Chemical Engineering Science, vol. 116, pp. 296–305, 2014.

P. Nomikos and J. F. MacGregor, “Monitoring batch processes using multiway principal component analysis,” AIChE Journal, vol. 40, no. 8, pp. 1361–1375, 1994.

P. Nomikos and J. F. MacGregor, “Multi-way partial least squares in monitoring batch processes,” Chemometrics and Intelligent Laboratory Systems, vol. 30, no. 1, pp. 97–108, 1995.

C. Zhao, F. Wang, Z. Mao, N. Lu, and M. Jia, “Adaptive monitoring based on independent component analysis for multiphase batch processes with limited modeling data,” Industrial and Engineering Chemistry Research, vol. 47, no. 9, pp. 3104–3113, 2008.

S. Wold, N. Kettaneh, H. Fridén, and A. Holmberg, “Modelling and diagnostics of batch processes and analogous kinetic experiments,” Chemometrics and Intelligent Laboratory Systems, vol. 44, no. 1-2, pp. 331–340, 1998.

C. Zhao, F. Wang, Z. Mao, N. Lu, and M. Jia, “Improved batch process monitoring and quality prediction based on multiphase statistical analysis,” Industrial and Engineering Chemistry Research, vol. 47, no. 3, pp. 835–849, 2008.

M. Kano and Y. Nakagawa, “Data-based process monitoring, process control, and quality improvement: recent developments and applications in steel industry,” Computers and Chemical Engineering, vol. 32, no. 1-2, pp. 12–24, 2008.

N. Lu and F. Gao, “Stage-based process analysis and quality prediction for batch processes,” Industrial and Engineering Chemistry Research, vol. 44, no. 10, pp. 3547–3555, 2005.

C. Ündey, S. Ertunç, and A. Çinar, “Online batch/fed-batch process performance monitoring, quality prediction, and variable-contribution analysis for diagnosis,” Industrial and Engineering Chemistry Research, vol. 42, no. 20, pp. 4645–4658, 2003.

K. P. Murphy, Machine Learning: A Probabilistic Perspective, MIT Press, 2012.

Z. Ge and Z. Song, “Nonlinear probabilistic monitoring based on the gaussian process latent variable model,” Industrial and Engineering Chemistry Research, vol. 49, no. 10, pp. 4792–4799, 2010.

B. Schölkopf, A. Smola, and K. R. Müller, “Nonlinear component analysis as a kernel eigenvalue problem,” Neural Computation, vol. 10, no. 5, pp. 1299–1319, 1998.

N. D. Lawrence, “Gaussian process latent variable models for visualisation of high dimensional data,” in Advances in Neural Information Processing Systems, vol. 16, pp. 329–336, 2004.

S. J. Qin, “Statistical process monitoring: Basics and beyond,” Journal of Chemometrics, vol. 17, no. 8-9, pp. 480–502, 2003.

G. E. Box, “Some theorems on quadratic forms applied in the study of analysis of variance problems. I. Effect of inequality of variance in the one-way classification,” The Annals of Mathematical Statistics, vol. 25, pp. 290–302, 1954.

G. H. Golub and C. F. van Loan, Matrix Computations, The Johns Hopkins University Press, Baltimore, Md, USA, 2012.

G. Birol, C. Ündey, and A. Çinar, “A modular simulation package for fed-batch fermentation: penicillin production,” Computers & Chemical Engineering, vol. 26, no. 11, pp. 1553–1565, 2002.

G. Birol, C. Ündey, S. J. Parulekar, and A. Çinar, “A morphologically structured model for penicillin production,” Biotechnology and Bioengineering, vol. 77, no. 5, pp. 538–552, 2002.

L. Zhou, Z. Song, J. Chen, Z. Ge, and Z. Li, “Process-quality monitoring using semi-supervised probability latent variable regression models,” in Proceedings of the World Congress, pp. 8272–8277, Cape Town, South Africa, 2014.