35

まず第一に、私はknitrと再現性のある分析の概念に非常に慣れていないことを認めなければなりませんが、現在のワークフロー（ワードドキュメントへの多くのコピー貼り付けを含む）を改善する可能性を見ることができます。

私はしばしばグループ（この例では病院）ごとに複数のレポートを作成する必要があり、各病院内には、結果を報告しているさまざまな病棟が存在する可能性があります。以前は、ループを使用してRですべてのプロットと分析を実行し、その後、コピー/貼り付け作業を開始しました。ただし、この投稿を読んだ後（Sweaveは多くのPDFを自動的に生成できますか？）、実際に多くの手順をスキップして、RからRnw/knitrを介してレポートすることができるかもしれないという希望を与えました。

ただし、試してみると、うまく機能していないものがあることがわかります（Rnw内のR環境は、渡そうとしているループ変数を認識していないようです??）。 

   ##  make my data
Hospital <- c(rep("A", 20), rep("B", 20))
Ward <- rep(c(rep("ICU", 10), rep("Medicine", 10)), 2)
Month <- rep(seq(1:10), 4)
Outcomes <- rnorm(40, 20, 5)
df <- data.frame(Hospital, Ward, Month, Outcomes)


##  Here is my current work flow-- produce all plots, but export as png and cut/paste
for(hosp in unique(df$Hospital)){
  subgroup <- df[ df$Hospital == hosp,]
  for(ward in unique(subgroup$Ward)){
    subgroup2 <- subgroup[subgroup$Ward == ward,]
    savename <- paste(hosp, ward)
    plot(subgroup2$Month, subgroup2$Outcomes, type="o", main=paste("Trend plot for", savename))
  }
}
# followed by much copy/pasting


##  Here is what I'm trying to go for using knitr 
library(knitr)
for (hosp in unique(df$Hospital)){
  knit("C:file.path\\testing_loops.Rnw", output=paste('report_', Hospital, '.tex', sep=""))
}

## With the following *Rnw file
## start *.Rnw Code
\documentclass[10pt]{article}
\usepackage[margin=1.15 in]{geometry}
<<loaddata, echo=FALSE, message=FALSE>>=
  Hospital <- c(rep("A", 20), rep("B", 20))
Ward <- rep(c(rep("ICU", 10), rep("Medicine", 10)), 2)
Month <- rep(seq(1:10), 4)
Outcomes <- rnorm(40, 20, 5)
df <- data.frame(Hospital, Ward, Month, Outcomes)
subgroup <- df[ df$Hospital == hosp,]
@

\begin{document}
<<setup, echo=FALSE >>=
  opts_chunk$set(fig.path = paste("test", hosp , sep=""))
@

Some infomative text about hospital \Sexpr{hosp}

<<plots, echo=FALSE >>=
  for(ward in unique(subgroup$Ward)){
    subgroup2 <- subgroup[subgroup$Ward == ward,]
    #     subgroup2 <- subgroup2[ order(subgroup2$Month),]
    savename <- paste(hosp, ward)
    plot(subgroup2$Month, subgroup2$Outcomes, type="o", main=paste("Trend plot for", savename))
  }
@
\end{document}


##  To be then turned into pdf with this
tools::texi2pdf("C:file.path\\report_A.tex", clean = TRUE, quiet = TRUE)

knit（）コードチャンクを実行しようとすると、次のエラーが発生します。

Error in file(con, "w") : invalid 'description' argument

また、*。texファイルが作成されるディレクトリを調べると、病院Aから2つのpdfプロットが作成され（Bの場合はなし）、pdfに編成する病院固有の*.texファイルがないことがわかります。あなたが提供できるどんな助けにも前もって感謝します！

4

3 に答える 3

16

ファイル内のデータを再定義する必要はありません。警告は、出力名を（ループインデックス）ではなく（病院の完全なベクトル）と.Rnw一緒に配置しているという事実から来ていると思います。Hospitalhosp

あなたの例に従うtestingloops.Rnwと、

\documentclass[10pt]{article}
\usepackage[margin=1.15 in]{geometry}
<<loaddata, echo=FALSE, message=FALSE>>=
subgroup <- df[ df$Hospital == hosp,]
@

\begin{document}
<<setup, echo=FALSE >>=
  opts_chunk$set(fig.path = paste("test", hosp , sep=""))
@

Some infomative text about hospital \Sexpr{hosp}

<<plots, echo=FALSE >>=
  for(ward in unique(subgroup$Ward)){
    subgroup2 <- subgroup[subgroup$Ward == ward,]
    #     subgroup2 <- subgroup2[ order(subgroup2$Month),]
    savename <- paste(hosp, ward)
    plot(subgroup2$Month, subgroup2$Outcomes, type="o", main=paste("Trend plot for", savename))
  }
@
\end{document}

ドライバーのRファイルは

##  make my data
Hospital <- c(rep("A", 20), rep("B", 20))
Ward <- rep(c(rep("ICU", 10), rep("Medicine", 10)), 2)
Month <- rep(seq(1:10), 4)
Outcomes <- rnorm(40, 20, 5)
df <- data.frame(Hospital, Ward, Month, Outcomes)

## knitr loop
library("knitr")
for (hosp in unique(df$Hospital)){
  knit2pdf("testingloops.Rnw", output=paste0('report_', hosp, '.tex'))
}
于 2013-03-13T22:34:59.150 に答える
11

