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

Question

我正在循环中对许多数据集进行线性拟合,并将结果绘制在pdf文件中.是否可以直接将摘要(fit)的输出保存在同一个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.

1. 我创建一个新的R形文件(使用菜单)

2. 在哪里插入2个chuncks(使用右边的块)

   <<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))
    @

1. 您可以使用编译PDF按钮生成pdf文件.

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