Research Article  Open Access
M. Borowski, S. Siebig, C. Wrede, M. Imhoff, "Reducing False Alarms of Intensive Care OnlineMonitoring Systems: An Evaluation of Two Signal Extraction Algorithms", Computational and Mathematical Methods in Medicine, vol. 2011, Article ID 143480, 11 pages, 2011. https://doi.org/10.1155/2011/143480
Reducing False Alarms of Intensive Care OnlineMonitoring Systems: An Evaluation of Two Signal Extraction Algorithms
Abstract
Onlinemonitoring systems in intensive care are affected by a high rate of false threshold alarms. These are caused by irrelevant noise and outliers in the measured time series data. The high false alarm rates can be lowered by separating relevant signals from noise and outliers online, in such a way that signal estimations, instead of raw measurements, are compared to the alarm limits. This paper presents a clinical validation study for two recently developed online signal filters. The filters are based on robust repeated median regression in moving windows of varying width. Validation is done offline using a large annotated reference database. The performance criteria are sensitivity and the proportion of false alarms suppressed by the signal filters.
1. Introduction
In intensive care, the condition of a patient is supervised by onlinemonitoring systems which measure several vital signs with a high sampling rate of up to one observation per second. These devices produce different types of alarms to alert the clinical staff. Most frequently these alarms are socalled simple threshold alarms, which are given when the measured values of a vital sign lie outside of specified alarm limits. The high rate of irrelevant threshold alarms, caused by artefacts and short fluctuations, is a wellknown problem and shown in several studies; see [1] or [2], for instance. High false alarm rates may lead to a dangerous desensitization of medical staff toward alarms and, thus, dramatically reduce the effective sensitivity of the entire alarm system [3–5]. Moreover, alarm limits may be set inadequately wide or alarms disabled completely. Obviously, onlinemonitoring systems are in need of improvement with respect to the high rates of false positive alarms.
Several approaches for improved alarm systems have been proposed. For instance, median filters can be used to eliminate noise and artefacts [6]; in [7], a method based on control charts is developed to detect the onset of changes in systolic blood pressure; in [8], a preprocessing algorithm is proposed which provides a basis for an online trend extraction methodology [9]. An overview of alarm algorithms in critical care monitoring is given in [10].
One approach to decreasing the false alarm rate is statistical signal extraction or filtering. Assume that the data consist of a true but unknown relevant signal overlaid with irrelevant noise and outliers. The signal can then be extracted online, meaning that the signal is extracted sequentially with each new incoming measurement. Then, the alarm limits can be compared to the online extracted signal instead of to the raw measured data, leading to fewer threshold alarms. Figure 1 illustrates this approach.
Time series from clinical online monitoring are not stationary but show level shifts, enduring and changing trends, a high level of noise, and are corrupted by (patches of) outliers. A suitable approach for online signal extraction, given such difficult data, is robust linear regression in a moving time window [11–13]. Then, the level of the regression line at the central or, alternatively, rightmost window position is used as the signal estimation. Using the level at the rightmost window position has the advantage that in an online application, the signal is estimated at the current time point, that is, without a time delay. We will, therefore, consider this type of moving window regression.
Since outliers appear frequently, robust regression methods should be applied. The robust Repeated Median (RM) regression [14] has proven to be a suitable candidate for the special situation of clinical online monitoring [11, 13]. It outperforms other robust regression methods with respect to robustness, computing time, and efficiency, in the sense of low bias and variability. However, as for any moving window technique, the RM signal estimation is strongly affected by the window size: large windows induce “smooth” signal estimations with little variability whereas small windows lead to signal estimations that are closer to the data. The adaptive online RM (aoRM) [15] chooses the window width for the RM automatically; as long as the data are “stable”, the window width gradually grows but when a structural change, for example, a level shift, occurs, the window width is set to a predetermined minimum value. The aoRM is enhanced to a filtering procedure for multivariate time series, namely, the adaptive online Trimmed Repeated MedianLeast Squares (aoTRMLS) filter [16]. This procedure factors in local crosscorrelations (e.g., systolic and diastolic arterial pressure are highly correlated) in order to improve the filtering outcome.
A good filtering procedure suppresses “many” irrelevant (i.e., false) alarms while suppressing only “a few” relevant (i.e., true) alarms. In this paper, we will investigate aoRM and aoTRMLS with respect to these criteria. The investigation is done using recorded monitoring data from an intensive care unit; all given alarms have been annotated retrospectively by an experienced physician as true or false. The aoRM and aoTRMLS filters are applied offline to the data in order to determine the number of suppressed false alarms and the number of correctly reproduced true alarms. The corresponding performance criteria are sensitivity and false alarm reduction rate.
In the next section, we explain the aoRM and aoTRMLS algorithms. Section 3 describes the study setting, focussing on data collection and annotation. In Section 4, we illustrate the considered performance criteria and the way we determine them. Section 5 presents the results of our analysis. A summary and an outlook are given in Section 6.
2. Signal Extraction Algorithms
The monitoring of vital signs (e.g., heart rate, oxygen saturation, and systolic, diastolic, and mean arterial pressure) leads to a variate time series , where is the time index, indicating each second, for instance. The online signal extraction approach is based on the assumption that the observed data can be decomposed into a true but unknown signal which is overlaid with noise and outliers Here, denotes the dimensional signal at time . The noise term is , where are errors with zero median and timedependent variances . The errors may be correlated, that is, possibly for any . An outlier generating mechanism that produces impulsive spiky noise is denoted by .
Both the aoRM and aoTRMLS can be used to extract the signal vector online: with each incoming new observation at time , a new signal estimation output is produced. The aoRM is developed for univariate time series, meaning that it must be applied separately to each univariate component of a variate time series. On the contrary, the aoTRMLS filter is developed for multivariate time series. It takes into account local dependences between variables and ensures robustness against outliers regarding the local covariance structure. Since aoTRMLS is based on aoRM, we introduce the univariate aoRM first.
2.1. The Adaptive Online Repeated Median Filter (aoRM)
In [15], it is assumed that the underlying signal of a univariate time series is approximately linear in a small moving window of length where and . Here, is the level of the straight line at the rightmost time point , and is the associated slope in the time window. The rightmost window time point corresponds to the current time point in the online case. Then, the RM estimates of and in (2) are The aoRM estimates the signal by the RM, after the window width has been adapted to the current data situation. The window width adaption is done automatically when a new observation comes in, that is, at each time point . The adapted window width at time is, therefore, denoted by . After is determined, the signal is estimated by RM regression in the time window .
In [15], it is demanded that . The minimum width guarantees robustness against a certain number of outliers while the maximum width limits the computing time. Both and must be set beforehand by the user. A flow chart for the complete aoRM algorithm is shown in Figure 2. The main step of the aoRM algorithm is the decision of whether or not the RM fit in the time window , that is, on the window sample , is adequate. This decision is made by means of a test procedure which is based on the fact that an RM regression results in an equal number of positive and negative residuals. If the positive and negative residuals are not balanced for the rightmost residuals within the time window of width , the RM fit is regarded as inadequate. The input parameter should be chosen such that . For more details regarding the test procedure and the choice of , see [15].
The aoRM filter estimates the signal at the rightmost or current time point, meaning that the signal is extracted without relevant time delay. (Its computing time depends on the chosen input parameters. The R package robfilter [17] provides a function of the aoRM; see Section 5.2. We applied this aoRM function using a 2.3 GHz computer with 2 GB RAM and obtained a mean computing time of 0.007 seconds for one iteration, resp., time point.) However, estimating the signal by the level at the right end of the regression line implies that the signal estimates possibly deviate distinctly from the data, especially when level shifts occur. Then, signal estimates “overshoot”; that is, they leave the range given by the window observations, see Figure 3(a). It shows a generated time series (dotted) with upwards and downwards level shifts at time points and and the corresponding aoRM signal estimation time series (solid). Around time and , aoRM signal estimations overshoot. Those overshoots are crucial in our context since a sudden change in the data may cause the signal estimations to cross an alarm limit although the measurements do not, as can be seen in the figure. That is, an aoRMbased alarm system could theoretically cause more false alarms than a system based on raw measurements.
(a)
(b)
In order to prevent signal estimations from overshooting and causing false alarms, it is suggested in [15] to restrict the signal estimation to a value within the range of the most recent observations . (Note that this is the subsample which is used for the decision of whether or not the RM fit is adequate.) Defining the minimum and maximum of the most recent observations by the restricttorange rule is The effect of this rule can be seen in Figure 3(b). The signal estimations around time and are “pulled back” to the measurements. They do not violate the alarm limits, and, therefore, unjustified alarms are prevented.
Figure 3 also shows that aoRM signal estimations trace changes in the data time series with a certain time delay. This inert reaction to sudden data changes is due to the RM’s robustness against outliers. It has a finite sample replacement breakdown point of 50% [18]. That is, in order for a patch of levelshifted observations to affect the RM, the patch must consist of more than half of the window observations. Since a sudden change in the data implies that the window width is set down to , the choice of is crucial: it defines the distinction between outlier patches and level shifts as well as the tracing delay for structural changes in the data, which is approximately time points. In Figure 3, , so the tracing delay is approximately 20 time points.
The aoRM can be used for filtering multivariate time series by applying it separately to each univariate component time series. However, it does not account for dependences between the variables. The aoTRMLS [16] is developed as a multivariate enhancement of the aoRM, which uses information given by the local covariance structure.
2.2. The Adaptive Online Trimmed Repeated MedianLeast Squares Filter (aoTRMLS)
The aoTRMLS arises out of a combination of the aoRM and the multivariate TRMLS [19]. Similarly to (2), it is assumed that each component of the underlying variate signal is approximately linear in a short moving time window where and . Here, is the level vector at the rightmost or current time point , and the vector of the slopes; see (2).
The aoTRMLS filter can be used to estimate and in (6). At each time , the aoTRMLS algorithm searches for an overall window width which is adequate for all of the variables. That is, it searches for the greatest such that the linear approximation (6) with is adequate. This search is done using the aoRM window width adaption principle. Thereafter, within the time window specified by , the signal vector is estimated by means of multivariate TRMLS regression which is explained later on. The aoTRMLS algorithm is as follows: (0)initialization: wait until observations , , are present, (1)apply the aoRM window width adaption procedure to each of the window samples in order to obtain appropriate individual window widths , and set the overall window width ; (2)perform TRMLS regression on the multivariate sample , where , and store the TRMLS signal estimation ,(3) set , and update the window: set and go to step (1).
After the overall window width is determined at step (1), the signal is estimated by means of TRMLS regression at step (2). An outline of the TRMLS regression algorithm is as follows (for simplicity, we set ).
First, RM regression is performed separately on each of the window samples . Then the RM residuals are regarded as a multivariate sample , where . The local error covariance matrix is then estimated on this residual sample using a robust estimator proposed in [20]. The estimate is utilised to detect residual vectors that are outliers regarding the local covariance structure, that is, residual vectors , , with where is an adequate upper bound. (For more details see [16, 19] or [20].) Then, observation vectors , which correspond to outlying residual vectors , are removed from the window sample. Finally, a multivariate LS regression is performed on the outlierfree window sample, and the levels at the right end of the LS regression lines build the signal estimation vector.
The aoTRMLS yields robust but also efficient signal estimations since the signal vector is finally estimated by means of LS regression. Just like aoRM, aoTRMLS estimates the signal vector at the right end point of the moving time window. Hence, the restricttorange rule (5) is also recommended for aoTRMLS. Furthermore, for data that exhibit a known correlation structure, it is suggested to apply the aoTRMLS filter not to the whole variate time series but separately to blocks that consist of highly positively correlated variables [16]. For instance, systolic, mean, and diastolic blood pressure are highly positively correlated and, therefore, can be combined into a correlation block. A block wise application of aoTRMLS improves the window width adaption, leading to smoother signal extraction time series.
We apply both aoRM and aoTRMLS retrospectively to recorded onlinemonitoring data from intensive care in order to evaluate their ability to suppress false and to reproduce true alarms. In the following, we describe the study setting, including data collection and annotation.
3. Study Setting
In accordance with the declaration of Helsinki, the study was approved by the Ethics Committee of the University of Regensburg. The data were collected at an intensive care unit at the University Hospital Regensburg [21]. Only adult patients with continuous monitoring of at least invasive arterial blood pressure, heart rate, and oxygen saturation were included into the study. The deployed monitoring system was an Infinity Monitor by Dräger Medical, Lübeck, Germany. Data acquisition from the Infinity Monitoring System took place using the special software eData by Dräger Medical, Lübeck, Germany. The data were recorded at a 250 Hz sampling rate, labeled with a time index, and stored. The reference data that we consider were extracted from these data at a sampling rate of one per second and stored in text files. This reference database includes (i)numerical measurements of the vital signs, (ii)all monitoring system alarms with corresponding time and alarm message (e.g., “heart rate lower limit violation”), (iii)the alarm limits that were set by the medical staff, (iv)information when alarms were deactivated (“alarms off” periods).
For more details about the database, see Section 5.1.
All alarm situations were retrospectively annotated by an experienced physician by means of graphical representations of the collected data from the monitoring system in combination with video recordings showing the patient and the screens of the monitoring system. Since the data contain the alarm time points, the physician could wind the video tape and watch the patient and the screens of the monitoring system at these alarm time points. Then, by means of a specially developed JavaScript program [22], each alarm situation was assessed in terms of whether it was clinically relevant and whether it was technically true.
In this study, an alarm is regarded as technically true if it is based on a correct measurement or if the monitoring system correctly recognizes a technical problem and gives a technical alarm. All other alarms are technically false. Furthermore, a situation is defined as alarm relevant or true, if it implicates a diagnostic or therapeutic decision or the correction of a technical problem. A technically true alarm is annotated as advisory if it does not require immediate action; that is, it is not alarm relevant but judged to be helpful. All other alarms are annotated as not alarm relevant or false, respectively. A technically false alarm is always a false alarm.
4. Performance Criteria
Our aim is to evaluate the “new” aoRM and aoTRMLS alarm systems and compare their performances to that of the “old” alarm system based on raw measurements. In this context, an alarm system is a diagnostic tool: an alarm corresponds to the diagnosis “alarm relevant situation”. Common performance criteria for diagnostic methods are sensitivity (SE) and specificity (SP). In our context, SE is the conditional probability of an alerting monitoring system given that the situation is actually alarm relevant, and SP is the conditional probability of a nonalerting monitoring system, given that the situation is actually not alarm relevant While SE assesses the performance with respect to the detection of alarm relevant situations, SP quantifies the liability of the alarm system to produce false alarms. Both SE and SP of a diagnostic method should be large; at best, both equal to 1.
The common approach to estimate an alarm system’s SE is to determine the ratio of detected alarm relevant situations and all alarm relevant situations. In our study, alarms given by the old system (positive alarms) were annotated as either true positive or false positive. Nongiven alarms, that is, negative alarms, including false negative alarms, do not occur. Conversely, this means that the old system detects all alarm relevant situations and, therefore, has 100% SE. However, this is not de facto but follows from the study design.
The SE of a new alarm system based on signal filtering can be estimated by the ratio of the number of alarm relevant situations that are detected, that is, reproduced correctly, and the number of all alarm relevant situations. A situation is regarded as detected by a new system if the signal estimations violate the alarm limits close to the time when the alarm was given by the old system. A detailed description is given later on.
Similar to SE, SP is usually estimated by the ratio of the number of true negative alarms and the number of not alarm relevant situations. However, due to the study design, there is no information about these numbers. Therefore, in [23], an alternative approach for estimating the SP of an onlinemonitoring alarm system is used: defining as the highest expected number of false positive alarms and FP as the actual number of false positive alarms, the difference is an estimate for the number of true negative alarms. This number is then divided by the worst case number to get an alternative estimation of SP: For SP estimation as proposed in [23], the analyst has to decide whether or not to count in alarms of the new system that occur during “alarms off” periods. In these periods, alarms were deactivated by the clinical staff—even if the measurements violated the alarm thresholds, no alarm was given. Since alarms are often deactivated as a reaction to threshold alarms, a new system would be favoured if it cannot produce alarms in these periods; see Figure 4. It shows a part of a time series of mean arterial blood pressure measurements (dotted) and the corresponding signal estimations by the simple RM with window width (solid). Due to the RM’s robustness against outliers, its signal estimations react to changes in the data with a tracing delay of approximately time points, as explained in Section 2. Since the staff deactivated all alarms at , only four seconds after the false alarm at occurred, the signal extraction time series violates the lower alarm limit within the “alarms off” period. If “alarms off” periods were not considered in the analysis, the false alarm at would be regarded as suppressed by the new system although its signal estimations violate the lower alarm limit. Thus, in order to not favour the new system, one must include alarms of the new system that would be given during “alarms off” periods. However, during those periods the old system cannot produce threshold alarms, but the new system can. Moreover, each alarm of the new system is regarded as false if it does not belong to a situation annotated as relevant. That is, within “alarms off” periods, each limit violation of a new system would be regarded as false alarm. In a nutshell, due to the given study design and data basis, for SP estimation the analyst must decide in favour of the old or new system.
Due to these considerations, we refrain from determining SP. Instead, we consider the ratio of false alarms, which are suppressed by the new system, to all false alarms, denominated as false alarm reduction rate (FARR) of the new system: Similarly to the detection of true alarms, a false alarm is regarded as suppressed by the new system if its signal estimations do not violate the alarm limits close to the time when the false alarm occurs.
Obviously, the old alarm system has a FARR of 0%; that is, the new system cannot be worse with respect to this criterion. Moreover, there might be alarms of the new system that are not regarded. Hence, alarms of the new system should rather be considered independently from false alarms of the old system.
However, the restricttorange rule (5) guarantees that the estimated signal only violates the upper (e.g.) alarm limit, if at least one of the most recent measurements also violates the upper limit. Thus, an aoRM/aoTRMLS alarm system cannot generate more alarms than the old system. However, the old system already includes a simple algorithm to suppress false alarms. It gives an alarm only if observations lie outside the alarm limits for a certain time span, the alarm validation time; see Table 1. For instance, an alarm caused by a too high systolic blood pressure is not given until four consecutive measurements are above the upper alarm limit whereas a single heart rate measurement below the lower alarm limit causes an alarm immediately.

