Compressed Anomaly Detection with Multiple
Mixed Observations

arXiv:1801.10264v1 [] 31 Jan 2018

Natalie Durgin, Rachel Grotheer, Chenxi Huang, Shuang Li, Anna Ma,
Deanna Needell and Jing Qin

Abstract We consider a collection of independent random variables that are identically distributed, except for a small subset which follows a different, anomalous
distribution. We study the problem of detecting which random variables in the collection are governed by the anomalous distribution. Recent work proposes to solve
this problem by conducting hypothesis tests based on mixed observations (e.g. linear
combinations) of the random variables. Recognizing the connection between taking
mixed observations and compressed sensing, we view the problem as recovering
the “support” (index set) of the anomalous random variables from multiple meaNatalie Durgin
Spiceworks
Austin, TX, 78746.
e-mail: njdurgin@gmail.com
Rachel Grotheer
Goucher College
Baltimore, MD 21204.
Chenxi Huang
Yale University
New Haven, CT 06511.
e-mail: chenxi.huang@yale.edu
Shuang Li
Colorado School of Mines
Golden, CO 80401.
Anna Ma
Claremont Graduate University
Claremont, CA 91711.
Deanna Needell
University of California, Los Angeles
Los Angeles, CA 90095.
Jing Qin
Montana State University
Bozeman, MT 59717.

1

2

Durgin, Grotheer, Huang, Li, Ma, Needell, Qin

surement vectors (MMVs). Many algorithms have been developed for recovering
jointly sparse signals and their support from MMVs. We establish the theoretical
and empirical effectiveness of these algorithms at detecting anomalies. We also extend the LASSO algorithm to an MMV version for our purpose. Further, we perform
experiments on synthetic data, consisting of samples from the random variables, to
explore the trade-off between the number of mixed observations per sample and the
number of samples required to detect anomalies.

1 Introduction
The problem of anomaly detection has been the focus of interest in many fields
of science and engineering, including network tomography, cognitive radio, and
radar [36, 40, 27, 4]. In this paper, we study the problem of identifying a small
number of anomalously distributed random variables within a much larger collection of independent and otherwise identically distributed random variables. We call
the random variables following the anomalous distribution anomalous random variables. A conventional approach to detecting these anomalous random variables is
to sample from each random variable individually and then apply hypothesis testing
techniques [29, 30, 31, 31].
A recent paper [13] proposes to perform hypothesis testing on mixed observations (e.g. linear combinations) of random variables instead of on samples from individual random variables. They call this technique compressed hypothesis testing.
Such an approach is motivated by the recent development of compressed sensing
[7, 15, 23, 21], a signal processing paradigm that shows a small number of random
linear measurements of a signal is sufficient for accurate reconstruction. Now a large
body of work in this area shows that optimization-based [17, 8, 16, 10, 37] and iterative [38, 33, 6] methods can reconstruct the signal accurately and efficiently when
the samples are taken via a sensing matrix satisfying certain incoherence properties
[8, 9]. Compressed sensing is also studied in a Bayesian framework, where signals
are assumed to obey some prior distribution [25, 41, 3].
The results presented in [13] show that the “mixed” measurement approach
achieves better detection accuracy from fewer samples when compared to the conventional “un-mixed” approach. However, compressed hypothesis testing requires
that the distributions of the random variables are known a priori, which may not be
available in practice. Further, as the authors pointed out, their proposed approach
requires conducting a large number of hypothesis tests, especially when the number
of random variables in the collection is large, rendering such an approach computationally prohibitive. Two efficient algorithms are proposed as alternatives in [13],
but no analytical study of their performance is provided.
We propose new methods for detecting anomalous random variables that require
minimal knowledge of the distributions, are computationally efficient, and whose
performance is easy to characterize. We begin by generalizing the compressed hypothesis testing method and posing our problem as a multiple measurement vector

Compressed Anomaly Detection with Multiple Mixed Observations

3

(MMV) problem [24, 2, 19, 18, 11, 14, 32, 5]. In the MMV compressed sensing setting, a collection of signals are recovered simultaneously, under the assumption that
they have some commonalities, such as sharing the same support. A related vein of
work involves signals that are smoothly varying, where the support may not be consistent but changes slowly over time [1, 22, 35]. While the compressed hypothesis
testing in [13] is certainly motivated by compressed sensing techniques, the authors
do not formally frame the anomaly detection problem in the compressed sensing setting. Also, they do not focus on compressed sensing algorithms that might eliminate
the need for prior knowledge of the distributions, and might lead to more efficient
detection for large collections of random variables.
In the following, we view the collection of random variables as a random vector
and aim to identify the indices of the anomalous random variables within the random
vector. We also draw an analogy between the collection of independent samples
from the random vector and an ensemble of signals where in practice these signals
often become available over time. More specifically, we consider a random vector,
X = (X1 , . . . , XN ), where the Xn ’s are independent random variables. We assume
that each Xn follows one of two distributions, D1 , D2 . We call D1 the prevalent
distribution, and D2 the anomalous distribution. We let N = {n ∈ N : 1 ≤ n ≤ N}
denote the index set of the random variables, Xn , and let K denote the index set of the
K random variables that follow the anomalous distribution. Let x(,t) ∈ RN denote
the independent realization of the random vector at time t. At each time-step t, we
obtain M mixed observations by applying the sensing matrix φt ∈ RM×N ,
yt = φt x(,t) , 1 ≤ t ≤ T,
with yt ∈ RM . Thus the goal of the anomaly detection problem in this setting is to
recover the index set K from the MMVs yt , t = 1, · · · , T .
The signals x(,t) in our formulation are not necessarily sparse and may have different supports since they are samples from a random vector and are changing over
time. Nevertheless, there is still a close connection between our formulation and that
for recovering the common sparse support of a collection of signals from MMVs.
The index set of the anomalous random variables, which corresponds to the index set
of the anomalies (realizations of anomalous random variables) in the signals x(,t) ,
is shared by all signals. This index set can thus be viewed as the common “support”
of the anomalies in the signals, which motivates us to consider the applicability of
many MMV algorithms designed for signal reconstruction. Further, the analytical
studies of many of these algorithms are readily available. We therefore investigate
which of these MMV algorithms can be applied or adapted to the anomaly detection
problem under consideration and analyze their performance in detection accuracy in
theory and through numerical experiments. We focus on algorithms presented in [2].

4

Durgin, Grotheer, Huang, Li, Ma, Needell, Qin

