4

我正在尝试使用库lattice和latticeExtra. R但是，没有正确表示垂直 y 轴上的值。不是选择数值变量的实际值，而是R绘制值的等级。也就是说，有值[375, 500, 625, 750, ..., 3000]并R绘制它们的等级[1,2,3,4,...23]并相应地选择比例。有人遇到过这样的问题吗？如何使用(0, 500, 1000, 1500, ...)垂直 y 刻度上的刻度来管理获得正确的表示？

到目前为止的程序代码如下：

df.dose <- read.table("data.csv", sep=",", header=TRUE)
library(lattice); library(latticeExtra)

useOuterStrips(dotplot(z ~ sample.size | as.factor(effect.size)*as.factor(true.dose),
               groups=as.factor(type), data=df.dose, as.table=TRUE))

（从下面的评论中添加）：另外，可以将误差线添加到图表中吗？我想到了以下（要添加到调用中），但它似乎不起作用。有可能吗？

up=z+se, lo=z-se, panel.groups=function(x,y,..., up, lo, subscripts){ 
   up <- up[subscripts]
   lo <- lo[subscripts]
   panel.segments(lo, as.numeric(y), up, as.numeric(y), ...)
}

这是我的数据：https ://www.dropbox.com/s/egy25cj00rhum40/data.csv

expand.grid补充：这是使用和的数据的相关部分dput：

df.dose <- expand.grid(effect.size=c(-.5, -.625, -0.75),
                       sample.size=c(40L, 60L, 80L),
                       true.dose=c(375L, 500L, 750L, 1125L),
                       type=c("dose", "categ", "FP2", "FP1"))
df.dose$z <- c(875L, 875L, 750L, 750L, 750L, 625L, 625L, 625L, 625L, 875L, 
875L, 750L, 1000L, 1000L, 1000L, 1125L, 1000L, 875L, 1000L, 1000L, 
875L, 1000L, 1000L, 875L, 1125L, 1000L, 1000L, 1250L, 1125L, 
1000L, 1250L, 1250L, 1125L, 1250L, 1000L, 1000L, 500L, 500L, 
500L, 500L, 500L, 500L, 500L, 500L, 500L, 625L, 625L, 625L, 625L, 
625L, 625L, 625L, 625L, 625L, 750L, 750L, 625L, 750L, 750L, 750L, 
750L, 750L, 750L, 875L, 875L, 750L, 750L, 875L, 875L, 875L, 875L, 
875L, 2500L, 1500L, 1125L, 2000L, 1000L, 1750L, 250L, 500L, 500L, 
1250L, 750L, 625L, 875L, 500L, 500L, 875L, 500L, 375L, 1250L, 
875L, 750L, 1000L, 625L, 625L, 875L, 500L, 500L, 1125L, 1000L, 
875L, 1125L, 875L, 625L, 1125L, 1000L, 625L, 2500L, 2125L, 2375L, 
2000L, 750L, 2625L, 250L, 625L, 250L, 875L, 875L, 500L, 625L, 
500L, 625L, 1000L, 500L, 375L, 1000L, 875L, 625L, 875L, 500L, 
500L, 875L, 500L, 500L, 1250L, 1125L, 875L, 1125L, 875L, 750L, 
1250L, 1000L, 625L)
4

3 回答 3

7