素晴らしい質問です！これは、あなたがあなたの質問で提供した他のビットで私のために働きます。私はあなたhospをただのに置き換えたことに注意してくださいx。私はあなたのRnwファイルを呼び出しましたtest.rnw

# input data
Hospital <- c(rep("A", 20), rep("B", 20))
Ward <- rep(c(rep("ICU", 10), rep("Medicine", 10)), 2)
Month <- rep(seq(1:10), 4)
Outcomes <- rnorm(40, 20, 5)
df <- data.frame(Hospital, Ward, Month, Outcomes)

# generate the tex files, one for each hospital in df
library(knitr)
lapply(unique(df$Hospital), function(x) 
       knit("C:\\emacs\\test.rnw", 
            output=paste('report_', x, '.tex', sep="")))

# generate PDFs from the tex files, one for each hospital in df
lapply(unique(df$Hospital), function(x)
       tools::texi2pdf(paste0("C:\\emacs\\", paste0('report_', x, '.tex')), 
                       clean = TRUE, quiet = TRUE))

私はあなたのループをlapply無名関数に置き換えましRた。

ここで、ファイル内のどこでを置き換えたかを確認hospできxますrnw。

\documentclass[10pt]{article}
\usepackage[margin=1.15 in]{geometry}
<<loaddata, echo=FALSE, message=FALSE>>=
  Hospital <- c(rep("A", 20), rep("B", 20))
Ward <- rep(c(rep("ICU", 10), rep("Medicine", 10)), 2)
Month <- rep(seq(1:10), 4)
Outcomes <- rnorm(40, 20, 5)
df <- data.frame(Hospital, Ward, Month, Outcomes)
subgroup <- df[ df$Hospital == x,]
@

\begin{document}
<<setup, echo=FALSE >>=
  opts_chunk$set(fig.path = paste("test", x , sep=""))
@

Some informative text about hospital \Sexpr{x}

<<plots, echo=FALSE >>=
  for(ward in unique(subgroup$Ward)){
    subgroup2 <- subgroup[subgroup$Ward == ward,]
    #     subgroup2 <- subgroup2[ order(subgroup2$Month),]
    savename <- paste(x, ward)
    plot(subgroup2$Month, subgroup2$Outcomes, type="o", main=paste("Trend plot for", savename))
  }
@
\end{document}

結果は、2つのtexファイル（report_A.tex、report_B.tex）、図用の4つのPDF（A1、A2、B1、B2）、およびレポート用の2つのPDF（report_A.pdf、report_B.pdf）であり、それぞれに図があります。それらの中で。それはあなたが求めていたものですか？

于 2013-03-13T21:59:52.067 に答える
1

この回答では、特定の例ではなく、「ループを使用して複数のPDFレポートを作成する」というより一般的な質問に回答するつもりです。これは、この傾向を初心者として追跡するのが非常に困難だったためです。私はそれを最終的に（htmlバージョンで）動作させることができたので、これは私の謙虚な解決策です。おそらくここに公開されているより良いものがいくつかありますが、私はまだそれらを完全に理解することはできません。

  1. デザインを使用してRMDファイルを作成し、working \ inputディレクトリ（Rstudio：file-> newfile-> Rマークダウン）に保存します。このファイルには、レポートでプロットを作成するために必要なすべての関数が含まれている必要があります（これらのコードチャンクの1つでそれらを宣言するだけです）。このファイルは、今後のすべてのレポートのテンプレートと考えてください。以前にデータをかみ砕いた後、データをその環境に渡すことについて心配する必要はありません。これについては（2）で説明します。理解しておくべき重要な問題は、すべての計算がパイプラインのさらに下流で行われることです（RMDファイルをレンダリングする時点で）。

  2. 異なるコントロールrファイルで使用する必要のあるループを作成します。私の場合、ディレクトリ内のすべてのファイルを繰り返し処理し、それらをデータフレームに取り込むループがあります。次に、それらをプロットするために、他のデータ変数とともに、それらのデータフレームをRMDに渡します。これがその方法です。

    run_on_all<-function(path_in="path:\\where\\your\\input\\and\\RMD\\is", path_out="path:\\where\\your\\output\\will\\be") setwd(path_in) ibrary(rmarkdown) library(knitr) list_of_file_names=list.files(path = getwd, pattern = "*.csv") #this gets a list of the input files names for (file_name in list_of_file_names) { data=read.csv(file_name) #read file into data frame report_name=paste(some_variable_name,".html",sep="") render("your_template.Rmd",output_file =report_name,output_dir =path_out,list(data,all other parameters you want to input into the RMD))} }

  3. 最も重要なコマンドは、render関数の呼び出しです。それはあなたが望むどんなparamenterでもRMD環境に投げることを可能にします。また、レポートの名前を変更したり、出力場所を変更したりすることもできます。さらに、これを呼び出すとレポートも生成されるため、すべてを1行で取得できます（RMDの呼び出しが関数内にある場合、入力した変数が欠落していることに気付く場合がありますが、レポートはまだ正しく公開されています）

まとめ

必要なファイルは2つあります。RMDファイルは、すべての追加レポートのテンプレートになり、制御ファイルになります。制御ファイルはデータを取得し、それを噛み砕き、噛み砕いたパラメーターをRMDに渡します（レンダリング関数を介して）。RMDはデータを取得し、いくつかの計算を実行し、それをプロットして、新しいファイルに公開します（これもrender関数によって）。お役に立てば幸いです。

于 2018-01-01T18:51:48.547 に答える