Factoring positive operators on reproducing kernel Hilbert spaces

We characterize when positive operators can be factored by analytic Toeplitz type operators. As a corollary, we give an operator theory characterization of those invariant subspaces of doubly commuting unilateral shifts, which are generated by a single inner function on the bidisk. The last result extends to shifts of arbitrary (countable) multiplicity.