Because of alarm validation times, aoRM/aoTRMLS signal estimations can violate an alarm limit (and hence produce an alarm) although the old system has not given an alarm beforehand. However, such cases are very rare and occur only if measurements fluctuate around an alarm limit. In those cases, an alarm is rather helpful: either the alarm limits are set too narrow, or the patient’s condition deteriorates.
Due to these considerations and taking into account the fact that SP does not allow for a fair comparison, we conclude that FARR is a more sensible performance criterion in our situation.
In the following, we explain how we estimate SE and determine the FARR of a new alarm system based on signal extractions. For simplicity, SE now denotes the estimated sensitivity.
As mentioned above, SE of a new alarm system can be determined by the ratio of detected true alarms and all true alarms. A true alarm can be regarded as detected if the signal estimations violate the alarm limits close to the time when the true alarm occurred. Hence, we consider all true alarm time points , , which correspond to true threshold alarms of the old system. If signal estimations violate the alarm limits in a certain time range in proximity to , the referring true alarm is regarded as detected by the new alarm system.
In contrast to SE, FARR is determined using all alarm time points which are annotated as false, that is, . If aoRM/aoTRMLS signal estimations do not violate the alarm limits in a certain time range in proximity to , the corresponding false alarm is regarded as suppressed by the new alarm system. Then, the FARR of the new system is the ratio of suppressed false alarms to all false alarms.
In order to specify the range around true alarm time points , we build true alarm intervals . The detection tolerance time was chosen by the physicians involved in the study. For ART.M, ART.S, and SpO_{2}, this time is , and for HR it is seconds. However, we decided to build true alarm intervals which merely include time points after an alarm time point, that is, , since aoRM and aoTRMLS were developed for suppressing false alarms. They do not forecast alarm situations but react to the data which is reflected in the tracing delay. Thus, including time points left from into the alarm interval makes little sense.
Similar to true alarms, the range after a false alarm is specified by a false alarm interval . The suppression tolerance time must be chosen such that it is greater than the tracing delay of the signal estimations. This prevents false alarms from being assessed as suppressed by the new system just because its signal estimations violate the alarm limits too late. The largest minimum window width is in our study, which corresponds to a tracing delay of approximately 45 seconds. Hence, we think that is an ample suppression tolerance time.
The determination of SE and FARR based on true and false alarm intervals does not work without further complications: we have to bear in mind that some false alarms appear close to true alarms so that the corresponding alarm intervals overlap. Thus, there might be inconsistencies since a new alarm could be regarded as detected true alarm and as not suppressed false alarm. Since the detection of true alarms is more crucial than the suppression of false alarms, we simply exclude those false alarms from the analysis whose false alarm intervals overlap with true alarm intervals. Hence, approximately 10% of false alarms are excluded from the analysis.
5. Evaluation of the aoRM and aoTRMLS Alarm Systems
In this section, we first analyse the database and introduce a strategy to handle advisory alarms; that is, alarms that were not alarm relevant but helpful. Then, the results of our analysis regarding SE and FARR of the new aoRM and aoTRMLS alarm systems are presented.
5.1. The Database and Reassessment of Advisory Alarms
Our reference database contains recorded onlinemonitoring data from 85 different cases between January 2006 and May 2008. A case stands for one disease episode of a patient, meaning that some patients correspond to several cases.
The overall monitoring time is 1245:52:28 hours, and the mean monitoring time for each case is 14:39:26 hours (minimum 0:55:01 hours, maximum 31:35:25 hours). The monitoring system generated a total number of 9290 alarms, of which 4825 (52%) were simple threshold alarms. A total number of 9290 alarms means a frequency of 78 alarms per hour. However, considering the 85 cases, thealarmper hour frequency was quite dissimilar. Furthermore, the alarm frequency is higher in the morning and the afternoon, which is probably related to nursing actions during the day.
Several vital parameters were monitored like the heart rate, pulse, blood oxygen saturation and temperature, respiratory rate, and systolic, mean, and diastolic artery blood pressure. However, only mean artery blood pressure, heart rate, and blood oxygen saturation were monitored in each case.
Since 90% of all threshold alarms are caused by systolic and mean artery blood pressure (ART.S and ART.M), heart rate (HR), and oxygen saturation (SpO_{2}), our analysis concentrates on these four vital signs. The numbers of true, false, and advisory threshold alarms regarding ART.S, ART.M, HR, and SpO_{2} are listed in Table 2.

