[step:Prove the exterior-algebra adjoint identity on basis forms]
For $1\le j\le n$, define exterior multiplication
\begin{align*}
\varepsilon_j:\Lambda^{p,q-1}&\to\Lambda^{p,q}\\
\alpha&\mapsto d\bar z_j\wedge\alpha,
\end{align*}
and contraction
\begin{align*}
\iota_j:\Lambda^{p,q}&\to\Lambda^{p,q-1}\\
\beta&\mapsto \iota_{\partial/\partial\bar z_j}\beta.
\end{align*}
The contraction $\iota_j$ is the graded derivation of degree $-1$ determined on coordinate one-forms, for $1\le \ell\le n$, by
\begin{align*}
\iota_j(dz_\ell)&=0,\\
\iota_j(d\bar z_\ell)&=
\begin{cases}
1, & \ell=j,\\
0, & \ell\ne j,
\end{cases}
\end{align*}
and, for homogeneous forms $\gamma$ and $\delta$ with total exterior degree $|\gamma|$, by
\begin{align*}
\iota_j(\gamma\wedge\delta)
=
\iota_j\gamma\wedge\delta
+(-1)^{|\gamma|}\gamma\wedge\iota_j\delta.
\end{align*}
We prove that $\iota_j$ is the adjoint of $\varepsilon_j$ with respect to the Euclidean Hermitian inner products:
\begin{align*}
(\varepsilon_j\alpha,\beta)_{\Lambda^{p,q}}
=
(\alpha,\iota_j\beta)_{\Lambda^{p,q-1}}.
\end{align*}
It suffices to check basis vectors. Fix $I,I'\in\mathcal I_p$, $J\in\mathcal I_{q-1}$, and $K\in\mathcal I_q$. If $j\in J$, then $d\bar z_j\wedge d\bar z_J=0$, so $\varepsilon_j e_{I,J}=0$. Also $K$ cannot be the increasing rearrangement of $\{j\}\cup J$ because this set has a repeated element, and hence the coefficient of $e_{I,J}$ in $\iota_j e_{I',K}$ is zero. Both inner products vanish.
Assume now that $j\notin J$. Let $K(j,J)\in\mathcal I_q$ be the increasing rearrangement of the set $\{j,j_1,\dots,j_{q-1}\}$. Define $\sigma(j,J)\in\{1,-1\}$ by
\begin{align*}
d\bar z_j\wedge d\bar z_J=\sigma(j,J)d\bar z_{K(j,J)}.
\end{align*}
Since $d\bar z_j$ must first pass across the $p$ holomorphic factors in $dz_I$, set
\begin{align*}
\tau(j,J)=(-1)^p\sigma(j,J).
\end{align*}
Then
\begin{align*}
\varepsilon_j e_{I,J}
&=d\bar z_j\wedge dz_I\wedge d\bar z_J\\
&=\tau(j,J)e_{I,K(j,J)}.
\end{align*}
By the defining rule for contraction with $\partial/\partial\bar z_j$,
\begin{align*}
\iota_j e_{I',K}
=
\begin{cases}
\tau(j,J)e_{I',J}, & K=K(j,J),\\
0, & K\ne K(j,J).
\end{cases}
\end{align*}
Therefore
\begin{align*}
(\varepsilon_j e_{I,J},e_{I',K})_{\Lambda^{p,q}}
&=
\begin{cases}
\tau(j,J), & I=I'\text{ and }K=K(j,J),\\
0, & \text{otherwise},
\end{cases}\\
(e_{I,J},\iota_j e_{I',K})_{\Lambda^{p,q-1}}
&=
\begin{cases}
\tau(j,J), & I=I'\text{ and }K=K(j,J),\\
0, & \text{otherwise}.
\end{cases}
\end{align*}
The two expressions agree on basis vectors, and sesquilinearity gives the identity for all $\alpha\in\Lambda^{p,q-1}$ and $\beta\in\Lambda^{p,q}$.
[/step]