How to use GCDkit in a batch mode II

In one of the earlier blogs, Jeff has presented my exercise showing effective use of batch mode borrowed from our joint workshop at University of Catania.  

Here I would like to carry on and show the new capabilities of GCDkit to use multiple datasets and to overplot them on a single diagram (or plate thereof).

Binary plot

# Loading and storing the 1st dataset

data(atacazo)

accessVar("atacazo")

assign1col("gray")

pokeDataset("atacazo",overwrite.warn=FALSE)

# Loading and storing the 2nd dataset

data(sazava)

accessVar("sazava")

# Switching back to the atacazo dataset

peekDataset("atacazo")

# Overplotting classification diagram

plotDiagram("TAS",F,F)

overplotDataset("sazava",pch="*",col="red",cex=2)

Spiderplot

# Overplotting on spiderplot

peekDataset("atacazo")

# Atacazo as grey field

spider(WR,"Boynton",0.1,1000,field=TRUE,fill.col=TRUE,

  shaded.col="khaki")

overplotDataset("sazava")

A plate (Harker plots)

# blatna dataset

data(blatna)

accessVar("blatna")

multiple("SiO2","MgO,CaO,Na2O,K2O")

plateXLim(c(53,70)) # Common x axis scaling

plate0YLim() # All y axes start from 0

plateCex(1.5)

plateCexLab(1.5)

plateCol("gray")

# Overplot sazava on top of it

overplotDataset("sazava",cex=2)

# Getting rid of the datasets other than the current (blatna)

purgeDatasets()