Method and system for detecting change in data streams

FIELD OF THE INVENTION

0001The present invention relates to the detection of changes in data streams. In particular, the present invention seeks to detect a change in a data stream created by a target event.

BACKGROUND OF THE INVENTION

0002Change detection systems find a wide variety of applications, including, but not limited to, fraud detection (for example, looking for changes in patterns on credit card usage), security systems (for example, detecting attacks on computer networks), process, fault and condition monitoring (for example, looking for changes in the pattern of vibration in vehicle engines), environmental monitoring systems (for example, identifying chemical spillage and pollution), and health monitoring (for example, to alert medical workers of sudden changes in the condition of patients). To detect changes in practice, entities and processes are typically monitored by taking regular measurements of critical parameters such as those parameters that are most likely to contain information about the changes of interest. Changes are usually identified by comparing the parameters to thresholds, designed to indicate a target event, with an alert generated if the thresholds are exceeded. Such thresholds are usually not fixed, but are functions of statistics extracted from a data stream containing information about the system being monitored, such as its mean, variance or percentiles.

0003The change detection performance of such systems is limited by several factors. Firstly, unless the thresholds are effectively set to extreme percentiles (such as the 99^th) the number of false alerts - those generated even when no significant change in the data stream has occurred - will be too high for many practical applications. For example, a classic application of the invention is to monitor calls in large telecommunications networks for unusual activity that may be indicative of fraud. Since it is not uncommon for such networks to support 100 million calls per day, a threshold based on the 99^th percentile would generate around 1 million false alerts per day (assuming that the percentile estimate is accurate, and that the call statistics are ergodic), which is over 1,000 times the number that can be processed by the fraud investigation teams employed by most telecommunications companies. In practice, gradual changes in the way in which the telecommunications network is used will cause the false alert rate to be even higher than the 1 percent that the 99^th percentile would prima facie imply. For example, the growth of internet usage has produced a gradual increase in the number of exceptionally long calls - those over two hours - which has been sustained over several years. The affects of such changes on the false alert rate of percentile-based algorithms cannot readily be alleviated by increasing the percentile used to, for example, the 999^th because not only does that increase the risk of missing the more subtle frauds, but more extreme percentiles also take longer to estimate. This time factor is important because a reasonable estimate of the percentile must be formed before the algorithm can produce useful fraud alerts, increasing the risk that major frauds are missed because they are committed before the algorithm is ready to detect them.

0004A more fundamental limitation of the use of thresholds is that not all changes that can conceivably occur in a data stream can be detected by them. Assume, for example, that a change detection algorithm is used to monitor the condition of the suspension of a car through the periodic measurement of the extension of a spring attached to one of the car's wheels. Ten seconds worth of simulated measurements are shown in the top graph shown in <figref idref="f0001">Figure 1</figref>. Assuming that the normal behaviour of the suspension is that there is wide variation in the spring length, as shown in the regions of the top graph of <figref idref="f0001">Figure 1</figref> that lie outside the dashed lines, it is possible that some modes of failure, such as periodic seizure, can cause the spring length to show less variation, as occurs between the dashed lines. This change in behaviour cannot be detected using thresholds because there is no threshold that can be placed on the measured parameter - that is, spring length - that would cause substantially more alerts to be generated when the suspension behaves pathologically than when it behaves normally.

0005In summary, threshold based systems can only detect changes where there is a substantial shift in the proportion of the probability mass of the contents of the data stream from below to above the threshold (of vice versa) and the choice of threshold is necessarily restricted to extreme percentiles to minimise the rate at which false alerts are generated. M. Iguchi and S. Goto, in their paper "Network Surveillance for Detecting Intrusions", 0-7803-9939-0/99, 1999IEEE, describe a method of capturing and analysing a network traffic pattern with "port profiling" by comparing a live traffic pattern with expected behaviour recorded in the port profile.

SUMMARY OF THE PRESENT INVENTION

0006The present invention seeks to provide an improved method of detecting changes in data streams, and the information they contain.

