new power function

gistfile1.txt

power <- function(B1 = B1, B2 = B2, v.old = v, tol = 10^-10, p = .15) {
  N <- length(v)
  B.p <- matrix(0, ncol=5242, nrow=5242)
  B.p <- sapply(1:5242, function(x) (1-p)/sum(B1[x,]) * B1[x,])

  v.old <- v
  
  v.new <- B.p %*% v.old + (p/N) * B2 %*% v.old
  v.new <- v.new / norm(v.new)
  i <- 0
  
  while( norm(v.new - v.old) > tol) {
    v.old <- v.new
    v.new <- B.p %*% v.old + (p/N) * B2 %*% v.old
    v.new <- v.new / norm(v.new)
    i <- i + 1
  }
  return(list(vector = v.new, 
              value = t(v.new) %*% (B.p + (p/N) * B2) %*% v.new, 
              topnode = which.max(v.new),
              iter = i
              ))
}

##results
$vector
huge vector of nonsense

$value
     [,1]
[1,]    1

$topnode
[1] 2775

$iter
[1] 132