Prais: Prognose mit Prais-Winsten Transformation

Wie erweitere ich R um eigene Funktionen oder Pakete? Welches Paket ist passend für meine Fragestellung?

Moderatoren: EDi, jogo

Antworten
FranziskaBönte
Beiträge: 6
Registriert: Mi Feb 03, 2021 10:06 am

Prais: Prognose mit Prais-Winsten Transformation

Beitrag von FranziskaBönte »

Liebe Community,

wieder einmal ein kleines R Problem...

Und zwar schätze ich via library(prais) den Prais-Winsten Schätzer bei Autokorrelation.
Jetzt möchte ich gerne eine Prognose erstellen. allerdings sagt mir R: "Fehler in UseMethod("predict") :
nicht anwendbare Methode für 'predict' auf Objekt der Klasse "prais"

Wie schätze ich dann die zukünftigen Werte?

Hier ein Auszug aus meinen Daten und der Quellcode passend dazu

Code: Alles auswählen

library(prais)
library(sandwich)

y = c(2413.10, 2421.26, 2434.95, 2429.89, 2431.36, 2421.47, 2431.23, 2420.21 ,2402.80, 2330.75, 2269.42, 2258.62, 
      2177.25, 2141.12)
x1= c(177.0424, 174.2006, 178.3835, 168.5238, 168.4668, 170.4619, 176.6676, 165.3332, 160.8694, 166.1567, 168.7098,
      171.1672, 168.3668,182.1042)
x2 = c( 2839,  3249,  3929,  3725,  3204,  3061,  2485,  2073, 15202,  4074,  1996,  1887,   420,  1022)
Daten = data.frame(y, x1, x2)

Modell = lm(y ~ x1+x2)
summary(Modell)

dwtest(lm(Modell), alternative = c("greater","two.sided", "less"),iterations=15, tol=1e-10, data  =list())

Modell_pw = prais_winsten(lm(Modell), data=Daten)
summary(Modell_pw)
NewData = data.frame(y= 2500, x1= 190, x2= 1500)

#Prognose des "einfachen" Modells
predict(Modell, newdata = NewData)

#Prognose des transformierten Modells
predict(Modell_pw, newdata = NewData)
Über einen Tipp würde ich mich sehr freuen.
Franziska
Antworten