0007According to a first aspect of the present invention there is provided a method of detecting changes in a data stream comprising the steps of: <ul id="ul0001" list-style="none" compact="compact"><li>maintaining a short term distribution that models the data stream, wherein the short term distribution is updated when the input is received;</li><li>maintaining a long term distribution that models the data stream, wherein the long term distribution is updated when the input data is received;</li><li>determining a difference between the short term distribution and the long term distribution; and</li><li>applying a statistical measure to the difference and generating an alert if the measure of the difference exceeds a threshold,</li></ul> characterised in that the method further comprises the step of: <ul id="ul0002" list-style="none" compact="compact"><li>returning at least one of the short term distribution and the long term distribution to a state just before it was updated to include the input that caused the alert when the alert is generated.</li></ul>

0008Preferably the short term distribution is a model of probability distributions that describe the data stream. Preferably the long term distribution is a model of probability distributions that describe the data stream.

0009Preferably the short term distribution is a recursively estimated weighted distribution of all the data received thus far. Preferably the long term distribution is a recursively estimated weighted distribution of all of the data received thus far.

0010Preferably when an alert is generated the short term distribution is returned to a state just before it was updated to include an input that caused the alert. Preferably the long term distribution is also returned to a state just before it was updated to include an input that caused the alert.

0011Preferably an alert is generated if the difference between the distributions exceeds an adaptive alert threshold. Preferably the adaptive alert threshold is determined as some function of the short and long term distributions, previous values of the threshold, the contents of the data stream, or a subset thereof. Preferably the adaptive alert threshold is formed from a predictability measure and a variability measure. Preferably the predictability measure is a moving average of the difference between the short term distribution and the long term distribution. Preferably the variability measure is a moving average of the absolute difference between the short term distribution and the long term distribution and the predictability measure. Preferably neither the predictability nor variability measures are updated when an alert is generated.

0012Preferably the alert generated includes information in the data stream and/or a function of information in the data stream. Typically, the alert includes a propensity measure as an indication of the severity of change. Preferably the propensity measure is calculated by dividing the difference between the measure of difference between the short term distribution and the long term distribution and the alert threshold by the variability measure.

0013Preferably the method includes a step of maintaining an estimate of the amounts by which the sensitivity to the alert threshold would be needed to have been adjusted in order to not have generated an alert that turned out not to be caused by an event of interest since the last time it was instructed to adapt.

0014Preferably the sensitivity adjustment estimate is increased by an additive constant each time an alert is generated and decays exponentially with each input received.

0015Preferably a lead period may be provided during which alerts cannot be generated and the short term distribution and the long term distribution are adapted to all inputs within that period.

0016Alerts may also be suppressed by inputs that are above a configurable lower percentile or below a configurable upper percentile. The lower percentile and upper percentile are both estimated from the long term probability distribution.

0017According to a second aspect of the present invention there is provided a system for detecting changes in a data stream comprising: <ul id="ul0003" list-style="none" compact="compact"><li>distribution maintenance engine configured to: <ul id="ul0004" list-style="none" compact="compact"><li>maintain a short term distribution that models data stream, wherein the short term distribution is updated when input data is received; and</li><li>maintain a long term distribution that models the data stream, wherein the long term distribution is updated when the input data is received;</li></ul></li><li>means for determining a difference between the short term distribution and the long term distribution, and</li><li>an alert generation engine configured to apply a statistical measure to the difference and to generate an alert if the measure of the difference exceeds a threshold,</li></ul> characterised in that the distribution maintenance engine is further configured to return at least one of the short term distribution and the long term distribution to a state just before it was updated to include the input that caused the alert when the alert is generated.

0018Preferably the system further comprises a threshold adaption engine configured to adaptively determine the threshold. Preferably the distribution maintenance engine is configured to return the short term distribution to a state just before is was updated, when an alert is generated. Preferably the distribution maintenance engine is also configured to return the long term distribution to a state just before it was updated, when an alert is generated.

BRIEF DESCRIPTION OF THE DRAWINGS

0019In order to facilitate a better understanding of the nature of the invention, preferred embodiments will now be described in greater detail, by way of example only, with reference to the accompanying drawings in which: <ul id="ul0005" list-style="none"><li><figref idref="f0001">Figure 1</figref> is an example of simulated measurements;</li><li><figref idref="f0001">Figure 2</figref> is a schematic representation of a preferred embodiment of a system for detecting changes in the properties of a data stream according to the present invention;</li><li><figref idref="f0002">Figure 3</figref> is a flow chart showing the preferred embodiment of a method of detecting changes in the properties of a data stream in accordance with the present invention; and</li><li><figref idref="f0001">Figure 4</figref> is a representation of examples of alerts generated in response to the change in simulated measurements of <figref idref="f0001">Figure 1</figref>.</li></ul>

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS

