code: def get_marginal_topic_distrib(doc_topic_distrib, doc_lengths):
      
      unnorm = (doc_topic_distrib.T * doc_lengths).sum(axis=1)
      return unnorm / unnorm.sum()
docs: Return marginal topic distribution p(T) (topic proportions) given the document-topic distribution (theta)
      `doc_topic_distrib` and the document lengths `doc_lengths`. The latter can be calculated with `get_doc_lengths()`.