1.1 Contributions
In this paper, by extending the definitions of two so-called joint sparsity models
(JSMs) from [2], we introduce two new signal models, JSM-2R and JSM-3R, for
the problem of anomaly detection. For JSM-2R and JSM-3R signals, we adapt several MMV signal reconstruction algorithms to anomaly detection. Additionally, we
develop a new algorithm for the JSM-2R model that extends the Least Absolute
Shrinkage and Selection Operator (LASSO) algorithm [12] to the MMV framework. We show theoretically and numerically that these algorithms accurately detect
the anomalous random variables. We also provide numerical results which demonstrate the trade-off between the number of time-steps, and the number of mixed
observations per time-step needed to detect the anomalies.

1.2 Organization
In Section 2, we introduce the models JSM-2R, JSM-3R and the four algorithms we
have repurposed from MMV signal recovery into MMV anomaly detection, as well
as our new LASSO algorithm. We also provide theoretical guarantees in this section. In Section 3, we explore the performance of these algorithms by conducting
numerical experiments for some strategic choices of the parameters involved. Finally, we conclude in Section 4. To help keep track of notation, we provide a handy
reference table in Section 4. We adopt the convention that random variables will be
upper case and their realizations will be lower case. All matrix entries will have two,
subscripted indices. The first index will indicate the row position, the second will
indicate the column position.

2 Method
In this section, we introduce two new signal models for the anomaly detection problem and describe five algorithms for detecting anomalous random variables under
these signal models. We also provide theoretical guarantees for the algorithms.
Recall that we consider the problem of detecting K anomalous random variables
from a collection of N random variables where K  N. The anomalous random
variables have a different probability distribution from that of the remaining N − K
random variables. We seek to identify the K anomalous random variables, from T
independent realizations of the N random variables. To emphasize our framing of
this random variable problem as a compressed sensing problem, we refer to the independent realizations as signals. These T signals have an important commonality:
they share the same indices of anomalous entries (realizations of anomalous random variables). Commonality among signals has already been explored in the field
of distributed compressed sensing for recovering signals that have specific corre-

Compressed Anomaly Detection with Multiple Mixed Observations

5

lation among them. Three joint sparsity models (JSMs) were introduced in [2] to
characterize different correlation structures. To utilize the commonality of the signals for anomaly detection, we propose two new signal models that are motivated
by two of the JSMs defined in [2], namely, JSM-2 and JSM-3. Since the signals
under consideration are realizations of random variables, we term the new models
JSM-2R and JSM-3R, respectively, where the appended “R” indicates the “random
variable” version of the existing JSMs.
Before we define the new models, We first briefly describe JSM-2 and JSM-3.
The JSM-2 signals are jointly sparse signals that share the same support (the indices
of non-zero entries). The JSM-3 signals consist of two components: a non-sparse
“common component” shared by all signals and a sparse “innovation component”
that is different for each signal. But the innovation components of the JSM-3 signals share the same support. We next extend these definitions to the signals in the
anomaly detection setting. The new JSM-2R and JSM-3R models are defined as
follows:
Definition 1 (JSM-2R and JSM-3R). Let the random variable Xn ∼ D1 if n ∈
/K
and Xn ∼ D2 if n ∈ K where K is the set of the anomalous indices. For a signal
ensemble x ∈ RN×T where each of its entries x(n,t) denotes the realization of Xn at
time t,
1. x is a JSM-2R signal ensemble when: |x(n,t) | is small if n ∈
/ K and |x(n,t) | is large
if n ∈ K;
I
I
2. x is a JSM-3R signal ensemble when: x(n,t) = xCn + x(n,t)
such that |x(n,t)
| is small
I
C
if n ∈
/ K and |x(n,t) | is large if n ∈ K. xn is a common component shared by all t,
I
and x(n,t)
is an innovation component that is different for different t.
The JSM-2R signal model assumes a small amplitude for variables generated
from the prevalent distribution and a large amplitude for variables generated from
the anomalous distribution. Such a model characterizes a scenario where anomalies
exhibit large spikes. This model relates to a sparse signal model where the support
of the sparse signal corresponds to the set of indices of the anomalous random variables. In fact, when D1 = N (0, σ 2 ) and D2 = N (µ, σ 2 ) with µ 6= 0, the JSM-2R
signal is a sparse signal with additive Gaussian noise.
Different from the JSM-2R signals, the JSM-3R signal model introduced above
does not have constraints on the amplitude of the signal entries x(n,t) . Rather the
signals at different time-steps are assumed to share an unknown common compoI
for signals at different
nent xCn while having a different innovation component x(n,t)
C
time-steps. Of note, the common component xn from the prevalent distribution may
or may not be the same as that from the anomalous distribution. Further, the innovaI
tion component x(n,t)
is assumed to follow the JSM-2R signal model. Such a model
characterizes a scenario where there exists a background signal that does not change
over time and the anomalies exhibit large spikes on top of the background signal.
Because of the common component, the JSM-3R signals no longer correspond to a
sparse signal model. Fig. 1 provides a visual illustration of the model nuances.

Durgin, Grotheer, Huang, Li, Ma, Needell, Qin

1 2 3 4 5 6 7 8 9 10

t

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20

JSM-2R

1 2 3 4 5 6 7 8 9 10

t

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20

JSM-3

n

JSM-2

n

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20

n

n

6

1 2 3 4 5 6 7 8 9 10

t

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20

JSM-3R

10.0

7.5

5.0

2.5

0.0

2.5
1 2 3 4 5 6 7 8 9 10

t

Fig. 1 Depiction of the existing joint sparsity models (JSM-2 and JSM-3) and the new models
developed for anomaly detection (JSM-2R and JSM-3R). The distributions used to generate this
example are the same as the ones used for the numerical experiments in Section 3, see Table 1. The
index set of the anomalies is K = {6, 10}.

2.1 Algorithms
We briefly describe the five algorithms we study in this paper, among which three
are for JSM-2R signals and two are for JSM-3R signals. Two of the algorithms for
JSM-2R signals were originally proposed for recovering JSM-2 signals, including
the one-step greedy algorithm (OSGA) and the multiple measurement vector simultaneous orthogonal matching pursuit (MMV-SOMP) algorithm. We further propose
a new MMV version of the LASSO algorithm for detecting anomalies for JSM-2R
signals and investigate its performance via numerical experiments. The two algorithms for JSM-3R were also proposed in [2] for recovering JSM-3 signals, including the Transpose Estimation of Common Component (TECC) algorithm and the
Alternating Common and Innovation Estimation (ACIE) algorithm.
For each of the presented algorithms, the goal is to identify the indices of
the anomalous random variables from the mixed measurements yt = φt x(,t) for
t = 1, 2, . . . , T. The number of anomalies K is assumed to be known a priori. We first describe three algorithms we applied to anomaly detection for JSM-2R
signals.

