98

我经常有一个主要的 R Markdown 文件或 knitr LaTeX 文件,其中我source有一些其他 R 文件(例如,用于数据处理)。但是,我认为在某些情况下,将这些源文件作为它们自己的可重现文档会是有益的(例如,一个 R Markdown 文件,它不仅包含用于数据处理的命令,而且还生成一个可重现的文档来解释数据处理决策)。

因此,我想source('myfile.rmd')在我的主 R Markdown 文件中有一个命令。这将提取和获取 .R 代码块中的所有 R 代码myfile.rmd。当然,这会产生错误。

以下命令有效:

```{r message=FALSE, results='hide'}
knit('myfile.rmd', tangle=TRUE)
source('myfile.R')
```

results='hide'如果需要输出,可以省略where 。即,knitr 从myfile.rmdinto输出 R 代码myfile.R

但是,它似乎并不完美:

  • 它会导致创建一个额外的文件
  • 如果需要控制显示,它需要出现在自己的代码块中。
  • 它不像简单那么优雅source(...)

因此我的问题是: 有没有更优雅的方式来获取 R Markdown 文件的 R 代码?

4

12 回答 12

37

看来您正在寻找单线。把这个放在你的 中.Rprofile怎么样?

ksource <- function(x, ...) {
  library(knitr)
  source(purl(x, output = tempfile()), ...)
}

但是,我不明白你为什么要source()在 Rmd 文件本身的代码。我的意思是knit()将运行本文档中的所有代码,如果您提取代码并以块的形式运行它,那么当您编写本文档时,所有代码将运行两次knit()(您自己在内部运行)。这两个任务应该是分开的。

如果你真的想运行所有代码,RStudio 让这变得相当简单:Ctrl + Shift + R. 它基本上调用purl()source()幕后。

于 2012-06-10T15:44:45.243 回答
22

将通用代码分解到一个单独的 R 文件中,然后将该 R 文件源到您想要的每个 Rmd 文件中。

例如,假设我需要制作两份报告,流感爆发和枪支与黄油分析。自然我会创建两个 Rmd 文档并完成它。

现在假设老板出现并想查看流感爆发与黄油价格的变化(控制 9 毫米弹药)。

  • 复制并粘贴代码以分析报告到新报告中对于代码重用等来说是一个坏主意。
  • 我希望它看起来不错。

我的解决方案是将项目分解为这些文件:

  • 流感
    • 流感数据导入.R
  • Guns_N_Butter.Rmd
    • guns_data_import.R
    • 黄油数据导入.R

在每个 Rmd 文件中,我都会有类似的内容:

```{r include=FALSE}
source('flu_data_import.R')
```

这里的问题是我们失去了可重复性。我的解决方案是创建一个公共子文档以包含到每个 Rmd 文件中。因此,在我创建的每个 Rmd 文件的末尾,我添加了以下内容:

```{r autodoc, child='autodoc.Rmd', eval=TRUE}
``` 

当然,还有 autodoc.Rmd:

Source Data & Code
----------------------------
<div id="accordion-start"></div>

