Sesión 11 (C) — Del ajuste al informe: edición y exportación de gráficos y resultados en Gretl

Las prácticas A y B de esta sesión se centraron en el contenido econométrico (ortogonalidad, \(R^2\)). Esta práctica es distinta: no introduce ningún concepto nuevo de la asignatura, sino lo necesario para preparar un informe o trabajo escrito con resultados de Gretl: personalizar y exportar gráficos con calidad suficiente, y exportar tablas de resultados en un formato que se pueda incluir directamente en un documento (en vez de copiar números a mano o pegar capturas de pantalla). Reutilizamos el ajuste price~rooms ya estimado en la práctica B de esta sesión.

Objetivo

1. Personalizar un diagrama de dispersión (título, etiquetas de los ejes, colores).
2. Exportar gráficos en distintos formatos (PNG, PDF, EPS) y entender cuándo conviene cada uno.
3. Exportar los resultados de una regresión en un formato apto para un informe (LaTeX/RTF), en vez de copiarlos a mano.
4. Exportar directamente a fichero un histograma y un diagrama de caja (ya vistos en la práctica A de la sesión 4).
5. Fijar una breve lista de comprobaciones de buenas prácticas para preparar un informe reproducible.

Comandos nuevos de esta práctica. Además de ols (práctica A de esta sesión) y de gnuplot, freq y boxplot (sesiones 4 y 7):

list PL = price rooms

da un nombre a una lista de variables, para poder pasarla a otros comandos como un solo argumento.

plot ... end plot

un bloque, en lugar de una sola línea, para construir un gráfico con opciones y con instrucciones propias de gnuplot (el programa que Gretl usa para dibujar). Se explica en la Actividad 2.

--output"fichero.png"=

opción de gnuplot, boxplot y end plot (y --plot en freq) que guarda el gráfico en un fichero en vez de abrirlo en una ventana; la extensión decide el formato.

eqnprint y tabprint

escriben el último modelo estimado en un fichero, como ecuación o como tabla, listo para un documento.

Esta práctica es autocontenida: reabrimos hprice2 y volvemos a estimar el modelo price~rooms de la práctica B de esta sesión.

Abrimos hprice2, ponemos etiquetas a las dos variables que usaremos y estimamos price sobre la constante y rooms. Los comandos que siguen (eqnprint, tabprint) actúan sobre el último modelo estimado, así que conviene estimarlo justo antes de usarlos.

en línea de comandos:

open hprice2.gdt --quiet
setinfo price -d "Precio mediano de la vivienda ($)" -n "Precio ($)"
setinfo rooms -d "Número medio de habitaciones" -n "Habitaciones"

ols price const rooms

• GUI: genere el gráfico Gráfico X-Y de price sobre rooms (como en la práctica B de esta sesión). Sobre la ventana del gráfico, botón derecho, Editar: se abre un panel con pestañas para el título, las etiquetas de los ejes, el estilo de los puntos/líneas y los colores. Cambie el título a algo descriptivo (p. ej. ``Precio mediano frente a número de habitaciones'') y añada unidades explícitas en las etiquetas de los ejes (p. ej. ``Precio ($)'', ``Habitaciones (número medio)'').

  El equivalente en línea de comandos, línea a línea: list PL = price rooms define la lista de las dos variables, la última de ellas en el eje horizontal; plot PL abre el bloque; options fit=none pide el gráfico sin recta ajustada; cada línea literal pasa a gnuplot una orden tal cual (aquí, el título y las etiquetas de los dos ejes); y end plot --output"…"= cierra el bloque y guarda el gráfico en el fichero indicado.

  el equivalente en línea de comandos usa un bloque ``literal'' de gnuplot, entre llaves, para personalizar el gráfico:

  list PL = price rooms
  plot PL
      options fit=none
      literal set title "Precio mediano frente al numero de habitaciones"
      literal set xlabel "Habitaciones (numero medio)"
      literal set ylabel "Precio ($)"
  end plot --output="price_rooms_personalizado.png"
  

La estructura con el bloque plot es:

list MiLista = variable_y variable_x
plot MiLista
    options fit=none          # opciones de gretl para el gráfico (fit, with-lines, etc.)
    literal set title "..."   # cualquier comando gnuplot "en crudo"
    literal set xlabel "..."
    literal set ylabel "..."