2.1.1 OSGA
The OSGA algorithm is a non-iterative greedy algorithm introduced in [2] to recover the support of JSM-2 signals based on inner products of the measurement
and columns of the sensing matrix (Algorithm 1). We show in Theorem 1 that with
some condition on the prevalent and anomalous distributions, the OSGA algorithm
is able to recover the anomaly indices under the JSM-2R model, using a small number of measurements per time-step. Although the OSGA algorithm is shown to work

Compressed Anomaly Detection with Multiple Mixed Observations

7

asymptotically, it may not perform well when only a small number of time-steps are
available. Empirical evidence has confirmed this conjecture when the OSGA algorithm is used to reconstruct JSM-2 signals [2]. Thus we further consider approaches
like matching pursuit [28, 34] for our problem. Next, we describe the MMV version
of orthogonal matching pursuit algorithm proposed in [2].
Algorithm 1 OSGA
1:
2:
3:
4:
5:
6:

Input: yt , φt , K.
b
Output: K.
for n = 1, 2, . . . , N do
T
Compute ξn = T1 ∑t=1
hyt , φt(,n) i2
end for
b = {n, for the K largest ξn }
return K

2.1.2 MMV-SOMP
The MMV-SOMP algorithm is an iterative greedy pursuit algorithm for recovery
of jointly sparse signals. SOMP was first proposed in [39] and was adapted to the
MMV framework in [2]. Since our focus is not on signal recovery but on detecting anomalous entries, we adapt this algorithm for our JSM-2R signal model. The
adapted algorithm is presented in Algorithm 2, which identifies the anomaly indices
one at a time. In each iteration, the column index of the sensing matrices that accounts for the largest residual across signals of all time-steps is selected. Then the
remaining columns of each sensing matrix (for each time-step) are orthogonalized.
The algorithm stops after K iterations where K is the number of anomalous random
variables. We show through numerical experiments in Section 3 that the adapted
MMV-SOMP algorithm performs better than the OSGA algorithm for a small number of time-steps.

2.1.3 MMV-LASSO
The LASSO algorithm aims to find a sparse solution to the regression problem by
constraining the L1 norm of the solution [12]. The LASSO algorithm was also considered in [13] as an efficient algorithm for anomaly detection from mixed observations. However, the authors of [13] considered the LASSO algorithm when using only one measurement at each time-step. In this paper, we further extend the
LASSO algorithm to a more general setting for MMV and term it the MMV-LASSO
algorithm. The MMV-LASSO algorithm is described in Algorithm 3. The measurements yt ∈ RM up to T time-steps are concatenated vertically to become a vector
y ∈ R(MT )×1 ; the sensing matrices φt ∈ RM×N are also concatenated vertically to become φ ∈ R(MT )×N . The concatenated measurements and sensing matrices are then

8

Durgin, Grotheer, Huang, Li, Ma, Needell, Qin

Algorithm 2 MMV-SOMP
1:
2:
3:
4:
5:

Input: yt , φt , K.
b
Output: K.
b = 0,
Initialize: K
/ residual rt0 = yt .
for k = 1, . . . , K do.
Select

|hrtk−1 , φt(,n) i|
.
kφt(,n) k2
t=1
T

nk = arg max ∑
n

6:
7:

b = [K,
b nk ].
Update K
Orthogonalize selected basis vector against previously selected vectors for all t, 1 ≤ t ≤ T :
γt0 = φt(,k) ,

if k = 1,

hφt(,nk ) , γtl i l
γt ,
γtk = φt(,nk ) −
kγtl k22
l=0
k−1

∑

8:

if k > 1.

Update the residual for all t, 1 ≤ t ≤ T ,
rtk = rtk−1 −

hrtk−1 , γtk i k
γt .
kγtk k22

9: end for
b
10: return K

fed to the regular LASSO algorithm, where the anomaly indices are found by taking
indices corresponding to the K largest amplitudes of the estimate. The LASSO problem, that is, Step 4 in Algorithm 3, can be tackled by various approaches [20, 26],
which is out of scope of this paper.
Algorithm 3 MMV-LASSO
Input: yt , φt , K.
b
Output: K.
T T
T
T T
Let y = [yT
1 , . . . , yT ] and φ = [φ1 , . . . , φT ]
Solve
1
x̂ = arg min ky − φ xk22 +λ kxk1
2
x
5: Let x̂n denote the n-th element of x̂
b = {n, for the K largest |x̂n |}
6: return K
1:
2:
3:
4:

We next describe two algorithms for anomaly detection for JSM-3R signals.

Compressed Anomaly Detection with Multiple Mixed Observations

9

2.1.4 TECC
The key difference between JSM-2R and JSM-3R signals is that JSM-3R signals
share a common component that is unknown. Thus the two algorithms for the JSM3R signals aim to first estimate the common component from the mixed measurement and subtract the contribution of this component from the measurement. The
TECC algorithm was proposed in [2] for recovering JSM-3 signals. We also adapt
the algorithm to focus only on detecting the anomalous indices of JSM-3R signals,
and the adapted algorithm can be found in Algorithm 4. The first step of the TECC
algorithm estimates the common component of the JSM-3R signals. Using this estimate, the contribution of the remaining innovation component to the measurement
can be estimated. Then algorithms for JSM-2R signals can be applied to identify the
anomaly indices. We show in Theorem 2 that the TECC algorithm is able to identify the anomalous variables under some conditions on the prevalent and anomalous
distributions. Similar to the OSGA algorithm, while Theorem 2 guarantees the success of the TECC algorithm in the asymptotic case as T goes to infinity, it may not
perform well for a small T . Next we describe an alternative algorithm also proposed
in [2] for cases with a small T .
Algorithm 4 TECC
Input: yt , φt , K.
b
Output: K.
T T
T
T T
Let y = [yT
1 , . . . , yT ] , and φ = [φ1 , . . . , φT ]
1
T
b
C
Calculate x = T M φ y.
5: Calculate ybt = yt − φt xbC .
b from ybt by Algorithm 1, 2 or 3
6: Estimate K
b
7: return K
1:
2:
3:
4:

