Closed bcipolli closed 8 years ago
I liked doing R- or L-only to both rather than RL to both precisely because for any lateralized whole-brain component will have a hard time getting a good match with artificially concatenated RL. By doing R- or L-only to both, some lateralized wb components will have a good match with one hemi but not the other, and we get to see that information by comparing R- to both and L-to both matrices.
Separately from that issue, I was thinking about whether we should be forcing one-to-one matches. I was actually playing with adding a histogram of all the dissimilarity scores to see whether we can find a good cut-off of scores to determine match vs non-match, because I'm not sure how meaningful it is to grab the best match when a component really doesn't have any good corresponding component. On the other hand, one wb component could split up in two separate sub-components in R- or L- only, in which case that component might have two different more or less equally good (but not perfect) matches.
So my thought was to either to come up with a threshold based on the distribution of the scores or just take the best 20 matches, for example, allowing for reusing of components...
Sounds good :+1: . I will also open an issue to rename both
to wb
;)
Again, this will be made easier by an interactive plot (e.g. #21 ). If you can make a PR with a histogram of similarity scores, I think that'd be cool!
I'm going to close this since the main.py in newRLmatch now let you choose unique (i.e. one-to-one matching) or allow reusing of matches with --force_match option.
@bcipolli:
@atsuch: