showClass("piper")
l <- list( Ca = c(43,10,73,26,32),
Mg = c(30,50,3,14,12),
Cl = c(24,10,12,30,43),
SO4 = c(24,10,12,30,43),
WaterType = c(2,2,1,2,3),
IDs = c("A","B","C","D","E") )
lp <- piper(l)
plot( lp, main="Piper-Hill Diagram of Water Quality" )
Run the code above in your browser using DataLab