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

2 votes301 views1 answer
Python buffer
Python
6 votes320 views3 answers
Python how to draw a square
Python
2 votes304 views3 answers
Sum elements array with step numy
Python
1 vote291 views1 answer
Extract image from pdf python
Python
0 vote326 views1 answer
Select text in a div selenium python
Python
0 vote288 views1 answer
Read://https_www.tumblr.com/?url=https://www.tumblr.com/login?redirect_to=%2Fneue_web%2Fiframe%2Freblog%2F629907744681590784%2FjQw7OUs8&3739a18c-0c68-43cc-a4cb-b8b99e9bfd72=a52e06db-92b6-4b86-b3c5-fa2ab267405c
Python
2 votes325 views1 answer
Python IRR calculation
Python
4 votes267 views1 answer
Telegram bot send file
Python
0 vote318 views1 answer
Pandas not is na
Python
1 vote306 views1 answer
Evaluate how much a python program memory
Python