Only 10% of all alarms were true whereas 43% were false and 47% advisory alarms. This low rate of relevant alarms matches other studies; see [10]. (All alarms means all alarms given by the four considered vital signs in the following.) The most alarms were produced by ART.S (42%). This vital sign exhibits a lower rate of false alarms but more advisory alarms than average. SpO_{2} produced 26% of all alarms. Its false alarm rate is highest, and it has the lowest proportion of advisory alarms. For HR, which produces 21% of all alarms, the opposite is true: most HR alarms were advisory. ART.M produced the least alarms (11%), of which 11% were true, 49% advisory, and 40% false alarms.
Advisory alarms are problematic, since our analysis is based on true/false decisions. One possibility is to simply exclude advisory alarms from the analysis. However, this approach is not really satisfying since nearly half the alarms were advisory, for heart rate even more than threequarters. Too many situations would be excluded, so the informative value of the analysis would be lowered substantially. Moreover, advisory alarms also involve an alarm sound and, therefore, need to be considered. We, therefore, decide to assess each individual advisory alarm as either true or false. The assessment is done regarding the alarm length, which is the time an alarm is active, that is, the time span for which the measurements lie outside the alarm limits. The idea is that short alarm limit violations do not exhibit clinical relevance and are interfering rather than helpful [24, 25]. According to the physicians involved in the study, advisory alarms of ART.S and ART.M that are shorter than 10 seconds can be regarded as irrelevant. An advisory HR or SpO_{2} alarm is regarded as irrelevant if it is shorter than 5 or 15 seconds, respectively. Using this strategy, the set of 2032 advisory alarms is split into 1201 true and 831 false alarms. That is, there are 1650 (38%) true and 2709 (62%) false alarms now. However, due to overlapping alarm intervals (see Section 4), we exclude 418 false alarms from the analysis so that we obtain 1650 (42%) true and 2291 (58%) false alarms (Table 3). The highest false alarm rate is still produced by SpO_{2} with 79% irrelevant alarms. ART.S exhibits a 55% and ART.M a 51% false alarm rate. HR has the lowest false alarm rate (39%).

