Solution (source code)

= Solution

Let
$$
v=p_1(1-p_1)+p_0(1-p_0),
\qquad I_j=\frac{m_j}{v},
$$
where $m_j$ is the cumulative sample size per arm. The <canonical joint distribution for group sequential test statistics> is
$$
\begin{pmatrix}Z_1\\Z_2\end{pmatrix}
\dot\sim N_2\!\left[
\begin{pmatrix}\delta\sqrt{I_1}\\\delta\sqrt{I_2}\end{pmatrix},
\begin{pmatrix}
1&\sqrt{I_1/I_2}\\
\sqrt{I_1/I_2}&1
\end{pmatrix}
\right].
$$
Thus $\operatorname{Corr}(Z_1,Z_2)=\sqrt{m_1/m_2}$. This correlation arises because the second statistic reuses all first-stage observations.