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

Lattice 将为每个条件变量创建一个条带。如果您只想要一条，请尝试对交互进行调节。例如

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

这给出了

您可以提前创建交互并自定义关卡，如果您愿意，可以使用您想要的任何描述。

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

Faceted Lattice Plots in R, e.g., wireframes: How to remove strips and add 1-Line subtitles

r

lattice