Zeitreihe
Verfasst: Mo Nov 21, 2016 2:12 pm
Liebes Forum,
schade um das alte Forum, aber hier sind neue Daten und Probleme .
Ich habe Pegelmessungen eines Sees über mehrere Jahre und möchte jahreszeitliche Schwankungen des Wasserstandes betrachten.
1. Wie erstelle ich in R eine Zeitreihe mit den festen Zeitpunkten der Messungen (Messungen erfolgten ca. einmal je Woche, aber unregelmäßig und mit Ausfällen)? Ich erhoffe mir über die Zeitreihe vergleichbare (interpolierte Werte?) zu festen Zeitpunkten sowie über den Trend das Fallen und das Steigen des Wasserstandes anzeigen lassen zu können.
2. Wie kann ich mir Min, Max und Mittelwerte für bestimmte Zeitfenster berechnen, z.B. Min, Max, Mean von Jan bis Mrz 2011 bzw. die Werte aller Jahre?
3. Gibt es eine elegantere Lösung für das Filtern nach Jahren/ erstellen eines Vektors je Jahr?
4. Warum erscheint bei mir kein Plot nach dem Befehl lines(pegel....)?
Fragen über Fragen....
Danke schon mal für Eure Hilfe.
Viele Grüße
Das ist der Code:
schade um das alte Forum, aber hier sind neue Daten und Probleme .
Ich habe Pegelmessungen eines Sees über mehrere Jahre und möchte jahreszeitliche Schwankungen des Wasserstandes betrachten.
1. Wie erstelle ich in R eine Zeitreihe mit den festen Zeitpunkten der Messungen (Messungen erfolgten ca. einmal je Woche, aber unregelmäßig und mit Ausfällen)? Ich erhoffe mir über die Zeitreihe vergleichbare (interpolierte Werte?) zu festen Zeitpunkten sowie über den Trend das Fallen und das Steigen des Wasserstandes anzeigen lassen zu können.
2. Wie kann ich mir Min, Max und Mittelwerte für bestimmte Zeitfenster berechnen, z.B. Min, Max, Mean von Jan bis Mrz 2011 bzw. die Werte aller Jahre?
3. Gibt es eine elegantere Lösung für das Filtern nach Jahren/ erstellen eines Vektors je Jahr?
4. Warum erscheint bei mir kein Plot nach dem Befehl lines(pegel....)?
Fragen über Fragen....
Danke schon mal für Eure Hilfe.
Viele Grüße
Das ist der Code:
Code: Alles auswählen
pegel.data <- read.table("Pegel.txt", header=TRUE,
sep="\t", na.strings="NA", dec=".", strip.white=TRUE)
pegel.data$Datum <- as.POSIXct(pegel.data$Datum, format = '%d.%m.%Y')
str(pegel.data)
plot(pegel.data$Datum,pegel.data$WSP, type="b")
pegel2008 <- pegel.data[pegel.data$Datum > as.POSIXct("2008-01-01")
& pegel.data$Datum < as.POSIXct("2008-12-31"),]
pegel2009 <- pegel.data[pegel.data$Datum > as.POSIXct("2009-01-01")
& pegel.data$Datum < as.POSIXct("2009-12-31"),]
plot(pegel2008$Datum,pegel2008$WSP, type="l", col="blue",ylim=c(24,26))
lines(pegel2009$Datum,pegel2009$WSP, type="l", col="blue",ylim=c(24,26))