I came across an equation that I don't understand how it was derived. The equation was calculating the expectation of the covariance matrix of an error.
The error is $\tilde\theta_N$: $$ \tilde\theta_N=R(N)^{-1}\sum_{t=1}^N\varphi(t)e(t) $$ where $$ \varphi(t) = \left[u(t-1)~~~u(t-2)~~~...~~~u(t-m)\right]^T $$ $u(t)$ is a time-series signal and the number $m$ is unknown at this time.
$e(t)$ is a white noise sequence with variance $\lambda$. According to the book,
$e(t)$ can be described as a sequence of independent random variables with zero mean values and variances $\lambda$.
$$ R(N) = \sum_{t=1}^N\varphi(t)\varphi^T(t) $$
The expectation of its covariance matrix is $P_N$: $$ P_N=E\tilde\theta_N\tilde\theta_N^T=ER(N)^{-1}\sum_{t,s=1}^N\varphi(t)e(t)e(s)\varphi^T(s)R(N)^{-1} $$ My question is: how did the equation for $P_N$ came along?
This is what I've got so far:
$$ P_N = E\tilde\theta_N\tilde\theta_N^T = ER(N)^{-1}\sum_{t=1}^N\varphi(t)e(t)\left[\sum_{s=1}^N\varphi(s)e(s)\right]^T\left(R^{-1}(N)\right)^T $$