Calculer itérativement le coefficient de dérive à partir de la fonction de dérive aléatoire dans R, compiler dans la liste

Mon objectif est de répertorier le coefficient de dérive d'une marche aléatoire avec fonction de prévision de dérive, appliqué à un ensemble de données historiques (ci-dessous). Spécifiquement j'essaie de rassembler le coefficient de dérive à partir de la marche aléatoire avec le modèle de dérive de la première année, puis cumulativement à la dernière, en enregistrant le coefficient chaque fois, signifiant itérativement ou chaque année supplémentaire (l'enregistrant dans une liste? approprié). Pour être clair, chaque nouvelle prévision de marche aléatoire inclut toutes les années précédentes.Calculer itérativement le coefficient de dérive à partir de la fonction de dérive aléatoire dans R, compiler dans la liste

Les données sont une liste de 241 niveaux de consommation, et je tente de discerner comment le coefficent de dérive changerait au cours de itérativement progressant de n = 1 à n = 241

Où par exemple la marche aléatoire avec le modèle de dérive est Y [t] = c + Y [t-1] + Z [t] où Z [t] est une erreur normale et c est le coefficient que je cherche. Mes tentatives actuelles de cette impliquent pour la fonction de boucle et extraire le coefficient c de la fonction RWF() du paquet « Prévisions » dans R.

Pour extraire, je fais en tant que telle

rwf(x, h = 1, drift = TRUE)$model[[1]]

qui extrait le coefficient de dérive. Le problème est que mes tentatives de sous-ensemble des données dans l'appel rwf ont échoué, et je ne crois pas, par essais et erreurs, que rwf() prenne en charge l'argument du sous-ensemble, comme le fait par exemple un modèle lm. Dans ce sens, mes tentatives de boucler la fonction ont également échoué.

Un exemple de code tel est

for (i in 1:5){print((rwf(x[1:i], h = 1, drift = TRUE))$model[[1]])}

qui me donne l'erreur suivante

Error in lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) : 
    0 (non-NA) cases 
In addition: Warning message: 
In is.na(rows) : is.na() applied to non-(list or vector) of type 'NULL'

Toute aide serait grandement appréciée. J'ai beaucoup lu pour de l'aide mais c'est la première fois que je pose une question.

Les données sont les suivantes

Source

2012-04-30 BryanC

FYI c'est 'forecast' avec un minuscule' f' – Chase

Vous avez besoin d'au moins deux points pour l'adapter à votre modèle. Voici comment j'aborder le problème après avoir lu vos données dans un data.frame nommé x:

library(forecast) 
drifts <- sapply(2:nrow(x), function(zz) rwf(x[1:zz,], drift = TRUE)$model$drift)

Je ne sais pas si c'est ce que vous attendiez ou non, mais voici une parcelle de vos valeurs de dérive : enter image description here

Source

2012-05-01 00:12:33 Chase

Calculer itérativement le coefficient de dérive à partir de la fonction de dérive aléatoire dans R, compiler dans la liste

Répondre

Questions connexes