我有一个数据框,例如下面名为dataset的数据框,其中一列代表一个因子,另外两个代表值。我还计算了每个因子级别的一个列的平均值。
dataset <- data.frame(year=rep(1971:1977,3),
value=rnorm(21,160,6),
region=as.factor(c(rep("A",7),rep("B",7),rep("C",7))))
meancomp <-by(dataset$value,dataset$region,mean)
meanset <- data.frame(region=as.factor(c("A","B","C")),
meanvalue=c(meancomp[["A"]],meancomp[["B"]],meancomp[["C"]]))
xyplot(value~year|region,
data=dataset,
panel=function(...) ({panel.xyplot(...}))
我想在每个面板中添加一个图例,其中包含为该组计算的平均值。我想从指定的面板函数内部执行此操作(在示例代码中,面板只调用默认面板)。我当然可以将数据帧均值传递给面板,但是,我如何知道每个面板中绘制的因子的值? (如果我有这个值,我可以从meanvale获得平均值)。是否有另一种更简单的方法来实现相同的结果? 提前谢谢。