Question

我对R很新，所以我希望你们可以帮我解决这个问题。我试图寻找有类似问题的其他帖子，但无法找到任何帖子。

这就是我的代表性数据集的构建方式：

library(lattice)
data(Theoph)

DF <- Theoph
DF$Time <- rep(c(0,0.25,0.5,1,2,4,5,7,9,12,24), times=12)
DF$group[DF$Subject==1|DF$Subject==2|DF$Subject==3] <- "A"
DF$group[DF$Subject==4|DF$Subject==5|DF$Subject==6] <- "B"
DF$group[DF$Subject==7|DF$Subject==8|DF$Subject==9] <- "C"
DF$group[DF$Subject==10|DF$Subject==11|DF$Subject==12] <- "D"

y <- DF$conc
x <- DF$Time
groups <- as.factor(DF$group)

我制作了一个xyplot，它显示了各个数据点以及使用此代码显示每组平均值（A，B，C和D）的行：

xyplot(y ~ x, type="p", groups=groups, auto.key=T,
panel = function(x, y, groups, subscripts, ...){
        panel.superpose(x, y, groups, subscripts,...,
        panel.groups=function(x,y,col,col.symbol,...){
            panel.xyplot(x, y, col=col.symbol,...)
            panel.average(x, y, col=col.symbol, lwd = 2, horizontal = FALSE)
        })
    })

我想在平均值中添加误差条（例如+/-标准偏差）。我知道有一个名为panel.arrows的函数，但是我没有成功地将它与panel.groups函数结合起来。

有没有办法做到这一点？或者我应该首先使用SD以及误差线的下限和上限制作单独的数据帧？

Answer 1

这可能不是你想要的，但它非常接近。在这里阅读： R - Lattice xyplot - How do you add error bars to groups and summary lines?

r-将错误条添加到xyplot中的不同组

1 个答案: