Extract columns 1, 3, 4
DataxzStein = DataxyzStein[[All, {1, 3, 4}]]
FittXZStein = 
  NonlinearModelFit[
   DataxzStein, {ModelxzStein}, {Subscript[T, 0], Subscript[T, 2], 
    Subscript[T, 4], Subscript[T, 6], Subscript[T, 8], Subscript[T, 10]}, {x, z}]
FittXZStein["RSquared"]
(* 0.9845 *)
FittXZStein["ParameterTable"]
