R 晶格中跨面板的条件层

我想在一个带有两个面板的格子

stripplot

中绘制两种不同条件下的单个主题的意思。我还想添加我已计算并存储在单独数据框中的受试者置信区间。我试图用latticeExtra的

层

函数覆盖这些置信区间。当我添加层时，两组间隔显示在两个面板上（如下面的代码和第一幅图像所示），或者如果我在

层

命令中将

[下标]

添加到x和y（如下面的第二个代码剪辑和图像所示），两组间隔仅显示在第一个面板上。如何在适当的面板上显示适当的间隔

library(latticeExtra)

raw_data <- data.frame(subject = rep(1:6, 4), cond1 = as.factor(rep(1:2, each = 12)), cond2 = rep(rep(c("A", "B"), each = 6), 2), response = c(2:7, 6:11, 3:8, 7:12))
summary_data <- data.frame(cond1 = as.factor(rep(1:2, each = 2)), cond2 = rep(c("A", "B"), times = 2), mean = aggregate(response ~ cond2 * cond1, raw_data, mean)$response, within_ci = c(0.57, 0.54, 0.6, 0.63))
summary_data$lci <- summary_data$mean - summary_data$within_ci
summary_data$uci <- summary_data$mean + summary_data$within_ci

subject_stripplot <- stripplot(response ~ cond1 | cond2, groups = subject, data = raw_data, 
  panel = function(x, y, ...) {
    panel.stripplot(x, y, type = "b", lty = 2, ...)
    panel.average(x, y, fun = mean, lwd = 2, col = "black", ...)    # plot line connecting means
  }
)
addWithinCI <- layer(panel.segments(x0 = cond1, y0 = lci, x1 = cond1, y1 = uci, subscripts = TRUE), data = summary_data, under = FALSE)
plot(subject_stripplot + addWithinCI)

库（latticeExtra）
原始数据一种可能的解决方案是
打印条带图（例如，在
png
或任何其他图形设备内），然后使用
网格焦点修改每个子面板

## display stripplot
print(subject_stripplot)

## loop over grops
for (i in c("A", "B")) {

  # subset of current group
  dat <- subset(summary_data, cond2 == i)

  # add intervals to current panel
  trellis.focus(name = "panel", column = ifelse(i == "A", 1, 2), row = 1)
  panel.segments(x0 = dat$cond1, y0 = dat$lci, 
                 x1 = dat$cond1, y1 = dat$uci, subscripts = TRUE)
  trellis.unfocus()
}

一种可能的解决方案是
打印
条带图（例如，在
png
或任何其他图形设备内），然后使用
grillis.focus
修改每个子面板

## display stripplot
print(subject_stripplot)

## loop over grops
for (i in c("A", "B")) {

  # subset of current group
  dat <- subset(summary_data, cond2 == i)

  # add intervals to current panel
  trellis.focus(name = "panel", column = ifelse(i == "A", 1, 2), row = 1)
  panel.segments(x0 = dat$cond1, y0 = dat$lci, 
                 x1 = dat$cond1, y1 = dat$uci, subscripts = TRUE)
  trellis.unfocus()
}

另一个不需要latticeExtra的解决方案（来自Duncan Mackay）：

summary_data$cond3另一个不需要latticeExtra的解决方案（来自Duncan Mackay）：

summary_数据$cond3
p_seg <- xyplot(lci ~ cond1 | cond2, data = summary_data, ylim = c(1, 13),
       panel = function(...) {
         # lower and upper y values
         y0 <- list(summary_data$lci[c(1, 3)], summary_data$lci[c(2, 4)])
         y1 <- list(summary_data$uci[c(1, 3)], summary_data$uci[c(2, 4)])
         # insert vertical lines depending on current panel
         panel.segments(x0 = 1:2, x1 = 1:2,
                        y0 = y0[[panel.number()]], 
                        y1 = y1[[panel.number()]])
       })

p_comb <- subject_stripplot + 
  as.layer(p_seg)

# print(p_comb)

summary_data$cond3 <- sapply(summary_data$cond2, pmatch, LETTERS)

mypanel <- function(x, y, ..., lci, uci, scond1, scond3, groups, type, lty){
pnl = panel.number()
panel.xyplot(x, y, ..., groups = groups, type = type, lty = lty)
panel.average(x, y, horizontal = FALSE, col = "black", lwd = 3)
panel.segments(x0 = scond1[scond3 == pnl],
               y0 = lci[scond3 == pnl],
               x1 = scond1[scond3 == pnl],
               y1 = uci[scond3 == pnl])
}
with(summary_data,
 stripplot(response ~ cond1 | cond2, data = raw_data,
           groups = subject,
           lci = lci,
           uci = uci,
           scond1 = summary_data$cond1,
           scond3 = cond3,
           type = "b",
           lty = 2,
           panel = mypanel)
)