cefoperazone sodium ( cef ) , ( 6r,7r)-7-[[(2r)-2-[[(4-ethyl-2,3-dioxopiperazin-1-yl)carbonyl]amino]-2-(4-hydroxyphenyl)acetyl]amino]-3-[[(1-methyl-1h - tetrazol-5-yl ) sulphanyl]methyl]-8-oxo-5-thia-1-azabicyclo[4.2.0]oct-2-ene-2-carboxylate ( figure 1(a ) ) , is a third - generation cephalosporin that acts by inhibiting bacterial cell wall biosynthesis . according to b. p. , 7-aminocephalosporanic acid ( 7-aca ) ( figure 1(b ) ) , ( 6r,7r)-3-[(acetyloxy)methyl]-7-amino-8-oxo-5-thia-1-azabicyclo[4.2.0]oct-2-ene-2-carboxylic acid , and 5-mercapto-1-methyl - tetrazole ( 5-mer ) ( figure 1(c ) ) , 1-methyl-1h - tetrazole-5-thiol , are deemed as potential impurities for cef . 7-aca is a significant and key intermediate needed for the synthesis of semisynthetic cephalosporin antibiotics in pharmaceutical industries . chemical compounds containing this core are known to be relatively stable to hydrolysis and tolerant to -lactamase , the enzyme that hydrolyzes cephalosporins . accordingly , 7-aca is a very important moiety for the activity of all cephalosporins . 5-mer is one of tetrazole derivatives that are used as intermediate of cephalosporin side chains . literature review presents several analytical methods for assay of cef in its pharmaceutical formulation including spectrophotometric methods for determination of cef [ 5 , 6 ] , near - infrared reflectance spectroscopy , and derivative uv spectrophotometry for determination of cef in combination with sulbactam . chromatographic methods were applied for analysis of cef and sulbactam [ 9 , 10 ] ; besides , an hplc method with -cyclodextrin stationary phase for determination of cef , ampicillin , and sulbactam was reported . cef and sulbactam were determined in plasma also by lc - ms / ms method . additionally , electrochemical behavior and voltammetric determination of cef [ 13 , 14 ] were reported . firstly , the presented chemometric models show the ability of multivariate models to analyze selectively cef in ternary mixtures with its two reported impurities using cost - effective and available instruments like uv spectrophotometer . second , the presented study aims to compare two methods of multivariate calibration , plsr and linear svr chemometric models , through assay of cef , 7-aca , and 5-mer mixtures indicating the advantages and limitations of each model . the selected models offer better trueness and precision for quantitative determination of cef in pharmaceutical formulation compared to the reported hplc method . a double beam uv - vis spectrophotometer ( shimadzu , japan ) model uv-1601 pc equipped with a quartz cell of 1  cm width and connected to ibm compatible computer was used . ( cas number 957 - 68 - 6 ) and 5-mer ( cas number 13183 - 79 - 4 ) were bought from sigma - aldrich through the egyptian international center for import and export ( eic , egypt ) . cefobid 0.5  gm vials ( batch number ( 2203 ) ) were produced by pfizer pharmaceutical industries co. and cefoperazone 1  gm vials ( batch number ( 1240395 ) ) were manufactured by sigmatec pharmaceutical industries co. there are 3 concentrations available in the market for vials : 0.5  gm , 1  gm , and 1.5  gm . all chemicals and solvents used in this study were of analytical grade , including methanol ( e. merck , germany ) , water for injection b. p. 2003 ( egypt otsuka pharmaceutical co. , s.a.e . , 10th of ramadan city , are ) , and dipotassium hydrogen phosphate ( k2hpo4 ) ( el - nasr pharmaceutical chemicals co. , abu - zaabal , cairo , egypt ) . ( a ) stock standard solutions of 1  mg  ml for each of cef , 7-aca , and 5-mer were prepared in 3  ml of 0.05  m k2hpo4 solution and volume was completed with pure methanol . 7-aca is soluble only in slightly alkaline solvent ; accordingly the solubility was achieved by addition of fixed small volume of 0.05  m k2hpo4 solution before adding methanol . this step was carried out with all stock solutions of cef , 7-aca , 5-mer , and pharmaceutical formulation as well . ( b ) working standard solution of 100  g  ml of cef was prepared in methanol . working standard solution ( 1 ) of 100  g  ml and working standard solution ( 2 ) of 10  g  ml for both 7-aca and 5-mer were prepared to allow preparation of final mixtures with very small concentrations of the impurities , up to 3% of cef calculated on molar basis . uv spectra of a set of standards of cef from 1 to 70  g  ml were recorded from 210 to 300  nm . cef exhibited linearity between 5 and 50  g  ml at its max is 229  nm . the superimposed spectra of 10  g  ml of cef , 7-aca , and 5-mer are shown in figure 2 . a 4-level , 3-factor calibration design was implemented using 4 concentration levels coded as + 2 , + 1 , 1 , and + 1 , where ( 1 ) is the central level for each of the components to be analyzed including the main drug ( cef ) and its two impurities ( 7-aca and 5-mer ) . the design aims to span the mixture space appropriately , where 4 mixtures for every compound at every concentration level exist , ending up with 16 mixtures for the training set . the central level selected for the design was 20  g  ml for cef and the concentration of each level for cef depended on its calibration range . concentration levels of impurities were based on involving impurities in an amount up to 3% of cef calculated on molar basis to span most of the probabilities in future analyses . 2d scores plot of the first two pcs of the concentration matrix was drawn to affirm the orthogonality , rotatability , and symmetry of the training set mixtures ( presented as circles ) as shown in figure 3 . mean centering of all types of data was the best preprocessing method to provide best results . to examine the validity and prediction capabilities of the introduced chemometric models , the independent test set mixtures were obtained by preparation of nine independent mixtures other than the training set mixtures but within the concentration space of the design as indicated in table 1 . the well position of the mixtures of both training set and test set mixtures is shown in figure 3 . for each of the two dosage forms , accurately weighed aliquot equivalent to 100  mg of cef was transferred into 100  ml volumetric flask . to prepare stock solution , 3  ml of 0.05  m k2hpo4 the solution was then diluted to prepare 100  g  ml working solution using methanol as solvent . lastly , 2  ml portion of the working solution was diluted to 10  ml with methanol . this experiment was replicated six times and the produced spectra were analyzed by the suggested models . codes for plsr ( pls1 algorithm , bootstrap , and grid search were written using matlab 7.5.0 ( r2007b ) ) . multivariate calibration models are chemometric tools that set a relation between the spectra in data matrix x and the concentrations in data vector c. multiple linear regression ( mlr ) , principal component regression ( pcr ) , and partial least squares regression ( plsr ) are among the common models used for pharmaceutical analysis . pcr and plsr methods are more preferable than mlr , because mlr needs further variable selection steps to perform optimally and to avoid multicollinearity . additionally , pcr and plsr can cope with a large number of spectral variables by decomposing the x data matrix into a relatively small number of scores . the scores matrix t then replaces the original x data matrix in the further steps . plsr is more developed than pcr , where the c data vector shares in construction of the scores as well [ 17 , 18 ] . furthermore , the compression to a small number of scores works as a useful filter for noise . the plsr model depends on the theory of existence of a linear relation between the x data matrix and the independent variables in concentration vector c . the data matrix x and the response vector c are decomposed using a given number of pls components ( latent variable lvs ) [ 2124 ] as follows:(1)x = tp+e , c = tq+f , where t and p are the scores and loadings for x and q is the loading vector for c. e and f are the residuals for x and c , respectively . plsr is considered one of the best in multivariate calibration because it is reported to perform better than mlr and pcr in several pharmaceutical applications . optimization of number of pls components for plsr model . for prediction of optimum number of pls components , bootstrap technique [ 25 , 26 ] was used . this technique is based on dividing the original training set to two - thirds ( bootstrap training set ) and one - third ( bootstrap test set ) . the plsr model is then applied through building a model with the bootstrap training set to predict concentrations in the bootstrap test set and calculating the prediction error through the following equation:(2)rmsep=1nn=1ncnc^na2,where n is the number of bootstrap test set samples , cn is the known concentration of sample n , and c^na is the corresponding predicted concentration at the defined number of pls components . equation ( 2 ) represents just one iteration out of 1000 iterations that were implemented in the presented study . the higher the number of iterations , the higher the probability of selecting all samples in both data sets ( training set and test set ) . finally a plot was established of the average of the 1000 root mean square error of prediction ( rmsep ) values for different number of pls components against the corresponding number of pls components to choose the optimum number of optimum pls components . the bootstrap training set was mean centered with every iteration . for a data set x  ( i  j ) of an output vector c , svr models aim to find a multivariate regression function f(x ) that depends on x to predict a desired response ( e.g. , concentration of a chemical compound ) from an object ( e.g. , a spectrum ) . svr model equations are illustrated in literature [ 2628 ] and the summary equation can be given as follows:(3)fx=i , j=1niixixi+b , where i and i are the lagrange multipliers that fit to the constraint 0  i , i  c. c is known as penalty error or regularization constant which determines the trade - off between model simplicity and training error . further illustrations of ( 3 ) and the parameters  and c are found in literature [ 27 , 29 , 30 ] . -insensitive loss function is an additional factor commonly applied for svr and will be used and optimized in this study [ 31 , 32 ] . linear svr model is used in this study , where the used spectral data exhibit linearity guaranteed through the well planned experimental design . finally , in prediction step , unknown c^ value can be calculated as given below : ( 4)c^=i=1iiixxj+b . optimum  and c values were calculated by using a grid search that depends on 4-fold cross validation to give the lowest root mean square error of cross validation ( rmsecv ) . primary range of values was set for  ( 0.011 ) and c ( 301000 ) . for each set of svr parameters , 4 samples ( n = 4 ) were taken out ; the linear established svr model was applied on the remaining 12 ( i  n ) samples . further , rmsecv was predicted for the n samples and finally , the average of rmsecv after removal of all samples was calculated as follows:(5)rmsecv=1ii=1icic^i2,where ci is the true concentration for sample n and c^i is the corresponding predicted concentration . bootstrap technique was applied to choose optimum number of pls components to build the best calibration model based on the training set . concerning svr , the grid search that gave the lowest rmsecv ( 5 ) resulted in the values (  = 0.02 and c = 280 ) . this study aims to introduce a comparative study between two chemometric methods known as plsr and linear svr via analysis of cef in presence of its reported impurities : 7-aca and 5-mer . the two multivariate models could handle the uv data and overcome the overlapping difficulty of the components ' spectra shown in figure 2 . both models were successfully able to determine the concentrations of cef in training set and test set manifested by high recovery % with low sd as presented in table 2 . the rmsep is a parameter used to evaluate the prediction abilities of the two models , table 2 . rmsep comparative plot between plsr and linear svr for prediction of test set samples is presented in figure 5 . there are many ways of comparison included in our study such as root mean square error of calibration ( rmsec ) , root mean square error of prediction ( rmsep ) , and calculation intensity and optimization steps . first , it was observed that linear svr gives the least rmsec ( 0.1960 ) compared to plsr ( 0.2306 ) indicating better trueness and the corresponding standard deviation is smaller also indicating better precision . second , the comparative bar plot in figure 5 indicates that linear svr gives the least rmsep ( 0.3386 ) compared to plsr ( 0.4457 ) indicating higher ability to process future samples and better generalization ability of linear svr when compared to plsr . concerning calculation intensity , plsr is computationally simpler than svr which is deemed to be more intense in calculation and time - consuming because of intense calculations for optimization . choice of optimum parameters ' values for svr model could be misleading , where parameters that give the lowest rmsecv can be used ( considered as overfitting as in plsr ) but still give better rmsep or vice versa . in this study , we apply 4-fold cross validation to optimize svr parameters to evade overfitting through prediction of small subsets of data instead of single sample ( as used in leave - one - out cross validation technique ) , so improving the robustness of model and its generalization ability . accordingly , with plsr , selecting too many components will lead to less robust model , that is , with less ability to predict future samples that have unknown signals . svr parameters are calculated by using k - fold cross validation to evade overfitting , showing better robustness and higher prediction ability for future samples and hence being a more general model . the probability of adherence of svr to overfitting is less common than plsr [ 35 , 36 ] . the implemented linear svr model in this study has another merit of optimizing only two parameters , unlike the nonlinear svr models that use kernels and need optimization of more parameters and hence are more time - consuming in optimization process . in conclusion , plsr model is considered as one of the best in multivariate calibration and usually used in quality control routine applications and industry . it proved to be simple in optimization and computation and gives comparable results to reference hplc methods in spite of processing simple uv data . however , svr is still considered as a more general model with higher predictive ability for future samples . the proposed chemometric methods were implemented for analysis of cef in cefobid 0.5  gr vial and cefoperazone 1  gr vial and satisfactory results with good recoveries were obtained . these results were statistically compared to the results obtained by applying the reported hplc method   using t- and f - tests . the calculated t- and f - values are less than the tabulated ones showing no significant difference between the two introduced models and the reference hplc method with regarding both trueness and precision , table 3 . in general , the goals of this paper are presenting two multivariate chemometric models , plsr and linear svr , for analysis of cef in presence of its reported impurities and making a modest comparison between the two models highlighting the advantages and limitations of each . concerning predictive ability , the linear svr proved to be better than plsr according to rmsep values indicating better generalization ability . the two chemometric methods were also applied for the pharmaceutical formulations and statistically compared to reference hplc method . the calculated t- and f - values were found to be less than tabulated ones showing no significant difference in respect to both trueness and precision . the presented advantages of the proposed models suggest their use for routine quality control analysis without interference of normally encountered pharmaceutical additives or impurities that could be present in minor ratios . additionally , the obtained results affirm the possibility of using modern chemometric approaches , especially linear svr , for assay of different pharmaceutical dosage forms using accessible cheap and simple instruments like uv spectrophotometer even in presence of large number of interfering components with extremely overlapped spectra .
a comparison between partial least squares regression and support vector regression chemometric models is introduced in this study . 
 the two models are implemented to analyze cefoperazone sodium in presence of its reported impurities , 7-aminocephalosporanic acid and 5-mercapto-1-methyl - tetrazole , in pure powders and in pharmaceutical formulations through processing uv spectroscopic data . 
 for best results , a 3-factor 4-level experimental design was used , resulting in a training set of 16 mixtures containing different ratios of interfering moieties . for method validation , 
 an independent test set consisting of 9 mixtures was used to test predictive ability of established models . 
 the introduced results show the capability of the two proposed models to analyze cefoperazone in presence of its impurities 7-aminocephalosporanic acid and 5-mercapto-1-methyl - tetrazole with high trueness and selectivity ( 101.87  0.708 and 101.43  0.536 for plsr and linear svr , resp . ) . 
 analysis results of drug products were statistically compared to a reported hplc method showing no significant difference in trueness and precision , indicating the capability of the suggested multivariate calibration models to be reliable and adequate for routine quality control analysis of drug product . 
 svr offers more accurate results with lower prediction error compared to plsr model ; however , plsr is easy to handle and fast to optimize .