Due to the reassessment of advisory alarms, we are able to perform the analysis regarding SE and FARR as described in Section 4. In the following, we explain how aoRM and aoTRMLS are applied to the data. Afterwards, we present the results.
5.2. Application of aoRM and aoTRMLS
For the offline application of aoRM and aoTRMLS to the recorded monitoring data, we use the open source software R, version 2.10.1. The R package robfilter [17] contains functions of aoRM and aoTRMLS, the function names are adore.filter (equates to aoRM) and madore.filter (equates to aoTRMLS).
Using R, one can handle a broad range of data storage formats, also including the text file format of the reference data. Each of the 85 cases refers to one data set, that is, one text file. Each data set text file is loaded in R in form of a data matrix with rows, , where each column contains the measurements of one vital parameter. The number of rows equals the monitoring time in seconds of case .
The adore.filter function is designed for univariate time series which correspond to vectors in R. Hence, the adore.filter function is applied column by column. In contrast, the madore.filter function is designed for multivariate time series and can, therefore, be applied to a whole data matrix. The signal estimation outputs of the adore.filter and madore.filter functions are stored in text files, so that they can be analysed regarding SE and FARR as explained in Section 4. This analysis is also done by means of R 2.10.1.
For an application of the filters in clinical practice, the treatment of missing values is an issue we are concerned with: due to technical problems, onlinemonitoring data can contain missing values at single points as well as long stretches of missing values. The functions adore.filter and madore.filter have similar strategies to deal with missing values. Their algorithms give an output only if enough observations for a reasonable signal estimation are present. Otherwise, the signal estimation output is a missing value. For more details, see the R help or [15, 16].
Since the monitored vital signs hold a block dependence structure (cf. Section 2.2), we apply aoTRMLS block wise to one block consisting of ART.S, ART.M, and diastolic artery blood pressure, and to one block consisting of HR and pulse. SpO_{2} has an exceptional position since it is a “block” on its own. In this case, the madore.filter algorithm applies the univariate adore.filter meaning that the madore.filter and adore.filter signal estimations are equal for SpO_{2}. Our analysis still concentrates solely on ART.M, ART.S, HR, and SpO_{2} although aoTRMLS is also applied to diastolic artery blood pressure and pulse rate measurements—these measurements provide additional information for the aoTRMLS filter.
As explained in Section 2, the minimum window width is crucial, since it determines (a) the tracing delay of aoRM/aoTRMLS signal estimations when reacting to sudden data changes and (b) the number of outliers the filtering procedures can resist. Therefore, we apply aoRM and aoTRMLS to the whole data set using . The width of the time window, which is used to assess the RM regression fit and as comparison sample for the restricttorange rule, is always . The maximum window width is always .
5.3. Results
We determine SE and FARR of the aoRM and aoTRMLS alarm systems separately for each of the four vital signs ART.M, ART.S, HR, and SpO_{2}. In Figure 5, we plot FARR against SE for each vital sign and each minimum window width . SE and FARR of aoRM are indicated by dots, SE and FARR of aoTRMLS by crosses. The number below (above) a dot (cross) indicates the used minimum window width . Note that for HR the plot range is whereas for the other vital signs it is . As can be seen in Figure 5, the greater the , the greater the FARR and the smaller the SE. According to the physicians involved in the study, we demand at least 95% SE, marked by the grey lines in Figure 5. This helps with finding the which induces the largest FARR under the restriction of at least 95% SE.
(a)
(b)
(c)
(d)
For ART.M, aoTRMLS with offers and . The aoRM filter with yields comparable results with and FARR . Regarding ART.S and demanding at least 95% SE, we obtain somewhat worse results with for aoRM () and for aoTRMLS (). In order to obtain at least 95% SE for HR, a smaller minimum window width must be chosen. Using aoRM with , approximately onethird of all HR false alarms can be suppressed with . Using the same minimum window width for aoTRMLS, we obtain but merely FARR ; using leads to but with . Since SpO_{2} builds a block on its own, the multivariate madore.filter algorithm applies the univariate adore.filter to this block. Thus, the aoRM and aoTRMLS signal estimations are equal for SpO_{2}. Using we obtain and .
6. Summary and Outlook
The aoRM and aoTRMLS signal filters can be used to extract signals online from nonstationary, noisy, and outlier contaminated onlinemonitoring time series from intensive care. By comparing signal estimates instead of raw measurements to upper and lower alarm limits, the number of threshold alarms can be reduced.
The evaluation of aoRM and aoTRMLS for an application in clinical practice is done based on recorded data from an intensive care unit. In this data record, each alarm has been annotated retrospectively as either true or false. In order to evaluate the aoRM and aoTRMLS alarm systems, we estimate their sensitivity (SE) but not their specificity. The reason is that for a specificity estimation one has to decide whether or not to regard limit violations of the new system within “alarms off” periods, meaning that either the new or old alarm system is favoured. Hence, we refrain from estimating specificity and determine the false alarm reduction rate (FARR) of the new system instead.
FARR and SE of the new system are determined by building false and true alarm intervals which consist of a certain number of time points after a false and true alarm, respectively. If the aoRM/aoTRMLS signal estimations violate the alarm limits within a true alarm interval, the according true alarm is regarded as detected. Analogically, a false alarm of the old system is regarded as suppressed by the new system, if its signal estimations do not violate the alarm limits within the false alarm interval.
An analysis of the database shows that more than 90% of all threshold alarms were caused by ART.S, ART.M, HR, and SpO_{2}. Hence, we concentrate our analysis on these four vital signs. Furthermore, almost half of all alarms were assessed as advisory; these alarms do not require immediate action but are helpful. Since our evaluation strategy is based on true/false decisions, we reassess advisory alarms as false or true.
The application of aoRM and aoTRMLS is performed retrospectively using the R functions adore.filter and madore.filter from the R package robfilter. Since the choice of the minimum window width has a great impact on the aoRM/aoTRMLS signal extraction, we apply the filters using . The function madore.filter (aoTRMLS) has been applied block wise with one block consisting of arterial blood pressures and one of heart rate and pulse; oxygen saturation builds an own block, meaning that the multivariate madore.filter algorithm applies the univariate adore.filter to this block.
We found that both filters are able to suppress around a quarter to a third of all false alarms while providing at least 95% sensitivity. Using a larger minimum window width, a greater FARR can be obtained at the cost of a lower SE. For instance, merely demanding at least 90% SE, an aoRM or aoTRMLS alarm system can reduce 31 to 38% of all threshold alarms.
The need for an improvement of the actual situation on intensive care units is obvious. Preprocessing the raw monitoring data by the aoRM or aoTRMLS filter is a good possibility to achieve this aim. The proposed filters could be implemented into the monitoring systems. Then, the practical performance of a filtering based alarm system can be compared to that of the “old” system in a test phase. However, one has to keep in mind the lower sensitivity of filtering based alarm systems. Therefore, the physician must have access to both the filtered and the raw measurements at the bedside.
The potential of the aoRM and aoTRMLS signal filters is not restricted to false alarm suppression. For instance, they can also be used to ease patient monitoring since smooth signal extractions are easier to interpret than noisy and outlier contaminated measurements. The aoRM and aoTRMLS signal extraction can also be beneficial in other fields, for example, for highfrequency measurements from industry or finance.
Furthermore, the aoRM and aoTRMLS can be enhanced or used as a basis for other filters. For instance, there might be improved principles for the window width adaption, possibly based on methods for the detection of structural breaks; for the aoTRMLS, an automatic and timedependent choice of the correlation blocks (see Section 2.2) would be an improvement, especially for multivariate time series with an unknown and possibly changing dependence structure.
Acknowledgments
The authors sincerely thank the referees for the valuable and helpful comments and gratefully acknowledge the support of the Deutsche Forschungsgemeinschaft (SFB 475, Reduction of complexity in multivariate data structures, and SFB 823, Statistical modelling of nonlinear dynamic processes).
References
 M. C. Chambrin, D. CalveloAros, A. Jaborska, C. Chopin, P. Ravaux, and B. Boniface, “Multicentric study at monitoring alarms in the adult intensive care unit (ICU): a descriptive analysis,” Intensive Care Medicine, vol. 25, no. 12, pp. 1360–1366, 1999. View at: Publisher Site  Google Scholar
 C. L. Tsien and J. C. Fackler, “Poor prognosis for existing monitors in the intensive care unit,” Critical Care Medicine, vol. 25, no. 4, pp. 614–619, 1997. View at: Google Scholar
 M. C. Chambrin, “Alarms in the intensive care unit: how can the number of false alarms be reduced?” Critical Care, vol. 5, no. 4, pp. 184–188, 2001. View at: Publisher Site  Google Scholar
 S. T. Lawless, “Crying wolf: false alarms in a pediatric intensive care unit,” Critical Care Medicine, vol. 22, no. 6, pp. 981–985, 1994. View at: Google Scholar
 Z. Reslan, Clinical alarm management and noise reduction in hospitals, M.S. thesis, University of Connecticut, 2007.
 A. Mäkivirta, E. Koski, A. Kari, and T. Sukuvaara, “The median filter as a preprocessor for a patient monitor limit alarm system in intensive care,” Computer Methods and Programs in Biomedicine, vol. 34, no. 23, pp. 139–144, 1991. View at: Google Scholar
 R. R. Kennedy, “A modified Trigg's tracking variable as an 'advisory' alarm during anaesthesia,” International Journal of Clinical Monitoring and Computing, vol. 12, no. 4, pp. 197–204, 1995. View at: Publisher Site  Google Scholar
 S. Charbonnier, G. Becq, and L. Biot, “Online segmentation algorithm for continuously monitored data in intensive care units,” IEEE Transactions on Biomedical Engineering, vol. 51, no. 3, pp. 484–492, 2004. View at: Publisher Site  Google Scholar
 S. Charbonnier and S. Gentil, “A trendbased alarm system to improve patient monitoring in intensive care units,” Control Engineering Practice, vol. 15, no. 9, pp. 1039–1050, 2007. View at: Publisher Site  Google Scholar
 M. Imhoff and S. Kuhls, “Alarm algorithms in critical monitoring,” Anesthesia and Analgesia, vol. 102, no. 5, pp. 1525–1537, 2006. View at: Publisher Site  Google Scholar
 P. L. Davies, R. Fried, and U. Gather, “Robust signal extraction for online monitoring data,” Journal of Statistical Planning and Inference, vol. 122, no. 12, pp. 65–78, 2004. View at: Publisher Site  Google Scholar
 U. Gather, R. Fried, and V. Lanius, “Robust detailpreserving signal extraction,” in Handbook of Time Series Analysis, M. Winterhalder, B. Schelter, and J. Timmer, Eds., pp. 6–32, Wiley, Berlin, Germany, 2005. View at: Google Scholar
 U. Gather, K. Schettlinger, and R. Fried, “Online signal extraction by robust linear regression,” Computational Statistics, vol. 21, no. 1, pp. 33–51, 2006. View at: Publisher Site  Google Scholar
 A. F. Siegel, “Robust regression using repeated medians,” Biometrika, vol. 69, no. 1, pp. 242–244, 1982. View at: Publisher Site  Google Scholar
 K. Schettlinger, R. Fried, and U. Gather, “Realtime signal processing by adaptive repeated median filters,” International Journal of Adaptive Control and Signal Processing, vol. 24, no. 5, pp. 346–362, 2010. View at: Publisher Site  Google Scholar
 M. Borowski, K. Schettlinger, and U. Gather, “Multivariate real time signal processing by a robust adaptive regression filter,” Communications in Statistics—Simulation and Computation, vol. 38, no. 2, pp. 426–440, 2009. View at: Google Scholar
 R. Fried, K. Schettlinger, and M. Borowski, robfilter: Robust Time Series Filters, R package version 2.6.1, 2010.
 P. Rousseeuw and A. Leroy, Robust Regression and Outlier Detection, Wiley, New York, NY, USA, 1987.
 V. Lanius and U. Gather, “Robust online signal extraction from multivariate time series,” Computational Statistics and Data Analysis, vol. 54, no. 4, pp. 966–975, 2010. View at: Publisher Site  Google Scholar
 R. A. Maronna and R. H. Zamar, “Robust estimates of location and dispersion for highdimensional datasets,” Technometrics, vol. 44, no. 4, pp. 307–317, 2002. View at: Google Scholar
 S. Siebig, S. Kuhls, M. Imhoff et al., “Collection of annotated data in a clinical validation study for alarm algorithms in intensive carea methodologic framework,” Journal of Critical Care, vol. 25, no. 1, pp. 128–135, 2010. View at: Publisher Site  Google Scholar
 S. Kuhls, S. Siebig, F. Stóßel, and M. Imhoff, “Entwicklung einer Eingabemaske für die Erfassung klinischer Annotationen,” Tech. Rep. 15/06, SFB 475, Technische Universität Dortmund, 2006. View at: Google Scholar
 S. Kuhls, Planung, Dürchfuhrung und Analyse einer klinischen Studie zur Bewertung und zum Vergleich von Alarmsystemen in der Intensivmedizin, Ph.D. dissertation, Fakultät Statistik, Technische Universität Dortmund, 2008.
 M. Görges, B. A. Markewitz, and D. R. Westenskow, “Improving alarm performance in the medical intensive care unit using delays and clinical context,” Anesthesia & Analgesia, vol. 108, no. 5, pp. 1546–1552, 2009. View at: Publisher Site  Google Scholar
 M. Imhoff and R. Fried, “The crying wolf: still crying?” Anesthesia and Analgesia, vol. 108, no. 5, pp. 1382–1383, 2009. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2011 M. Borowski 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.