Question

我在循环中对许多数据集进行线性拟合，并将结果绘制在pdf文件中。是否可以直接将摘要（拟合）的输出保存在同一个pdf文件中，而不是通过控制台观察大约100个数据集的摘要？

  LMmodel <- y ~ x
  fit <- lm(LMmodel, data = Dataset)

  pdf(file = OutputFile, width = 10, height = 6, paper = "a4r")

  xLim = range(x)
  yLim = range(y)

  plot(x, y, type = "p", xlim = xLim, ylim = yLim,
       main = plotTitle, xlab = "x [m]", ylab = "y [dB]",
       pch = 20, cex = .9)
  regLine(fit, col=palette()[2], lwd=2, lty=1)
  grid(lwd = 1.5)

  plot(density(residuals(fit)), main = "Density Plot of the Residuals"))

  dev.off()
  graphics.off()
  return(summary(fit))

Answer 1

我真的建议Knitr与Rstudio生成报告。

在这里，我使用您的代码通过3个简单的步骤生成PDF。我假设你已经安装了Rstudio。

我创建了一个新的R形文件（使用菜单）

我插入2个chuncks（使用右边的 Chunks ）

<<myplot,echo=FALSE,fig=TRUE>>=
library(car)
x <- rnorm(n=20,mean=30,sd=20)
y <- rnorm(n=20,mean=180,sd=10)
Dataset <- data.frame(x=x,y=y)
LMmodel <- y ~ x
fit <- lm(LMmodel, data = Dataset)
xLim = range(x)
yLim = range(y)
plot(x, y, type = "p", xlim = xLim, ylim = yLim,
    main = "plotTitle", xlab = "x [m]", ylab = "y [dB]",
    pch = 20, cex = .9)
regLine(fit, col=palette()[2], lwd=2, lty=1)
grid(lwd = 1.5)
plot(density(residuals(fit)), main = "Density Plot of the Residuals")
@

摘要是：

    <<mysummary>>=
    print(summary(fit))
    @

使用编译PDF 按钮生成pdf文件。

您可以在摘要和图表之间插入所需内容，以构建复杂的报告。

在pdf文件或类似文件中写下fit的摘要

1 个答案: