Hallo ok ok verstehe, also ich probiere es nochmal.
Ich habe eine Datenreihe von zwei Jahren (732) Werte, hier liegt wenn ich es richtig verstanden habe bereits das erste Problem, weil ich nur eine y Koordinate habe. Also wäre meine erste Frage, wie weise ich direkt beim einlesen der Datei dem jeweiligen Wert 1:732 zu?
Code: Alles auswählen
> str(q1)
num [1:732] 0.12 0.53 0.49 0.4 0.3 0.43 0.46 0.38 0.27 0.23 ...
Nach dem Import habe ich direkt den Plot erstellt und ein paar optische Anpassungen durchgeführt hierbei scheint es aber zu funktionieren, weil mein Plot der Ganglinie richtig aussieht.
Im nächsten Schritt habe ich die 95%quantil berechnet, die in meinem Fall etwa 0.7 beträgt. Dann wollte ich die Werte im Plot markieren die einen Peak darstellen, der größer als die 0.7 ist. Mein erste Idee war über die erste Ableitung, aber hierfür habe ich für meine Zeitreihe keine passende Funktion in R gefunden. Also habe ich mir einfach die Werte die größer als 0.7 sind ausgeben lassen. Hier tritt aber das nächste Problem auf, nämlich ist nicht jeder Wert der größer als 0.7 ist ein Peak. Hier also die nächste Frage gibt es eine Möglichkeit bzw. Eine Funktion mit der ich die erwünschten Hochpunkte (y>0.7) auswählen kann?
Hier der Code soweit ich mir sicher bin, das ich weiß was ich tue
Code: Alles auswählen
q1 <- scan("file1.txt", skip=20)
q1
q1976=q1[1:366]
q1976
q1977=q1[367:732]
q1977
plot(1:366, q1976[1:366], type="l", xlab="Day of the year (1976)", ylab="Discharge (m³/s)", main="Gauge Ibm / Moosache")
quantile(q1976, probs=0.95)
abline(h=0.77)
Die dritte und letzte Frage ist jetztwie ich diese dann am besten markiere mit der Points funktion ist es mir nicht gelungen den y Koordinaten die richtige Koordinaten zuzuordnen. Beziehungsweise habe ich eine Fehlermeldung bekommen, ich hoffe das hilft jetzt beim Verständnis von meinem Problem. Vielen Dank für die Zeit, ich habe jetzt mal meine Daten im Anhang dazugegeben, ich hoffe das hilft!
Liebe Grüße!
e: hab jetzt eine Lösungsmöglichkeit gefunden, um alle Punkte größer 0.7 zu markieren, allerdings mit einer if Schleife also falls jemand noch andere Ideen hat, bitte her damit