notebook.community

Edit and run 





In [1]:



    


import networkx as nx
import numpy as np



    


















    






Couldn't import dot_parser, loading of dot files will not be possible.

















In [13]:



    


file = "/Users/mark/seriationct/experiments/seriationct-14/rawnetworkmodels/seriationct-14-networkmodel/seriationct-14-002.gml"



    











In [14]:



    


g_cur = nx.read_gml(file)



    











In [11]:



    


def calc_migration_matrix(g_cur, migration_fraction):
    g_mat = nx.to_numpy_matrix(g_cur)

    # get the column totals
    rtot = np.sum(g_mat, axis = 1)
    scaled = (g_mat / rtot) * migration_fraction
    diag = np.eye(np.shape(g_mat)[0]) * (1.0 - migration_fraction)
    g_mat_scaled = diag + scaled

    print g_mat
    
    return g_mat_scaled.tolist()



    











In [15]:



    


calc_migration_matrix(g_cur, 0.2)



    


















    






[[  0.  10.  10.  10.  10.  10.  10.  10.  10.  10.   0.   0.   0.   0.
    0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.
    0.   0.]
 [ 10.   0.  10.  10.  10.  10.  10.  10.  10.  10.   0.   0.   0.   0.
    0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.
    0.   0.]
 [ 10.  10.   0.  10.  10.  10.  10.  10.  10.  10.   0.   0.   0.   0.
    0.   0.   0.   0.   1.   0.   0.   0.   0.   0.   0.   0.   0.   0.
    0.   0.]
 [ 10.  10.  10.   0.  10.  10.  10.  10.  10.  10.   0.   0.   0.   0.
    0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.
    0.   0.]
 [ 10.  10.  10.  10.   0.  10.  10.  10.  10.  10.   0.   0.   0.   1.
    0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.
    0.   0.]
 [ 10.  10.  10.  10.  10.   0.  10.  10.  10.  10.   0.   0.   0.   0.
    0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.
    0.   0.]
 [ 10.  10.  10.  10.  10.  10.   0.  10.  10.  10.   0.   0.   0.   0.
    0.   0.   0.   0.   0.   0.   0.   0.   1.   0.   0.   0.   0.   0.
    0.   0.]
 [ 10.  10.  10.  10.  10.  10.  10.   0.  10.  10.   0.   0.   0.   0.
    0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.
    0.   0.]
 [ 10.  10.  10.  10.  10.  10.  10.  10.   0.  10.   0.   0.   0.   0.
    0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   1.   0.
    0.   0.]
 [ 10.  10.  10.  10.  10.  10.  10.  10.  10.   0.   0.   0.   0.   0.
    0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.
    0.   0.]
 [  0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.  10.  10.  10.
   10.  10.  10.  10.  10.  10.   0.   0.   0.   0.   0.   0.   0.   0.
    0.   0.]
 [  0.   0.   0.   0.   0.   0.   0.   0.   0.   0.  10.   0.  10.  10.
   10.  10.  10.  10.  10.  10.   0.   0.   0.   0.   0.   0.   0.   0.
    0.   0.]
 [  0.   0.   0.   0.   0.   0.   0.   0.   0.   0.  10.  10.   0.  10.
   10.  10.  10.  10.  10.  10.   0.   0.   0.   0.   0.   0.   0.   0.
    1.   0.]
 [  0.   0.   0.   0.   1.   0.   0.   0.   0.   0.  10.  10.  10.   0.
   10.  10.  10.  10.  10.  10.   0.   0.   0.   0.   0.   0.   0.   0.
    0.   0.]
 [  0.   0.   0.   0.   0.   0.   0.   0.   0.   0.  10.  10.  10.  10.
    0.  10.  10.  10.  10.  10.   0.   0.   1.   0.   0.   0.   0.   0.
    0.   0.]
 [  0.   0.   0.   0.   0.   0.   0.   0.   0.   0.  10.  10.  10.  10.
   10.   0.  10.  10.  10.  10.   0.   0.   0.   0.   0.   0.   0.   0.
    0.   0.]
 [  0.   0.   0.   0.   0.   0.   0.   0.   0.   0.  10.  10.  10.  10.
   10.  10.   0.  10.  10.  10.   0.   0.   0.   0.   0.   0.   0.   0.
    0.   0.]
 [  0.   0.   0.   0.   0.   0.   0.   0.   0.   0.  10.  10.  10.  10.
   10.  10.  10.   0.  10.  10.   0.   0.   0.   0.   0.   0.   0.   0.
    0.   0.]
 [  0.   0.   1.   0.   0.   0.   0.   0.   0.   0.  10.  10.  10.  10.
   10.  10.  10.  10.   0.  10.   0.   0.   0.   0.   0.   0.   0.   0.
    0.   0.]
 [  0.   0.   0.   0.   0.   0.   0.   0.   0.   0.  10.  10.  10.  10.
   10.  10.  10.  10.  10.   0.   0.   0.   0.   0.   0.   0.   0.   0.
    0.   0.]
 [  0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.
    0.   0.   0.   0.   0.   0.   0.  10.  10.  10.  10.  10.  10.  10.
   10.  10.]
 [  0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.
    0.   0.   0.   0.   0.   0.  10.   0.  10.  10.  10.  10.  10.  10.
   10.  10.]
 [  0.   0.   0.   0.   0.   0.   1.   0.   0.   0.   0.   0.   0.   0.
    1.   0.   0.   0.   0.   0.  10.  10.   0.  10.  10.  10.  10.  10.
   10.  10.]
 [  0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.
    0.   0.   0.   0.   0.   0.  10.  10.  10.   0.  10.  10.  10.  10.
   10.  10.]
 [  0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.
    0.   0.   0.   0.   0.   0.  10.  10.  10.  10.   0.  10.  10.  10.
   10.  10.]
 [  0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.
    0.   0.   0.   0.   0.   0.  10.  10.  10.  10.  10.   0.  10.  10.
   10.  10.]
 [  0.   0.   0.   0.   0.   0.   0.   0.   1.   0.   0.   0.   0.   0.
    0.   0.   0.   0.   0.   0.  10.  10.  10.  10.  10.  10.   0.  10.
   10.  10.]
 [  0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.
    0.   0.   0.   0.   0.   0.  10.  10.  10.  10.  10.  10.  10.   0.
   10.  10.]
 [  0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   1.   0.
    0.   0.   0.   0.   0.   0.  10.  10.  10.  10.  10.  10.  10.  10.
    0.  10.]
 [  0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.
    0.   0.   0.   0.   0.   0.  10.  10.  10.  10.  10.  10.  10.  10.
   10.   0.]]










    
