<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.2 20190208//EN" "https://jats.nlm.nih.gov/publishing/1.2/JATS-journalpublishing1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:mml="http://www.w3.org/1998/Math/MathML" article-type="research-article" dtd-version="1.2" xml:lang="en">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">IJDST</journal-id>
<journal-title-group><journal-title>International journal of Dairy Science and Technology</journal-title></journal-title-group>
<publisher><publisher-name>Thomson &amp; Ryberg Publications</publisher-name></publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.62715/ijdst.2026.v1.i1.pp21-26</article-id>
<article-id pub-id-type="publisher-id">1e631d5e37</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Methodology Article</subject></subj-group></article-categories>
<title-group><article-title>Rapid Detection of Milk Adulterants Using Portable Near-Infrared Spectroscopy and Chemometric Calibration</article-title></title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes"><name><surname>Chukwu</surname><given-names>Aiko</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1">*</xref></contrib>
<contrib contrib-type="author"><name><surname>Gutierrez</surname><given-names>Vikram</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib>
<contrib contrib-type="author"><name><surname>Harutyunyan</surname><given-names>Agnieszka</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref></contrib>
<aff id="aff1"><label>1</label>Department of Agri, Tagus Academy of Agri, Croatia</aff>
<aff id="aff2"><label>2</label>Department of Agri, Thames School of Agri, Mexico</aff>
<aff id="aff3"><label>3</label>Department of Agri, Periyar Academy of Agri, Czech Republic</aff>
</contrib-group>
<author-notes><corresp id="cor1">* Correspondence: <email>aiko.chukwu@tagus.edu</email></corresp></author-notes>
<pub-date pub-type="epub" iso-8601-date="2026-08-04"><day>04</day><month>08</month><year>2026</year></pub-date>
<volume>1</volume>
<issue>1</issue>
<fpage>21</fpage>
<lpage>26</lpage>
<history>
<date date-type="received" iso-8601-date="2026-04-11"><day>11</day><month>04</month><year>2026</year></date>
<date date-type="accepted" iso-8601-date="2026-06-17"><day>17</day><month>06</month><year>2026</year></date>
</history>
<permissions><copyright-statement>© 2026 The Author(s)</copyright-statement><copyright-year>2026</copyright-year><license license-type="open-access" xlink:href="https://creativecommons.org/licenses/by/4.0/"><license-p>This is an open access article distributed under the Creative Commons Attribution 4.0 International License (CC BY 4.0).</license-p></license></permissions>
<self-uri xlink:href="https://trpub.org/article.php?id=1e631d5e37"/>
<abstract><p>Handheld near-infrared spectrometers now cost a fraction of what bench instruments do and can be carried to a collection centre or a roadside checkpoint, which makes them attractive for screening raw milk for added water, urea, starch, whey and similar adulterants. A spectrometer alone detects nothing, however. Its usefulness rests entirely on the calibration that turns a spectrum into a decision, and published feasibility studies often leave out the steps that determine whether a calibration will survive contact with real milk. This article proposes a complete procedure for building, validating and maintaining such a calibration. It covers the definition of the screening question, the design of a reference sample set that spans natural variation in authentic milk as well as the adulterants of concern, controlled acquisition of spectra, pre-processing, exploratory analysis, the choice between class-modelling, discriminant and regression approaches, validation with independent samples, and the monitoring and transfer of models once instruments are in the field. The procedure is drawn from the chemometric literature and from international guidance on near-infrared analysis of milk. It is presented as a proposal and has not been validated experimentally here, and the limits of what near-infrared screening can be expected to detect are stated plainly.</p></abstract>
<kwd-group kwd-group-type="author">
<kwd>Milk adulteration</kwd>
<kwd>Near-infrared spectroscopy</kwd>
<kwd>Portable spectrometers</kwd>
<kwd>Chemometrics</kwd>
<kwd>Partial least squares</kwd>
<kwd>Food fraud</kwd>
</kwd-group>
</article-meta>
</front>
<body>
<sec id="sec1"><title>1. Introduction</title>
<p>Milk is among the foods most often named in records of economically motivated adulteration (Moore et al., 2012). The commonest fraud is the oldest: water is added to increase volume, and other substances are then added to hide the dilution by restoring density, apparent protein, sweetness or keeping quality. Reviews list urea, starch, sucrose, salt, whey, vegetable oil, detergents, neutralisers and preservatives such as formalin and hydrogen peroxide among the materials found (Poonia et al., 2017). The melamine contamination of infant formula in China in 2008 showed how serious the consequences can be when a nitrogen-rich compound is used to defeat a protein test based on total nitrogen, and the burden of milk fraud falls most heavily where control systems are weakest (Handford et al., 2016).</p>
<p>Reference methods for each adulterant exist, but they are slow, need a laboratory and test for one substance at a time. Near-infrared spectroscopy offers a different bargain: a measurement that takes seconds, needs no reagents and responds to the overall composition of the sample. Miniaturised instruments have brought that capability to the point of collection (Crocombe, 2018). Feasibility studies have shown that infrared spectra, combined with multivariate analysis, can distinguish adulterated from authentic milk under laboratory conditions (Kasemsumran et al., 2007; Santos et al., 2013).</p>
<p>The gap between such demonstrations and a dependable screening tool is methodological. This article proposes a procedure for closing it. No experiment is reported: the procedure is assembled from established chemometric practice and has yet to be validated on milk by the route described.</p>
</sec>
<sec id="sec2"><title>2. Background</title>
<p>Near-infrared radiation, conventionally 780 to 2500 nm, is absorbed by overtones and combinations of the fundamental vibrations of C-H, N-H and O-H bonds. The bands are broad and overlapping, and in milk they are dominated by water, whose strong absorptions near 1450 and 1940 nm leave relatively narrow windows in which fat, protein and lactose can be seen. Milk also scatters light strongly, because fat globules and casein micelles are comparable in size to the wavelengths used, so a spectrum carries physical as well as chemical information. International guidance on near-infrared analysis of milk products sets out the consequences for sampling, calibration and validation (ISO and IDF, 2006).</p>
<p>Portable instruments add constraints of their own. Many cover only part of the range, often the short-wave region up to about 1700 nm, with lower resolution and signal-to-noise ratio than bench spectrometers, and units of the same model differ from one another (Crocombe, 2018). None of this prevents useful work, but it means that an adulterant present at a low concentration contributes a small signal on top of large natural variation, and that only a carefully built multivariate model can separate the two.</p>
</sec>
<sec id="sec3"><title>3. The Proposed Procedure</title>
<p>The procedure is arranged as a cycle of six stages, shown in Figure 1, because a calibration is never finished: monitoring in use feeds back into the sample set and the model. The stages are described below.</p>
<fig id="fig1"><label>Figure 1</label><caption><p>The proposed calibration life cycle for portable near-infrared screening of milk</p></caption><graphic xlink:href="IJDST_V1I1_2026_04_Figure1.png" mimetype="image" mime-subtype="png"/></fig>
<sec id="sec3-2"><title>3.1. Stage 1: Define the screening question</title>
<p>The first decision concerns what the device is to answer. One formulation is targeted: is a named adulterant present above a stated level? The other is non-targeted: does this sample look like authentic milk? The two lead to different models and different sample sets. The developer should also fix the consequence of a positive result, normally that the lot is held and a confirmatory reference test is run, and the lowest adulteration level that is economically worthwhile for a fraudster, since there is little point in pursuing detection far below it.</p>
</sec>
<sec id="sec3-3"><title>3.2. Stage 2: Assemble reference samples</title>
<p>A calibration can only recognise variation it has seen. The authentic set should therefore span the sources of natural variability in the milk the device will meet: breeds, stages of lactation, seasons, feeding systems, farms, and morning and evening milkings. Authenticity must be assured, in practice by supervised milking or collection. Adulterated samples are prepared by spiking authentic milk at graded levels around the level of interest, using many different base milks so that the model cannot learn the base instead of the adulterant. Combined adulterations, such as water with urea or with sucrose, should be included because they are what fraud looks like in practice. Reference values, where regression is intended, come from the gravimetric preparation and from standard methods.</p>
</sec>
<sec id="sec3-4"><title>3.3. Stage 3: Acquire spectra under controlled conditions</title>
<p>Spectra of milk are sensitive to temperature, which shifts the water bands, and to the state of the fat. Samples should be brought to a fixed temperature and mixed gently before measurement, in a cell or probe of fixed path length, with the measurement mode (transmission, transflectance or diffuse reflectance) kept constant. Each sample is scanned in replicate with repositioning, a reference or background is recorded at set intervals, and the order of measurement is randomised so that drift does not line up with adulterant level. Instrument serial number, operator, date and ambient conditions are logged, since they will be needed later to explain model failures.</p>
</sec>
<sec id="sec3-5"><title>3.4. Stage 4: Pre-process, explore and model</title>
<p>Pre-processing removes variation unrelated to composition. Standard normal variate transformation (Barnes et al., 1989) and multiplicative scatter correction (Geladi et al., 1985) reduce the effect of scatter, and Savitzky-Golay derivatives remove baseline offsets and sharpen overlapping bands (Savitzky and Golay, 1964). Rinnan et al. (2009) review these techniques and caution that every additional step can remove information as well as noise; the simplest treatment that works should be preferred. A principal component analysis of the pre-processed spectra is then examined for outliers, for clustering by day or instrument, and for any visible trend with adulteration.</p>
<p>Three families of model are available, and Table 1 compares them. Partial least squares regression relates spectra to a concentration and suits a targeted, quantitative question (Wold et al., 2001). Partial least squares discriminant analysis assigns a sample to one of several predefined classes (Barker and Rayens, 2003). Class-modelling methods, of which soft independent modelling of class analogy is the best known, describe the authentic class alone and reject whatever falls outside it. Rodionova et al. (2016) argue that authentication is properly a class-modelling problem, because the set of possible adulterants is open-ended and a discriminant model is forced to place an unfamiliar sample in one of its known classes. The procedure proposed here follows that reasoning: a one-class model of authentic milk as the primary screen, supplemented by regression models for the few adulterants of greatest local concern.</p>
<table-wrap id="tbl1"><label>Table 1</label><caption><p>Modelling approaches for near-infrared screening of milk and their suitability</p></caption><table frame="hsides" rules="groups"><thead><tr><th>Approach</th><th>Question answered</th><th>Strength</th><th>Limitation</th></tr></thead><tbody><tr><td>One-class modelling (e.g. SIMCA)</td><td>Is the sample consistent with authentic milk?</td><td>Can flag adulterants not in the training set</td><td>Needs a broad authentic set; no identification</td></tr><tr><td>PLS discriminant analysis</td><td>Which predefined class does it belong to?</td><td>Good separation of known adulterants</td><td>Unknown adulterants are forced into a known class</td></tr><tr><td>PLS regression</td><td>How much of a named adulterant is present?</td><td>Gives a level to compare with a threshold</td><td>One model per adulterant; needs reference values</td></tr><tr><td>Principal component analysis</td><td>What structure do the spectra show?</td><td>Reveals outliers, drift and batch effects</td><td>Exploratory only; makes no decision</td></tr><tr><td>Combined scheme (proposed)</td><td>Authentic or not, then what and how much?</td><td>Screens broadly, then characterises</td><td>More models to validate and maintain</td></tr></tbody></table></table-wrap>
</sec>
<sec id="sec3-6"><title>3.5. Stage 5: Validate with independent samples</title>
<p>Validation decides whether the model can be believed. The samples are divided into a calibration set and a test set before modelling begins, and the test set is not touched until the model is final. The split should keep all spectra of a given base milk on the same side; otherwise replicates leak across and the error is flattered. The Kennard-Stone algorithm offers one way to choose a calibration set that covers the spectral space evenly (Kennard and Stone, 1969). Cross-validation within the calibration set guides the number of latent variables, but Westerhuis et al. (2008) showed how readily discriminant models overfit and recommend permutation tests and a strictly separate test set as safeguards.</p>
<p>Performance should be expressed in terms the user can act on. For a classification screen these are sensitivity and specificity at each adulteration level, reported with the number of independent base milks behind them; for regression, the root mean square error of prediction, bias, and an estimate of the lowest level reliably detected. A final test on milk collected at a different time and place from the calibration samples is the strongest evidence, because it exposes the model to variation the developer did not choose.</p>
</sec>
<sec id="sec3-7"><title>3.6. Stage 6: Deploy, monitor and update</title>
<p>A model built on one instrument rarely transfers unchanged to another. Calibration transfer methods map the response of a second instrument onto the first using a small set of samples measured on both; direct and piecewise direct standardisation are the classical approaches (Wang et al., 1991), and the wider family is reviewed by Feudale et al. (2002). In the field, a stable check sample is measured daily and plotted on a control chart, operators follow a written measuring routine, and a share of screened lots, negatives as well as positives, is sent for reference analysis. Samples that the model handles badly are added to the calibration set at the next revision, which closes the cycle.</p>
</sec>
</sec>
<sec id="sec4"><title>4. Fitting the Procedure to Practice</title>
<p>The effort involved is front-loaded. Collecting authentic milk across seasons takes a year, and the developer should resist releasing a model built on a single season&#x27;s samples. A dairy company, a cooperative union or a regulator is better placed than a single collection centre to carry that cost, and a shared calibration maintained centrally and distributed to many devices is the natural arrangement. Operators need training in sample handling far more than in chemometrics, since most field errors arise from temperature, mixing and dirty optics.</p>
<p>Results must be interpreted within a decision rule fixed in advance. A positive screen is grounds for holding a lot and testing it by a reference method; it is not proof of fraud and should not, by itself, lead to a penalty. A negative screen means only that the sample resembles authentic milk within the sensitivity of the method.</p>
</sec>
<sec id="sec5"><title>5. Limitations and What Remains to Be Shown</title>
<p>The procedure is a proposal and has not been tested here. What a portable instrument can detect in liquid milk by following it is an open, empirical question that depends on the adulterant, the device and the variability of the local milk. Near-infrared spectroscopy is a bulk technique. Dilution with water and additions that change composition by an appreciable fraction are plausible targets; contaminants that are harmful at trace levels, of which melamine is the best-known example, are not realistic targets for direct screening of liquid milk with this technology, and claims to the contrary should be examined against the validation design that supports them.</p>
<p>Further matters would need to be established before routine use: how many authentic samples are needed for a one-class model to be stable across seasons, how often models require updating, how well transfer methods cope with the unit-to-unit differences of low-cost devices, and whether the rate of false alarms is low enough for collection staff to keep trusting the instrument. Each can be answered by a study that follows the stages set out above and publishes its independent test results in full.</p>
</sec>
<sec id="sec6"><title>6. Conclusion</title>
<p>Portable near-infrared spectrometers put a rapid, reagent-free measurement in the hands of those who receive milk, but the measurement becomes a screening test only through calibration. The procedure proposed in this article treats that calibration as a managed life cycle: a clearly stated question, a reference set that captures natural variation, disciplined acquisition, restrained pre-processing, a one-class model backed by targeted regressions, validation on independent milk, and continuing control after deployment. Followed carefully, it should yield models whose limits are known, which is the property a screening method most needs.</p>
</sec>
<sec id="sec7"><title>Conflict of interest</title>
<p>To be declared by the author(s) before publication.</p>
</sec>
<sec id="sec8"><title>Funding</title>
<p>To be declared by the author(s) before publication.</p>
</sec>
</body>
<back>
<ref-list><title>References</title>
<ref id="ref1"><label>1</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Barker, M. and Rayens, W.</string-name></person-group> (<year>2003</year>). <article-title>Partial least squares for discrimination</article-title>. <source>Journal of Chemometrics</source>, <volume>17(3)</volume>, <fpage>166-173</fpage> <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1002/cem.785">https://doi.org/10.1002/cem.785</ext-link></mixed-citation></ref>
<ref id="ref2"><label>2</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Barnes, R.J., Dhanoa, M.S. and Lister, S.J.</string-name></person-group> (<year>1989</year>). <article-title>Standard normal variate transformation and de-trending of near-infrared diffuse reflectance spectra</article-title>. <source>Applied Spectroscopy</source>, <volume>43(5)</volume>, <fpage>772-777</fpage> <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1366/0003702894202201">https://doi.org/10.1366/0003702894202201</ext-link></mixed-citation></ref>
<ref id="ref3"><label>3</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Crocombe, R.A.</string-name></person-group> (<year>2018</year>). <article-title>Portable spectroscopy</article-title>. <source>Applied Spectroscopy</source>, <volume>72(12)</volume>, <fpage>1701-1751</fpage></mixed-citation></ref>
<ref id="ref4"><label>4</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Feudale, R.N., Woody, N.A., Tan, H., Myles, A.J., Brown, S.D. and Ferré, J.</string-name></person-group> (<year>2002</year>). <article-title>Transfer of multivariate calibration models: A review</article-title>. <source>Chemometrics and Intelligent Laboratory Systems</source>, <volume>64(2)</volume>, <fpage>181-192</fpage></mixed-citation></ref>
<ref id="ref5"><label>5</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Geladi, P., MacDougall, D. and Martens, H.</string-name></person-group> (<year>1985</year>). <article-title>Linearization and scatter-correction for near-infrared reflectance spectra of meat</article-title>. <source>Applied Spectroscopy</source>, <volume>39(3)</volume>, <fpage>491-500</fpage></mixed-citation></ref>
<ref id="ref6"><label>6</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Handford, C.E., Campbell, K. and Elliott, C.T.</string-name></person-group> (<year>2016</year>). <article-title>Impacts of milk fraud on food safety and nutrition with special emphasis on developing countries</article-title>. <source>Comprehensive Reviews in Food Science and Food Safety</source>, <volume>15(1)</volume>, <fpage>130-142</fpage></mixed-citation></ref>
<ref id="ref7"><label>7</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>ISO and IDF</string-name></person-group> (<year>2006</year>). <article-title>ISO 21543:2006 | IDF 201:2006, Milk products: Guidelines for the application of near infrared spectrometry</article-title>. <source>International Organization for Standardization, Geneva</source></mixed-citation></ref>
<ref id="ref8"><label>8</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Kasemsumran, S., Thanapase, W. and Kiatsoonthon, A.</string-name></person-group> (<year>2007</year>). <article-title>Feasibility of near-infrared spectroscopy to detect and to quantify adulterants in cow milk</article-title>. <source>Analytical Sciences</source>, <volume>23(7)</volume>, <fpage>907-910</fpage></mixed-citation></ref>
<ref id="ref9"><label>9</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Kennard, R.W. and Stone, L.A.</string-name></person-group> (<year>1969</year>). <article-title>Computer aided design of experiments</article-title>. <source>Technometrics</source>, <volume>11(1)</volume>, <fpage>137-148</fpage></mixed-citation></ref>
<ref id="ref10"><label>10</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Moore, J.C., Spink, J. and Lipp, M.</string-name></person-group> (<year>2012</year>). <article-title>Development and application of a database of food ingredient fraud and economically motivated adulteration from 1980 to 2010</article-title>. <source>Journal of Food Science</source>, <volume>77(4)</volume>, <fpage>R118-R126</fpage> <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1111/j.1750-3841.2012.02657.x">https://doi.org/10.1111/j.1750-3841.2012.02657.x</ext-link></mixed-citation></ref>
<ref id="ref11"><label>11</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Poonia, A., Jha, A., Sharma, R., Singh, H.B., Rai, A.K. and Sharma, N.</string-name></person-group> (<year>2017</year>). <article-title>Detection of adulteration in milk: A review</article-title>. <source>International Journal of Dairy Technology</source>, <volume>70(1)</volume>, <fpage>23-42</fpage></mixed-citation></ref>
<ref id="ref12"><label>12</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Rinnan, Å., van den Berg, F. and Engelsen, S.B.</string-name></person-group> (<year>2009</year>). <article-title>Review of the most common pre-processing techniques for near-infrared spectra</article-title>. <source>TrAC Trends in Analytical Chemistry</source>, <volume>28(10)</volume>, <fpage>1201-1222</fpage> <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/j.trac.2009.07.007">https://doi.org/10.1016/j.trac.2009.07.007</ext-link></mixed-citation></ref>
<ref id="ref13"><label>13</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Rodionova, O.Y., Titova, A.V. and Pomerantsev, A.L.</string-name></person-group> (<year>2016</year>). <article-title>Discriminant analysis is an inappropriate method of authentication</article-title>. <source>TrAC Trends in Analytical Chemistry</source>, <volume>78</volume>, <fpage>17-22</fpage></mixed-citation></ref>
<ref id="ref14"><label>14</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Santos, P.M., Pereira-Filho, E.R. and Rodriguez-Saona, L.E.</string-name></person-group> (<year>2013</year>). <article-title>Rapid detection and quantification of milk adulteration using infrared microspectroscopy and chemometrics analysis</article-title>. <source>Food Chemistry</source>, <volume>138(1)</volume>, <fpage>19-24</fpage></mixed-citation></ref>
<ref id="ref15"><label>15</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Savitzky, A. and Golay, M.J.E.</string-name></person-group> (<year>1964</year>). <article-title>Smoothing and differentiation of data by simplified least squares procedures</article-title>. <source>Analytical Chemistry</source>, <volume>36(8)</volume>, <fpage>1627-1639</fpage> <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1021/ac60214a047">https://doi.org/10.1021/ac60214a047</ext-link></mixed-citation></ref>
<ref id="ref16"><label>16</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Wang, Y., Veltkamp, D.J. and Kowalski, B.R.</string-name></person-group> (<year>1991</year>). <article-title>Multivariate instrument standardization</article-title>. <source>Analytical Chemistry</source>, <volume>63(23)</volume>, <fpage>2750-2756</fpage></mixed-citation></ref>
<ref id="ref17"><label>17</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Westerhuis, J.A., Hoefsloot, H.C.J., Smit, S., Vis, D.J., Smilde, A.K., van Velzen, E.J.J., van Duijnhoven, J.P.M. and van Dorsten, F.A.</string-name></person-group> (<year>2008</year>). <article-title>Assessment of PLSDA cross validation</article-title>. <source>Metabolomics</source>, <volume>4(1)</volume>, <fpage>81-89</fpage></mixed-citation></ref>
<ref id="ref18"><label>18</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Wold, S., Sjöström, M. and Eriksson, L.</string-name></person-group> (<year>2001</year>). <article-title>PLS-regression: A basic tool of chemometrics</article-title>. <source>Chemometrics and Intelligent Laboratory Systems</source>, <volume>58(2)</volume>, <fpage>109-130</fpage> <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/S0169-7439(01)00155-1">https://doi.org/10.1016/S0169-7439(01)00155-1</ext-link></mixed-citation></ref>
</ref-list>
</back>
</article>