We investigated the reliability and applicability of so-calledmagnetic-force linear responsemethods to calculate spin-spin interaction strengths from first principles. We examined the dependence on the numerical parameters, including the number of basis orbitals and their cutoff radii within nonorthogonal LCPAO (linear combination of pseudoatomic orbitals) formalism. It is shown that the parameter dependence and the ambiguity caused by these choices are small enough in comparison to the other computation approach and experiments. Furthermore, we tried to pursue the broader applicability of this technique. We showed that magnetic-force theory can provide a reasonable prediction, especially for the case of strongly localized moments, even when the ground-state configuration is unknown or the total-energy value is not accessible. The formalism is extended for LCPAO to carry the orbital resolution from which the matrix form of the magnetic coupling constant is calculated. From the applications to Fe-based superconductors including LaFeAsO, NaFeAs, BaFe2As2, and FeTe, the distinctive characteristics of orbital-resolved interactions are clearly noticed in between single-stripe pnictides and double-stripe chalcogenides.