2.1.5 ACIE
The ACIE algorithm is an extension of the TECC algorithm, also introduced in [2],
based on the observation that the initial estimate of the common component may not
be sufficiently accurate for subsequent steps. Instead of one-time estimation in the
TECC algorithm, the ACIE algorithm iteratively refines the estimates of the common component and the innovation components. The ACIE algorithm can also be
easily adapted for the JSM-3R signals for anomaly detection. In the ACIE algorithm
described in Algorithm 5, we first obtain an initial estimate of the anomaly index
b using the TECC algorithm. Then for each iteration, we build a basis Bt for RM
set K
where M is the number of measurements at each time-step: Bt = [φt,Kb , qt ], where
b and qt
φ b is the subset of the basis vectors in φt corresponding to the indices in K
t,K

has orthonormal columns that spans the orthogonal complement of φt,Kb . Then we

10

Durgin, Grotheer, Huang, Li, Ma, Needell, Qin

can project the measurements onto qt to obtain the part of the measurement caused
by signals not in K:
yet = qt T yt ,
φet = qt T φt .

(1)
(2)

Then yet and φet are used to refine the estimate of the common component. After
subtracting the contribution of this estimated common component, algorithms such
as OSGA and MMV-SOMP described above can be applied to detect the anomalies.
Algorithm 5 ACIE
1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:

Input: yt , φt , K, L (iteration counter).
b
Output: K.
T T
Let y = [yT
1 , . . . , yT ]
b from Algorithm 4
Obtain an initial estimate of K
for l = 1, 2, . . . , L do
Update yet and φet according to Equations (1) and (2) for all t, 1 ≤ t ≤ T
T e
eT
eT
e†
eT e −1 eT
eT
Update xeC = φe† ye, where ye = [e
yT
T ] , φ = [φ1 , . . . , φT ] and φ = (φ φ ) φ
1 ,···,y
end for
Calculate ybt = yt − φt xeC
b from ybt by Algorithm 1, 2 or 3
Estimate K
b
return K

2.2 Theoretical Guarantees
In this section we show theoretically that Algorithm 1 and Algorithm 4 (coupled
with Algorithm 1 in step 6) can detect anomalies for the JSM-2R and JSM-3R settings, respectively.
Recall that Algorithm 1 is designed for JSM-2R signals where variables generated from the prevalent distribution are much smaller in amplitude than those
from the anomalous distribution. The following theorem shows that for JSM-2R
signals, the OSGA algorithm is able to identify the indices of the anomalous variables asymptotically, with very few measurements at each time-step.
Theorem 1. [Adapted from [2] Theorem 8] Let the M × N sensing matrix, φt , contain entries that are i.i.d ∼ N (0, 1) at each time-step t. Suppose the random variables, Xn , are distributed with D1 = N (0, σ12 ) if n ∈
/ K and D2 = N (µ2 , σ22 ) if
2
2
2
n ∈ K. Assuming µ2 +σ2 > σ1 , then with M ≥ 1 measurements per time-step, OSGA
recovers K with probability approaching one as T → ∞.
Proof. Before diving into the proof of Theorem 1, we first observe that the signals
correspond to the JSM-2R signals: with a zero mean and a potentially small variance

Compressed Anomaly Detection with Multiple Mixed Observations

11

σ12 for the prevalent distribution D1 , the signal entry x(n,t) , n ∈
/ K (i.e. the realization
of Xn at the time-step t) is expected to have small amplitude. In contrast, with a
non-zero mean µ2 and a similar or possibly larger variance σ22 for the anomalous
distribution D2 , the amplitude of x(n,t) , n ∈ K is expected to be much larger.
We assume, for convenience and without loss of generality, that the anomalous random variables are indexed by, K = {1, 2, . . . , K}, and the prevalent random variables are indexed by N\K = {K + 1, . . . , N}. Consider that the test statistic
T
hyt , φt(,n) i2 is the sample mean of the random variable hY, Φ(,n) i2 , so
ξn = T1 ∑t=1
by the Law of Large Numbers,
lim ξn = E[hY, Φ(,n) i2 ].

T →∞

We select an arbitrary index n from each of the anomalous random variable index
set and the prevalent random variable index set, and compute E[hY, Φ(,n) i2 ] in each
case. As the final step, we compare the expected values of the two ξn and establish
that they are distinguishable under very general conditions. Without loss of generality, we select n = K + 1 for the “prevalent case” and n = 1 for the “anomalous
case”. Note that [2] refers to these cases respectively as the “bad statistics” and the
“good statistics” in their setting. For them, “bad” reflects an incorrect estimate of
the sparse support and “good” reflects a correct estimate of the sparse support.
Prevalent Case: Substituting ΦX for Y in hY, Φ(,K+1) i and rearranging we obtain
hY, Φ(,K+1) i = ∑Nn=1 Xn hΦ(,n) , Φ(,K+1) i. We can then write,

! 
2

N

2

E[hY, Φ(,K+1) i ] = E 

∑ Xn hΦ(,n) , Φ(,K+1) i



n=1

"

#

N
2

2

∑ (Xn ) hΦ(,n) , Φ(,K+1) i

=E

n=1




N


+ E∑

N

∑ Xn Xl hΦ(,l) , Φ(,K+1) ihΦ(,n) , Φ(,K+1) i


n=1 l=1
l6=n

N

=

∑ E[(Xn )2 ]E[hΦ(,n) , Φ(,K+1) i2 ]

n=1

N

+

N

∑ ∑ E[Xn ]E[Xl ]E[hΦ(,l) , Φ(,K+1) ihΦ(,n) , Φ(,K+1) i].

n=1 l=1
l6=n

The last step follows from the independence of Φ and X and the independence of
the Xn ’s from each other. We claim that the cross-terms above sum to zero. To see
this, we set Φ(,l) = a, Φ(,K+1) = b and Φ(,n) = c, where the entries of the vectors
a, b, c are all i.i.d. N (0, 1). We note that if l, K + 1, and n are mutually distinct, then
a, b, c are mutually independent. In this case we have,

12

Durgin, Grotheer, Huang, Li, Ma, Needell, Qin

E[ha, bihc, bi] = E[aT bcT b]
= E[aT ]E[bcT b]
= 0.
Since the cross-terms assume l 6= n, we consider the cases when either n = K + 1
or l = K + 1. In the case where n = K + 1 we have,
E[ha, bihb, bi] = E[aT bbT b]
= E[aT ]E[bbT b]
= 0.
Similarly, in the case where l = K + 1 we have,
E[hb, bihc, bi] = E[bT bcT b]
= E[cT bbT b]
= E[cT ]E[bbT b]
= 0.
Thus, all cross-terms vanish so returning to our original goal we may claim,
N

E[hY, Φ(,K+1) i2 ] =

