El proceso habitual para hacer un informe (o unas transparencias) en el que aparezcan gráficos o tablas resumen de algún análisis estadístico consiste en :
Escribir el texto en un programa (Word, Powerpoint, Prezi, etc.)
Realizar los cálculos estadísticos y gráficos en otro programa (R, Stata, Eviews, etc.)
Pegar los gráficos y tablas en el documento de texto.
Este proceso tiene ciertas desventajas: dificulta la investigación reproducible y puede ser tedioso de rehacer si por ejemplo cambian ligeramente los datos, etc.
En R es posible realizar todo el informe, tanto la escritura del texto como la realización de los cálculos y gráficos, en un único documento. Hay varias formas de hacerlo, pero nos centraremos en los documentos RMarkdown (.Rmd).
Los documentos RMarkdown (.Rmd) facilitan mucho la realización de informes y transparencias ya que permiten combinar texto, código y resultados de la evaluación del código en un único documento. Si por ejemplo cambian los datos sólo habría que cambiar la ruta a los nuevos datos y el informe (gráficos, tablas etc…) se volvería a generar automáticamente con los nuevos datos.
La razón para hacer los análisis reproducibles no sólo es cumplir con los estándares científicos, que también, sino también hay un interés personal para el analista. Para entenderlo puedes ver el siguiente video de 1’44’’ que muestra una de las principales ventajas de usar documentos .Rmd.
Hace poco David Keyes pregunto en Twitter a la comunidad R cuales eran los principales beneficios de usar ficheros .Rmd para tus investigaciones/informes:
Help me out #rstats hive mind. I've been asked to explain the benefits of RMarkdown. I'm developing a list on a continuum from simple (i.e. do what you currently do, just better) to complex (i.e. completely transform how you work). Here 'tis:
Con las respuestas al tweet, David elaboró este post.
Ya hemos trabajado con algunos documentos Rmarkdown (.Rmd); de hecho, todos los tutoriales del curso (incluido éste que estás viendo) se han elaborado usando ficheros .Rmd. Los tutoriales se escriben en ficheros RMarkdown, pero estos ficheros .Rmd se convierten a html, que suele ser el formato final en el que se muestran los tutoriales. En RStudio el proceso de conversión de .Rmd a .html consiste simplemente en pinchar un botón; además, simplemente cambiando una linea podemos convertir los documentos .Rmd a una gran variedad de formatos: html, pdf, word, ioslides, beamer, etc…
Si quieres ver algunos ejemplos de la gran variedad de formatos a los que puedes transformar un documento .Rmd ve a la siguiente galería. Para iniciarte en el universo Rmarkdown puedes consultar R Markdown: The Definitive Guide.
Proceso para convertir los .Rmd a otro(s) formatos [OPCIONAL]
Como trabajamos con RStudio, en la práctica, procesar los ficheros .Rmd consistirá solamente en pinchar en el icono Knit1. Muy fácil!!
No es necesario, pero quizá os interese saber cómo se procesan realmente los ficheros .Rmd para acabar convirtiéndose en html, pdf, etc.
La respuesta es que se ocupa de ello el paquete rmarkdown que llama otro paquete de R, knitr y a un programa llamado pandoc. Si quieres más detalles puedesir aquí
En palabras: knitr se ocupa de ejecutar todos los trozos con código R que haya en el fichero .Rmd, después de ejecutar el código, pegará los resultados de la evaluación del código (gráficos, tablas etc…) junto con el texto en un documento intermedio (con extensión .md), para después transferir, con la ayuda del paquete rmarkdown, este documento .md a pandoc que se encargará de traducirlo al formato elegido (html, pdf, …)2. Knitr sabe diferenciar el texto del código R porque éste se señaliza con unas marcas.
Visualmente:
2 Creando .Rmd’s en RStudio
En la práctica, RStudio facilita mucho la creación de documentos Rmarkdown. Para generar un documento .Rmd basta con seguir la siguiente ruta de menús: File > New File > R Markdown ...
Se abrirá una ventana que nos solicitará un título y un autor para nuestro .Rmd, así como el formato de salida.
Cuando aceptemos nos generará un documento/plantilla para nuestro .Rmd. Si queremos procesarlo o “knitearlo” tendremos que hacer click en el icono Knit
3. ¿Qué son los documentos .Rmd?
Son simplemente ficheros de texto (se pueden escribir en cualquier editor de texto, por ejemplo Notepad); PERO facilitan mucho la tarea de generar informes o transparencias con contenido estadístico, ya que permiten mezclar en un mismo documento texto y código R.
El código R (así como los resultados de la evaluación del código) se mostrarán automáticamente (gracias a knitr) en el documento final; de esta forma, se facilita mucho la realización de informes y transparencias ya que evita el tener que ir copiando los resultados (tablas, gráficos etc…) en el informe.
Veamos un ejemplo con un documento Rmarkdown muy sencillo:
---
title: "El título de mi informe"
author: "Yo mismo"
date: "Marzo de 2017"
output: html_document
---
En este informe haremos un resumen (o summary) del conjunto de datos iris.
```{r}
summary(iris)
```
El trozo de arriba es un chunk de código R. Cuando compile el documento, knitr ejecutará el código y mostrará los resultados en el documento final.
Con esto finalizamos este informe.
Como veis, los documentos Rmarkdown tienen 3 partes o elementos:
encabezamiento
trozos de código R
texto
Luego hablaremos de ellas. Antes veamos como quedaría el documento .Rmd tras ser procesado por knitr (“kniteado”). Tras pasar por knitr y pandoc se generará un documento .html que se verá así:
Como veis, el chunk con código R summary(iris) se ha ejecutado y se ha mostrado tanto el código como el resultado de su evaluación en unas cajas de texto.
4. Partes de los ficheros .Rmd
Los documentos Rmarkdown tienen 3 partes o elementos:
Encabezamiento o YAML header
Trozos de código R(chunks)
Texto (escrito en markdown)
Veámoslas una a una.
Encabezamiento (YAML header)
Abajo tenéis un ejemplo sencillo de un YAML header:
---
title: "El título de mi informe"
author: "Yo mismo"
date: "Marzo de 2017"
output: html_document
---
Como podéis imaginar, el encabezamiento se coloca al principio del documento y comienza y acaba con una marca de 3 guiones: ---
En el encabezamiento se introducen elementos básicos del documento como el título, el autor, fecha, y el formato de salida del documento. En el ejemplo hemos elegido como formato de salida html; si preferimos pdf, habría que sustituir output: html_document por output: pdf_document
El YAML header puede incluir otros elementos para personalizar un poco más el documento final o output. Abajo tenéis un ejemplo más complejo de un YAML header. Es el que utilizo para hacer los tutoriales del curso.
En este post y en esta vignette del paquete ymlthis tienes más opciones/posibilidades a especificar en el YAML.
Chunks (o código R)
Los trozos de código R o chunks permiten hacer análisis estadísticos y mostrar los resultados en el documento final.
Para que knitr distinga las instrucciones de R del texto normal tenemos que poner las instrucciones de R dentro de unas marcas o identificadores: ```{r} al principio y ``` al final.
Por ejemplo:
```{r}
summary(iris)
```
Knitr interpreta ese trozo de texto como instrucciones de R porque van dentro de las marcas, y hará que R las ejecute y muestre los resultados en el documento final.
Los chunks pueden tienen diversas opciones que permiten una mayor flexibilidad en como se muestra el código y los resultados en el documento final. Las opciones más usadas son:
echo
eval
Por ejemplo, si introducimos este texto en nuestro fichero .Rmd:
En este caso, se mostrará el chunk(echo = TRUE) y también se evaluará (eval = TRUE) y, por tanto, se mostrarán los resultados en el documento final. Se verá algo como:
summary(iris)
#> Sepal.Length Sepal.Width Petal.Length Petal.Width
#> Min. :4.300 Min. :2.000 Min. :1.000 Min. :0.100
#> 1st Qu.:5.100 1st Qu.:2.800 1st Qu.:1.600 1st Qu.:0.300
#> Median :5.800 Median :3.000 Median :4.350 Median :1.300
#> Mean :5.843 Mean :3.057 Mean :3.758 Mean :1.199
#> 3rd Qu.:6.400 3rd Qu.:3.300 3rd Qu.:5.100 3rd Qu.:1.800
#> Max. :7.900 Max. :4.400 Max. :6.900 Max. :2.500
#> Species
#> setosa :50
#> versicolor:50
#> virginica :50
#>
#>
#>
Mientras que si en el .Rmd escribimos lo siguiente:
NO se mostrará el código (echo = FALSE), SI se evaluará (eval = FALSE), PERO como results = "hide" NO se mostrarán los resultados en documento final.
Hay más opciones sobre los chunks que nos permiten una mayor flexibilidad sobre como mostrar los resultados y el código; pero si quieres ver todas las opciones tendrás que ir a la página web de knitr o al cheat sheet sobre Rmarkdown.
Una opción útil es include = FALSE; en este caso, el chunk se ejecutará, pero ni se mostrará en el documento final ni se mostrarán los resultados de la ejecución del código. Esta opción es muy útil para los chunks que se utilizan para hacer el “setup”.
Generalmente los documentos .Rmd tienen un primer chunk (chunk de setup) donde se fijan opciones globales para los chunks, fijra opciones globales, incluso se puede utilizar para cargar paquetes básicos como el tidyverse. Por ejemplo este suele ser el primer chunk en mis documentos .RMd:
La parte principal de un informe suele ser texto (narratives). En un fichero .Rmd, todo lo que no sea encabezamiento o chunks será interpretado por knit como texto y lo mostrará tal cual; es decir, como texto.
Aquí podríamos acabar nuestro tutoría sobre como escribir texto en un fichero .Rmd; pero generalmente en un texto queremos resaltar ciertas palabras con negrita, o ponerlas en cursiva, o poner un titulo de sección y de sub-secciones. Todo esto lo tendremos que hacer utilizando Markdown. Markdown es un lenguaje de marcas ligero y muy sencillo de aprender (lo básico se aprende en unos 10 minutos) pero muy utilizado. Es posible que alguno de vosotros haya utilizado una variante de Markdown al escribir en negrita en Whatsapp.
El texto de un documento .Rmd es “simplemente” texto PERO está escrito en Markdown. Ahora lo veremos, pero antes vamos a ver un poco más sobre qué es Markdown.
¿Qué es Markdown? [OPCIONAL]
Mardown es un lenguaje de marcado ligero ideado en 2004 por Jhon Grueber y Aaron Swartz. Hay diversas variantes de Markdown, aquí peuedes leer sobre un intento de estandarización.
Podemos pensar que Markdown es un método de escritura3: evidentemente sirve para escribir. La ventaja de escribir en Markdown es que es un lenguaje muy fácil de aprender y que como está basado en un formato de texto plano, es y será compatible con la mayoría de plataformas.
La mayoría de vosotros escribís en Word. Es muy fácil escribir en Word pero un archivo word solo es posible leerlo en el programa WORD de Microsoft. Si intentas abrir un documento .doc en Notepad, éste será completamente ilegible; sin embargo, la mayoría de plataformas y servicios web saben interpretar y mostrar correctamente un documento escrito en Markdown.
Lo que escribas en Rmarkdown se mostrará tal cual en el documento final, pero lo más habitual es que quieras dar un poco de formato el texto: negritas, cursivas, listas, enlaces de internet, etc…
Todos estos formatos (negrita, …) se introducen en rmarkdown con marcas; por ejemplo si quieres que una palabra se resalte en negritas tienes que escribirla enmarcada en **: **esto se mostraría en negrita**
Para aprender las principales reglas de Rmarkdown podéis usar un editor on-line de Markdown y probar a escribir algo.
Markdown is a way to write content for the web. It’s written in what nerds like to call “plaintext”, which is exactly the sort of text you’re used to writing and seeing. Plaintext is just the regular alphabet, with a few familiar symbols. Unlike cumbersome word processing applications, text written in Markdown can be easily shared between computers, mobile phones, and people. It’s quickly becoming the writing standard for academics, scientists, writers, and many more. Websites like GitHub and reddit use Markdown to style their comments.
Aquí tienes algunas reglas de Markdown
5. Más cosas de RMarkdown
En la página web de Markdown, concretamente aquí nos avisan de lo siguiente:
Markdown is not a replacement for HTML, or even close to it. Its syntax is very small, corresponding only to a very small subset of HTML tags … The idea for Markdown is to make it easy to read, write, and edit prose. HTML is a publishing format; Markdown is a writing format. Markdown was not designed to solve everything.
A pesar de que, como nos avisó Jhon Grueber, el desarrollador de Markdown, Markdown no está diseñado para resolver todas las necesidades de un escritor/científico, en Rmarkdown se pueden introducir también elementos como:
Ecuaciones
Se pueden introducir formulas matemáticas escritas en Látex. Para formulas en linea se usa la marca $ y para formulas independientes se usa $$.
Para fórmulas en linea (o dentro del texto) se utiliza la marca $ al principio y final de la formula. Por ejemplo $\sum_{i=1}^n X_i$ se mistaría así: \(\sum_{i=1}^n X_i\). Ves, la formula está dentro del texto, en una linea del texto.
Para presentar una ecuación independiente (en una linea independiente), se usa la marca $$ al principio y final de la formula. Si escribes $$E = mc^{2}$$, se mostrará en una linea independiente tal que así:
\[E = mc^{2}\]
¿Que pasa, que no sabes escribir formulas o ecuaciones en Latex? Yo tampoco mucho, pero puedes utilizar un programa como Lyx, o mucho más fácil, puedes utilizar algún editor online de Latex, por ejemplo: este o este. En este libro tienen algunos ejemplos de ecuaciones en Latex.
Recientemente, el paquete equatiomatic permite obtener fácilmente la ecuación de un modelo:
# remotes::install_github("datalorax/equatiomatic")
library(equatiomatic)
# Fit a simple model
mod1 <- lm(mpg ~ cyl + disp, mtcars)
# Give the results to extract_eq
extract_eq(mod1)
Para poner hyperlinks se puede escribir lo siguiente: <http://www.wikipedia.es> y se mostrará así: http://www.wikipedia.es
Pero es mejor ponerlo así: [enlace a la Wikipedia](http://www.wikipedia.es) y se mostrará así: enlace a la Wikipedia.
Si quieres que el enlace se abra en el navegador en una página nueva has de añadir {target="_blank"}
Notas al pie de página
Para poner notas al pie has de poner [^1] y luego al final del documento poner [^1]: esto es una nota al pie., y se verá esto:4
O alternativamente poner, en el sitio del texto donde quieras insertar una nota al pie, esta marca: ^[el texto que quiero que se lea en la nota al pie] y se mostrará como una nota al pie.
Tablas
Habrá un tutorial especifico para tablas, pero para mostrar una tabla, lo más básico y sencillo es utilizar la función knitr::kable():
```{r}
knitr::kable(summary(iris))
```
Sepal.Length
Sepal.Width
Petal.Length
Petal.Width
Species
Min. :4.300
Min. :2.000
Min. :1.000
Min. :0.100
setosa :50
1st Qu.:5.100
1st Qu.:2.800
1st Qu.:1.600
1st Qu.:0.300
versicolor:50
Median :5.800
Median :3.000
Median :4.350
Median :1.300
virginica :50
Mean :5.843
Mean :3.057
Mean :3.758
Mean :1.199
NA
3rd Qu.:6.400
3rd Qu.:3.300
3rd Qu.:5.100
3rd Qu.:1.800
NA
Max. :7.900
Max. :4.400
Max. :6.900
Max. :2.500
NA
6. Aún más cosas
Hay más cosas, como referencias bibliográficas, pero ya será para el segundo curso de R o lo tendrás que aprender por tu cuenta en R Markdown: The Definitive Guide.
A pesar de que Markdown permite formatear el texto, en ciertos sentidos es limitado; pero si quieres aún más flexibilidad en el output, tendrás que aprender html y mejor html+CSS.
Algunos ejemplos:
Si quieres introducir un párrafo en otro color tendrás que hacerlo en html, tendrás que escribir: <FONT COLOR="Red">Esto se mostrará en ROJO!!</FONT>
Esto se mostrará en ROJO!!
si quieres centrar un párrafo, tendrás que escribir: <CENTER>Este párrafo irá centrado</CENTER>
Este párrafo irá centrado
Puedes insertar en tu documento (sólo si el output es html) una página web completa, para ello has de escribir: <iframe src="http://www.eldiario.es/" height="400" width="800"></iframe>
Un vídeo. Has de escribir: <iframe width="560" height="315" src="https://www.youtube.com/embed/ACv9zaBa1A4" frameborder="0" allowfullscreen></iframe>
Aunque también hay un paquete, al menos, que facilita la inserción de videos en documentos .Rmd: es el paquete vembedr.
Aún más cosas (II)
El paquete demoR ayuda a presentar el código R en documentos Rmd. Por ejemplo, permite marcar/highligt algunas partes de una sentencia o código R. Para aprender puedes ir a la viñeta del paquete. Por ejemplo, en el próximo chunk voy a marcar en amarillo el operador pipe %>%
Podemos usar iconos en nuestros Rmds. Por ejemplo: + =
Podemos usar cajas de colores para resaltar un trozo de texto por ejemplo para poner conclusiones. Esto lo aprendí aquí. En este otro post nos enseñan a hacer esos cuadros realmente bonitos.
Conclusiones:
- This is my first conclusion - This is my second conclusion
si queremos que los chunks puedan copiarse en el portapapeles, podemos usar el paquete klippy
Hemos hablado muy poco de investigación reproducible pero, es evidente que usar documentos .Rmd y trabajar con Rprojects facilita la investigación reproducible. No la garantiza del todo, para ello habría que ver/usar más herramientas (docker, packrat, github, …). No lo vamos a hacer, pero al menos señalar dos prácticas que facilitan que tus análisis se acerquen a ser reproducibles:
Cuando estás haciendo un análisis y quieres compartirlo hay que estar seguro de los packages que se necesitan cargar para replicarlo, por eso es bueno hacer el análisis desde una sesión nueva/fresca de R y cargar los paquetes al principio del script. Para ello, puede ser de utilidad saber que paquetes tienes cargados en un momento dado, y puedes saberlo con: (.packages()).
A pesar de que yo no lo suelo hacer (👎), es recomendable introducir al final de los ficheros .Rmd la siguiente instrucción: sessionInfo(). De esta forma, proporcionarás información sobre que ordenador, sistema operativo y versión de R utilizaste en tu análisis, así como de las opciones locales de tu sistema (idioma, etc …) y de los paquetes que tienes cargados en memoria. Por ejemplo, ahora mismo mi sessionInfo es:
Otra forma de poner los detalles de la sesión es utilizar sessioninfo::session_info(); además lo usamos junto a details::details() que genera un desplegable para ver (o no) el outpt).