# Basic HITS Algorithm for the Stanford MMDS Course

HITS.py

# Basic HITS Algorithm for the Stanford MMDS Course
import numpy as np

# L = [[ 0.,1.,1.,1.,0.], [1.,0.,0.,1.,0.], [0.,0.,0.,0.,1.], [0.,1.,1.,0.,0.], [0.,0.,0.,0.,0.]]
L = [[0,1,1,0],[1,0,0,0],[0,0,0,1],[0,0,1,0]]

Lt = np.transpose(L)

h = np.matrix([1.,1.,1.,1.]).T

for i in range(1,3):
  a = np.dot(Lt,h)  
  a = np.dot(a,1/a.max())
  #print a
  h = np.dot(L,a)
  h = np.dot(h,1/h.max())
  #print h

print h #Hubiness