0020Referring to <figref idref="f0001">Figure 2</figref> there is shown a system 10 for detecting changes in the properties of a data stream 12 to produce alerts 22 that a target event may be indicated within the data stream 12. The system includes a distribution maintenance engine 14 for maintaining a short term distribution 24 and a long term distribution 26; a means 16 for applying a statistical measure to the difference between the short term distribution 24 and long term distribution 26; an alert generation engine for generating alerts 22 and a threshold adaption engine 20 for adapting an alert generating threshold to certain conditions.

0021In the present invention, the system 10 may be implemented in the form of a computer or a network of computers programmed to perform each of the stages of processing the data from the data stream. For example, a single computer could be programmed to run the system 10 or a dedicated computer may be programmed to run each of the components of the system 10.

0022Referring now to <figref idref="f0001">Figures 2</figref> and <figref idref="f0002">3</figref>, the short term distribution 24 and long term distribution 26 are recursively estimated weighted distributions of all input data 12 received. They require no explicit storage of the contents of the data stream 12. The short term distribution 24 is a model that weights more recent information in the data stream 12 more highly than the model of the data in the long term distribution, and thus responds more rapidly to changes.

0023A preferred form of the method of the invention is shown as 30 in <figref idref="f0002">Figure 3</figref>. The process 30 is undertaken each time the system 10 receives 31 an input. The distribution maintenance engine 14 maintains the short 24 and long 26 term distribution models updating them when an input 12 is received 31 from the data stream.

0024The optimal periods required for the short 24 and long 26 term models depend on the particular application of the invention and the timescales of the variations that occur within it. For example, when applied to the problem of telecommunications fraud detection, effective periods for the short 24 and long 26 term models have been found to be around 100 and 250 calls respectively. The invention overcomes the limitations of simple threshold-based change detection systems by explicitly modelling the distributions 24 and 26 that describe the data stream 12 and its contents. The short term 24 and long 26 term probability distribution models can each be joint or conditional of several data streams 12, thus permitting the detection of changes in the relationship between the streams and the information they carry.

0025Additionally, conditional short term and long term probability distribution models allow changes in patterns of sequences of information relating to one or more data streams 12 to be detected. The use of distributions 24 and 26 also means that the system 10 is able to handle arbitrary combinations of continuous and discrete data/information, and hence numeric and categorical information. For example, while searching for telecominunications fraud, it may be useful monitor the cost of calls, and to categorise them as local, national, international, etc.

0026It is noted that the term probability distribution is usually used when modelling discrete values like integers and densities are usually used when modelling continuous values like real numbers. The present invention can be used for either and throughout this specification the term distribution is intended to encompass densities as well, unless this is clearly not intended.

0027The different rates of response of the short term 24 and long term 26 distributions to changes in a data stream 12 means a difference between the distributions 24 and 26 that can be measured 33 by suitable means 16 using any of a number of standard statistical measures, such as χ² or the Kullback-Leibler divergence. A measure of the difference between the short 24 and long 26 term distributions is passed to the alert generation engine 18. The alert generation engine 18 checks 34 whether the difference exceeds an adaptive alert threshold, which is maintained by the threshold adaptation engine 20. The adaptive alert threshold is formed from two components: a predictability measure, which is a moving average of the difference between the short 24 and long 26 term distributions; and a variability measure, which is a moving average of the absolute difference between the difference between the short 24 and long term 26 distributions and the predictability measure. Other means of adapting the threshold, and the statistics from which it is formed, can also be considered. For example, in order to improve robustness, the predictability measure may not be updated when the measured difference between the short 24 and long 26 term distributions exceeds some function of the variability measure.

