Pregunta

Tengo un zoológico con datos diarios que se ven así:> head(almorol)
1973-10-02 1973-10-03 1973-10-04 1973-10-05 1973-10-06 1973-10-07
183.9 208.2 153.7 84.8 52.5 35.5

Quiero trazar totales anuales y un promedio móvil, así que lo hice:

y<-apply.yearly(almorol, FUN=sum)
plot(y, main="Annual totals - Tagus (Almorol)",ylab="Q (m3/s)")
lines(rollapply(y, 10, mean, na.rm=TRUE), col="red", lwd=2)

Lo que funciona bien, pero debido a que los datos no son continuos, un gráfico de línea no es correcto. Si lo hago con puntos es difícil de seguir, así que quería hacer una placa. Sin embargo cuando lo hago barplot(y) El X AX aparece en el formato YY-MM-DD en lugar de solo años y no puedo agregar una línea con el promedio móvil (no da un error que simplemente no aparece).

¡Gracias por ayudar!

editar: > y
1973-12-31 1974-12-31 1975-12-31 1976-12-31
19224.20 103766.30 72180.90 55939.80
1977-12-31 1978-12-31 1979-12-31 1980-12-31
215905.00 231014.21 319481.02 58979.84
1981-12-31 1982-12-31 1983-12-31 1984-12-31
32931.17 67989.06 83920.62 99431.75
1985-12-31 1986-12-31 1987-12-31 1988-12-31
161357.10 82910.87 101154.81 147541.80
1989-12-31 1990-12-31 1991-12-31 1992-12-31
137684.21 134974.39 89039.02 20774.72
1993-12-31 1994-12-31 1995-12-31 1996-12-31
58916.95 55187.38 52361.77 271064.34
1997-12-31 1998-12-31 1999-12-31 2000-12-31
220510.88 125116.62 42170.95 103915.99
2001-12-31 2002-12-31 2003-12-31 2004-12-31
244513.59 66811.15 167458.93 67223.66
2005-12-31 2006-12-31 2007-12-31 2008-12-31
8906.31 116874.33 79105.30 40142.61
2009-12-31 2010-12-31 2011-05-31
43835.34 161491.77 73093.06

¿Fue útil?

Solución

Las coordenadas X en una placa no están relacionadas con la escala original de los datos. Puede ver el valor de retorno de la función de la placa de barras para ver cuál es el rango del eje x (la razón probable para no ver la línea es trazada completamente fuera de la pantalla).

La función UpdateUSR en el paquete TeachingDemos se puede usar para cambiar la escala de coordenadas para que coincida con los datos que desea agregar.

Otra alternativa es usar la función de trazado con type='h' (y mira lwd y lend opciones en ?par) Para crear su propia trama como la trama usando las coordenadas de interés, luego agregue su línea de referencia.

Licenciado bajo: CC-BY-SA con atribución
No afiliado a StackOverflow
scroll top