Observed-link-flow-based OD matrix estimation with stochastic factors