code: def get_relative_abundance(biomfile):
      
      biomf = biom.load_table(biomfile)
      norm_biomf = biomf.norm(inplace=False)
      rel_abd = {}
      for sid in norm_biomf.ids():
          rel_abd[sid] = {}
          for otuid in norm_biomf.ids("observation"):
              otuname = oc.otu_name(norm_biomf.metadata(otuid, axis="observation")["taxonomy"])
              otuname = " ".join(otuname.split("_"))
              abd = norm_biomf.get_value_by_ids(otuid, sid)
              rel_abd[sid][otuname] = abd
      ast_rel_abd = bc.arcsine_sqrt_transform(rel_abd)
      return ast_rel_abd
docs: Return arcsine transformed relative abundance from a BIOM format file.
  
      :type biomfile: BIOM format file
      :param biomfile: BIOM format file used to obtain relative abundances for each OTU in
                       a SampleID, which are used as node sizes in network plots.
  
      :type return: Dictionary of dictionaries.
      :return: Dictionary keyed on SampleID whose value is a dictionarykeyed on OTU Name
               whose value is the arc sine tranfsormed relative abundance value for that
               SampleID-OTU Name pair.