R中的刻面格子图，例如线框：如何删除条带并添加1行字幕

Question

我使用这种功能（来自虹膜数据集）

model_test <- lm(Sepal.Length ~( Petal.Length + Sepal.Width + Petal.Width +Species)^2,
data=iris)
gg<-expand.grid(Petal.Length=0:6,Species=levels(Species))
vv<-expand.grid(Sepal.Width=0:4,Petal.Width=1:4)

pd<-do.call(rbind,Map(function(Petal.Length,Species,Sepal.Width,Petal.Width){
              nd <- cbind(vv, Petal.Length=Petal.Length,Species=Species,
                               Sepal.Width=Sepal.Width, Petal.Width=Petal.Width)
              cbind(nd, pred=predict(model_test, nd, type="response"))},
              Petal.Length=iris$Petal.Length,Species=iris$Species,
               Sepal.Width=iris$Sepal.Width,Petal.Width=iris$Petal.Width))

wireframe(pred~Sepal.Width+Petal.Width|Species*Petal.Length,
           pd, drape=FALSE,scale=list(arrows=FALSE),subset=(Species=="setosa"),
           layout = c(3, 3))

我得到这个情节

我的问题： 如果我在模型中包含更多因素，（彩色）条带会使图形非常小。如何删除绿色和浅红色条带并添加1行标题 例如Species==setosa & Petal.Length==0，....，Species==setosa & Petal.Length==6。

Answer 1

莱迪思将为每个条件变量创建一个条带。如果您只想要一个条带，请尝试调整交互。例如

wireframe(pred~Sepal.Width+Petal.Width|interaction(Species,Petal.Length),
           pd, drape=FALSE,scale=list(arrows=FALSE), subset=(Species=="setosa"), layout=c(3,3))

给出了

您可以提前创建交互，并根据需要自定义级别。