0028As with the short 24 and long 26 term distributions, the predictability and variability measures (and hence the alert threshold) are updated 35 every time the invention receives an input 12, provided that an alert 22 is not generated. The adaptive alert threshold is typically formed from the sum of the predictability measure and some multiple of the variability measure, where the multiple is a sensitivity parameter that is used to control the invention's sensitivity (note that increasing the sensitivity parameter reduces the invention's sensitivity). Other functions may be used to adapt the alert threshold. For example, setting the threshold to equal the predictability plus a multiple of the variability squared may give improved performance in some applications.

0029If the difference between the short 24 and long 26 term distributions exceeds the adaptive alert threshold, then the alert generation engine 18 generates 36 an alert signal 22. The alert signal 22 can contain any of the information available within the system, information presented to the system in the data stream 12, and functions of either and/or both. Typically, however, it contains a propensity measure, which indicates the severity of the change, and detailed information about its exact nature, which can be obtained by examining the differences between the short term 24 and long term 26 probability distributions. One of the simplest forms of propensity measure can be derived by dividing the difference between the measure of the difference between the distributions 24 and 26 and the alert threshold by the variability measure. More sophisticated schemes scale the propensity according to the number of alerts recently generated, and non-linearly processes it so that it always lies between zero and one, in order to improve its interpretability. Detailed information about the exact nature of a change in the data stream can be extracted by measuring the contributions of points on the short 24 and long 26 term distributions to the difference that was measured between them. For example, assuming that the χ² measure was used by means 16, the contribution of a point on the short 24 and long 26 term distributions to the measured difference would be given by the χ² measure between the distributions 24 and 26 at the point, divided by the total difference for all points that are sampled. In some circumstances, such as when the distributions 24 and 26 are represented by histograms, the contributions are the same for all points within each histogram bin, and hence only need to be computed once within each. The alert 22 signal may also contain additional information, such as a representation of the long term distribution 26, which provides a description of the statistics of the data stream 12 just prior to the change that generated the alert.

0030The distribution maintenance engine 14 also unwinds 37 an update that caused an alert 22 to be generated, in order to restore the distributions 24 and 26 to their pre-updated state. Restoring the distributions in this way stops the invention adapting to abnormal changes in behaviour. This prevents it from becoming used to any abnormal behaviour that may be sustained in the short term, guaranteeing that it continues to generate alerts 22 throughout, and ensuring that it does not generate alerts 22 when normal behaviour is resumed.

0031It is not uncommon for changes in data streams to occur that are of sufficient magnitude to cause any change detection algorithm to generate alerts, but which are not target events of interest to the user of the system 10. For example, when the system 10 is used to monitor the behaviour of a subscriber to a telecommunications service, it may generate alerts if they suddenly start to make international calls, perhaps because a friend or relative has moved to another country. When these alerts are investigated by the telecommunications service provider, they will be found to be innocuous, and the provider will expect to be able to instruct the change detection algorithm that no significant change has occurred, and that it should stop generating alerts for the new behaviour and should adapt to it. The system 10 provides such a facility by maintaining an estimate of the amount by which its sensitivity would have needed to have been adjusted in order not to have generated any alerts since the last time it was instructed to adapt. One simple but effective way to produce such an estimate (which works well when the alert threshold is a linear function of the sensitivity parameter, the difference between the short 24 and long 26 term distributions, and the predictability and variability measures) is to maintain a sensitivity adjustment estimate that is increased by an additive constant each time an alert 22 is generated, and decays exponentially with each input the invention receives. The additive constant and the rate of decay can be determined empirically to give a good trade-off between over-adjustment and the resulting under-sensitivity, and under-adjustment and inadequate desensitisation. When the system 10 receives an adapt signal from the user, it increases its sensitivity parameter by an amount equal to the largest value taken by the sensitivity adjustment estimate since the last time it was instructed to adapt.

0032Two additional features are supported by the system 10. The first feature allows the method 30 to be applied to data streams 12 that exhibit systematic (i.e. non-random) variation over time. Such variation is frequently encountered when searching for signs of fraud in telecommunications networks. For example, businesses make many calls during weekdays, but few at weekends; conversely, domestic users may make more calls at weekends when some network operators offer lower tariffs. The naïve application of change detection algorithms to these types of data streams can result in large numbers of false alerts as the algorithms respond to the changes in behaviour that occur at the start and end of each weekend. To overcome this problem, the system 10 may support the specification of a lead in period, during which alerts 22 cannot be generated, thereby adapting to all inputs 12. At the end of the lead in period, the number of inputs that were passed to the distribution maintenance engine 14 is recorded and used as the period for updates to the predictability and variability measures, if the number of inputs is larger than the default (which is typically 250 for telecommunications applications). This forces the threshold adaption engine 20 to adjust the alert threshold more slowly than it would by default, with the result that systematic variations in the input stream 12 similar to those that occurred during the lead in period do not generate alerts 22.

