在具有多个面板的晶格xyplot图中创建标准误差线

Question

我有一个数据集，其中我绘制了2个地点随时间推移的4种治疗方法及其标准误的平均值。标准错误条没有正确分配给它们各自的均值-它们将同时显示在两个面板上-请您指教吗？参见示例：

d <- data.frame(site=rep(1:2,each=12),time=rep(1:3,8),trt=rep(rep(1:4,each=3),2)) 
d$mn <- rnorm(24,4,1)
d$se <- rnorm(24,2,1)
d$ul <- d$mn+d$se  # create y value for standard error upper limit

my.panel <- function(x,y, ...){
    panel.xyplot(x, y, ...)
    panel.arrows(x, y, x, d$ul, length = 0.1,
      angle = 90)
      }

xyplot(mn ~ time|site,data=d,
    group = trt,
    type=c('p','l'),
    panel = my.panel
)