end plot --output="nombre_fichero.png"   # o .pdf, .eps; el --output va en end plot

Tres reglas del bloque plot:

plot (a diferencia de gnuplot) solo admite un argumento: una lista, una serie o una matriz; de ahí que haya que definir la list antes. Dentro del bloque, cada línea empieza por una palabra clave reservada (options, literal, printf…). El destino del fichero (--output=) es una opción de end plot, igual que --output= es opción del comando gnuplot clásico.

• GUI: sobre la ventana del gráfico, botón derecho, Guardar como, y elija el formato (PNG, PDF, EPS, SVG según disponibilidad).

  o bien, en línea de comandos, repitiendo el mismo gráfico con distinta extensión. Son tres veces el mismo gnuplot; solo cambia la extensión del fichero de --output, y con ella el formato en que se guarda:

  gnuplot price rooms --output="price_rooms.png"
  gnuplot price rooms --output="price_rooms.pdf"
  gnuplot price rooms --output="price_rooms.eps"
  

Discusión: para un documento LaTeX, se prefiere normalmente un formato vectorial (PDF o EPS) frente a uno de mapa de bits (PNG): un PDF/EPS se puede escalar sin perder nitidez (líneas y texto siguen nítidos a cualquier tamaño), mientras que un PNG, al ampliarse, se pixela. El PNG es preferible sobre todo para páginas web o cuando el propio gráfico contiene un número muy grande de puntos (una nube de dispersión muy densa en PDF puede generar ficheros muy grandes).

• GUI: en la ventana del modelo estimado (price~rooms de la Actividad 1), menú LaTeX -> Tabular o LaTeX -> Ecuación, con las opciones de Ver, Copiar o Guardar. También, desde Archivo -> Guardar como, puede exportar en formato RTF (para Word) directamente.

  o bien, en línea de comandos, Gretl ofrece los comandos eqnprint (la ecuación estimada) y tabprint (una tabla completa de coeficientes y errores estándar) para exportar directamente a LaTeX. El bloque reestima el modelo con --quiet, para asegurarse de que es el último, y después escribe la ecuación en un fichero y la tabla en otro; el nombre del fichero va en --output y su extensión decide el formato (.tex para LaTeX, .rtf para un procesador de textos):

  ols price const rooms --quiet
  eqnprint --output="ecuacion_price_rooms.tex"
  tabprint --output="tabla_price_rooms.tex"
  

(Consulte help tabprint y help eqnprint en la consola de Gretl si las opciones exactas difieren en su versión. En tabprint el formato lo decide la extensión del fichero de --output: .tex para LaTeX, .rtf para RTF y .csv para valores separados por comas; eqnprint está documentado para LaTeX.)

Ventaja frente a copiar a mano: estos comandos generan directamente el código (LaTeX o RTF) listo para pegar en el informe, con los números exactos que calculó Gretl. Se evitan errores de transcripción y el trabajo de dar formato a una tabla a mano.

Esto es lo que ha quedado en ecuacion_price_rooms.tex y que puede pegar, tal cual, en su documento:

%%% the following needs the amsmath LaTeX package

\begin{gather}
\widehat{\rm price} = 
-\underset{(2651,5)}{34796,2}
+\underset{(419,34)}{9119,55}\,\mbox{rooms}
 \notag \\
T = 506 \quad \bar{R}^2 = 0,4831 \quad F(1,504) = 472,95 \quad \hat{\sigma} = 6620,9 \notag \\
\centerline{(Desviaciones típicas entre paréntesis)} \notag
\end{gather}

El fichero tabla_price_rooms.tex contiene, de forma análoga, el entorno tabular completo con los coeficientes, sus desviaciones típicas y el resto de columnas. Nótese que ambos ficheros contienen los números con todas sus cifras y con la coma decimal: no es necesario transcribir nada.

Recordando la práctica A de la sesión 4 (histograma y diagrama de caja de price). Son los mismos freq y boxplot de entonces; la diferencia es que, en vez de abrir una ventana, guardan el gráfico en el fichero indicado, freq con la opción --plot y boxplot con --output:

en línea de comandos:

freq price --plot="hist_price_informe.png"
boxplot price --output="boxplot_price_informe.png"

