Skip to content

Riemannian Measure μ ? #24

Description

@Noeloikeau

Hello,

Thank you for the wonderful library! I have a question regarding how to obtain the Riemannian Measure / Volume Element μ, also called the leading-left eigenvector of the transition probability matrix μP = μ.

In the original NLSA paper page 3, column 2, paragraph 2 and algorithm page 4, line 40, μ is given explicitly.

Am I correct in thinking that this quantity is calculated as the reciprocal of row_sum on line 136 of diffusion_map.py during the _left_normalize function, and subsequently discarded?

Thank you for any assistance.

Activity

  1. ehthiede commented on Apr 27, 2024

    @ehthiede
    Member

    Hi,

    Realizing I never responded to this: terribly sorry. I think you are correct in your assessment.

    For fixed bandwith kernels I think you are correct. Then the generator is given by (D_p^{-1} K - I) / e, where p is the row-sum of K. It is easy to see that p is then a left-eigenvector of the generator. For variable-bandwidth kernels we later left-multiply by (1/h(x)^2) where h(x) is the bandwidth function. Then, to construct the density you would need to take the output of _left_normalize and divide by h(x)^2, I believe. Hope that helps, and sorry again for the delay.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

      Milestone

      No milestone

      Relationships

      None yet

      Development

      No branches or pull requests

      Issue actions