De esta respuesta en Stats Stackexchange https://stats.stackexchange.com/a/353715/274404 No estoy de acuerdo con la última parte.
\begin{align}AB^TBA^T &= \left(\left(AB^TBA^T\right)^T\right)^T \\ &= \left(BA^TAB^T\right)^T \\ &= B\left(BA^TA\right)^T \\ &= BAA^TB^T \end{align}
Si considero que B es una matriz cuadrada y A es la matriz identidad, implicaría que $ B^T B= B B^T$ . ¿Me he perdido algo?