6个模型的一张图

Question

我在一个单独的图中绘制了六个模型。

我的问题是:如何使用R将它们组合到一个图形绘图中？

我使用的代码是：

fitemax <- fitMod(dose, resp, data=biom, model="emax",bnds = c(0.00, 1))
plot(fitemax)

fitlinearlog <- fitMod(dose, resp, data=biom, model="linlog",bnds = c(0.00, 1))
plot(fitlinearlog)

fitlinear <- fitMod(dose, resp, data=biom, model="linear",bnds = c(0.00, 1))
plot(fitlinear)

fitquadratic <- fitMod(dose, resp, data=biom, model="quadratic",bnds = c(0.00, 1))
plot(fitquadratic)

fitexponential <- fitMod(dose, resp, data=biom, model="exponential",bnds = c(0.00,1))
plot(fitexponential)

fitlogistic <- fitMod(dose, resp, data=biom, model="logistic",defBnds(MaxEff, 
logistic = matrix(c(0.0, 0.2, 0.4, 0.8)*MaxEff, 2)))
plot(fitlogistic)

DoseFinding包中R中的数据仓库

Answer 1

使用par(mfrow=c(2,3))将以下绘图排列在一个2*3网格中。

如果您想要精细控制，请继续阅读here(布局)，here(ggplot+gridExtra)

png(filename="C:\\Users\\datafireball.com\\Documents\\R\\stackoverflow_7144118.png")
par(mfrow=c(3,2))
plot(rnorm(10), rnorm(10))
plot(rnorm(10), rnorm(10))
plot(rnorm(10), rnorm(10))
plot(rnorm(10), rnorm(10))
plot(rnorm(10), rnorm(10))
plot(rnorm(10), rnorm(10))
dev.off()

您可以删除第一行和最后一行，以便可以将其打印到突出的输出。

​

​

更新:在您的例子中，看起来不会工作，因为我不认为它实际上是在调用基本的plot方法，相反，从fitMod方法返回的对象实际上是一个名为"trellis“的类型，它属于网格包。如果你想了解更多关于trellis的信息，请阅读here。但是，如果您只想知道如何使其工作，我可以使用来自gridExtra的grid.arrange方法使其工作。

library(DoseFinding)
library(gridExtra)
data(biom)
# here, the bnds argument has been ignored so the default value from defBnds will be applied. 
fitemax <- fitMod(dose, resp, data=biom, model="emax")
p1 <- plot(fitemax)
fitlinearlog <- fitMod(dose, resp, data=biom, model="linlog")
p2 <- plot(fitlinearlog)
fitlinear <- fitMod(dose, resp, data=biom, model="linear")
p3 <- plot(fitlinear)
fitquadratic <- fitMod(dose, resp, data=biom, model="quadratic")
p4 <- plot(fitquadratic)
fitexponential <- fitMod(dose, resp, data=biom, model="exponential")
p5 <- plot(fitexponential)
fitlogistic <- fitMod(dose, resp, data=biom, model="logistic")
p6 <- plot(fitlogistic)
grid.arrange(p1, p2, p3, p4, p5, p6)
# Message: Need bounds in "bnds" for nonlinear models, using default bounds from "defBnds".

​

​

这是您想要的输出吗？