The influence of various slip system combination models on crystal plasticity finite element simulation of NiTi shape memory alloy subjected to uniaxial compression deformation is investigated according to three combinations of slip systems, including combination of {010}〈100〉 and {110}〈111〉 slip modes, combination of {110}〈100〉 and {110}〈111〉 slip modes and combination of {110}〈100〉, {010}〈100〉 and {110}〈111〉 slip modes, which consist of 18, 18 and 24 slip systems, respectively. By means of simulating mechanical response, strain distribution, stress distribution and Schmid factor, it can be found that in terms of simulation accuracy, combination of {110}〈100〉 and {110}〈111〉 slip modes is in good agreement with combination of {110}〈100〉, {010}〈100〉 and {110}〈111〉 slip modes. The contribution of {110}〈100〉 slip mode to plastic strain is primary in plastic deformation of NiTi shape memory alloy, whereas {010}〈100〉 slip mode, which makes small contribution to plastic deformation, can be regarded as the unfavorable slip mode. In the case of large plastic strain, the {010}〈100〉 slip mode contributes to the formation of (001) [01¯0] texture component, while {110}〈100〉 and {110}〈111〉 slip modes facilitate the formation of γ-fibre (〈111〉) texture. Keywords: Shape memory alloy, Plastic deformation, Crystal plasticity, Finite element method, Texture