你需要做z一个因素：dotplot(factor(z) ~ ...

此外，您可能希望情节中有一些抖动以防止重叠；尝试添加jitter.x=TRUEor jitter.y=TRUE，或两者兼而有之。

从您在下面的评论和再次查看数据来看，我认为您以错误的方式绘制点图。我认为您希望这些线条适用于样本大小，而不是适用于z's. 如果您真的想要z在垂直轴上，则需要添加horizontal=TRUE. 您还可以交换水平轴和垂直轴上的内容。

useOuterStrips(dotplot(z ~ factor(sample.size) | 
                             as.factor(effect.size)*as.factor(true.dose),
                  groups=as.factor(type), data=df.dose,  
                  as.table=TRUE, horizontal=FALSE, jitter.x=TRUE))

要添加错误栏，它有点复杂，因为您在面板中有组，因此您需要使用一个panel.groups函数；此外，为了使线条不重叠，您可能希望将它们从一侧抖动到另一侧，这最好在自定义panel函数中完成。

df.dose$se <- 200
df.dose$type <- factor(df.dose$type)
df.dose$sample.size <- factor(df.dose$sample.size)

panel.groups.mydotplot <- function(x, y, subscripts, up, lo, 
                                   col=NA, col.line=NA, ...) {
  panel.points(x, y, ...)
  panel.segments(x, lo[subscripts], x, up[subscripts], col=col.line, ...)
}
panel.mydotplot <- function(x, y, subscripts, groups, ..., jitter=0.1) {
  jitter <- seq(-1,1,len=nlevels(groups))*jitter
  xx <- as.numeric(x) + jitter[as.numeric(groups[subscripts])]
  panel.dotplot(x, y, groups=groups, subscripts=subscripts, pch=NA, ...)
  panel.superpose(xx, y, groups=groups, subscripts=subscripts,  
                  panel.groups=panel.groups.mydotplot, ...)
}
pp <- dotplot(z ~ sample.size | as.factor(effect.size)*as.factor(true.dose),
              groups=type, data=df.dose, as.table=TRUE, horizontal=FALSE,
              up=df.dose$z + df.dose$se, lo=df.dose$z - df.dose$se,
              panel=panel.mydotplot, auto.key=list(space="right"))
useOuterStrips(pp)

在此处输入图像描述

于 2012-10-11T14:11:49.490 回答
3

我不确定我是否理解问题并且您要求提供格子解决方案，但我认为使用 ggplot2 完成此操作可能会有所帮助：

ggplot(data=df.dose, aes(x=sample.size, y=as.factor(z), colour=type)) +
    geom_point() + facet_grid(true.dose~effect.size)

产量： 在此处输入图像描述

或者我们可以通过以下方式释放天平：

ggplot(data=df.dose, aes(x=sample.size, y=as.factor(z), colour=type)) +
    geom_point() + facet_grid(true.dose~effect.size, scales="free")

产量：

在此处输入图像描述

于 2012-10-11T14:51:12.150 回答
3

您也可以使用xYplot包Hmisc来实现类似于@Aaron 的解决方案，尽管要获得与他相同的抖动可能有点棘手：

a <- xYplot(Cbind(z, z-se, z+se) ~ sample.size | as.factor(effect.size) * as.factor(true.dose),
            groups=as.factor(type), data=df.dose, as.table=TRUE, auto.key=list(space="top"))
useOuterStrips(a)

在此处输入图像描述

但真的是翔实的情节吗？它是否很好地显示了您的数据效果，突出了您的比较？它是否探索数据中的任何趋势？为了更好地查看您想要绘制的所有因素，我将首先在您的组之间建立线条连接，以便更好地查看不同sample.size.

key.variety <- list(space = "top", 
                    text = list(levels(df.dose$type)),
                    points = list(pch = 0:3, col = "black"))
a <- xyplot(z ~ as.factor(sample.size) | as.factor(effect.size)*as.factor(true.dose),
            df.dose, type = "o", as.table=TRUE, groups = type, key = key.variety, 
            lty = 1, pch = 0:3, col.line = "darkgrey", col.symbol = "black")
useOuterStrips(a)

在此处输入图像描述

但是那里隐藏着一些东西，由于数据的密度，仍然存在太多的噪音。让我们摆脱effect.size并绘制回归线，尽管使用这么少的数据点可能是一种罪过。

a <- xyplot(z ~ as.factor(sample.size) | as.factor(type)*as.factor(true.dose), 
            data=df.dose, as.table=TRUE, 
            panel = function(x, y){
               panel.xyplot(x, y, jitter.x = T, col=1);
               panel.lmline(x, y, col=1, lwd=1.5);
            })
useOuterStrips(a)

在此处输入图像描述

我知道我可能没有说服您，但有时最好从太多因素中卸载图表以更好地查看数据。有时，如果您将因素分开显示，它可能更容易在视觉上访问。

于 2012-10-16T19:59:17.683 回答