Realización de boxplots de datos horarios en R

Question

Tengo mediciones de la temperatura del aire en dos lugares fijos medidos a intervalos de una hora. El código siguiente es un conjunto aleatorio de números utilizados para representar el formato de mis datos:

set.seed(1)
RandData <- rnorm(8760*2,sd=10)
Locations <- rep(c('UK','France'),each=8760)

Date = seq(from=as.POSIXct("1991-01-01 00:00"), 
              to=as.POSIXct("1991-12-31 23:00"), length=8760)

Final <- data.frame(Loc = Locations,
                    Doy = as.numeric(format(Date,format = "%j")),
                    Tod = as.numeric(format(Date,format = "%H")),
                    Temp = RandData)

Puedo trazar la variación de la temperatura en función del día del año con el siguiente código:

require(lattice)
xyplot(Temp~Doy | Loc, data = Final, col = "black", type = "l")

Esto mostraría el patrón anual de los datos. Sin embargo, lo que me gustaría hacer es producir boxplots de la variación de la temperatura para diferentes momentos del día. Así, para el ejemplo anterior me gustaría tener dos figuras, una para cada país y cada figura debería estar compuesta por gráficos de caja que muestren la variación de la temperatura a las 00:00, 01:00... y así sucesivamente, refiriéndose a Final$Tod. ¿Cómo se puede conseguir esto?

Muchas gracias por su ayuda.

Answer 1

¿Algo así?

library(ggplot2)
ggplot(Final, aes(x = as.factor(Tod), y = Temp)) + geom_boxplot()  + facet_wrap(~ Loc)

Answer 2

library(robustbase)
adjbox(Final$Temp[Final$Loc=="UK"]~Final$Tod[Final$Loc=="UK"])

Los boxplots son una herramienta de visualización, así que te daré un consejo visual. Lo que tienes es esencialmente datos funcionales por lo que usted por razones de visualización) utilizar una herramienta de diagrama de caja que reconozca eso. Pruebe el boxplot funcional en el paquete paquete fda.