As alternative visualisation to the level plot, the following figures are bubble plots of the same dataset. The figures seem to show better visualisation of the element concentration in the groundwater and river water.

1. load library and data

library("lattice")    # for plotting
library("dplyr")      # for data manipulation
library("gridExtra")  # for multiple plot
library("sp")         # for spatial analysis
data <- read.csv("BandungData.csv", header=TRUE)
coordinates(data) <- ~ x + y # assigning x and y columns as coordinates

2. Making bubble plot objects

We used the following code to build bubble plot object.

bubbleNO3 <- bubble(data, zcol="NO3", 
          xlab="X coord", ylab="Y coord",
          main="Bubble plot NO3",
          scales=list(tck=0.5))

3. Plotting

grid.arrange(bubbleNO3, bubbleNO2, 
             bubbleSO4, bubbleCl, 
             ncol=2)

plot of chunk Fig_plotriversig

grid.arrange(bubbleSiO2, bubbleFe, 
             bubbleMg, bubbleMn, 
             ncol=2)

plot of chunk Fig_plotriversig

4. Remarks