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

1 vote209 views1 answer
Show all numpy array
Python
0 vote254 views1 answer
Matplotlib remove ticks and lines
Python
1 vote211 views2 answers
Python leetcode
Python
0 vote198 views1 answer
Tensorboard 2.1.0 has requirement grpcio>=1.24.3, but you'll have grpcio 1.15.0 which is incompatible
Python
0 vote219 views2 answers
Python error: command 'x86_64-linux-gnu-gcc' failed with exit status 1
Python
0 vote202 views1 answer
Python long numbers as string
Python
0 vote216 views1 answer
Lenet 5 keras
Python
29 votes213 views5 answers
Python date and time
Python
1 vote248 views2 answers
Json.stringify equivalent in python
Python
0 vote202 views1 answer
How remove spaces python input
Python