code: def raw_abundance(biomf, sampleIDs=None, sample_abd=True):
      
      results = defaultdict(int)
      if sampleIDs is None:
          sampleIDs = biomf.ids()
      else:
          try:
              for sid in sampleIDs:
                  assert sid in biomf.ids()
          except AssertionError:
              raise ValueError(
                  "\nError while calculating raw total abundances: The sampleIDs provided "
                  "do not match the sampleIDs in biom file. Please double check the "
                  "sampleIDs provided.\n")
      otuIDs = biomf.ids(axis="observation")
  
      for sampleID in sampleIDs:
          for otuID in otuIDs:
              abd = biomf.get_value_by_ids(otuID, sampleID)
              if sample_abd:
                  results[sampleID] += abd
              else:
                  results[otuID] += abd
      return results
docs: Calculate the total number of sequences in each OTU or SampleID.
  
      :type biomf: A BIOM file.
      :param biomf: OTU table format.
  
      :type sampleIDs: List
      :param sampleIDs: A list of column id's from BIOM format OTU table. By default, the
                        list has been set to None.
  
      :type sample_abd: Boolean
      :param sample_abd: A boolean operator to provide output for OTUID's or SampleID's. By
                         default, the output will be provided for SampleID's.
  
      :rtype: dict
      :return: Returns a dictionary keyed on either OTUID's or SampleIDs and their
               respective abundance as values.