First hitting time python

from numpy import *

hit_idx = (0,4)

# Define a graph by edge list
edges = [[0,1],[1,2],[2,3],[2,4]]

# Create adj. matrix
A = zeros((5,5))
A[zip(*edges)] = 1
# Undirected condition
A += A.T

# Make the final state an absorbing condition
A[hit_idx[1],:] = 0
A[hit_idx[1],hit_idx[1]] = 1

# Make a proper Markov matrix by row normalizing
A = (A.T/A.sum(axis=1)).T

B = A.copy()
Z = []
for n in xrange(100):
    Z.append( B[hit_idx] )
    B = dot(B,A)

from pylab import *
plot(Z)
xlabel("steps")
ylabel("hit probability")
show()

Source: stackoverflow.com

Similar codes

0 vote314 views1 answer
Xpath start-with python
Python
13 votes321 views2 answers
Python remove directory not empty
Python
0 vote273 views1 answer
Get method to create a set of counters in python
Python
0 vote351 views10 answers
EMAIL_BACKEND = 'django.core.mail.backends.console.EmailBackend'
Python Django
1 vote335 views1 answer
Django don't redirect on submission
Python Django
0 vote320 views1 answer
Python pandas series to title case
Python
0 vote311 views1 answer
Opencv waitkey example
Python
110 votes338 views6 answers
Python while
Python
0 vote318 views1 answer
Extract all bounding boxes using OpenCV Python
Python
2 votes322 views2 answers
Correlation plot python seaborn
Python