∑ E[(Xn )2 ]E[hΦ(,n) , Φ(,K+1) i2 ]

n=1
K

=

∑ E[(Xn )2 ]E[hΦ(,n) , Φ(,K+1) i2 ] + E[(XK+1 )2 ]E[kΦ(,K+1) k4 ]

n=1

N

+

∑

E[(Xn )2 ]E[hΦ(,n) , Φ(,K+1) i2 ].

n=K+2

Examining each expected value individually, we recall that for n ∈ {1, . . . , K} = K
the Xn were distributed with D2 and thus E[(Xn )2 ] = E[Xn ]2 + Var(Xn ) = µ22 +
σ22 . Recalling that the rest of the Xn are distributed with D1 which has µ = 0,
we have that E[(Xn )2 ] = σ12 in the subsequent cases. In [2] they establish that
E[kΦ(,K+1) k4 ] = M(M + 2) and E[hΦ(,n) , Φ(,K+1) i2 ] = M, and we may use these
results without further argument because we make the same assumptions about Φ.
Finally, substituting the expected values we have just calculated, we have that as T
grows large, the statistic ξn when n ∈
/ K converges to
2
2
2
E[hY, Φ(,K+1) i ] = K(µ2 + σ2 )M + σ12 M(M + 2) + (N − K − 1)σ12 M
= M[K(µ22 + σ22 ) + (M + 1 + N − K)σ12 ].

(3)

Compressed Anomaly Detection with Multiple Mixed Observations

13

Anomalous Case: With n = 1, we proceed as in the previous case,

! 
2

N

2

E[hY, Φ(,1) i ] = E 

∑ Xn hΦ(,n) , Φ(,1) i



n=1
N

=

∑ E[(Xn )2 ]E[hΦ(,n) , Φ(,1) i2 ],

n=1

K

= E[(X(1,t) )2 ]E[kΦ(,1) k4 ] +

∑ E[(Xn )2 ]E[hΦ(,n) , Φ(,1) i2 ]

n=2

N

+

∑

E[(Xn )2 ]E[hΦ(,n) , Φ(,1) i2 ]

n=K+1

