Verfügbare Informationen zu "Zettel 5"
Qualität des Beitrags: Beteiligte Poster: Nils Forum: das Problem R aus dem Unterforum: Vergleichen Antworten: 1 Forum gestartet am: Freitag 18.05.2007 Sprache: deutsch Link zum Originaltopic: Zettel 5 Letzte Antwort: vor 16 Jahren, 11 Monaten, 2 Tagen, 10 Stunden, 39 Minuten
Alle Beiträge und Antworten zu "Zettel 5"
Re: Zettel 5
Nils - 24.05.2007, 19:14Zettel 5
also ich habs so:
Code: #-Aufgabe 22. a)
rm(list=ls(all=TRUE))
#"windows(record=T)" wird verwendet, da diese Aufgabe
#mehrere Grafiken beinhaltet und man so mittels
#PgUp und PgDn zwischen diesen blättern kann.
windows(record=T)
#Die Datei "peake.csv" wird eingelesen, als "Daten"
#gespeichert und durch attach werden die Werte
#verfügbar gemacht.
Daten=read.csv("peake.csv",header=TRUE)
attach(Daten)
#Die Werte von AREA werden in dem Vektor x und die
#Werte von SPECIES im Vektor y gespeichert.
#Anschließend wird das Streuungsdiagramm der beiden
#Vektoren geplotet.
x=AREA
y=SPECIES
plot(x,y,las=1,main="Streuungsdiagramm von Area und Spezies",
xlab="Area",ylab="Spezies")
#Nun wird das Diagramm noch einmal geplotet, allerdings
#werden diesem dieses Mal 5 "lowess"-Glättungskurven
#mit verschiedenen Glättungsparametern hinzugefügt.
plot(x,y,las=1,
main="Glättung eines Streuungsdiagrammes \n mittels lowess",
xlab="Area",ylab="Spezies")
lines(lowess(x,y,f=1/8),col="red")
lines(lowess(x,y,f=1/4),col="green")
lines(lowess(x,y,f=1/3),col="yellow")
lines(lowess(x,y,f=1/2),col="skyblue")
lines(lowess(x,y,f=1),col="purple")
#Dann wird eine Legende für die Graphen erstellt.
#"bottomright" ist die Position, also unten rechts,
#die Klammer hinter c beinhaltet die Parameter,
#wie sie in der Legende ausgegeben werden sollen,
#die Klammer hinter dem col=c gibt die Farben der
#Markierungen in der Legende an, pch bestimmt das
#Aussehen besagter Markierungen, mit title wird die
#Überschrift bestimmt und mit bg die Hintergrundfarbe.
legend("bottomright",c("1/8","1/4","1/3","1/2","1"),
col=c("red","green","yellow","skyblue","purple"),
pch=c(16,16,16,16,16),title="Glättungsparameter:",bg="#EEEEEE")
#Das Streuungsdiagramm wird wieder geplotet und es werden
#5 "smooth.spline"-Glättungskurven mit verschiedenen
#Glättungsparametern hinzugefügt.
plot(x,y,las=1,
main="Glättung eines Streuungsdiagrammes \n mittels smooth.spline",
xlab="Area",ylab="Spezies")
lines(smooth.spline(x,y,df=2),col="red")
lines(smooth.spline(x,y,df=4),col="green")
lines(smooth.spline(x,y,df=9),col="yellow")
lines(smooth.spline(x,y,df=17),col="skyblue")
lines(smooth.spline(x,y,df=25),col="purple")
#Die Legende funktioniert wie im "lowess"-Fall.
legend("bottomright",c("2","4","9","17","25"),
col=c("red","green","yellow","skyblue","purple"),
pch=c(16,16,16,16,16),title="Glättungsparameter:",bg="#EEEEEE")
#Ein letztes Mal wird das Streuungsdiagramm geplotet,
#dieses Mal mit 5 "loess-smooth"-Glättungskurven mit
#verschiedenen Glättungsparametern
plot(x,y,las=1,
main="Glättung eines Streuungsdiagrammes \n mittels loess.smooth",
xlab="Area",ylab="Spezies")
lines(loess.smooth(x,y,span=1/4),col="red")
lines(loess.smooth(x,y,span=1/3),col="green")
lines(loess.smooth(x,y,span=1/2),col="yellow")
lines(loess.smooth(x,y,span=1),col="skyblue")
lines(loess.smooth(x,y,span=4),col="purple")
#Die Legende funktioniert wie in den beiden vorigen
#Fällen.
legend("bottomright",c("1/4","1/3","1/2","1","4"),
col=c("red","green","yellow","skyblue","purple"),
pch=c(16,16,16,16,16),title="Glättungsparameter:",bg="#EEEEEE")
Code: #-Aufgabe 22. b)
rm(list=ls(all=TRUE))
#"windows(record=T)" wird verwendet, da diese Aufgabe
#mehrere Grafiken beinhaltet und man so mittels
#PgUp und PgDn zwischen diesen blättern kann.
windows(record=T)
#Die Excel-Datei "hunte.xls" wurde auf Fehler hin
#untersucht und im *.csv-Format gespeichert.
#Diese *.csv-Datei wird nun eingelesen, als
#"Daten" gespeichert und durch attach werden die Werte
#verfügbar gemacht.
Daten=read.csv2("hunte.csv",header=TRUE)
Daten
#Es werden Zeitreihen für Niedrigwasser (nw) und
#Hochwasser (hw) erstellt. Der Beginn liegt bei 1996
#und die Abstände sollen 12 Monate betragen. Diese
#Informationen genügen R, um beim ploten dieser
#Zeitreihen die Achsen richtig zu beschriften.
nw=ts(Daten$Niedrigwasser,start=1996,freq=12)
hw=ts(Daten$Hochwasser,start=1996,freq=12)
#Mittels ts.plot werden nun beide Zeitreihen in einem
#Bild geplotet.
ts.plot(nw,hw,gpars=list(xlab="Jahr",ylab="mittlere Tide in cm",las=1,
col=c("red","green"),main="Hoch- und Niedrigwasser der Hunte"))
#Dann wird eine Legende für die Graphen erstellt.
#"right" ist die Position, also rechts,
#c("Hochwasser","Niedrigwasser")sind die in der Legende beschrieben
#Zeitreihen, col=c("green","red") die Farbe der beiden
#Markierungen in der Legende, pch=c(16,16) bestimmt das
#Aussehen besagter Markierungen, mit title wird die
#Überschrift bestimmt und mit bg die Hintergrundfarbe.
legend("right",c("Hochwasser","Niedrigwasser"),
col=c("green","red"),pch=c(16,16),title="Legende:",bg="#EEEEEE")
#Die gleitenden Durchschnitte der Ordnungen 3 und 12 werden
#mittels filter berechnet und in den Variablen nw3 und nw12
#gespeichert.
nw3 = filter(nw, rep(1/3, 3))
nw12 = filter(nw, rep(1/12, 12))
#Die Zeitreihe des Niedrigwassers wird geplotet und dieses noch
#die beiden zuvor gespeicherten gleitenden Durchschnitte
#hinzugefügt. Zudem wird eine Legende erstellt.
plot(nw,col="blue",las=1,xlab="Jahr",ylab="mittlere Tide in cm",
main="Niedrigwasser der Hunte und gleitender \n Durchschnitt der Ordnungen 3 und 12")
lines(nw12,col="red")
lines(nw3,col="green")
legend("topleft",c("Niedrigwasser","Durchschnitt (Ordnung: 3)",
"Durchschnitt (Ordnung: 12)"),col=c("blue","green","red"),pch=c(16,16,16),
title="Legende:",bg="#EEEEEE")
#Die Zeitreihe wird wieder mit Legende geplotet, dieses Mal
#werden ihr aber 2 "lowess"-Glättungskurven hinzugefügt.
plot(nw,col="blue",las=1,xlab="Jahr",ylab="mittlere Tide in cm",
main="Niedrigwasser der Hunte und 'lowess'- \n Durchschnitt mit verschiedenen Glättungsparametern")
lines(lowess(nw,f=1/8),col="green")
lines(lowess(nw,f=1/2),col="red")
legend("topleft",c("Niedrigwasser","Parameter: 1/8",
"Parameter: 1/2"),col=c("blue","green","red"),pch=c(16,16,16),
title="Legende:",bg="#EEEEEE")
#Wie zuvor, dieses Mal jedoch mit 2 "smooth.spline"-Glättungskurven.
plot(nw,col="blue",las=1,xlab="Jahr",ylab="mittlere Tide in cm",
main="Niedrigwasser der Hunte und 'smooth.spline'- \n Durchschnitt mit verschiedenen Glättungsparametern")
lines(smooth.spline(nw,df=17),col="green")
lines(smooth.spline(nw,df=4),col="red")
legend("topleft",c("Niedrigwasser","Parameter: 17",
"Parameter: 4"),col=c("blue","green","red"),pch=c(16,16,16),
title="Legende:",bg="#EEEEEE")
#Für die "loess.smooth"-Glättungskurven kann man leider
#nicht so vorgehen, da R nötige y-Werte fehlen.
plot(nw,col="blue",las=1,xlab="Jahr",ylab="mittlere Tide in cm",
main="Niedrigwasser der Hunte und 'loess.smooth'- \n Durchschnitt mit verschiedenen Glättungsparametern")
lines(loess.smooth(nw,span=1/3),col="green")
lines(loess.smooth(nw,span=1),col="red")
legend("topleft",c("Niedrigwasser","Parameter: 1/3",
"Parameter: 1"),col=c("blue","green","red"),pch=c(16,16,16),
title="Legende:",bg="#EEEEEE")
Jo, loess.smooth klappt bei mir bei 22b nicht.
Sollte das richtig sein und ihr wollts übernehmen, ändert das bitte irgendwie etwas ab, geht bei R ja eigentlich ganz gut.
Mit folgendem Code, können Sie den Beitrag ganz bequem auf ihrer Homepage verlinken
Weitere Beiträge aus dem Forum das Problem R
36b - gepostet von ClaudiaB am Donnerstag 14.06.2007
Sammeln für Klausur - gepostet von Shorty am Freitag 25.05.2007
Ähnliche Beiträge wie "Zettel 5"
Zettel Nr 6 Aufgabe 23a - ahnungslos (Mittwoch 30.05.2007)
Mahte Zettel? - sebo69 (Samstag 07.10.2006)
Zettel :: Irgendjemand - Irgendjemand - Philipp2005 (Sonntag 17.06.2007)
Physik-Zettel - july (Mittwoch 17.10.2007)
Zettel :: Kira - Louisa - Philipp2005 (Sonntag 17.06.2007)
Deutschhausaufgaben zum 15.05.07 / Zettel - Matze (Dienstag 15.05.2007)
Zettel 5 19B (kein r) - daniel (Donnerstag 24.05.2007)
>>>> Der Lehmann Zettel<<<<< - ZANDER (Mittwoch 18.10.2006)
Zettel vom Leisse die ich kopieren soll - Zealer (Montag 15.05.2006)
Zettel-Jagd - Ganralf (Dienstag 19.09.2006)