Closed #119 by adding node_quad_census()

DESCRIPTION

@@ -29,6 +29,7 @@ Imports:
     igraph,
     magrittr,
     network,
+    oaqc,
     purrr,
     RColorBrewer,
     readxl,

NAMESPACE

@@ -132,6 +132,7 @@ export(node_components)
 export(node_constraint)
 export(node_degree)
 export(node_eigenvector)
+export(node_quad_census)
 export(node_smallworld)
 export(node_tie_census)
 export(node_triad_census)
@@ -240,6 +241,7 @@ importFrom(network,as.matrix.network.adjacency)
 importFrom(network,as.matrix.network.incidence)
 importFrom(network,as.network)
 importFrom(network,is.network)
+importFrom(oaqc,oaqc)
 importFrom(purrr,map)
 importFrom(readxl,read_xlsx)
 importFrom(rlang,.data)

R/census.R

@@ -81,6 +81,32 @@ node_triad_census <- function(object){
   out # This line says the function returns the output
 }
 
+#' @rdname census
+#' @importFrom oaqc oaqc
+#' @examples 
+#' (quad_cen <- node_quad_census(southern_women))
+#' @export
+node_quad_census <- function(object){
+  graph <- object %>% as_tidygraph() %>% 
+    activate(edges) %>% as.data.frame()
+  out <- oaqc::oaqc(graph)[[1]]
+  out <- out[-1,]
+  rownames(out) <- node_names(object)
+  colnames(out) <- c("E4", # co-K4
+                     "I41","I40", # co-diamond
+                     "H4", # co-C4
+                     "L42","L41","L40", # co-paw
+                     "D42","D40", # co-claw
+                     "U42","U41", # P4
+                     "Y43","Y41", # claw
+                     "P43","P42","P41", # paw
+                     "04", # C4
+                     "Z42","Z43", # diamond
+                     "X4") # K4
+  if(is_twomode(object)) out <- out[,-c(8,9,14,15,16,18,19,20)]
+  out
+}
+
 #' @rdname census
 #' @examples 
 #' group_tie_census(task_eg, cutree(cluster_structural_equivalence(task_eg), 4))