We investigate roles of electron correlation effects in the determination of the gj factors of the 4s 2 S 1/2 , 4p 2 P 1/2 , 4p 2 P 3/2 , 3d 2 D 3/2 , and 3d 2 D 5/2 states, representing to different parities and angular momenta, of the Ca + ion. Correlation contributions are highlighted with respect to the mean-field values evaluated using the Dirac-Hartree-Fock method, relativistic second order manybody theory, and relativistic coupled-cluster (RCC) theory with the singles and doubles approximation considering only the linear terms and also accounting for all the non-linear terms. This shows that it is difficult to achieve reasonably accurate results employing an approximated perturbative approach. We also find that contributions through the non-linear terms and higher-level excitations such as triple excitations, estimated perturbatively in the RCC method, are found to be crucial to attain precise values of the gj factors in the considered states of Ca + ion.