• Reproducibilidad: conserve el guion completo (.inp) usado para generar cada gráfico y cada tabla, no solo el resultado final; si usa datos simulados, fije siempre set seed.
• Formato vectorial para figuras en LaTeX: prefiera PDF/EPS a PNG cuando el documento final sea un PDF compilado con LaTeX.
• Ejes con unidades: toda figura debe indicar en sus ejes qué se mide y en qué unidades (dólares, número de habitaciones, etc.), no solo el nombre abreviado de la variable en Gretl.
• Tablas generadas, no transcritas a mano: use tabprint/eqnprint (o la vía GUI equivalente) en lugar de copiar números de la pantalla, para evitar errores de transcripción.
• Cite la fuente de los datos: indique siempre el origen del dataset (aquí, hprice2, Harrison y Rubinfeld, 1978).

1. ¿Por qué para un documento LaTeX suele preferirse exportar un gráfico en PDF o EPS en lugar de PNG?
2. ¿Qué ventaja tiene usar tabprint / eqnprint frente a copiar manualmente los números de la salida de Gretl a un documento?
3. ¿Por qué es importante incluir las unidades en las etiquetas de los ejes de un gráfico que se va a incluir en un informe?
4. ¿Por qué conviene conservar el guion completo (.inp) de Gretl, y no solo los gráficos y resultados finales, al preparar un informe reproducible?
5. Si tuviera que incluir en un informe la comparación entre los modelos price~rooms y price~nox de la práctica B de esta sesión, ¿qué elementos (gráficos, tablas, estadísticos) incluiría, y por qué?

Para profundizar:

• Cottrell, A. y Lucchetti, R. (2023). Gretl User's Guide. Secciones sobre el comando gnuplot (incluidos los bloques ``literal'' de personalización) y sobre tabprint/eqnprint para exportar resultados.

1. PDF/EPS frente a PNG en LaTeX. Un formato vectorial (PDF, EPS) describe el gráfico mediante instrucciones geométricas (líneas, curvas, texto), de modo que se puede escalar a cualquier tamaño sin perder nitidez; un PNG es un mapa de bits (una cuadrícula de píxeles de tamaño fijo), y al ampliarlo dentro del documento (o al imprimirlo a alta resolución) se pixela y pierde calidad, especialmente notorio en el texto de los ejes.
2. Ventaja de tabprint/eqnprint frente a copiar a mano. Copiar números manualmente introduce el riesgo de errores de transcripción (un dígito mal copiado, un signo perdido) y obliga a dar formato a la tabla desde cero cada vez que cambia la regresión. Generar la tabla directamente desde Gretl garantiza que los números que aparecen en el informe coinciden con los calculados, y que si se vuelve a estimar el modelo con datos actualizados basta con volver a ejecutar el mismo comando para regenerar la tabla, sin repetir ningún paso manual.
3. Unidades en las etiquetas de los ejes. Sin unidades explícitas, un lector no puede saber si ``price'' está en dólares, en miles de dólares, o en otra escala, ni si ``rooms'' es un conteo o un promedio; esto puede llevar a interpretaciones erróneas de la magnitud de los efectos mostrados en el gráfico. Las etiquetas con unidades hacen el gráfico autosuficiente, comprensible sin necesidad de leer el texto que lo acompaña.
4. Conservar el guion completo (.inp). Un informe reproducible debe permitir que otra persona (o uno mismo, meses después) obtenga los mismos resultados a partir de los mismos datos, sin depender de pasos manuales no documentados (clics de menú, ajustes manuales no registrados en una ventana gráfica). El guion completo documenta, de forma exacta y ejecutable, cada transformación de los datos y cada estimación realizada; los gráficos y tablas finales son solo su consecuencia, no una fuente independiente de verdad.
5. Elementos para comparar price~rooms y price~nox. Como mínimo: (i) una tabla con \(\hat\beta_1,\hat\beta_2,R^2\) de ambos modelos (generada con tabprint, no copiada a mano); (ii) los dos diagramas de dispersión con la recta ajustada, con ejes debidamente etiquetados y unidades explícitas, para que el lector aprecie visualmente la diferencia de ajuste; (iii) una breve discusión textual de los signos obtenidos y de por qué un \(R^2\) mayor no debe leerse como ``una relación causal más fuerte'', sino solo como ``una nube de puntos más alineada con la recta'' (lección 8). No haría falta incluir capturas de pantalla de la consola de Gretl.