Out[15]:








[[0.8,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0],
 [0.022222222222222223,
  0.8,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0],
 [0.02197802197802198,
  0.02197802197802198,
  0.8,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0021978021978021982,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0],
 [0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.8,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0],
 [0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.8,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.0,
  0.0,
  0.0,
  0.0021978021978021982,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0],
 [0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.8,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0],
 [0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.8,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0021978021978021982,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0],
 [0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.8,
  0.022222222222222223,
  0.022222222222222223,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0],
 [0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.8,
  0.02197802197802198,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0021978021978021982,
  0.0,
  0.0,
  0.0],
 [0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.8,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0],
 [0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.8,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0],
 [0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.022222222222222223,
  0.8,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0],
 [0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.02197802197802198,
  0.02197802197802198,
  0.8,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0021978021978021982,
  0.0],
 [0.0,
  0.0,
  0.0,
  0.0,
  0.0021978021978021982,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.8,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0],
 [0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.8,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.0,
  0.0,
  0.0021978021978021982,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0],
 [0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.8,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0],
 [0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.8,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0],
 [0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.8,
  0.022222222222222223,
  0.022222222222222223,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0],
 [0.0,
  0.0,
  0.0021978021978021982,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.8,
  0.02197802197802198,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0],
 [0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.8,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0],
 [0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.8,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223],
 [0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.022222222222222223,
  0.8,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223],
 [0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.002173913043478261,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.002173913043478261,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.021739130434782608,
  0.021739130434782608,
  0.8,
  0.021739130434782608,
  0.021739130434782608,
  0.021739130434782608,
  0.021739130434782608,
  0.021739130434782608,
  0.021739130434782608,
  0.021739130434782608],
 [0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.8,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223],
 [0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.8,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223],
 [0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.8,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223],
 [0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0021978021978021982,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.8,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198],
 [0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.8,
  0.022222222222222223,
  0.022222222222222223],
 [0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0021978021978021982,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.02197802197802198,
  0.8,
  0.02197802197802198],
 [0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.0,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.022222222222222223,
  0.8]]

















In [ ]:

Content source: mmadsen/experiment-seriationct

Similar notebooks:

notebook.community | gallery | about