I use this kind of function (from iris data set)
model_test <- lm(Sepal.Length ~( Petal.Length + Sepal.Width + Petal.Width +Species)^2,
data=iris)
gg<-expa
Lattice is going to create a strip for each conditioning variable. If you want just one strip, try conditioning on the interaction. For example
wireframe(pred~Sepal.Width+Petal.Width|interaction(Species,Petal.Length),
pd, drape=FALSE,scale=list(arrows=FALSE), subset=(Species=="setosa"), layout=c(3,3))
which gives
You can create the interaction ahead of time and customize the levels if you like with whatever description you want.