def raw_abundance(biomf, sampleIDs=None, sample_abd=True):
results = defaultdict(int)<EOL>if sampleIDs is None:<EOL><INDENT>sampleIDs = biomf.ids()<EOL><DEDENT>else:<EOL><INDENT>try:<EOL><INDENT>for sid in sampleIDs:<EOL><INDENT>assert sid in biomf.ids()<EOL><DEDENT><DEDENT>except AssertionError:<EOL><INDENT>raise ValueError(<EOL>"<STR_LIT>"<EOL>"<STR_LIT>"<EOL>"<STR_LIT>")<EOL><DEDENT><DEDENT>otuIDs = biomf.ids(axis="<STR_LIT>")<EOL>for sampleID in sampleIDs:<EOL><INDENT>for otuID in otuIDs:<EOL><INDENT>abd = biomf.get_value_by_ids(otuID, sampleID)<EOL>if sample_abd:<EOL><INDENT>results[sampleID] += abd<EOL><DEDENT>else:<EOL><INDENT>results[otuID] += abd<EOL><DEDENT><DEDENT><DEDENT>return results<EOL>
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.