Este código de R corresponde a la figura 6.13 del libro Análisis de datos con el programa estadístico R: una introducción aplicada de Salas-Eljatib (2021).

Última actualización de esta web: 13 noviembre, 2022

Cargando los datos de ejemplo

Ocuparemos la dataframe annualppCities del paquete datana.

> library(datana)
> data(annualppCities)
> pp <- annualppCities
> tail(pp) #ultimas seis filas
> nrow(pp) #numero de filas de la dataframe
> #que ciudades existen en la dataframe?
> unique(pp$city)

Esta dataframe contiene la informacion de precipitaciones en algunas ciudades de Chile desde hace 70 anhos. Los datos han sido compilados por el Centro de Resilencia y del Clima (CR2).

El grafico

Una serie de tiempo muestra la trayectoria de una variable en el tiempo, y donde cada linea representa una misma unidad muestral.

Para realizar esta figura del libro es necesario cargar el paquete lattice, el cual esta instalado por defecto en R, por lo tanto basta con cargarlo a la presente sesion mediante

> require(lattice)

Ahora procedemos a realizar el grafico.

> #solo dejaremos a cuatro ciudades
> pp <- subset(pp, city!="Talca" & city!="Puerto Montt")
> unique(pp$city)
[1] "Chillan"  "Temuco"   "Valdivia" "Santiago"
> col.stgo <- "black";col.chillan <- "brown";
> col.temuco <- "green";col.valdi <- "blue";
> col.h <- c(col.stgo,col.chillan,col.temuco,col.valdi)
> lty.stgo <- 1;lty.chillan <- 3;lty.temuco <- 4;lty.valdi <- 5;
> lty.h <- c(lty.stgo,lty.chillan,lty.temuco,lty.valdi)
> #re-definiendo los niveles de los factores
> pp$city <- factor(pp$city,levels=c("Santiago","Chillan","Temuco","Valdivia"))
> #el grafico de interes
> xyplot(annual.pp ~ year, group=city, 
+        xlab=list(label="A\u00f1o", cex=1.0),
+        ylab=list(label="Precipitaci\u00f3n anual (mm)", cex=1.0),
+        data=pp, t="b", 
+        auto.key=list(space="top", columns=4, 
+                      lines=TRUE, points=FALSE))