diff --git a/Schrick-Noah_CG-Analysis.R b/Schrick-Noah_CG-Analysis.R index c923326..c6dc7f3 100644 --- a/Schrick-Noah_CG-Analysis.R +++ b/Schrick-Noah_CG-Analysis.R @@ -4,6 +4,7 @@ # Noah L. Schrick - 1492657 library(igraph) +library(centiserve) ################## Read in the previously generated networks ################## setwd(dirname(rstudioapi::getActiveDocumentContext()$path)) @@ -21,3 +22,16 @@ car.katz <- katz.cent(car) hipaa.katz <- katz.cent(hipaa) pci.katz <- katz.cent(pci) +### Page Rank +car.pr <- page.rank(car) +hipaa.pr <- page.rank(hipaa) +pci.pr <- page.rank(pci) + +### K-path +car.kpe <- geokpath(car, V(car), "out") + + +############# Other- Tmp work +min_cut(car,"0", "2490") +min_cut(hipaa,"0","2320") +min_cut(pci,"0","60") diff --git a/centralities.R b/centralities.R index f81e6fc..067d14c 100644 --- a/centralities.R +++ b/centralities.R @@ -7,6 +7,9 @@ katz.cent <- function(A, alpha=NULL, beta=NULL){ #NULL sets the default value lam.dom <- eigen(A)$values[1] #dom eigenvec if (is.null(alpha)){ alpha <- 0.9 * (1/lam.dom) #Set alpha to 90% of max allowed + if (is.complex(alpha)){ + alpha <- Re(alpha) + } } n <- nrow(A) @@ -18,4 +21,4 @@ katz.cent <- function(A, alpha=NULL, beta=NULL){ #NULL sets the default value scores <- solve(diag(n) - alpha*A,beta) return(scores) -} \ No newline at end of file +}