J'ai un complot de série chronologique dans le paquet ggplot2 et j'ai effectué la moyenne mobile et je voudrais ajouter le résultat de la moyenne mobile à l'intrigue de la série chronologique. Exemple d'ensemble de données (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 Code appliqué pour la présentation de séries chronologiques: Échantillon de la courbe moyenne mobile Échantillon des résultats escomptés Le défi est que les données de séries chronologiques obtenues à partir d'un jeu de données qui comprend des horodatages et la température, mais les données de moyenne mobile incluent juste la colonne moyenne et non les horodatages et l'ajustement de ces deux peut entraîner des incohérences. N'ont pas une fonction intégrée pour calculer les moyennes mobiles. En utilisant la fonction de filtrage, on peut écrire une fonction courte pour les moyennes mobiles: On peut alors utiliser la fonction sur n'importe quelle donnée: mav (data), ou mav (data, 11) si on veut spécifier un nombre différent de points de données Que les 5 tracés par défaut comme prévu: plot (mav (data)). En plus du nombre de points de données sur lesquels la moyenne, nous pouvons également changer l'argument des côtés des fonctions de filtre: sides2 utilise les deux côtés, sides1 utilise des valeurs passées seulement. (S) (s) (s) Séries temporelles: Début 1 Fin 4 Fréquence 1 1 NA 4,333333 5,000000 NA Ici, j'essayais de faire une moyenne mobile qui prenait en compte Compte les 3 derniers chiffres, donc je m'attendais à obtenir seulement deux numéros de 8211 4,333333 et 5 8211 et si il y avait des valeurs NA, je pensais qu'ils étaient au début de la séquence. En fait, c'est ce que le paramètre 8216sides8217 contrôle: les côtés pour les filtres de convolution seulement. Si les côtés 1 les coefficients de filtre sont pour des valeurs passées seulement si les côtés 2 ils sont centrés autour du lag 0. Dans ce cas la longueur du filtre devrait être impair, mais si elle est pair, plus du filtre est vers l'avant dans le temps que vers l'arrière. Donc, dans notre fonction 8216mav8217, la moyenne mobile ressemble aux deux côtés de la valeur courante plutôt que juste aux valeurs passées. Nous pouvons tweak cela pour obtenir le comportement que nous voulons: gt bibliothèque (zoo) gt rollmean (c 4,5,4,6), 3) 1 4.333333 5.000000 J'ai également réalisé que je peux énumérer toutes les fonctions dans un paquet avec le 8216ls8217 Fonction donc I8217ll être balayage zoo8217s liste des fonctions la prochaine fois que j'ai besoin de faire quelque chose série de temps liés 8211 there8217ll probablement déjà être une fonction pour elle gt ls (quotpackage: zooquot) 1 quotas. Datequot quotas. Date. numericquot quotas. Date. tsquot 4 Quotas. Date. yearmonquot quotas. Date. yearqtrquot quotas. yearmonquot 7 quotas. yearmon. defaultquot quotas. yearqtrquot quotas. yearqtr. defaultquot 10 quotas. zooquot quotas. zoo. defaultquot quotas. zooregquot 13 quotas. zooreg. defaultquot quotautoplot. zooquot quotcbind. zooquot 16 quotcoredataquot quotcoredata. defaultquot quotcoredatalt-quot 19 quotfacetfreequot quotformat. yearqtrquot quotfortify. zooquot 22 quotfrequencylt-quot quotifelse. zooquot quotindexquot 25 quotindexlt-quot quotindex2charquot quotis. regularquot 28 quotis. zooquot quotmake. par. listquot quotMATCHquot 31 quotMATCH. defaultquot quotMATCH. timesquot Quotmedian. zooquot 34 quotmerge. zooquot quotna. aggregatequot quotna. aggregate. defaultquot 37 quotna. approxquot quotna. approx. defaultquot quotna. fillquot 40 quotna. fill. defaultquot quotna. locfquot quotna. locf. defaultquot 43 quotna. splinequot quotna. spline. defaultquot Quotna. StructTSquot 46 quotna. trimquot quotna. trim. defaultquot quotna. trim. tsquot 49 quotORDERquot quotORDER. defaultquot quotpanel. lines. itsquot 52 quotpanel. lines. tisquot quotpanel. lines. tsquot quotpanel. lines. zooquot 55 quotpanel. plot. customquotquotpanel. plot. defaultquot quotpanel. points. itsquot 58 quotpanel. points. tisquot quotpanel. points. tsquot quotpanel. points. zooquot 61 quotpanel. polygon. itsquot quotpanel. polygon. tisquot quotpanel. polygon. tsquot 64 quotpanel. polygon. zooquot quotpanel. rect Est. quelque. secondaire. secondaire. secondaire. secondes. secondes. secondes. secondes. secondes. secondes. quotpanel. text. tsquot 76 quotpanel. text. zooquot quotplot. zooquot quotquantile. zooquot 79 quotrbind. zooquot quotread. zooquot quotrev. zooquot 82 quotrollapplyquot quotrollapplyrquot quotrollmaxquot 85 quotrollmax. defaultquot quotrollmaxrquot quotrollmeanquot 88 quotrollmean. defaultquot quotrollmeanrquot quotrollmedianquot 91 quotrollmedian. defaultquot quotrollmedianrquot quotrollsumquot 94 quotrollsum. defaultquot quotrollsumrquot quotscalexyearmonquot 97 quotscalexyearqtrquot quotscaleyyearmonquot quotscaleyyearqtrquot 100 quotSys. yearmonquot quotSys. yearqtrquot quottimelt-quot 103 quotwrite. zooquot quotxblocksquot quotxblocks. defaultquot 106 quotxtfrm. zooquot quotyearmonquot quotyearmontransquot 109 quotyearqtrquot quotyearqtrtransquot quotzooquot 112 quotzooregquot Be Sociable, Share
No comments:
Post a Comment