We perform the global analysis of polarized semi-inclusive deep inelastic scattering (SIDIS), pion-induced polarized Drell-Yan (DY) process, and W-+/-/Z boson production data and extract the Sivers function for u, d, s and for sea quarks. We use the framework of transverse momentum dependent factorization at next-to-next-to-next-to leading order ((NLO)-L-3) accuracy. The Qiu-Sterman function is determined in a model independent way from the extracted Sivers function. We also evaluate the significance of the predicted sign change of the Sivers function in DY process with respect to SIDIS.