```{r sourcedata, echo=FALSE, results='asis', warnings=FALSE}

if(!exists(autodoc.skip.df)) {
  autodoc.skip.df <- list()
}

#Generate the following table:
for (i in ls(.GlobalEnv)) {
  if(!i %in% autodoc.skip.df) {
    itm <- tryCatch(get(i), error=function(e) NA )
    if(typeof(itm)=="list") {
      if(is.data.frame(itm)) {
        cat(sprintf("### %s\n", i))
        print(xtable(itm), type="html", include.rownames=FALSE, html.table.attributes=sprintf("class='exportable' id='%s'", i))
      }
    }
  }
}
```
### Source Code
```{r allsource, echo=FALSE, results='asis', warning=FALSE, cache=FALSE}
fns <- unique(c(compact(llply(.data=llply(.data=ls(all.names=TRUE), .fun=function(x) {a<-get(x); c(normalizePath(getSrcDirectory(a)),getSrcFilename(a))}), .fun=function(x) { if(length(x)>0) { x } } )), llply(names(sourced), function(x) c(normalizePath(dirname(x)), basename(x)))))

for (itm in fns) {
  cat(sprintf("#### %s\n", itm[2]))
  cat("\n```{r eval=FALSE}\n")
  cat(paste(tryCatch(readLines(file.path(itm[1], itm[2])), error=function(e) sprintf("Could not read source file named %s", file.path(itm[1], itm[2]))), sep="\n", collapse="\n"))
  cat("\n```\n")
}
```
<div id="accordion-stop"></div>
<script type="text/javascript">
```{r jqueryinclude, echo=FALSE, results='asis', warning=FALSE}
cat(readLines(url("http://code.jquery.com/jquery-1.9.1.min.js")), sep="\n")
```
</script>
<script type="text/javascript">
```{r tablesorterinclude, echo=FALSE, results='asis', warning=FALSE}
cat(readLines(url("http://tablesorter.com/__jquery.tablesorter.js")), sep="\n")
```
</script>
<script type="text/javascript">
```{r jqueryuiinclude, echo=FALSE, results='asis', warning=FALSE}
cat(readLines(url("http://code.jquery.com/ui/1.10.2/jquery-ui.min.js")), sep="\n")
```
</script>
<script type="text/javascript">
```{r table2csvinclude, echo=FALSE, results='asis', warning=FALSE}
cat(readLines(file.path(jspath, "table2csv.js")), sep="\n")
```
</script>
<script type="text/javascript">
  $(document).ready(function() {
  $('tr').has('th').wrap('<thead></thead>');
  $('table').each(function() { $('thead', this).prependTo(this); } );
  $('table').addClass('tablesorter');$('table').tablesorter();});
  //need to put this before the accordion stuff because the panels being hidden makes table2csv return null data
  $('table.exportable').each(function() {$(this).after('<a download="' + $(this).attr('id') + '.csv" href="data:application/csv;charset=utf-8,'+encodeURIComponent($(this).table2CSV({delivery:'value'}))+'">Download '+$(this).attr('id')+'</a>')});
  $('#accordion-start').nextUntil('#accordion-stop').wrapAll("<div id='accordion'></div>");
  $('#accordion > h3').each(function() { $(this).nextUntil('h3').wrapAll("<div>"); });
  $( '#accordion' ).accordion({ heightStyle: "content", collapsible: true, active: false });
</script>

注意,这是为 Rmd -> html 工作流设计的。如果您使用乳胶或其他任何东西,这将是一个丑陋的混乱。此 Rmd 文档在全局环境中查找所有 source() 文件,并在文档末尾包含它们的源代码。它包括 jquery ui、tablesorter,并将文档设置为使用手风琴样式来显示/隐藏源文件。这是一项正在进行的工作,但请随时根据自己的用途进行调整。

不是单行,我知道。希望它至少能给你一些想法:)

于 2013-05-31T13:40:15.723 回答
8

试试 knitr 的 purl 函数:

source(knitr::purl("myfile.rmd", quiet=TRUE))

于 2020-02-03T10:24:36.010 回答
4

也许一个人应该开始思考不同的事情。我的问题如下:在 .R 文件中的 .Rmd 块中编写您通常拥有的每个代码。而对于你用来编织的 Rmd 文档,即一个 html,你只剩下

```{R Chunkname, Chunkoptions}  
source(file.R)  
```

这样,您可能会创建一堆 .R 文件,并且您失去了使用 ctrl+alt+n(或 +c,但通常这不起作用)处理所有代码“一个接一个块”的优势。但是,我读了 Gandrud 先生关于可重复研究的书,并意识到,他绝对使用 knitr 和 .Rmd 文件来创建 html 文件。主要分析本身是一个 .R 文件。如果您开始在内部进行整个分析,我认为 .Rmd 文档会迅速变得太大。

于 2015-04-29T07:07:30.363 回答
3

如果您只是在代码之后,我认为这些方面的东西应该可以工作:

  1. 阅读markdown/R文件readLines
  2. 用于查找代码块,例如grep搜索以开头的行<<<
  3. 获取包含原始行的对象的子集以仅获取代码
  4. 使用转储到临时文件writeLines
  5. 将此文件源到您的 R 会话中

将其包装在一个函数中应该可以满足您的需求。

于 2012-06-10T12:16:09.777 回答
3

以下黑客对我来说效果很好:

library(readr)
library(stringr)
source_rmd <- function(file_path) {
  stopifnot(is.character(file_path) && length(file_path) == 1)
  .tmpfile <- tempfile(fileext = ".R")
  .con <- file(.tmpfile) 
  on.exit(close(.con))
  full_rmd <- read_file(file_path)
  codes <- str_match_all(string = full_rmd, pattern = "```(?s)\\{r[^{}]*\\}\\s*\\n(.*?)```")
  stopifnot(length(codes) == 1 && ncol(codes[[1]]) == 2)
  codes <- paste(codes[[1]][, 2], collapse = "\n")
  writeLines(codes, .con)
  flush(.con)
  cat(sprintf("R code extracted to tempfile: %s\nSourcing tempfile...", .tmpfile))
  source(.tmpfile)
}
于 2016-10-19T20:40:46.670 回答
3

我使用以下自定义函数

source_rmd <- function(rmd_file){
  knitr::knit(rmd_file, output = tempfile())
}

source_rmd("munge_script.Rmd")
于 2019-08-30T22:11:22.970 回答
1

我建议将主要分析和计算代码保留在 .R 文件中,并根据需要在 .Rmd 文件中导入块。我已经解释了这里的过程。

于 2016-05-19T22:44:03.373 回答
1

sys.source("./your_script_file_name.R", envir = knitr::knit_global())

在调用 your_script_file_name.R 中包含的函数之前放置此命令。

如果您已经创建了项目,则在 your_script_file_name.R 之前添加“./”以显示文件的方向。

您可以查看此链接了解更多详情:https ://bookdown.org/yihui/rmarkdown-cookbook/source-script.html

于 2020-06-08T19:38:21.087 回答
1

我使用这个单线:

```{r optional_chunklabel_for_yourfile_rmd, child = 'yourfile.Rmd'}
```

请参阅: 我的 .Rmd 文件变得非常冗长。是否可以将它和 source() 拆分为主要 .Rmd 的较小部分?

于 2020-06-12T14:55:46.380 回答
0

这对我有用

source("myfile.r", echo = TRUE, keep.source = TRUE)
于 2020-02-26T16:57:42.823 回答
0

我想说没有一种更优雅的方式来获取 Rmarkdown 文件。Rmd 的精神是报告是可重复的,并且充其量是自包含的。但是,添加到 OP 的原始解决方案中,以下方法避免了在磁盘上永久创建中间文件。它还做出了一些额外的努力来确保块输出不会出现在渲染文档中:

knit_loc <- tempfile(fileext = ".R")
knitr::knit("myfile.rmd",
            output = knit_loc,
            quiet = TRUE,
            tangle = TRUE)
invisible(capture.output(source(knit_loc, verbose = FALSE)))

我还要补充一点,如果子 Markdown 依赖项在您的 R 环境之外(例如,将文件写入磁盘、下载一些外部资源、与 Web api 交互等),那么knit()我会选择rmarkdown::render()

rmarkdown::render("myfile.rmd")
于 2021-12-12T08:24:38.483 回答