Mixing netplot with base graphics

The hardware and bandwidth for this mirror is donated by dogado GmbH, the Webhosting and Full Service-Cloud Provider. Check out our Wordpress Tutorial.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]dogado.de.

George G. Vega Yon

2023-05-22

Working with grid graphics provides flexibility that is harder to get when using base graphics. Nevertheless, a lot of users may still be more comfortable with base graphics, so this document shows a few examples of how to mix both of them using the gridBase package.

First, we load igraph to simulate a small network and complement that with uninteresting correlated data:

library(igraph)

# Simulation parameters
set.seed(12)
n <- 100

# Points
x0  <- rnorm(n)
x1  <- x0*.4 + rnorm(n)

# network
net <- sample_smallworld(1, n, 4, .1)

And here, we can create the netplot object. We use the magrittr package to make the workflow nicer:

library(netplot)

library(magrittr)

np <- nplot(net, bg.col = "gray70", vertex.nsides = "square") %>%
  set_edge_gpar("line", col = "gray50", alpha = .7) %>%
  set_vertex_gpar("frame", fill = "black", col = "black")

Imagine that we would like to include a plot of the network in the following plot

plot(x0, x1)

We will use the gridBase package. gridBase allows some interaction between the base and grid graphics systems. In particular, we will use the baseViewports() function, which returns a list with viewports that can be mapped to match the current base plot

library(gridBase)
library(grid)

# Getting the current viewports
vp <- baseViewports()

Once we have the viewports, we can work with these to add the current plot to the base plot. In this case, we will include the network that we just created, np, and this will be placed on the top left corner of the plot region and will occupy 1/4 of the area:

# Adding the plot to the current view
pushViewport(vp$plot)

# Sub viewport, top-left corner using 1/4 of the area
pushViewport(viewport(.25, .75, width = .5, height = .5)) 

# Drawing, we could have used `print(np, newpage=FALSE)`
grid.draw(np)

Or suppose that you want to frame it into a rectangle to the top-right corner of the entire device that uses 3/10 x 3/10 of the area, then:

plot(x0, x1)
# vp <- baseViewports() 
# pushViewport(vp$figure)
pushViewport(viewport(1, 1, just = c("right", "top"), width = .3, height = .3))
grid.draw(np)

Notice that we didn’t need to extract the viewports of the base plot. A simple viewport did the job since we are using the entire device. Now, imagine that you are feeling playful and would like to draw it a bit tilted, say -20 degrees:

plot(x0, x1)
pushViewport(
  viewport(x = .8, y = .8, width = .3, height = .3, angle = -20)
  )
grid.draw(np)

These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.
Health stats visible at Monitor.