Nothing
## ---- include = FALSE---------------------------------------------------------
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>"
)
## ----setup, message=FALSE-----------------------------------------------------
library(Rwclust)
library(igraph)
## ----load_example_data--------------------------------------------------------
data(example1, package="Rwclust")
head(example1)
## ----plot_example_data, fig.align='center', fig.width=7, fig.height=5---------
labels <- c(1,1,1,1,2,2,2,2,3,3,3,3,3,3,3,3,3,4,4,4,4)
G <- igraph::graph_from_edgelist(as.matrix(example1[, c(1,2)]), directed=FALSE)
G <- igraph::set_edge_attr(G, "weight", value=example1$weight)
plot(G, edge.label=E(G)$weight, vertex.color=labels, layout=layout_with_fr)
## ----sharpen_edge_weights-----------------------------------------------------
result <- rwclust(example1, iter=6, k=3)
## ----plot_edge_weights, fig.align='center', fig.width=7, fig.height=5---------
G_sharpened <- igraph::graph_from_edgelist(as.matrix(example1[, c(1,2)]), directed=FALSE)
E(G_sharpened)$weights <- round(weights(result),0)
plot(G_sharpened, edge.label=E(G_sharpened)$weights, vertex.color=labels, layout=layout_with_fr)
## ----plot_histogram, fig.align='center', fig.width=7, fig.height=5------------
plot(result, cutoff=25, breaks=20)
## ----compute_components, fig.align='center', fig.width=7, fig.height=5--------
# delete edges with weights below the threshold
edges_to_keep <- which(weights(result) > 25)
example1_c <- example1[edges_to_keep, ]
example1_c$weight <- weights(result)[edges_to_keep]
G_c <- igraph::graph_from_edgelist(as.matrix(example1_c[, c(1,2)]), directed=FALSE)
# compute the connected components
clusters <- igraph::components(G_c)$membership
plot(G_c, vertex.color=clusters)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.