Ich versuche die Funktion variogramST aus dem R-Paket gstat zu verwenden, um ein räumlich-zeitliches Variogramm zu berechnen.R gstat spatio-temporales Variogramm Kriging
Es gibt 12 Jahre Daten mit 20'000 Datenpunkten an unregelmäßigen Punkten in Raum und Zeit (kein vollständiges Gitter oder Teilgitter). Ich muss die STIDF aus dem Spacetime-Paket für einen unregelmäßigen Datensatz verwenden. Ich möchte ein temporales Semivariogramm mit Referenzpunkten bei 0, 90, 180, 270 Tagen, bis zu einigen Jahren usw. Leider treten sowohl Rechen- als auch Speicherprobleme auf. Wenn der Befehl
samplevariogram<-variogramST(formula=formula_gstat,data=STIDF1)
ohne weitere Argumente ausgeführt wird, wird die semiovariogram Berücksichtigung nur sehr kurze Zeiträume in Bezug auf die Referenzpunkte für die Semivariogramms, das scheint nicht in geeigneter Weise die inhärente Datenstruktur zu erfassen.
Es gibt mehr Argumente für diese Funktion zur Verfügung des Benutzers, aber ich bin nicht sicher, wie man sie richtig parametrisiert: tlag, tunit, twindow. Insbesondere frage ich mich, wie sie interagieren und wie ich mein oben beschriebenes Ziel erreiche. Also habe ich versucht, den folgenden Code
samplevariogram<-variogramST(formula=formula_gstat,data=STIDF1,tlag= ...., tunit=... , twindow= ...)
Die folgende Code Ergebnisse sind nicht durch Fragen Arbeitsspeicher in meinem 32Gbyte RAM Computer:
samplevariogram<-variogramST(formula=formula_gstat,data=STIDF1,tlag=90*(0:20), tunit="days")
aber könnte vielleicht fehlerhaft sein, sonst. Darüber hinaus scheint die letztere Codezeile auch hinsichtlich der Rechenzeit nicht durchführbar.
Kann jemand die VariogrammST-Funktion von der gstat-Verpackung korrekt angeben, mit dem Ziel auf die gewünschten Zeitintervalle?
Dank
Haben Sie versuchen, Alternativen mit einer kleineren (random) Teilmenge der Daten aus, in um ein Gefühl für das zu bekommen, was vor sich geht? Haben Sie versucht, Ihren Zeitschritt auf 90 Tage zu setzen, z. Ersetzen Sie die Zeit durch 1, 2, 3, usw., die Tag 0, 90, 180 usw. reflektieren? –
Wie richte ich die Zeitschritte (von 90 Tagen) korrekt ein? Wie wäre es mit "Beispielvariogramm <-VariogrammST (Formel = Formel_RENT_gstat, Daten = STIDF1, tunit =" Tage ", twindow = 90, tlags = 0: 6)"? – Tintin
Oder muss ich das POSIXct-Format verwenden, um den Zeitindex beim Erstellen des STIDF-Objekts aus dem Spacetime-Paket anzugeben? Meinst du, ich sollte meine Zeitdaten auf ein kleineres Intervall abbilden, für das ein passendes "tunit" -Argument zur Verfügung steht? – Tintin