fix bug: modindices() failed when parameterization = "theta"