= (µ22 + σ22 )M(M + 2) + (K − 2)(µ22 + σ22 )M + (N − K)σ12 M
= M[(M + 1 + K)(µ22 + σ22 ) + (N − K)σ12 ].
Combining the results of (3) and (4), we have
(
M[(M + 1 + K)(µ22 + σ22 ) + (N − K)σ12 ]
lim ξn =
T →∞
M[K(µ22 + σ22 ) + (M + 1 + N − K)σ12 ]

(4)

n∈K
n∈
/ K.

The difference in the two expectations is thus,
M(M + 1)(µ22 + σ22 − σ12 ).
For any M ≥ 1 and µ22 + σ22 > σ12 , the expected value of ξn in the “anomalous case”
is strictly larger than the expected value of ξn in the “prevalent case”. Therefore,
as T increases, OSGA can distinguish between the two expected values of ξn with
overwhelming probability.
t
u
The next theorem shows that asymptotically, Algorithm 4 is able to detect anomalous variables with very few measurements at each time-step, for JSM-3R signals.
Recall that JSM-3R signals have an unknown common component shared by signals at all time-steps, while each signal has a different innovation component that
follows the JSM-2R model. The following theorem and proof assume that Algorithm 1 is implemented for step 6 of Algorithm 4. Once results like Theorem 1 exist
for Algorithms 2 and 3, then any JSM-2R algorithm could be used in step 6, and
Theorem 2 would still hold.
Theorem 2 (Adapted from [2] Theorem 10). Let the M × N sensing matrix φt at
each time-step t contain entries that are i.i.d. ∼ N (0, 1). For random variables Xn
that are distributed with D1 = N (µ1 , σ12 ) if n ∈
/ K and D2 = N (µ2 , σ22 ) if n ∈ K, if
2
2
σ2 > σ1 and with M ≥ 1, TECC algorithm (with OSGA) recovers K with probability
approaching one as T → ∞.

14

Durgin, Grotheer, Huang, Li, Ma, Needell, Qin

Proof. We first note that the signals in Theorem 2 correspond to the JSM-3R signals:
I where xI are
for n ∈
/ K, the signal entries x(n,t) can be written as x(n,t) = µ1 + xn,t
n,t
2
i.i.d. ∼ N (0, σ1 ). With zero-mean and a potentially small variance, the amplitude
I ,n∈
of xn,t
/ K is expected to be small. For n ∈ K, the signal entries x(n,t) can be written
I where xI are i.i.d. ∼ N (0, σ 2 ). With a larger variance σ 2 , the
as x(n,t) = µ2 + xn,t
n,t
2
2
I
amplitude of xn,t , n ∈ K is expected to be much larger.
By the common component estimation from Algorithm 4, we have:
1 T
φ y
xbC =
TM
1 1 T T
=
∑ φt yt
M T t=1
1
=
M

1
T

T

∑ φtT φt x(,t)

!
.

t=1

Note that this is 1/M times the sample mean of the random variable ΦT ΦX. Letting
IN denote the N × N identity matrix, we note that since Φ has independent N (0, 1)
entries then E[ΦT Φ] = MIN . Since Φ is fully independent of X,
1
1
E[ΦT ΦX] = E[ΦT Φ]E[X] = IN E[X] = E[X].
M
M
Invoking the Law of Large Numbers, we have
lim xbC = E[X].

T →∞

Let Xb = X − xbC , then as T → ∞, Xbn is distributed as N (0, σ12 ) if n ∈
/ K and
b it follows from Theorem
N (0, σ22 ) if n ∈ K. Since Yb = Y −ΦxbC = Φ(X − xbC ) = ΦX,
1 that with M ≥ 1 and σ22 > σ12 , the TECC with OSGA algorithm recovers K with
probability approaching one as T → ∞.
t
u

3 Experiments
In this section, we evaluate numerically the performance of Algorithms 1, 2, 3,
4 and 5 for anomaly detection. More specifically, we examine the success rate of
determining the anomalous index set K from the signal matrix x ∈ RN×T , whose
columns are signals obtained at each time-step and share the same anomalous indices. The performance is assessed under various settings, by varying the number
of anomalies, the number of columns in x (i.e. the time-steps) and the number of
mixed measurement M at each time-step. Our focus is on the trade-off between the
number of measurements M and the number of time-steps T required to identify K
for varying numbers of anomalies.

Compressed Anomaly Detection with Multiple Mixed Observations

15

In all experiments, the measurement matrices φt ∈ RM×N comprise independent,
N (0, 1) entries and the measurement vectors yt ∈ RM are calculated by yt = φt x(,t)
for t = 1, . . . , T . To obtain an estimate of an algorithm’s recovery success rate with
high confidence, instead of using a fixed number of random trials across the different parameter combinations, we adaptively determine the necessary number of
trials with a Jeffreys interval, a Bayesian two-tailed binomial proportion confidence
interval. When the 95% confidence interval around the true success rate shrinks to a
width smaller than 0.1, we report the current proportion of successes as the recovery
accuracy for the algorithm. The signals (i.e. x(,t) ) are generated under two models
corresponding to the JSM-2R and JSM-3R signal definitions introduced in Section
2. Algorithms 1, 2 and 3 are applied to the JSM-2R signals while Algorithms 4 and
5 are applied to the JSM-3R signals.
The experiments are summarized in Table 1. The JSM-2R experiments assume a
mean zero for the prevalent distribution and a much larger mean for the anomalous
distribution while letting the variance be small. As shown in the previous section,
signals generated from these distributions satisfy the definitions of JSM-2R. For
JSM-3R experiments, we explore two settings: First, the prevalent and anomalous
distributions are assumed to have different means; second, the two distributions have
the same mean. Recall from the previous section, we show that the means of the
distributions are the common components for the JSM-3 signals generated from
these distributions. Note that the algorithms for JSM-3R signals have no knowledge
of the mean of the prevalent or anomalous distributions.
Table 1
Signal model
JSM-2R
JSM-3R
JSM-3R

D1

D2

Algorithms
OSGA, MMV-SOMP,
N (0, 1) N (7, 1)
MMV-LASSO
N (7, 1) N (0, 10)
TECC, ACIE
N (7, 1) N (7, 10)
TECC, ACIE

We chose the distributions in Table 1 for our numerical simulations to remain
consistent with [13]. We observe, in the JSM-2R experiments, that the distributions
N (0, 1) and N (7, 1) have their means separated by three standard deviations each,
with one additional standard deviation in between for good measure. This ensures
that the distributions are statistically distinct from each other. We have not explored
how the detection accuracy is affected as we vary the proportion of overlap in the
two distributions.

3.1 JSM-2R
We now present the results of recovering the anomalous index set for the JSM-2R
signals. The signal length is fixed at N = 100 and results of K = 1, 5, and 10 anoma-

16

Durgin, Grotheer, Huang, Li, Ma, Needell, Qin

lies are presented. For each K value, K random variables follow the distribution
N (7, 1) and the other N − K random variables follow another distribution N (0, 1).
The goal is to recover the index set K of these K random variables. Fig. 2 shows
the success rate of identifying K for the three K values using the OSGA Algorithm.
Each dot in the figure denotes the success rate for a specific M (number of measurements per time-step) and a specific T (number of time-steps) estimated from a
number of trials, and the value is indicated by the color (see the colorbar). Both M
and T take values from 1 to 100. Fig. 3 and 4 plot the success rate for MMV-SOMP
and MMV-LASSO algorithms respectively.
For all three algorithms, the success rate of anomaly identification increases as
the number of measurements M increases and/or as the number of time-steps T
increases. A 100% success of identification is obtained with a sufficiently large
number of measurements and time-steps. There are some important differences in
performance among the three algorithms.
Firstly, for the OSGA and the MMV-SOMP algorithms, with a sufficiently large
number of time-steps, the minimum number of measurements at each time-step required for anomaly detection increases with the number of anomalies present. The
MMV-LASSO performance seems less affected by varying the number of anomalies than the performance of the other two algorithms. Secondly, comparing Fig. 2
and 3 reveals that MMV-SOMP requires fewer time-steps than the OSGA algorithm to reach 100% success for a given number of measurements. Thirdly, the
MMV-LASSO algorithm requires significantly fewer measurements and time-steps
for 100% success compared with OSGA and MMV-SOMP. Finally, there is asymmetry between the effect of increasing the number of measurements versus that of
increasing the number of time-steps on the performance of OSGA and MMV-SOMP.
For these two algorithms, increasing the number of measurements is more effective
than increasing the number of time-steps for improving the performance. No obvious asymmetry of recovery performance is found for the MMV-LASSO algorithm.
The near symmetry phenomenon of the MMV-LASSO is expected since doubling
either M or T doubles the number of rows in the matrix φ in Algorithm 3, providing
similar amounts of information for the algorithm.

K=5

K=10
10

20

20

20

30

30

30

40

40

40

50

50

50

m

10

m

m

K=1
10

60

60

60

70

70

70

80

80

80

90

90

90

100

10 20 30 40 50 60 70 80 90 100

t

100

10 20 30 40 50 60 70 80 90 100

t

100

1.0
0.8
0.6
0.4
0.2

10 20 30 40 50 60 70 80 90 100

0.0

t

Fig. 2 The recovery phase transition for the OSGA algorithm with K = 1, 5, and 10 anomalous
random variables.

Compressed Anomaly Detection with Multiple Mixed Observations

K=5

K=10

10

10

20

20

20

30

30

30

40

40

40

50

50

50

m

10

m

m

K=1

60

60

60

70

70

70

80

80

80

90

90

90

100

10 20 30 40 50 60 70 80 90 100

100

t

17

10 20 30 40 50 60 70 80 90 100

100

t

1.0
0.8
0.6
0.4
0.2

10 20 30 40 50 60 70 80 90 100

0.0

t

Fig. 3 The recovery phase transition for the MMV-SOMP algorithm with K = 1, 5, and 10 anomalous random variables.

K=5

K=10
10

20

20

20

30

30

30

40

40

40

50

50

50

m

10

m

m

K=1
10

60

60

60

70

70

70

80

80

80

90

90

90

100

10 20 30 40 50 60 70 80 90 100

t

100

10 20 30 40 50 60 70 80 90 100

t

100

1.0
0.8
0.6
0.4
0.2

10 20 30 40 50 60 70 80 90 100

0.0

t

Fig. 4 The recovery phase transition for the MMV-LASSO algorithm with K = 1, 5, and 10 anomalous random variables.

3.2 JSM-3R
We next present the results of recovering the anomalous index set for the JSM3R signals. Similar to JSM-2R signals, the length of the signal is set to N = 100
and the number of anomalies takes values of K = 1, 5, and 10. Unlike the JSM2R signals, the N − K random variables now follow the distribution N (7, 1) while
the K anomalous random variables follow the distribution N (0, 10) or N (7, 10).
In order for a fair comparison between the algorithms, we implement the OSGA
algorithm for both step 6 of the TECC algorithm and step 10 of the ACIE algorithm.
The iteration L in the ACIE algorithm is set to L = 5. The performance of the TECC
and ACIE algorithms for varying numbers of measurements M and time-steps T
when the anomalous distribution follows N (0, 10) is presented in Figs. 5 and 6,
where both M and T range from 1 to 100. The performance for the setting where
the anomalous variables are distributed as N (7, 10) is similar to Figs. 5 and 6 and
is thus omitted.
With a sufficiently large number of measurements and time-steps, both algorithms are able to achieve 100% success in recovery of the anomalous index set. For
a fixed number of time-steps, the minimum number of measurements required for
identification increases as the number of anomalies increases for both algorithms.
There is improvement in performance of the ACIE algorithm over the TECC al-

18

Durgin, Grotheer, Huang, Li, Ma, Needell, Qin

gorithm. The ACIE algorithm requires fewer time-steps to reach 100% recovery
success, for a given number of measurements; similarly, it requires fewer measurements for 100% recovery success with a given number of time-steps.

K=5

K=10
10

20

20

20

30

30

30

40

40

40

50

50

50

m

10

m

m

K=1
10

60

60

60

70

70

70

80

80

80

90

90

90

100

10 20 30 40 50 60 70 80 90 100

100

t

10 20 30 40 50 60 70 80 90 100

100

t

1.0
0.8
0.6
0.4
0.2

10 20 30 40 50 60 70 80 90 100

0.0

t

Fig. 5 The recovery phase transition for the TECC algorithm with K = 1, 5, and 10 anomalous
random variables. Here the prevalent distribution is N (7, 1) and the anomalous distribution is
N (0, 10).

K=5

K=10
10

20

20

20

30

30

30

40

40

40

50

50

50

m

10

m

m

K=1
10

60

60

60

70

70

70

80

80

80

90

90

90

100

10 20 30 40 50 60 70 80 90 100

t

100

10 20 30 40 50 60 70 80 90 100

t

100

1.0
0.8
0.6
0.4
0.2

10 20 30 40 50 60 70 80 90 100

0.0

t

Fig. 6 The recovery phase transition for the ACIE algorithm with K = 1, 5, and 10 anomalous
random variables. Here the prevalent distribution is N (7, 1) and the anomalous distribution is
N (0, 10).

4 Conclusion
In this paper, we formally posed the problem of detecting anomalously distributed
random variables as an MMV problem, by drawing an analogy between samples
of the random variables and ensembles of signals. We further established two signal models characterizing possible correlation structures among signals that contain
anomalous entries. Based on the new signal models, we showed through theoretical
and numerical analysis that many of the MMV algorithms for sparse signal recovery can be adapted to the anomaly detection problem. For two of the algorithms, we

Compressed Anomaly Detection with Multiple Mixed Observations

19

provided theoretical guarantees of anomaly detection in the asymptotic case. Our
experimental results on synthetic data show good performance for signals conforming to either model, when a sufficiently large number of time-steps is available.
While these algorithms succeed in detecting anomalies, there is still room for
optimizing performance. Currently these algorithms require storing the sensing matrices at each time-step in memory. In future work, we would like to explore optimal
ways to design sensing matrices to reduce the memory burden. Having provided
asymptotic anomaly detection guarantees for two algorithms, we are further interested in providing such guarantees for all the algorithms presented. Additionally, we
are interested in characterizing the performance bounds for each algorithm’s finite
sample case.
Acknowledgements The initial research for this effort was conducted at the Research Collaboration Workshop for Women in Data Science and Mathematics, July 17-21 held at ICERM. Funding
for the workshop was provided by ICERM, AWM and DIMACS (NSF grant CCF-1144502). SL
was supported by NSF CAREER grant CCF−1149225. DN was partially supported by the Alfred
P. Sloan Foundation, NSF CAREER #1348721, and NSF BIGDATA #1740325. JQ was supported
by the faculty start-up fund of Montana State University.

20

Durgin, Grotheer, Huang, Li, Ma, Needell, Qin

Appendix
Here we provide a summary of notation for reference.
N Number of random variables
N Set of random variables indices, {n ∈ N : 1 ≤ n ≤ N}
K Number of anomalous random variables
K Set of anomalous random variable indices, K ⊂ N, |K|= K
M Number of measurements per time-step
m Measurement index, 1 ≤ m ≤ M
T Number of time-steps measured
t Time-step index, 1 ≤ t ≤ T
D1 Prevalent distribution
D2 Anomalous distribution
X Random vector comprising independent random variables X1 , . . . , XN
x N × T -dimensional matrix of independent realizations of X for all T time-steps
Φ M × N-dimensional sensing matrix, i.i.d. ∼ N (0, 1) entries
φt M × N-dimensional realization of Φ at time t
φ (M · T ) × N-dimensional vertical concatenation of the φt , [φ1T , . . . , φTT ]T
yt M-dimensional result of measuring the signal, φt · x(,t) , at time t
T T
y (M · T )-dimensional vertical concatenation of the yt , [yT
1 , . . . , yT ]
Y M-dimensional random vector defined by ΦX
JSM Joint Sparsity Model, introduced in [2]
JSM-2 Signals are nonzero only on a common set of indices
JSM-3 Signals consist of common non-sparse component and a sparse innovation
JSM-2R “Random variable” version of JSM-2
JSM-3R “Random variable” version of JSM-3
OSGA One-step greedy algorithm
MMV Multiple Measurement Vector
MMV-LASSO MMV Least Absolute Shrinkage and Selection Operator
MMV-SOMP MMV Simultaneous Orthogonal Matching Pursuit
TECC Transpose Estimation of Common Component
ACIE Alternating Common and Innovation Estimation
We adopt the convention that random variables will be upper case and their realizations will be lower case. All matrix entries will have two, subscripted indices. The
first index will indicate the row position, the second will indicate the column position. We will indicate row and column vectors by substituting  for the respective
index.

Compressed Anomaly Detection with Multiple Mixed Observations

21

References
[1] Angelosante D, Giannakis GB, Grossi E (2009) Compressed sensing of timevarying signals. In: Digital Signal Processing, 2009 16th International Conference on, IEEE, pp 1–8
[2] Baron D, Wakin MB, Duarte MF, Sarvotham S, Baraniuk RG (2006) Distributed compressed sensing. IEEE Trans Info Theory 52(12):5406–5425
[3] Baron D, Sarvotham S, Baraniuk RG (2010) Bayesian compressive sensing
via belief propagation. IEEE Trans Signal Process 58(1):269–280
[4] Basseville M, Nikiforov IV, et al (1993) Detection of abrupt changes: theory
and application, vol 104. Prentice Hall Englewood Cliffs
[5] Berg E, Friedlander MP (2009) Joint-sparse recovery from multiple measurements. Arxiv preprint arXiv:09042051
[6] Blumensath T, Davies ME (2009) Iterative hard thresholding for compressed
sensing. Appl Comput Harmon Anal 27(3):265–274
[7] Candès EJ (2006) Compressive sampling. In: Proc. Int. Congress of Mathematicians, Madrid, Spain, vol 3, pp 1433–1452
[8] Candès EJ, Tao T (2005) Decoding by linear programming. IEEE Trans Info
Theory 51:4203–4215
[9] Candès EJ, Tao T (2006) Near optimal signal recovery from random projections: Universal encoding strategies? IEEE Trans Info Theory 52(12):5406–
5425
[10] Candès EJ, Romberg J, Tao T (2006) Stable signal recovery from incomplete
and inaccurate measurements. Commun Pur Appl Math 59(8):1207–1223
[11] Chen J, Huo X (2006) Theoretical results on sparse representations of multiple
measurement vectors. IEEE Trans Signal Process 54(12):4634–4643
[12] Chen SS, Donoho DL, Saunders MA (2001) Atomic decomposition by basis
pursuit. SIAM Rev 43(1):129–159
[13] Cho M, Xu W, Lai L (2016) Compressed Hypothesis Testing: To Mix or Not
to Mix? arXiv preprint arXiv:160907528
[14] Cotter S, Rao B, Engan K, Kreutz-Delgado K (2005) Sparse solutions to linear inverse problems with multiple measurement vectors. IEEE Trans Signal
Process 53(7):2477–2488
[15] Donoho D (2006) Compressed sensing. IEEE Trans Info Theory 52(4):1289–
1306
[16] Donoho DL (2006) Compressed sensing. IEEE Trans Info Theory 52(4):1289–
1306
[17] Donoho DL, Huo X (2001) Uncertainty principles and ideal atomic decomposition. IEEE Trans Info Theory 47:2845–2862
[18] Duarte M, Wakin M, Baron D, Baraniuk R (2006) Universal distributed sensing via random projections. In: Proc. Inf. Process. Sensor Networks (IPSN)
[19] Duarte MF, Wakin MB, Baron D, Sarvotham S, Baraniuk RG (2013) Measurement bounds for sparse signal ensembles via graphical models. IEEE Trans
Info Theory 59(7):4280–4289

22

Durgin, Grotheer, Huang, Li, Ma, Needell, Qin

[20] Efron B, Hastie T, Johnstone I, Tibshirani R (2004) Least angle regression.
Annals of Statistics 32:407–499
[21] Eldar YC, Kutyniok G (2012) Compressed sensing: theory and applications.
Cambridge University Press
[22] Filos J, Karseras E, Dai W, Yan S (2013) Tracking dynamic sparse signals with
hierarchical kalman filters: a case study. In: Digital Signal Processing (DSP),
2013 18th International Conference on, IEEE, pp 1–6
[23] Foucart S, Rauhut H (2013) A mathematical introduction to compressive sensing. Birkhäuser Basel
[24] Haupt J, Nowak R (2006) Signal reconstruction from noisy random projections. IEEE Trans Info Theory 52(9):4036–4048
[25] Ji S, Xue Y, Carin L (2008) Bayesian compressive sensing. IEEE Trans Signal
Process 56(6):2346–2356
[26] Keerthi SS, Shevade S (2007) A fast tracking algorithm for generalized
lars/lasso. IEEE Transactions on Neural Networks 18(6):1826–1830
[27] Lai L, Fan Y, Poor HV (2008) Quickest detection in cognitive radio: A sequential change detection framework. In: Global Telecommunications Conference,
2008. IEEE GLOBECOM 2008. IEEE, IEEE, pp 1–5
[28] Mallat SG, Zhang Z (1993) Matching pursuits with time-frequency dictionaries. IEEE Transactions on Signal Processing 41(12):3397–3415
[29] Malloy M, Nowak R (2011) On the limits of sequential testing in high dimensions. In: 2011 Conference Record of the Forty Fifth Asilomar Conference on
Signals, Systems and Computers (ASILOMAR), pp 1245–1249
[30] Malloy M, Nowak R (2011) Sequential analysis in high-dimensional multiple testing and sparse recovery. In: 2011 IEEE International Symposium on
Information Theory Proceedings, pp 2661–2665
[31] Malloy ML, Tang G, Nowak RD (2012) Quickest search for a rare distribution. In: 2012 46th Annual Conference on Information Sciences and Systems
(CISS), pp 1–6
[32] Mishali M, Eldar YC (2009) Reduce and boost: recovering arbitrary sets of
jointly sparse vectors. IEEE Trans Signal Process 56(10):4692–4702
[33] Needell D, Vershynin R (2007) Signal recovery from incomplete and inaccurate measurements via Regularized Orthogonal Matching Pursuit. IEEE J Sel
Top Signa 4:310–316
[34] Pati YC, Rezaiifar R, Krishnaprasad PS (1993) Orthogonal matching pursuit:
recursive function approximation with applications to wavelet decomposition.
In: Proceedings of 27th Asilomar Conference on Signals, Systems and Computers, pp 40–44 vol.1
[35] Patterson S, Eldar YC, Keidar I (2014) Distributed compressed sensing for
static and time-varying networks. IEEE Trans Signal Process 62(19):4931–
4946
[36] Poor HV, Hadjiliadis O (2009) Quickest detection, vol 40. Cambridge University Press Cambridge
[37] Tropp JA (2006) Just relax: Convex programming methods for subset selection
and sparse approximation. IEEE Trans Info Theory 52(3):1030–1051

Compressed Anomaly Detection with Multiple Mixed Observations

23

[38] Tropp JA, Gilbert AC (2007) Signal recovery from random measurements via
Orthogonal Matching Pursuit. IEEE Trans Info Theory 53(12):4655–4666
[39] Tropp JA, Gilbert AC, Strauss M (2005) Simultaneous sparse approximation
via greedy pursuit. In: ICASSP
[40] Xia Y, Tse D (2006) Inference of link delay in communication networks. IEEE
J Sel Areas in Commun 24(12):2235–2248
[41] Yu G, Sapiro G (2011) Statistical compressed sensing of gaussian mixture
models. IEEE Trans Signal Process 59(12):5842–5858

