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).

Cargando los datos de ejemplo

Ocuparemos la dataframe annualppCities del paquete datana.

library(datana)
data(annualppCities)
pp <- annualppCities
tail(pp) #ultimas seis filas
            city year annual.pp
335      Chillan 2019     628.4
336 Puerto Montt 2019    1278.4
337     Santiago 2019      47.2
338        Talca 2019     320.0
339       Temuco 2019     802.3
340     Valdivia 2019    1071.3
nrow(pp) #numero de filas de la dataframe
[1] 340
#que ciudades existen en la dataframe?
unique(pp$city)
[1] "Chillan"      "Temuco"       "Valdivia"     "Santiago"     "Puerto Montt"
[6] "Talca"       

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)
Loading required package: lattice

Attaching package: 'lattice'
The following object is masked _by_ '.GlobalEnv':

    parallel

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)

pp$city <- factor(pp$city,levels=c("Santiago","Chillan","Temuco","Valdivia"))

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))