格子中的“ coplot”替代品

Question

我正在嘗試復制S中的統計模型中給出的示例。 這是書中給出的代碼。

data(ethanol)
attach(ethanol)
E.intervals <- co.intervals(E, number = 9, overlap = 1/4)
coplot(NOx ~ C | E, given.values = E.intervals,
                    panel = function(x,y) panel.smooth(x,y,degree = 1, span = 1))

該代碼僅在S ，因為panel.smooth函數沒有degree參數。 所以我用

coplot(NOx ~ C | E, given.values = E.intervals,
                    panel = function(x,y,...) panel.smooth(x,y, iter = 1, span = 1))

我想使用lattice圖形包，因為它具有與S panel.smooth類似的panel.loess函數，但是我沒有找到與coplot類似的coplot 。 我認為這是xyplot 。 如果是的話，您可以建議我給定的參數等效給given.values嗎？

Answer 1

這可能是完成的方式，雖然不盡相同，但非常接近。

require(gridExtra)
require(lattice)
data(ethanol)
E.intervals <- equal.count(E, number = 9, overlap = 1/4)
plot1 <- plot(E.intervals)
plot2 <- xyplot(NOx ~ C | E.intervals, data = ethanol, span = 1, layout = c(3,3),
                panel = function(x,y,span){
                panel.xyplot(x, y);
                panel.loess(x, y, span);
               })
grid.arrange(plot1,plot2, nrow=2)

差異是上部圖的高度，應該小些。