head

‎"oɯsıɯ ol ǝɹdɯǝıs sɐƃɐɥ ou ,soʇuıʇsıp sopɐʇlnsǝɹ sɐɔsnq ıS"

Estadísticos descriptivos en R: Cálculo y Exportación a LaTex

Como es muy sabido, todo trabajo empírico lleva como parte inherente un resumen de estadísticos descriptivos de los datos que se utilizarán en el modelo para la investigación, por lo cual es importante conocer cuál es la función de R que nos proporciona tal información. Así mismo, un buen trabajo culmina con una buena presentación en un informe final o un paper, con lo cual se deberá utilizar , yo personalmente utilizo Scientific Workplace (SW) en lugar de , dado que SW es mucho más amigable para utilizarlo, se puede obtener un versión de SW en esta web. Una completa descripción sobre SW se encuentra en su página oficial, también se pueden conseguir los mismos resultados de compilación en si se utiliza el editor MikTex cuyas instrucciones de instalación están aquí. Una vez entendido esto pasamos a lo práctico con R.

Primero se tiene que descargar el paquete que contiene  la función basicStats() que es la que utilizaremos, el paquete se llama fBasics.

Instrucciones para descargar y activar paquetes en R

Una vez que se ha iniciado el programa nos dirigiremos sobre esta ruta sobre el menú principal
Paquetes → Instalar Paquetes → Seleccionar CRAN mirror → Ubicar paquete por su nombre (fBasics) y clic en OK

El CRAN mirror que yo utilizo es USA(TX1). Luego para cargar (activar) el paquete se escribe el comando library(nombre del paquete), en nuestro caso, library(fBasics)

O si se prefiere se pude hacer directamente con uso de comandos así:

Seleccionado el CRAN mirror con chooseCRANmirror() y luego buscar el paquete con utils:::menuInstallPkgs() y buscar el paquete en la larga lista que parece, pero como ya sabemos el nombre del paquete que usaremos, es más viable utilizar directamente la función install.packages(“fBasics”) y ya lo tendremos, luego para cargarlo lo haremos usando library(fBasics). Más paquetes y su correspondiente descripción aquí.

Una vez descargado y activado el paquete podemos hacer el cálculo de los estadísticos descriptivos, para ello utilizaremos los datos de la base Canada alojada en el paquete vars, para acceder a esos datos no es necesario descargar el paquete, accederemos directamente a ella con la función data("nombre de la base", package="paquete donde está alojada").


Estadísticos descriptivos con R

La función para esto, se ha dicho anteriormente, es basicStats(), así pues el resultado de aplicarla a la base Canada es
> data("Canada", package="vars")
> basicStats(Canada)
                       e         prod           rw          U
nobs           84.000000    84.000000    84.000000  84.000000
NAs             0.000000     0.000000     0.000000   0.000000
Minimum       928.563336   401.306688   386.136109   6.700000
Maximum       961.765710   418.002980   470.011666  12.770000
1. Quartile   935.353242   404.750046   423.898467   7.782500
3. Quartile   950.027595   410.699869   461.095270  10.607500
Mean          944.257255   407.820911   440.750835   9.320833
Median        945.951797   406.458163   444.441407   9.450000
Sum         79317.609422 34256.956522 37023.070173 782.950000
SE Mean         0.999189     0.460056     2.539019   0.175361
LCL Mean      942.269909   406.905879   435.700829   8.972047
UCL Mean      946.244601   408.735943   445.800842   9.669619
Variance       83.863773    17.778692   541.515791   2.583123
Stdev           9.157717     4.216479    23.270492   1.607210
Skewness       -0.075910     0.625394    -0.557143   0.177289
Kurtosis       -0.979191    -0.558561    -0.810872  -1.123477


Para una descripción detallada de las variables de esta base vea esto en un extracto del documento o en el documento completo.

Exportando a

Para exportar la tabla de estadísticos descriptivos de R a es necesario instalar y cargar el paquete Hmisc, el procedimiento es exactamente igual al realizado anteriormente, por las dudas doy los comandos:
install.packages(“Hmisc”)
library(Hmisc)


La función latex() es la que nos convertirá la tabla de estadísticos descriptivos a códigos y la exportará directamente al directorio en el que se esté utilizando R, por defecto este es la carpeta donde está instalado en Archivos de Programa, yo recomiendo crear una carpeta en Mis Documentos y cambiar el directorio de trabajo para que las exportaciones se guarden allí, esto se hace sencillamente yendo a Archivo en el menú principal de R y seleccionar Cambiar dir… y allí escogemos el que deseamos y ya sabremos que en ese directorio encontraremos el resultado de la exportación.
Volvemos a calcular los descriptivos y los asignaremos al nombre “descriptivos” esto con el objetivo de asignarle un nombre al archivo que se exportará y también será el nombre de la tabla, los resultados son:

descriptivos<-basicStats(Canada)
latex(descritivos)

El correspondiente código es generado por R esta aquí

Tal y como se aprecia en los códigos generados, la Tabla se ha nombrado como “descriptivos” y también ese será el nombre de la primera columna, veamos mejor como queda después de compilarlo en pdf usando SW o en MikTex. La compilación la he hecho usando SW.

descriptivos

Esta función también nos permite controlar en número de decimales que deseamos que aparezcan en la tabla final. Para obtener la misma tabla anterior pero esta vez con 3 decimales le indicamos a R la siguiente opción dentro la función
latex(x, digits= Número de dígitos deseados)

 
latex(descriptivos, digits=3)  

y una vez compilado nos queda:
descriptivos2
Esta tiene una mejor presentación.

Usando Scientific Workplace para compilar en pdf los códigos producidos por R

Para convertir esos horribles códigos de a verdaderas y bonitas tablas usamos SW. Primero, ubicamos el archivo creado por R con terminación tex, lo abrimos en un blog de notas y seleccionamos todo su contenido y lo copiamos, luego vamos a SW y seguimos la siguiente ruta Insert → Typeset Object → TeX Field… y en el cuadro que abre pegamos el código, luego vamos al menú principal y hacemos Typeset → Preview PDF… ya lo tendremos, el documento pdf compilado se abre por defecto si no pues lo tendremos que buscar en la carpeta donde hayamos guardado el archivo de SW.

Otra forma de obtener resultados en R y exportarlos a LATEX se puede hacer con la función xtable en el paquete xtable, para ver ejemplos de esto en este post

Espero que haya sido útil. Espero que dejes tu comentario.
Jilber Urbina.

1 comentario:

  1. Buen aporte, os habéis dado al clavo con lo que buscaba y por supuesto también os habéis ganado un seguidor. Sugiero revisar el siguiente vídeo "capacidad de los procesos" https://www.youtube.com/watch?v=BNFipw7lG0I personalmente me sirvió y mucho.

    ResponderEliminar