Description Usage Arguments Details Author(s) References See Also Examples
Creates a plot for each element of two lists or each column of each dataset against each of the second.
1 2 3 4 5 6  pairwisePlot(X,Y,...)
## Default S3 method:
pairwisePlot(X,Y=X,...,
xlab=deparse(substitute(X)),ylab=deparse(substitute(Y)),
nm=c(length(Y),length(X)),panel=plot,
add.line=FALSE, line.col=2,add.robust=FALSE,rob.col=4)

X 
a list, a data.frame, or a matrix representing the first set of things to be displayed. 
Y 
a list, a data.frame, or a matrix representing the second set of things to be displayed. 
... 
furter parameters to the panel function 
xlab 
The sequence of labels for the elements of X. Alternatively the labels can be given as colnames or names of X. This option takes precedence if specified. 
ylab 
The sequence of labels for the elements of Y. Alternatively the labels can be given as colnames or names of Y. This option takes precedence if specified. 
nm 
the parameter to be used in the call

panel 
The panel function to plot the individual panels.
If the panel function admits a formula interface, it is called
as 
add.line 
logical, to control the addition of a regression line in each panel 
line.col 
in case the regression line is added, which color should be used? defaults to red. 
add.robust 
logical, to control the addition of a robust regression line
in each panel. Ignored if covariable is a factor. This is nowadays
based on 
rob.col 
in case the robust regression line is added, which color should be used? Defaults to blue. 
This is a lightweight convenience function to plot several
aspects of one dataset against several aspects of another dataset. It
is far more straightforward than e.g. the pairs
function and
does not do any internal computation rather than organizing the names.
Of course, the rows of the two data sets must be the same.
The current implementation may display a warning about the function
panel
dispatching methods for generic plot
. It can be
ignored without harm.
Optionally, classical and/or robust regression lines can be drawn, though only for nonfactor covariables.
It may be convenient to use par
capabilities to fit the device
characteristics to the plot, in particular arguments mar
and oma
.
K.Gerald v.d. Boogaart http://www.stat.boogaart.de, Raimon TolosanaDelgado
Boogaart, K.G. v.d. , R. Tolosana (2008) Mixing Compositions and Other scales, Proceedings of CodaWork 08.
http://ima.udg.es/Activitats/CoDaWork03/
http://ima.udg.es/Activitats/CoDaWork05/
http://ima.udg.es/Activitats/CoDaWork08/
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39  X < rnorm(100)
Y < rnorm.acomp(100,acomp(c(A=1,B=1,C=1)),0.1*diag(3))+acomp(t(outer(c(0.2,0.3,0.4),X,"^")))
pairs(cbind(ilr(Y),X),panel=function(x,y,...) {points(x,y,...);abline(lm(y~x))})
pairs(cbind(balance(Y,~A/B/C),X),
panel=function(x,y,...) {points(x,y,...);abline(lm(y~x))})
pairwisePlot(balance(Y,~A/B/C),X)
pairwisePlot(X,balance(Y,~A/B/C),
panel=function(x,y,...) {plot(x,y,...);abline(lm(y~x))})
pairwisePlot(X,balance01(Y,~A/B/C))
# A function to extract a portion representation of subcompsitions
# with two elements:
subComps < function(X,...,all=list(...)) {
X < oneOrDataset(X)
nams < sapply(all,function(x) paste(x[[2]],x[[3]],sep=","))
val < sapply(all,function(x){
a = X[,match(as.character(x[[2]]),colnames(X)) ]
b = X[,match(as.character(x[[2]]),colnames(X)) ]
c = X[,match(as.character(x[[3]]),colnames(X)) ]
return(a/(b+c))
})
colnames(val)<nams
val
}
pairwisePlot(X,subComps(Y,A~B,A~C,B~C))
## using Hydrochemical data set as illustration of mixed possibilities
data(Hydrochem)
xc = acomp(Hydrochem[,c("Ca","Mg","Na","K")])
fk = Hydrochem$River
pH = log10(Hydrochem$H)
covars = data.frame(pH, River=fk)
pairwisePlot(clr(xc), pH)
pairwisePlot(clr(xc), pH, col=fk)
pairwisePlot(pH, ilr(xc), add.line=TRUE)
pairwisePlot(covars, ilr(xc), add.line=TRUE, line.col="magenta")
pairwisePlot(clr(xc), covars, add.robust=TRUE)

Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.