0033The alert generation engine 18 also allows alerts 22 generated by inputs 12 that are above a configurable lower percentile or below a configurable upper percentile (both estimated from the long term probability distribution 26) to be suppressed. This feature is useful when examining records of call costs for fraudulent behaviour which requires alerts to be produced for increases in the number of medium and high cost calls, which are strongly indicative of fraud, but not for increases in the number of low cost calls, which usually are not. The means by which the percentile estimates are obtained depend on the representations of the distributions.

0034For example, pseudo-code for calculating a percentile from a histogram-based representation is given below: <pre listing-type="program-listing"> current_bin = 1; accumulated _long_term_probability = 0; while (accumulated_long_term_probability + long_term_probability_of (current_bin) &lt; required_percentile ) accumulated_long_term_probability=accumulated_long_term_probability+long_term_probability_of (current_bin); current_bin = current_bin + 1; end bin_that_contains_percentile = current_bin; percentile = min_of (bin_that_contains_percentile) + (max_of (bin_that_contains_percentile) - min_of (bin_that_contains_percentile) ) *(required_percentile - accumulated_long term_probability) / (long_term_probability_of (bin_that_contains_percentile) );</pre>

0035Most terms in the pseudo-code are self-explanatory, except for the min_of(.) and max_of(.) functions, which return the smallest and largest values of the input associated with the bin specified in their arguments. The bins are specified by contiguous integers starting at one, and it is assumed that they are ordered according to the ranges of the input that they are associated with - that is, bin one is associated with the smallest values of the input. The algorithm first finds the bin that contains the percentile by looking for the bin associated with the smallest values of the input for which the sum of its long term probability and those of all bins associated with smaller values of the input are greater than the required percentile. The specific location of the percentile within the bin is then calculated by assuming that the probability mass associated with the bin is distributed uniformly within it.

0036As an example of how the invention is applied in practice, consider an instance of the system 10 that models the distribution of telephone calls across ten different categories. These categories could indicate different groups of 'hot' destinations (telephone numbers frequently called by fraudsters). For example, one category could include all telephone numbers in the Caribbean and Pakistan, while another may include all numbers in Brazil and Argentina. To illustrate the effect of a single call on the invention, Table 1 shows the probabilities associated with each of ten categories in the short term 24 and long term 26 distributions within the invention. <tables id="tabl0001" num="0001"><table frame="all"><title>Table 1</title><tgroup cols="4"><colspec colnum="1" colname="col1" colwidth="18mm" /><colspec colnum="2" colname="col2" colwidth="36mm" /><colspec colnum="3" colname="col3" colwidth="36mm" /><colspec colnum="4" colname="col4" colwidth="42mm" /><thead><row><entry align="center" valign="middle">Category</entry><entry align="center" valign="middle">Short term probabilities</entry><entry align="center" valign="middle">Long term probabilities</entry><entry align="center" valign="middle">Per category χ² divergence</entry></row></thead><tbody><row><entry align="center">1</entry><entry align="char" char="." charoff="5">0.1326</entry><entry align="char" char="." charoff="5">0.1226</entry><entry align="char" char="." charoff="5">0.0008</entry></row><row><entry align="center">2</entry><entry align="char" char="." charoff="5">0.0794</entry><entry align="char" char="." charoff="5">0.1012</entry><entry align="char" char="." charoff="5">0.0047</entry></row><row><entry align="center">3</entry><entry align="char" char="." charoff="5">0.0032</entry><entry align="char" char="." charoff="5">0.0394</entry><entry align="char" char="." charoff="5">0.0332</entry></row><row><entry align="center">4</entry><entry align="char" char="." charoff="5">0.1429</entry><entry align="char" char="." charoff="5">0.1311</entry><entry align="char" char="." charoff="5">0.0011</entry></row><row><entry align="center">5</entry><entry align="char" char="." charoff="5">0.0774</entry><entry align="char" char="." charoff="5">0.0979</entry><entry align="char" char="." charoff="5">0.0043</entry></row><row><entry align="center">6</entry><entry align="char" char="." charoff="5">0.1071</entry><entry align="char" char="." charoff="5">0.0882</entry><entry align="char" char="." charoff="5">0.0040</entry></row><row><entry align="center">7</entry><entry align="char" char="." charoff="5">0.1378</entry><entry align="char" char="." charoff="5">0.1247</entry><entry align="char" char="." charoff="5">0.0014</entry></row><row><entry align="center">8</entry><entry align="char" char="." charoff="5">0.1604</entry><entry align="char" char="." charoff="5">0.1613</entry><entry align="char" char="." charoff="5">0.0000</entry></row><row><entry align="center">9</entry><entry align="char" char="." charoff="5">0.1285</entry><entry align="char" char="." charoff="5">0.1034</entry><entry align="char" char="." charoff="5">0.0061</entry></row><row><entry align="center">10</entry><entry align="char" char="." charoff="5">0.0307</entry><entry align="char" char="." charoff="5">0.0302</entry><entry align="char" char="." charoff="5">0.0000</entry></row><row><entry namest="col1" nameend="col3" align="right">Total χ² divergence</entry><entry align="char" char="." charoff="5">0.0555</entry></row></tbody></tgroup></table></tables>

