3

ループ内の多くのデータセットで線形フィッティングを実行し、結果を pdf ファイルにプロットしています。コンソールから約 100 個のデータセットの要約を観察する代わりに、summary(fit) の出力を同じ pdf ファイルに直接保存することは可能ですか?

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

1 に答える 1

5

レポートを生成するには、 KnitrとRstudioを使用することを強くお勧めします。

ここでは、コードを使用して、3 つの簡単な手順で pdf を生成します。Rstudio がインストールされていると仮定します。

  1. 新しい R スウィーブ ファイルを作成します (メニューを使用)

  2. 2 つのチャンクを挿入する場所 (右側のチャンクを使用)

    <<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 ファイルを生成します。

要約とプロットの間に必要なものを挿入して、複雑なレポートを作成できます。

于 2013-01-07T12:41:04.277 に答える