code: def pBIC(embedding: TSNEEmbedding) -> float:
      if not hasattr(embedding.affinities, "perplexity"):
          raise TypeError(
              "The embedding affinity matrix has no attribute `perplexity`"
          )
      n_samples = embedding.shape[0]
  
      return (
          2 * embedding.kl_divergence
          + np.log(n_samples)
          * embedding.affinities.perplexity
          / n_samples
      )
  
label_desc: Use image node_id 1 for calling a global function with example usage: pBIC(embedding) and returns: unknown