网络诊断的概念最早由Vardi在1996年提出，现今的研究主要分为：
所谓网络起始节点至终端节点(SD)流量强度估计主要是想要估计网络内各条SD路径的封包流量。其主要概念为假设我们能观测到网络内各节点互相传送封包的路径，以下简称连结和各条连结的封包流量，由各条连结所观测到的结果来估计各条SD路径的网络流量。
主要问题架构与符号定义如下：
假设网络模型内有
n
{\displaystyle n}
r
{\displaystyle r}
c
{\displaystyle c}
r
{\displaystyle r}
c
{\displaystyle c}
令
X
j
(
k
)
{\displaystyle X_{j}^{\left(k\right)}}
j
{\displaystyle j}
k
{\displaystyle k}
X
j
(
k
)
∼
P
o
i
s
s
o
n
(
λ
j
)
,
j
=
1
,
.
.
.
,
c
,
k
=
1
,
.
.
.
,
K
{\displaystyle X_{j}^{\left(k\right)}\sim Poisson(\lambda _{j}),\;j=1,...,c,\;k=1,...,K}
Y
(
k
)
=
A
X
(
k
)
,
k
=
1
,
.
.
.
,
K
{\displaystyle {\boldsymbol {Y^{\left(k\right)}=AX^{\left(k\right)}}},k=1,...,K}
其中，
Y
(
k
)
=
(
Y
1
(
k
)
,
.
.
.
,
Y
r
(
k
)
)
′
{\displaystyle {\boldsymbol {Y^{\left(k\right)}}}=\left(Y_{1}^{\left(k\right)},...,Y_{r}^{\left(k\right)}\right)'}
K
{\displaystyle K}
A
{\displaystyle {\boldsymbol {A}}}
X
(
k
)
=
(
X
1
(
k
)
,
.
.
.
,
X
r
(
k
)
)
′
{\displaystyle {\boldsymbol {X^{\left(k\right)}}}=\left(X_{1}^{\left(k\right)},...,X_{r}^{\left(k\right)}\right)'}
K
{\displaystyle K}
我们希望利用观测到的Y向量去估计X向量中的参数值
λ
=
(
λ
1
,
.
.
.
,
λ
c
)
′
{\displaystyle \mathrm {\lambda =\left(\lambda _{1},...,\lambda _{c}\right)'} }
最大概似法 参数可能解中让概似函数最大者
借由EM算法找最大概似估计量
动差法 贝氏方法 假设网络模型有2条连结、3条SD路径、1期的SD流量，即
r
=
2
,
c
=
3
,
K
=
1
{\displaystyle r=2,c=3,K=1\ }
X
i
∼
P
o
i
s
s
o
n
(
λ
i
)
,
i
=
1
,
2
,
3
{\displaystyle X_{i}\sim Poisson(\lambda _{i}),\;i=1,2,3}
X
1
+
X
2
=
1
,
X
1
+
X
3
=
2
{\displaystyle X_{1}+X_{2}=1,X_{1}+X_{3}=2}
Y
=
(
1
2
)
,
A
=
(
1
1
0
1
0
1
)
,
X
=
(
X
1
X
2
X
3
)
{\displaystyle Y={\binom {1}{2}},A={\begin{pmatrix}1&1&0\\1&0&1\end{pmatrix}},X={\begin{pmatrix}X_{1}\\X_{2}\\X_{3}\end{pmatrix}}}
因此模型可表示成
(
1
2
)
=
(
1
1
0
1
0
1
)
⋅
(
X
1
X
2
X
3
)
{\displaystyle {\binom {1}{2}}={\begin{pmatrix}1&1&0\\1&0&1\end{pmatrix}}\cdot {\begin{pmatrix}X_{1}\\X_{2}\\X_{3}\end{pmatrix}}}
以最大概似法求最佳解为例，
先将所有可能的参数解找出。因为封包流量必须为正整数，因此只有以下两组解：
X
=
(
1
0
1
)
′
{\displaystyle X={\begin{pmatrix}1&0&1\end{pmatrix}}'}
X
=
(
0
1
2
)
′
{\displaystyle X={\begin{pmatrix}0&1&2\end{pmatrix}}'}
将可能的参数解代入概似函数
L
(
λ
)
=
(
λ
1
λ
3
+
λ
2
λ
3
2
2
)
e
x
p
(
−
λ
1
−
λ
2
−
λ
3
)
{\displaystyle L(\lambda )=(\lambda _{1}\lambda _{3}+{\frac {\lambda _{2}\lambda _{3}^{2}}{2}})exp(-\lambda _{1}-\lambda _{2}-\lambda _{3})}
找出让概似函数最大的参数解，即为最佳参数解
F
i
n
d
λ
:
m
a
x
λ
L
(
λ
)
{\displaystyle Find\;\lambda :max_{\lambda }L(\lambda )}
因为
L
(
(
1
0
1
)
′
)
{\displaystyle L({\begin{pmatrix}1&0&1\end{pmatrix}}')}
L
(
(
0
1
2
)
′
)
{\displaystyle L({\begin{pmatrix}0&1&2\end{pmatrix}}')}
最佳解为
λ
=
(
1
0
1
)
′
{\displaystyle {\lambda }={\begin{pmatrix}1&0&1\end{pmatrix}}'}
先将概似函数整理成期望值的型态
λ
=
E
λ
[
X
|
Y
=
A
X
]
{\displaystyle \mathbf {\lambda } =E_{\lambda }[X|Y=AX]}
选定起始值
λ
(
0
)
{\displaystyle \lambda ^{\left(0\right)}}
λ
(
n
+
1
)
=
E
[
X
|
Y
,
λ
(
n
)
]
{\displaystyle \lambda ^{\left(n+1\right)}=E[X|Y,\lambda ^{\left(n\right)}]}
利用EM算法求解的缺点是当网络模型较大时，在计算上比较复杂；即使当期数够多时，EM算法仅能提高估计上的准确性并无法解决计算上的复杂。
假设当各条连结流量观测的期数够多时，根据中央极限定理
Y
¯
=
1
K
∑
k
=
1
K
Y
(
k
)
d
i
s
t
.
→
N
r
(
A
λ
,
K
(
−
1
)
A
Λ
A
′
)
,
Λ
=
d
i
a
g
(
λ
)
{\displaystyle {\bar {Y}}={\frac {1}{K}}\sum _{k=1}^{K}Y^{\left(k\right)}\;{dist.\to }\;N_{r}(A\lambda ,K^{\left(-1\right)}A\Lambda A'),\;\Lambda =diag(\lambda )}
利用动差法，令样本平均数等于母体平均数，样本共变异数等于母体共变异数，即
Y
¯
=
A
λ
{\displaystyle {\bar {Y}}=A\lambda }
S
=
K
(
−
1
)
A
Λ
A
′
{\displaystyle S=K^{\left(-1\right)}A\Lambda A'}
由上述式子即可估计出参数解
因此Vardi提出动差法可解决计算上的困难，也可以利用产生动差等式解决参数解不唯一的情况。
所谓网络连结层级参数推估问题主要是想要推论网络连结的特性，例如节点传输之间资讯遗失率或延迟分配等。其主要概念为假设已知网络的形式，包含节点、路径等，一般常见为树状，以及假设已知网络特性的模型，搜集端点所测量的结果来找出有最大几率产生观察结果的网络参数。
若考虑封包遗失率下，其主要问题架构与符号定义如下：
假设有一个树状网络定义为
T
=
{
V
,
E
}
{\displaystyle T=\left\{{V},{E}\right\}}
表示该网络有V个节点(包含起始节点0、终端节点R及中间节点I)、E条连结。令
X
=
(
X
k
)
k
∈
V
{\displaystyle X=(X_{k})_{k\in V}}
其中
X
k
{\displaystyle X_{k}}
k
{\displaystyle k}
X
k
=
0
{\displaystyle X_{k}=0}
k
{\displaystyle k}
X
k
=
1
{\displaystyle X_{k}=1}
k
{\displaystyle k}
此外，若
X
i
=
1
{\displaystyle X_{i}=1}
X
j
=
1
{\displaystyle X_{j}=1}
i
{\displaystyle i}
j
{\displaystyle j}
α
i
{\displaystyle \alpha \,_{i}}
举例来说，
上图为一个树状网络
T
=
{
V
=
{
0
,
1
,
2
,
3
,
4
,
5
,
6
,
7
}
,
E
=
{
(
0
,
1
)
,
(
1
,
2
)
,
(
1
,
3
)
,
(
2
,
4
)
,
(
2
,
5
)
,
(
3
,
6
)
,
(
3
,
7
)
}
}
{\displaystyle T=\left\{{V}=\left\{0,1,2,3,4,5,6,7\right\},{E}=\left\{(0,1),(1,2),(1,3),(2,4),(2,5),(3,6),(3,7)\right\}\right\}}
数字表示节点，起始节点为0，中间节点为1、 2、 3，而终端节点为4、5、6、7，
α
i
{\displaystyle \alpha \,_{i}}
i
{\displaystyle i}
令封包传送结果
X
(
R
)
=
(
X
k
)
k
∈
R
{\displaystyle X_{(R)}=(X_{k})_{k\in R}}
p
(
x
;
α
)
=
P
α
(
X
(
R
)
=
x
)
{\displaystyle p\,(x;\alpha )=P_{\alpha }(X_{(R)}=x)}
并假设发送了
n
{\displaystyle n}
n
(
x
)
{\displaystyle n(x)}
x
{\displaystyle x}
n
{\displaystyle n}
x
1
{\displaystyle x^{1}}
x
2
{\displaystyle x^{2}}
x
n
{\displaystyle x^{n}}
p
(
x
1
,
x
2
,
⋯
,
x
n
;
α
)
{\displaystyle p(x^{1},x^{2},\cdots ,x^{n};\alpha )}
=
∏
m
=
1
n
p
(
x
m
;
α
)
{\displaystyle =\prod _{m=1}^{n}p(x^{m};\alpha )}
=
∏
x
∈
Ω
p
(
x
;
α
)
n
(
x
)
{\displaystyle =\prod _{x\in \Omega }p(x;\alpha )^{n(x)}}
因此问题的目标即为估计
α
{\displaystyle \alpha }
从起始节点传送封包，并观察终端节点封包通过情况。传送封包主要有两种情况，一种为一次只传送到一个接收的端点，称为单一传送 ；另一种为封包传送到特定的一些接收端点，称为多重传送 。然而这两种传送方式较没有弹性，且无法使用不同的流量或不同时间下观察网域，因此Xi et al. (2006)及Lawrence et al. (2006)针对弹性传送(flexicast) 封包的情况作探讨。
此种观察封包传送情况来对网络做推论产生了统计反向问题 ，即利用观察结果来诊断连结中的分配或特征。有许多统计方法可解决此类推论问题，Castro et al. (2004)提到像是降低复杂性的阶层统计模型(Complexity-Reducing Hierarchical Statistical Models)、动差或最大概似法为主的估计、EM及马可夫链蒙地卡罗(Markov Chain Monte Carlo, MCMC)演算方法等已被使用；且认为而使用统计方法来解决此问题的领域仍具有发展性，而未来应有更多现存的统计方法可加以应用。
以下兹列举一种问题情况：“针对多重传送为主的网络来推论该网络的封包遗失率”来说明网络连结参数中的遗失率推估问题。估计封包遗失率为Cáceres et al.(1999)首先研究，在假设连结遗失为独立的伯努利分配下 ，利用最大概似法来估计多重传送的树状网络中连结遗失率；他们亦证明此估计量具备强烈一致性 ，并透过最大概似估计量之渐近常态性来推导出这些估计的比率会收敛到真正的比率。
以最大概似法求估计之连结遗失率方法如下：首先计算对数概似函数，
L
(
α
)
=
log
p
(
x
1
,
x
2
,
⋯
,
x
n
;
α
)
=
∑
[
n
(
x
)
l
o
g
p
(
x
;
α
)
]
{\displaystyle L(\alpha )=\log \,p(x^{1},x^{2},\cdots ,x^{n};\alpha )=\sum [{n(x)}log\,p(x;\alpha )]}
则
α
{\displaystyle \alpha }
α
^
=
a
r
g
m
a
x
α
∈
[
0
,
1
]
L
(
α
)
{\displaystyle {\hat {\alpha }}=argmax_{\alpha \in [0,1]}L(\alpha )}
另外，Cáceres et al.(1999)亦利用终端节点接收封包几率来估计
α
{\displaystyle \alpha }
R
(
k
)
{\displaystyle R(k)}
k
{\displaystyle k}
Ω
(
k
)
{\displaystyle \Omega (k)}
R
(
k
)
{\displaystyle R(k)}
γ
k
=
P
[
Ω
(
k
)
]
{\displaystyle \gamma _{k}=P[\Omega (k)]}
γ
k
{\displaystyle \gamma _{k}}
Σ
[
n
(
x
)
n
]
{\displaystyle \Sigma \left[{\frac {n(x)}{n}}\right]}
k
=
f
(
j
)
{\displaystyle k=f\,(j)}
k
{\displaystyle k}
j
{\displaystyle j}
f
n
(
j
)
=
f
(
f
n
−
1
(
j
)
)
{\displaystyle f\,^{n}\,(j)=f(f\,^{n-1}\,(j))}
n
{\displaystyle n}
l
(
k
)
{\displaystyle l(k)}
k
{\displaystyle k}
β
k
=
P
[
Ω
(
k
)
|
X
f
(
k
)
=
1
]
{\displaystyle \beta _{k}=P\,[\Omega (k)|X_{f(k)}=1]}
表示给定从第
k
{\displaystyle k}
γ
k
{\displaystyle \gamma _{k}}
α
{\displaystyle \alpha }
γ
k
=
β
k
∏
i
=
1
l
(
k
)
α
f
i
(
k
)
{\displaystyle \gamma _{k}=\beta _{k}\prod _{i=1}^{l(k)}\alpha _{f\,^{i}(k)}}
即将通过第k条连结所在从起始到终端节点的所有
α
k
{\displaystyle \alpha _{k}}
γ
k
{\displaystyle \gamma _{k}}
以两层的树状网络为例：
令通过此网络终端节点的可能情况集合为
Ω
=
{
00
,
01
,
10
,
11
}
{\displaystyle \Omega =\left\{00,01,10,11\right\}}
其中
00表示封包皆没有通过节点2跟3，
10表示封包通过节点2但没有通过节点3，
01表示封包通过节点3但没有通过节点2，
11表示封包节点2跟3皆有通过。 可计算
γ
i
{\displaystyle \gamma _{i}}
γ
1
{\displaystyle \gamma _{1}}
Ω
{\displaystyle \Omega }
γ
2
{\displaystyle \gamma _{2}}
Ω
{\displaystyle \Omega }
γ
3
{\displaystyle \gamma _{3}}
Ω
{\displaystyle \Omega }
则
γ
^
1
=
P
^
(
11
)
+
P
^
(
10
)
+
P
^
(
01
)
{\displaystyle {\hat {\gamma }}_{1}={\hat {P}}(11)+{\hat {P}}(10)+{\hat {P}}(01)}
γ
^
2
=
P
^
(
11
)
+
P
^
(
10
)
{\displaystyle {\hat {\gamma }}_{2}={\hat {P}}(11)+{\hat {P}}(10)}
γ
^
3
=
P
^
(
11
)
+
P
^
(
01
)
{\displaystyle {\hat {\gamma }}_{3}={\hat {P}}(11)+{\hat {P}}(01)}
γ
k
{\displaystyle \gamma _{k}}
α
{\displaystyle \alpha }
α
^
1
=
γ
^
2
γ
^
3
γ
^
2
+
γ
^
3
−
γ
^
1
=
(
P
^
(
01
)
+
P
^
(
11
)
)
(
P
^
(
10
)
+
P
^
(
11
)
)
P
^
(
11
)
{\displaystyle {\hat {\alpha }}_{1}={\dfrac {{\hat {\gamma }}_{2}{\hat {\gamma }}_{3}}{{\hat {\gamma }}_{2}+{\hat {\gamma }}_{3}-{\hat {\gamma }}_{1}}}={\dfrac {({\hat {P}}(01)+{\hat {P}}(11))({\hat {P}}(10)+{\hat {P}}(11))}{{\hat {P}}(11)}}}
α
^
2
=
γ
^
2
+
γ
^
3
−
γ
^
1
γ
^
3
=
P
^
(
11
)
P
^
(
01
)
+
P
^
(
11
)
{\displaystyle {\hat {\alpha }}_{2}={\dfrac {{\hat {\gamma }}_{2}+{\hat {\gamma }}_{3}-{\hat {\gamma }}_{1}}{{\hat {\gamma }}_{3}}}={\dfrac {{\hat {P}}(11)}{{\hat {P}}(01)+{\hat {P}}(11)}}}
α
^
3
=
γ
^
2
+
γ
^
3
−
γ
^
1
γ
^
2
=
P
^
(
11
)
P
^
(
10
)
+
P
^
(
11
)
{\displaystyle {\hat {\alpha }}_{3}={\dfrac {{\hat {\gamma }}_{2}+{\hat {\gamma }}_{3}-{\hat {\gamma }}_{1}}{{\hat {\gamma }}_{2}}}={\dfrac {{\hat {P}}(11)}{{\hat {P}}(10)+{\hat {P}}(11)}}}
![{\displaystyle \mathbf {\lambda } =E_{\lambda }[X|Y=AX]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2c269172d4cb693cc8ec56d7b4c424fdf22ef737)
![{\displaystyle \lambda ^{\left(n+1\right)}=E[X|Y,\lambda ^{\left(n\right)}]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1a8d5ada39135a816470402b37803a66e76a5851)
![{\displaystyle L(\alpha )=\log \,p(x^{1},x^{2},\cdots ,x^{n};\alpha )=\sum [{n(x)}log\,p(x;\alpha )]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3790bbbf76209841dfcbc0a883030a242df0df25)
![{\displaystyle {\hat {\alpha }}=argmax_{\alpha \in [0,1]}L(\alpha )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/90c49ebe29cd6a46b4caa620a59bf8ba473c7808)
![{\displaystyle \gamma _{k}=P[\Omega (k)]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f464e949ae78aa2a5b80ef6969ed40500b817b2d)
![{\displaystyle \Sigma \left[{\frac {n(x)}{n}}\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/968fbf934fd23c9c97f5ffb8e0a2bb51fdb30700)
![{\displaystyle \beta _{k}=P\,[\Omega (k)|X_{f(k)}=1]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f06e8f24d17092f43494e9b698d018208cfcde6a)
![{\displaystyle Q({\boldsymbol {\theta }}|{\boldsymbol {\theta }}^{(m-1)})=E_{\mathbf {Z} \mid \mathbf {X} ,{\boldsymbol {\theta }}^{(m-1)}}[logf(\mathbf {x} ,\mathbf {z} \mid {\boldsymbol {\theta }})]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/583d123e0bc060dbb21e94ff0b18f31b64dba0e5)
![{\displaystyle {\hat {\boldsymbol {\theta }}}=argmax_{{\boldsymbol {\theta }}\in {\boldsymbol {\Theta }}}E_{\mathbf {Z} \mid \mathbf {X} ,{\boldsymbol {\theta }}^{(m-1)}}[logf(\mathbf {x} ,\mathbf {z} \mid {\boldsymbol {\theta }})]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d8698c9e6f12100fb8e985eda232ac06f7cdfa19)