r - 直方图和 qqplot 与 par 相邻 - 不显示在 PDF 中

Question

我有这个我绑定(bind)的配对样本

paired <- cbind(c(before, after))

然后我构建了一个直方图和一个 Q-Q 图

hist(paired, freq=FALSE, ylim=range(0,0.7,by=0.2))
curve(dnorm(x, mean=mean(paired), sd=sd(paired)), add=TRUE)
qqnorm(paired)
qqline(paired)
par(mfrow = c(1,2))

如果我看“图”，它非常好，它显示这两个完美地相邻。但是，当我想编译成 PDF 时，它只给我 PDF 中的直方图？有人可以告诉我我做错了什么吗？

Answer 1

您的原始代码有问题

par() 在所有绘图完成后，您放置它时不起作用。您将在 2 个不同的图形设备上获得两个图。如果将它们保存到 .pdf 文件中，它们将出现在同一个文件中，但会出现在两个不同的页面上。

解决方案

# You should put `par(mfrow = c(1,2))` **before** plotting;
par(mfrow = c(1,2))  ## split the plotting region in to 1 row 2 columns
hist(paired, freq=FALSE, ylim=range(0,0.7,by=0.2))
curve(dnorm(x, mean=mean(paired), sd=sd(paired)), add=TRUE)
qqnorm(paired)
qqline(paired)

另外:您可以使用 pdf() 直接在 .pdf 文件上绘图。请参阅 ?pdf。但记得在打开文件之前使用 def.off() 关闭它。

pdf("test.pdf", height = 6, width = 6)   ## open a graphical device, in this case, a .pdf file
## code above for producing plot
dev.off()  ## close active graphical device, in this case, you .pdf file

生成的 .pdf 文件位于您当前的工作目录中。如果您不知道它在哪里，请在 R 中执行 getwd()。