We present the new skein invariants of classical links, H [ H ] , K [ K ] and D [ D ] , based on the invariants of links, H, K and D, denoting the regular isotopy version of the Homflypt polynomial, the Kauffman polynomial and the Dubrovnik polynomial. The invariants are obtained by abstracting the skein relation of the corresponding invariant and making a new skein algorithm comprising two computational levels: first producing unlinked knotted components, then evaluating the resulting knots. The invariants in this paper, were revealed through the skein theoretic definition of the invariants Θ d related to the Yokonuma–Hecke algebras and their 3-variable generalization Θ , which generalizes the Homflypt polynomial. H [ H ] is the regular isotopy counterpart of Θ . The invariants K [ K ] and D [ D ] are new generalizations of the Kauffman and the Dubrovnik polynomials. We sketch skein theoretic proofs of the well-definedness and topological properties of these invariants. The invariants of this paper are reformulated into summations of the generating invariants (H, K, D) on sublinks of the given link L, obtained by partitioning L into collections of sublinks. The first such reformulation was achieved by W.B.R. Lickorish for the invariant Θ and we generalize it to the Kauffman and Dubrovnik polynomial cases. State sum models are formulated for all the invariants. These state summation models are based on our skein template algorithm which formalizes the skein theoretic process as an analogue of a statistical mechanics partition function. Relationships with statistical mechanics models are articulated. Finally, we discuss physical situations where a multi-leveled course of action is taken naturally.