def nonzeroCoef(beta, bystep = False): result = scipy.absolute(beta) > 0 if not bystep: result = scipy.any(result, axis = 1) return(result) # end of nonzeroCoef # =========================================