R 6个模型的一个图

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

我的问题是：如何使用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)

fitemax使用（mfrow=c（2,3））将以下绘图安排在一个2*3的网格中

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

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


Update:在您的例子中，par（mfrow）似乎不起作用，因为我认为它实际上并没有调用base plot方法，相反，来自
fitMod
方法的返回对象实际上是一个名为“grillis”的类型，它属于lattice包。如果您想了解更多关于格架的信息，请阅读。但是，如果您只是想知道如何让它工作，我使用
grid.arrange
方法从
gridExtra
获得了它

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".

库（DoseFinding）
图书馆（gridExtra）
数据（生物医学）
#这里，bnds参数已被忽略，因此将应用defBnds的默认值。
fitemax理想的输出是什么？这个“一个情节”会是什么样子？您提供的代码基本上是无用的，因为您没有提供任何数据，所以我们可以实际运行它，所以我们不知道这些图是什么样子。请创建一个真正的答案，这不适用于我的数据，数据存储在剂量查找包下的r中。数据是（biom）。我尝试了所有方法，但仍然单独显示。如果这个答案不适用于您的特定数据集，那么您应该提供一个可复制的（！）与最重要的特性相似的数据示例；fitlogistic plot（fitlogistic）plot（fitlogistic）中的错误：未找到对象“fitlogistic”>@SAMA，很明显，您正试图对
bnd
参数传递一些限制，如何使用一个特定的问题已经超出了这个问题的范围，我建议你开始另一个线程，让有专业知识的人给出答案。然而，根据我所知。bnd
bnd
将拾取一些默认值，这些值实际执行起来非常不合时宜。正如你从情节中看到的那样。@B.M.W，我尝试了上面的命令，得到了这个消息fitemax
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".