def topic_jsds(self, print_len=10):
        topic_indices = np.arange(self.model.K)

        doc_tops = self.theta.T
        M = np.sum(doc_tops, axis=0) / len(doc_tops)
        pjsd = np.sum(np.array([(D_i * np.log(D_i / M)) / len(doc_tops) 
                                    for D_i in doc_tops]), axis=0)

        k_arr = enum_sort(pjsd).view(LabeledColumn)
        k_arr.col_header = 'Topic Partial JSD'
        k_arr.subcol_headers = ['Index', 'pJSD']
        k_arr.col_len = 10
        return k_arr[::-1]
Returns the partial N-way JSD of each topic, where N is the number of
dataset/ETHPy150Open inpho/vsm/vsm/viewer/ldacgsviewer.py/LdaCgsViewer.topic_jsds