Group Abstract Group Abstract

Message Boards Message Boards

[WSS20] Implementing Mutual Information

Posted 5 years ago

Karman vortex street and an estimation of mutual information

POSTED BY: Luigi Brancati
7 Replies
POSTED BY: Clarisse Wagner
Posted 5 years ago
POSTED BY: Syd Geraghty

enter image description here -- you have earned Featured Contributor Badge enter image description here

Your exceptional post has been selected for our editorial column Staff Picks http://wolfr.am/StaffPicks and Your Profile is now distinguished by a Featured Contributor Badge and is displayed on the Featured Contributor Board. Thank you!

POSTED BY: EDITORIAL BOARD
Posted 5 years ago
POSTED BY: Luigi Brancati
Posted 5 years ago

Thanks so much, Luigi. I think I am able to get it running. But to make certain, may ask you if your Community post syntax:

KraskovI1[list, list+T, 4], with T ranged 0,20

is equivalent to your new ResourceFunction:

ResourceFunction["MutualInformation"][list,20]

A smaller question, I wonder if there's a small index-shift problem? If I do

ListLinePlot[output from your MutualInformation] 

the plot will show the zero-shift MI at position 1 (not zero) in the plot (since the zero shift is the first element in the MI list). The plots in the post were probably shifted by one element from actual value (e.g, for minimum in MI). I suppose this could be fixed if I do my plot shifted so 1st element in MI list is plotted at origin, at 0.

Thanks so much. You have probably moved on from this project, but I am finding it potentially very useful in my Physics research. And I would like to acknowledge your work if/when I ever publish what I'm looking at.

-Jim Clarage, University of St. Thomas Houston, Texas

POSTED BY: Updating Name
Posted 5 years ago

Where is the code for your function "KraskovI1" ? It is missing from your post and notebook. If possible I'd like to use this function, or some function to calculate the Mutual Information of a time series.

Thanks.

POSTED BY: Jim Clarage
Posted 5 years ago

Hi, thanks for the interest! My function is currently under review for submission in the Wolfram Function Repository, but you can find a preliminary (not reviewed) version here. It doesn't explicitely support timeseries yet, but I plan on adding that soon.

(Ignore the comment of the reviewer in the examples section, I just resubmitted it and couldn't delete that)

POSTED BY: Luigi Brancati
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
Attachments
Remove
or Discard