函數繪製對（）變量對對角線的公共垂直軸變量

Question

我試圖定義一個函數來產生一對圖的對角線上的散點圖。訣竅是我想要在panel.splot函數中指定垂直軸。最終我認爲這裏使用的是在對角線上可視化時間序列。對於一個簡單的例子，雖然我想用虹膜數據集來嘗試，並使用Species作爲垂直軸變量。下面的陰謀產生一個漂亮的情節：

pairs(iris[,c(1:4)]) 

所以現在我需要定義一個對角線面板的功能。這裏是我的努力迄今：

panel.splot <- function(x, yvar,...) 
{ 
    usr <- par("usr"); on.exit(par(usr)) 
    par(usr = c(0, 2, usr[3:4])) 
    plot(yvar,x, xlim=c(min(x),max(x))) 
} 

然而，當我嘗試運行它，我得到的是我不知道如何解釋錯誤消息。

pairs(iris[,c(1:4)],diag.panel=panel.splot(x=x, yvar="Species")) 

Error in plot.window(...) : invalid 'xlim' value 
In addition: Warning messages: 
1: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by coercion 
2: In xy.coords(x, y, xlabel, ylabel, log) : NAs introduced by coercion 
3: In min(x) : no non-missing arguments to min; returning Inf 
4: In max(x) : no non-missing arguments to max; returning -Inf 

我一直沒能找到另外一個例子。許多其他的功能來創建不同類型的情節，但沒有任何事情做到這一點。

爲了清楚起見，這是情節我想象會沿着對角線pairs()通話產生的類型：

plot(iris$Species,iris$Petal.Width) 

Answer 1

儘管你的例子情節，我解釋這個稍有不同（錯誤地）由於使用物種作爲縱軸變量。休伯特給出了答案，但認爲它可能會增加一些東西。

使用pairs()

panel.splot <- function(Y, Adjust=0.2){ 
        function(x, y=Y, adjust=Adjust) 
        { 
        usr <- par("usr"); on.exit(par(usr)) 
        xs <- range(x) 
        ys <- if(is.factor(y)) c(1, nlevels(y)) else range(y) 
        par(usr = c(xs[1], xs[2], ys[1], ys[2]) + c(-adjust, adjust)) 
        points(x, y) 
        } 
       } 

pairs(iris[, 1:4], diag.panel=panel.splot(Y=iris$Species)) 

將會產生

要做到同樣的事情ggpairs()你也可以定義用戶功能

library(GGally) 
library(ggplot2)  

# Define diagonal function 
diag_fun <- function(data, mapping, ...){ 
     ggplot(data = data, mapping = mapping) + 
     geom_point(...) 
     } 

ggpairs(iris, columns=1:4, 
     diag = list(continuous = diag_fun, mapping=aes(y=Species))) 

其中給出

---

更新您的意見...

是，對角項，該變量爲x映射通過。您可以使用coord_flip()或反轉函數中的映射。下面的函數執行此操作，並使用annotate()（ps。）添加變量名稱。你應該能夠換出_point GEOM的發言權，_boxplot

diag_fun <- function(data, mapping, xnudge=0.95, ynudge=1, pts=lst(), ...){ 

    # create range for annotate 
    lbl <- as.character(mapping$x) 
    xmax <- max(data[, as.character(mapping$x)], na.rm=TRUE) 
    ymax <- mean(seq(nlevels(data[,as.character(mapping$y)]))) 

    # reverse mapping so no need for coord_flip() 
    tmp <- mapping$y 
    mapping$y <- mapping$x 
    mapping$x <- tmp 

    ggplot(data=data, mapping=mapping) + 
     do.call(geom_point, pts) + 
     annotate("text", x=ynudge*ymax, y=xnudge*xmax, label=lbl, ...) 
     } 

所以在ggpairs與默認

ggpairs(iris, columns=1:4, diag = list(continuous = diag_fun, mapping=aes(y=Species))) 

使用它，你可以使用wrap通過進一步論證

ggpairs(iris, columns=1:4, 
     diag = list(continuous = 
         wrap(diag_fun, 
          size=10, col="red", # make changes to the text label 
          pts=list(colour="blue")), # make changes to geom_point 
        mapping=aes(y=Species))) 

Answer 2

繼RAWR建議，我設法讓你的情節：

panel.splot <- function(y) 
{ 
     function(x,yv=y) 
     { 
       usr <- par("usr"); on.exit(par(usr)) 
       par(usr = c(0, 2, usr[3:4]), new = TRUE) 
       plot(yv,x) 
     } 
} 
pairs(iris[,1:4],diag.panel=panel.splot(iris[,5])) 

你仍然必須處理的文本大小