bmtool 0.7.6__py3-none-any.whl → 0.7.8__py3-none-any.whl

bmtool/util/util.py

@@ -1163,9 +1163,43 @@ def percent_connections(
         num_sources = s_list[source_id_type].value_counts().sort_index().loc[source_id]
         num_targets = t_list[target_id_type].value_counts().sort_index().loc[target_id]
 
-        total = round(total_cons / (num_sources * num_targets) * 100, 2)
-        uni = round(num_uni / (num_sources * num_targets) * 100, 2)
-        bi = round(num_bi / (num_sources * num_targets) * 100, 2)
+        # Check if this is a recurrent network (same source and target population)
+        # For recurrent networks (e.g., FSI->FSI), we need special handling because:
+        # - Each pair can have at most 2 directed connections (bidirectional)
+        # - We want to count unique pairs, not directed connections
+        # - The denominator should be n*(n-1)/2, not n*n
+        is_recurrent = source_id == target_id
+        
+        if is_recurrent:
+            # For recurrent networks, calculate connectivity based on unique undirected pairs
+            # This avoids double-counting reciprocal connections and uses correct denominator
+            pair_counts = {}
+            for _, row in cons.iterrows():
+                sid = row['source_node_id']
+                tid = row['target_node_id']
+                if sid != tid:  # Exclude self-connections
+                    # Use symmetric pair key to count connections per unique pair
+                    pair_key = (min(sid, tid), max(sid, tid))
+                    if pair_key not in pair_counts:
+                        pair_counts[pair_key] = 0
+                    pair_counts[pair_key] += 1
+            
+            # Count pairs with exactly 1 connection (unidirectional) vs 2 connections (bidirectional)
+            num_uni = sum(1 for count in pair_counts.values() if count == 1)
+            num_bi = sum(1 for count in pair_counts.values() if count == 2)
+            
+            # Total possible unique pairs (excluding self-connections)
+            total_possible = num_sources * (num_sources - 1) / 2
+            total = round((num_uni + num_bi) / total_possible * 100, 2)
+            uni = round(num_uni / total_possible * 100, 2)
+            bi = round(num_bi / total_possible * 100, 2)
+        else:
+            # For non-recurrent networks, use the original calculation
+            # Each connection is unique (no double-counting issues)
+            total = round(total_cons / (num_sources * num_targets) * 100, 2)
+            uni = round(num_uni / (num_sources * num_targets) * 100, 2)
+            bi = round(num_bi / (num_sources * num_targets) * 100, 2)
+        
         if method == "total":
             return total
         if method == "uni":
@@ -1309,9 +1343,10 @@ def gap_junction_connections(
         s_list = kwargs["source_nodes"]
 
         cons = edges[(edges[source_id_type] == source_id) & (edges[target_id_type] == target_id)]
-        # add functionality that shows only the one's with gap_junctions
+        # print(cons)
+
         try:
-            cons = cons[cons["is_gap_junction"]]
+            cons = cons[cons["is_gap_junction"]==True]
         except:
             raise Exception("no gap junctions found to drop from connections")
 

{bmtool-0.7.6.dist-info → bmtool-0.7.8.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: bmtool
-Version: 0.7.6
+Version: 0.7.8
 Summary: BMTool
 Home-page: https://github.com/cyneuro/bmtool
 Download-URL: 

{bmtool-0.7.6.dist-info → bmtool-0.7.8.dist-info}/RECORD

@@ -1,34 +1,34 @@
 bmtool/SLURM.py,sha256=UBfITY1MtYo95nyKglgGSqAC9Ds8PBvlHczsiNMFxvc,20573
 bmtool/__init__.py,sha256=r_8fXc-2uj1DndCdhB4jME51r1pn6ESTD5zRc355BrU,134
 bmtool/__main__.py,sha256=uCEqPwRxIuNRUASKhsvh4S8Nkp4dqnvfXTMUv-wWWRU,665
-bmtool/connectors.py,sha256=4SJRqBXM145_-nFycGrnlfjSyaoarOr0QkHl00jWq4I,74384
+bmtool/connectors.py,sha256=tkCh9oVsIgaV5zzBzTFTG2Y6YVAoIHv2Pm-qxZAbca8,93124
 bmtool/graphs.py,sha256=gBTzI6c2BBK49dWGcfWh9c56TAooyn-KaiEy0Im1HcI,6717
 bmtool/manage.py,sha256=lsgRejp02P-x6QpA7SXcyXdalPhRmypoviIA2uAitQs,608
 bmtool/plot_commands.py,sha256=Dxm_RaT4CtHnfsltTtUopJ4KVbfhxtktEB_b7bFEXII,12716
-bmtool/singlecell.py,sha256=xqdLM2TjjnL8YyTy-c3WR6mElTv3E4zkZgkfUh4S5X0,47161
-bmtool/synapses.py,sha256=GbnPPdXotl7MTeju_QTBwZfFxYFNCq0fJijntzfttds,106071
+bmtool/singlecell.py,sha256=KFCh9qN3bvOk_B8WArBHlp1fBDb4m3yCZQOijkhCbpA,67120
+bmtool/synapses.py,sha256=gMJFG-3XwCzbLX4o_xObQ3oVs1CU_5L6WgxDQupiPGI,157967
 bmtool/analysis/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-bmtool/analysis/entrainment.py,sha256=NQloQtVpEWjDzmkZwMWVcm3hSjErHBZfQl1mrBVoIE8,25321
-bmtool/analysis/lfp.py,sha256=S2JvxkjcK3-EH93wCrhqNSFY6cX7fOq74pz64ibHKrc,26556
+bmtool/analysis/entrainment.py,sha256=l_loxU3GEJX5cDjgMEj3vpYkA5G78r3DqsLXH6o6HS0,29542
+bmtool/analysis/lfp.py,sha256=b3uajzE_hyqJGhMM7dIBYerHusJ58XLQ_hA09uwXEUY,26570
 bmtool/analysis/netcon_reports.py,sha256=VnPZNKPaQA7oh1q9cIatsqQudm4cOtzNtbGPXoiDCD0,2909
 bmtool/analysis/spikes.py,sha256=3n-xmyEZ7w6CKEND7-aKOAvdDg0lwDuPI5sMdOuPwa0,24637
 bmtool/bmplot/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-bmtool/bmplot/connections.py,sha256=xACgnkLAZrEBpaxa9Ob9MX-B1DO782UAeMvih9DleDk,59400
-bmtool/bmplot/entrainment.py,sha256=BrBMerqyiG2YWAO_OEFv7OJf3yeFz3l9jUt4NamluLc,32837
-bmtool/bmplot/lfp.py,sha256=7JLozQQJ19ty0ZNyfhkuJAr_K8_pVP9C0flVJd_YXaY,2027
-bmtool/bmplot/netcon_reports.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-bmtool/bmplot/spikes.py,sha256=odzCSMbFRHp9qthSGQ0WzMWUwNQ7R1Z6gLT6VPF_o5Q,15326
+bmtool/bmplot/connections.py,sha256=qUtkq8Sy3SP6IwVWKx6h95d74qy_fMGtYPw7xr5bnvw,73261
+bmtool/bmplot/entrainment.py,sha256=q6eyftK3TUaLVzukZ1gBWTt-JoQq_69K138TjBxJghU,37046
+bmtool/bmplot/lfp.py,sha256=tHRPqAe8npwiTIQSzBHfZl_PPxqcTDZpIPdn6eHwRbI,7253
+bmtool/bmplot/netcon_reports.py,sha256=WkqQHqTpNxgDutcSgAJDtGhzp1LXO47XOiXzTT-6Pfk,22
+bmtool/bmplot/spikes.py,sha256=gxPfxe0o4Jzo__TK6HgIm1zJoddEL8wDU3IhxJkeA54,21707
 bmtool/debug/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 bmtool/debug/commands.py,sha256=VV00f6q5gzZI503vUPeG40ABLLen0bw_k4-EX-H5WZE,580
 bmtool/debug/debug.py,sha256=9yUFvA4_Bl-x9s29quIEG3pY-S8hNJF3RKBfRBHCl28,208
 bmtool/util/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 bmtool/util/commands.py,sha256=Nn-R-4e9g8ZhSPZvTkr38xeKRPfEMANB9Lugppj82UI,68564
-bmtool/util/util.py,sha256=HbVgPpvDBSkvXv5gAz7ad1V5GKL26q_LyJGt6ehpGQE,75929
+bmtool/util/util.py,sha256=oXrY-Cc4wxPvgJVg2TqaAPua_gAl1Kj1nqrLgsqsXS8,77815
 bmtool/util/neuron/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 bmtool/util/neuron/celltuner.py,sha256=lokRLUM1rsdSYBYrNbLBBo39j14mm8TBNVNRnSlhHCk,94868
-bmtool-0.7.6.dist-info/licenses/LICENSE,sha256=qrXg2jj6kz5d0EnN11hllcQt2fcWVNumx0xNbV05nyM,1068
-bmtool-0.7.6.dist-info/METADATA,sha256=iSIT7erBJsvbZ6vJE1LKOlQPF1ne1BiluVL-81Z8l4U,3595
-bmtool-0.7.6.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-bmtool-0.7.6.dist-info/entry_points.txt,sha256=0-BHZ6nUnh0twWw9SXNTiRmKjDnb1VO2DfG_-oprhAc,45
-bmtool-0.7.6.dist-info/top_level.txt,sha256=gpd2Sj-L9tWbuJEd5E8C8S8XkNm5yUE76klUYcM-eWM,7
-bmtool-0.7.6.dist-info/RECORD,,
+bmtool-0.7.8.dist-info/licenses/LICENSE,sha256=qrXg2jj6kz5d0EnN11hllcQt2fcWVNumx0xNbV05nyM,1068
+bmtool-0.7.8.dist-info/METADATA,sha256=hkf6T0q5HcIzDBSstste4hoAzJigGRTfPuGoEDP-V84,3595
+bmtool-0.7.8.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+bmtool-0.7.8.dist-info/entry_points.txt,sha256=0-BHZ6nUnh0twWw9SXNTiRmKjDnb1VO2DfG_-oprhAc,45
+bmtool-0.7.8.dist-info/top_level.txt,sha256=gpd2Sj-L9tWbuJEd5E8C8S8XkNm5yUE76klUYcM-eWM,7
+bmtool-0.7.8.dist-info/RECORD,,