code: def max(self, array, role = None):
          
          return self.reduce(array, reducer = np.maximum, neutral_element = - np.infty, role = role)
docs: Return the maximum value of ``array`` for the entity members.
  
              ``array`` must have the dimension of the number of persons in the simulation
  
              If ``role`` is provided, only the entity member with the given role are taken into account.
  
              Example:
  
              >>> salaries = household.members('salary', '2018-01')  # e.g. [2000, 1500, 0, 0, 0]
              >>> household.max(salaries)
              >>> array([2000])