Parcela Media Móvil En R
Tengo una trama de la serie de tiempo en el paquete ggplot2 y he realizado el promedio móvil y me gustaría añadir el resultado de la media móvil a la trama de series de tiempo. Ejemplo de conjunto de datos (p31): ambtemp dt -1.14 2007-09-29 00:01:57 -1.12 2007-09-29 00:03:57 -1.33 2007-09-29 00:05:57 -1.44 2007 -09-29 00:07:57 -1.54 2007-09-29 00:09:57 -1.29 2007-09-29 00:11:57 Código aplicado para la presentación de la serie de tiempo: Muestra del gráfico del promedio móvil Muestra de resultados esperados The El reto es que los datos de series de tiempo obtenidos de un conjunto de datos que incluye marcas de tiempo y temperatura, pero los datos del promedio móvil incluyen sólo la columna promedio y no las marcas de tiempo y el ajuste de estos dos puede causar inconsistencia. Moving Promedios en R A mi leal saber, No tienen una función incorporada para calcular promedios móviles. Usando la función de filtro, sin embargo, podemos escribir una función corta para medias móviles: Podemos utilizar la función en cualquier dato: mav (data), o mav (data, 11) si queremos especificar un número diferente de puntos de datos Que el predeterminado 5 trazado funciona como se espera: plot (mav (datos)). Además del número de puntos de datos sobre los cuales se puede hacer un promedio, también podemos cambiar el argumento de las funciones del filtro: sides2 usa ambos lados, sides1 usa sólo valores pasados. Comparte esta: Navegación de los artículos Comentario de navegación Comentario de navegaciónR Promedios móviles en ggplot2 Gabor Grothendieck Es probable que desee utilizar un paquete de series de tiempo para ello. Hay instalaciones de trazado específicamente dirigidas a series de tiempo en zoológico, xts, quantmod, timeSeries y latticeExtra. Ilustramos con el zoológico que tiene gráficos clásicos y métodos de gráficos de celosía: devAskNewPage (TRUE) biblioteca (zoo) set. seed (123) z lt - zoo (rnorm (100), Sys. Date () - 100: 0) plot (cbind (Z, rollmean (z, 10)), pantalla 1, col 1: 2) biblioteca (celosía) xyplot (cbind (z, rollmean (z, 10) Y el en Dic 10, 2009 a las 7:59 pm Probablemente quiera usar un paquete de series de tiempo para esto. Hay instalaciones de trazado específicamente dirigidas a series de tiempo en zoológico, xts, quantmod, timeSeries y latticeExtra. Ilustramos con el zoológico que tiene gráficos clásicos y métodos de gráficos de celosía: devAskNewPage (TRUE) biblioteca (zoo) set. seed (123) z lt - zoo (rnorm (100), Sys. Date () - 100: 0) plot (cbind (Z, rollmean (z, 10)), pantalla 1, col 1: 2) biblioteca (celosía) xyplot (cbind (z, rollmean (z, 10) Y el liso en diversos paneles omite la pantalla 1. Vea plot. zoo. Xyplot. zoo. Rollmean y las tres viñetas que vienen con el zoológico. El jue, 10 de diciembre 2009 a las 2:15 PM, fruminator escribió: Tenga algunos datos de la serie de tiempo almacenados en un data. frame, y estoy complotando con ggplot2 (que es totalmente impresionante). He explorado la documentación y los archivos de la lista de correo, y no puedo ver ninguna forma de trazar un 39smoother39 que es sólo el K-paso media móvil. Por ejemplo, imagino que tenía un data. frame llamado 39sleep39 con 39date39 como la fecha (de as. Date ()) y 39hours39 como las horas que dormí esa noche, me encantaría hacer algo como: qplot (fecha, horas, datos Sueño) statsmooth (método 39movingaverage39, k 7) ¿existe tal cosa existe. Si no, sé que el paquete es extensible, así que cualquier dirección en cuanto a cómo hacerlo para hacerlo sería muy apreciada. R-help at r-project. org lista de correo stat. ethz. ch/mailman/listinfo/r-help POR FAVOR, lea la guía de publicación R-project. org/posting-guide. html y proporcione comentarios, mínimos, autónomos , Código reproducible.
Comments
Post a Comment