[R] Error : unused arguments in pairs()

Sarah Goslee sarah.goslee at gmail.com
Tue Jan 6 15:55:39 CET 2009


You didn't do what Prof. Ripley suggested - adding a ... argument.
Here's a crude version of what you want; I'm sure there's a more elegant
solution for passing the needed data to the panel function.

panel.cor <- function(x, y, digits=2, prefix="", splitvar, col.cor, ...)

{
   usr <- par("usr"); on.exit(par(usr))
   par(usr = c(0, 1, 0, 1))
   r <- abs(cor(x, y))
   if(!missing(splitvar)) {
      r <- c(r, abs(sapply(lapply(split(cbind.data.frame(x, y),
splitvar), cor), function(x)x[1,2])))
   }
   txt <- format(c(r, 0.123456789), digits=digits)[1:4]
   txt <- paste(prefix, txt, sep="")
   if(missing(col.cor)) col.cor <- c("black", "red", "green3", "blue")
   for(i in 1:length(txt)) {
      text(0.5, (1/(length(txt)+1))*i, txt[i], col = col.cor[i])
   }
 }

 pairs(iris[1:4], main = "Anderson's Iris Data -- 3 species",
     pch = 21, bg = c("red", "green3", "blue")[unclass(iris$Species)],
lower.panel=panel.cor, splitvar=iris$Species)


On Mon, Jan 5, 2009 at 5:32 PM, herwig <bachmannherwig at hotmail.com> wrote:
>
> Dear Prof. Ripley,
>
> thanks for your reply.
> Unfortunately I still was not able to solve the problem.
> I tried it with the Iris data and you can find the code below:
>
>>  panel.cor <- function(x, y, digits=2, prefix="", cex.cor)
>
> {
>      usr <- par("usr"); on.exit(par(usr))
>    par(usr = c(0, 1, 0, 1))
>    r <- abs(cor(x, y))
>    txt <- format(c(r, 0.123456789), digits=digits)[1]
>    txt <- paste(prefix, txt, sep="")
>    if(missing(cex.cor)) cex.cor <- 0.8/strwidth(txt)
>    text(0.5, 0.5, txt, cex = cex.cor * r)
>  }
>
>  pairs(iris[1:4], main = "Anderson's Iris Data -- 3 species",
>      pch = 21, bg = c("red", "green3", "blue")[unclass(iris$Species)],
> lower.panel=panel.cor )
>
> What I would ideally like to have is that the upper panel shows the plot
> with the iris$Species color coded and the lower panel showing the
> correlation coefficients  for the total data in the plot, but also for the
> each iris$Species with the same color code as in the upper panels.
> So the lower panel should therefore actually show 4 correlation coefficients
> in each panel (color coded); It does not need to be sized according to the
> value of the correlation coefficient as shown in the example above.
> Any help with this would be appreciated a lot,
>
> Kind regards,
>
> Herwig
>
>
>
-- 
Sarah Goslee
http://www.functionaldiversity.org




More information about the R-help mailing list