ich habe gelesen, dass wen man erahnt oder testet, dass ein OLS-Modell nicht richtig gut passt, z.B. wegen Heteroskedastizität (oder auch Autokorrelation?), dann kann man es mit einem GLS-Modell probieren. Nachdem mir dieser Befehl (gls() aus dem package nlme) schon einige Nerven gekostet hat (die Residuen sind keine Zeitreihe und man kann es nicht so wie ein lm-Modell behandeln), wollte ich einen BP-Test machen (also Breusch–Pagan-Test) – um wieder auf Heteroskedastizität zu testen. Ich finde leider gar keinen Ansatz mein GLS-Modell mit einem solchen zu testen. Hat da jemand Erfahrung? Das muss doch gehen?
Danke im Voraus.
Hier noch meine Fehlermeldung für ein Modell GLS_0:
Code: Alles auswählen
> bptest(GLS_0)
Error in as.data.frame.default(data) :
cannot coerce class "c("glsStruct", "modelStruct")" to a data.frame