search

gflow / GFlow

Software for modeling circuit theory-based connectivity
GNU General Public License v3.0
22 stars 5 forks source link

Issue calculating pairwise interactions #25

Closed famarques closed 6 years ago

famarques commented 6 years ago

Hello, i´m using GFlow to calculate connectivity beetwen points in europe-wide lenght. I´ve changed the execute_example.sh flags to calculate pairwise interactions beetween all nodes. I tried a lot of changes on the flags but in all GFlow stops calculating before all pairwise nodes been calculated. In annex the files that i´m using and the execute_example.sh file modified with my configurations.

I´m doing something wrong?

Thank you,

GFlow_data.zip

pbleonard commented 6 years ago

@famarques Can you be more specific about

GFlow stops calculating before all pairwise nodes been calculated.

I attempted to execute your example and noticed many of the points appear to be invalid.

famarques commented 6 years ago

Thanks for response, i´m trying to calculate all interactions between points like circuitscape, but the program stops and write the results before all pairs are calculated. For example in the species of those nodes i´ve points in UK that are not showed in the results.

Yes, but It happened also before it another species that i´ve tried to calculate, and have all or almost all points valid.

The points and resistance raster: resistance 1

The results: myo_dau_cnc 1

pbleonard commented 6 years ago

If you chose not to include the converge_at flag then all pairwise will be calculated unless there are > 1 million pairs (in this case we are aware of a bug). gflow will not calculate any pairwise connection with infinite resistance or with pairs that are not constrained inside the resistance extent. Can you tell me how you're making your points text file?

famarques commented 6 years ago

maybe it is because i have more than a million pairs. i´m making the points with a asc file using the script posted by you.

pbleonard commented 6 years ago

Great. Can you test with a subset of your pairs?

famarques commented 6 years ago

I will test it, now im working in another thing, but when i have time i will try it

pbleonard commented 6 years ago

@eduffy fixed the bug here (should already be included in release v. 0.1.7):