0037Assuming that a call is made to a destination that falls in the first category, the system 10 updates 32 the probabilities in its short term 24 and long term 26 distributions to the values shown in Table 2. <tables id="tabl0002" num="0002"><table frame="all"><title>Table 2</title><tgroup cols="4"><colspec colnum="1" colname="col1" colwidth="18mm" /><colspec colnum="2" colname="col2" colwidth="36mm" /><colspec colnum="3" colname="col3" colwidth="36mm" /><colspec colnum="4" colname="col4" colwidth="42mm" /><thead><row><entry valign="top">Category</entry><entry valign="top">Short term probabilities</entry><entry valign="top">Long term probabilities</entry><entry align="center" valign="top">Per category χ² divergence</entry></row></thead><tbody><row><entry align="center">1</entry><entry align="char" char="." charoff="5">0.1412</entry><entry align="char" char="." charoff="5">0.1261</entry><entry align="char" char="." charoff="5">0.0018</entry></row><row><entry align="center">2</entry><entry align="char" char="." charoff="5">0.0786</entry><entry align="char" char="." charoff="5">0.1008</entry><entry align="char" char="." charoff="5">0.0049</entry></row><row><entry align="center">3</entry><entry align="char" char="." charoff="5">0.0032</entry><entry align="char" char="." charoff="5">0.0392</entry><entry align="char" char="." charoff="5">0.0331</entry></row><row><entry align="center">4</entry><entry align="char" char="." charoff="5">0.1415</entry><entry align="char" char="." charoff="5">0.1306</entry><entry align="char" char="." charoff="5">0.0009</entry></row><row><entry align="center">5</entry><entry align="char" char="." charoff="5">0.0766</entry><entry align="char" char="." charoff="5">0.0975</entry><entry align="char" char="." charoff="5">0.0045</entry></row><row><entry align="center">6</entry><entry align="char" char="." charoff="5">0.1060</entry><entry align="char" char="." charoff="5">0.0879</entry><entry align="char" char="." charoff="5">0.0038</entry></row><row><entry align="center">7</entry><entry align="char" char="." charoff="5">0.1364</entry><entry align="char" char="." charoff="5">0.1242</entry><entry align="char" char="." charoff="5">0.0012</entry></row><row><entry align="center">8</entry><entry align="char" char="." charoff="5">0.1588</entry><entry align="char" char="." charoff="5">0.1606</entry><entry align="char" char="." charoff="5">0.0000</entry></row><row><entry align="center">9</entry><entry align="char" char="." charoff="5">0.1272</entry><entry align="char" char="." charoff="5">0.1030</entry><entry align="char" char="." charoff="5">0.0057</entry></row><row><entry align="center">10</entry><entry align="char" char="." charoff="5">0.0304</entry><entry align="char" char="." charoff="5">0.0301</entry><entry align="char" char="." charoff="5">0.0000</entry></row><row><entry namest="col1" nameend="col3" align="right">Total χ² divergence</entry><entry align="char" char="." charoff="5">0.0558</entry></row></tbody></tgroup></table></tables> These have been computed assuming that the period of the short term model 24 is 100 calls, and that of the long term model 26 is 250 calls, and hence have increased the probability of category one by 1/00 and 1/250 in the short 24 and long 26 term distributions respectively, and renormalizing each distribution. The new χ² value is computed 33 in the standard way from the sum of the squared differences between the short 24 and long 26 term probabilities of each category, divided by the long term probabilities 26, and is given at the bottom Table 2.

0038Assuming that the predictability measure is 0.0612, the variability measure is 0.0048, and the sensitivity parameter is 5.0 (which has been determined empirically as previously described), the alert threshold is 0.0612+5.0x0.0048 = 0.0852. Since the new value of χ² is less than the alert threshold, no alert is generated, and the predictability and variability measures are updated 35 to be (1-α)×0.0612+α×(0.0558-0.0612) = 0.0609 and (1-α)×0.0048+α×|0.0558-0.0609| ≈ 0.0048 respectively (35), where α is the reciprocal of the period of the long term distribution (i.e. 1/250) and |.| returns the magnitude of its argument.

0039If the new χ² value had been larger than the alert threshold an alert would have been generated 36, the predictability and variability measures would not have been updated and the probabilities in the short 24 and long 26 term probability distributions would be restored 37 to their previous values by subtracting 1/(1+100) and 1/(1+250) from the probabilities for category one and renormalizing. It should be emphasised that this example is a specific realisation of the invention and most of the details are peculiar to the realisation and not fundamental to the invention itself.

0040A simulated application of the invention is shown in <figref idref="f0001">Figure 4</figref>. Each continuous block of alerts is shown as a vertical black bar. The present invention produces a large number of change alerts in response to the change in simulated measurements and it is clear that the invention produces no change alerts outside the region of the pathological behaviour, and many within it, thus giving good discrimination. Note that although the change in the data stream shown in <figref idref="f0001">Figure 4</figref> can be detected by placing a threshold on the variance of the spring length measurements, there may be other modes of failure or combinations thereof that cause changes that cannot be detected in ways that can easily be specified a priority.

0041The skilled addressee will realise that the present invention has the following advantages: <ul id="ul0006" list-style="bullet"><li>The present invention can discriminate between behaviours that cannot be distinguished using thresholds and summary statistics due to the behaviour of the date stream being modelled by short term and long term probability distributions;</li><li>The distributions and densities are estimated recursively, and hence the invention requires no explicit storage of information from the data stream;</li><li>The invention can model multivariate and conditional properties of a data stream, making it possible to detect changes in interactions between different streams, and changes in the patterns of sequences that occur within a stream;</li><li>The invention measures the rate of change of behaviour by measuring the divergence between the short term and long term distributions;</li><li>The invention maintains statistics concerning the typical rate of change in the data stream that are used to compute an adaptive threshold that, in turn, controls the generation of alerts;</li><li>The invention generates alerts when the measured divergence exceeds an adaptive threshold;</li><li>The alerts generated by the invention contain detailed information about its cause - i.e. the way in which the stream's behaviour has changed;</li><li>The invention estimates the amount by which its internal sensitivity would need to be adjusted in order to produce an alert rate of around one percent; and</li><li>The invention can be forced to adjust its internal sensitivity by an externally specified amount, or according to its own internal estimate.</li></ul>

0042The skilled addressee will realise that the present invention is applicable to detecting changes in the properties of data streams originating from various sources. Even though the example provided of the present invention was for telecommunications fraud detection, the present invention may be suitable for numerous other applications where data contains information about some entity or process. Examples include but are not limited to: credit card fraud detection system; computer network security; process, fault and condition monitoring; environmental monitoring systems; and health monitoring systems. Such modifications and variations as would be readily apparent to the skilled addressee to enable the present invention to detect changes in the properties of data streams from various sources are intended to fall within the scope of the present invention, the nature of which is to be determined from the foregoing description and appended claims.