Monatsmittel

Allgemeine Statistik mit R, die Test-Methode ist noch nicht bekannt, ich habe noch keinen Plan!

Moderatoren: EDi, jogo

saskiasal47

Re: Monatsmittel

Beitrag von saskiasal47 »

Code: Alles auswählen

Daten<-read.table(file="Brocken_Gruppe_02.txt",header=T,dec=".",sep="",strip.white=T)
hatte header=F gesetzt gehabt

Code: Alles auswählen

 $ Jahr                  : int  2015 2015 2015 2015 2015 2015 2015 2015 2015 2015 ...
 $ Monat                 : int  1 1 1 1 1 1 1 1 1 1 ...
 $ Tag                   : int  1 2 3 4 5 6 7 8 9 10 ...
 $ LUFTTEMPERATUR        : num  -2 -2 -4.1 -3.6 -4.8 -4.7 -3.2 -2.3 -0.8 2.2 ...
 $ LUFTTEMPERATUR_MAXIMUM: num  -0.4 -0.5 -3.6 -3.2 -3.3 -1.5 -1.1 0.3 2.3 6.1 ...
 $ LUFTTEMPERATUR_MINIMUM: num  -3.2 -3.7 -4.8 -4 -6.4 -5.5 -5.1 -5.4 -2.7 -4.6 ...
 $ NIEDERSCHLAGSHOEHE    : num  0.1 3.6 6.1 12.4 0.7 1.2 4.2 21.8 25 36.6 ...
 $ SONNENSCHEINDAUER     : num  6.717 0 0 0 0.317 ...
>
jogo
Beiträge: 2086
Registriert: Fr Okt 07, 2016 8:25 am

Re: Monatsmittel

Beitrag von jogo »

sehr schön. Das sieht doch schon viel besser aus.
zwei kleine Hinweise:
1. TRUE bzw. FALSE möglichst nicht abkürzen.
2. viewtopic.php?f=20&t=29

Die Mittelwerte kannst Du jetzt so berechnen:

Code: Alles auswählen

tapply(Daten$LUFTTEMPERATUR, Daten$Monat, FUN=mean)
Gruß, Jörg
p.s.:
zu 1.:

Code: Alles auswählen

TRUE==T
FALSE==F
T <- 0
F <- 1
TRUE==T
FALSE==F
saskiasal47

Re: Monatsmittel

Beitrag von saskiasal47 »

ich habe mir mit dem Befehl:
Daten.Januar<-subset(d,d[,1]==1)

erst Mal die Daten für den Januar extrahiert, geht das nicht, dass ich da iwie den mean Befehl mit einbaue und das ganze noch in eine for Schleife packe für die drauffolgenden Monate?
jogo
Beiträge: 2086
Registriert: Fr Okt 07, 2016 8:25 am

Re: Monatsmittel

Beitrag von jogo »

natürlich geht auch das:

Code: Alles auswählen

MoMi <- numeric(12) # Vektor für Monatsmittel initialisieren
for (i in 1:12) {
  MoMi[i] <- mean(subset(Daten, Monat==i)$LUFTTEMPERATUR)
}
Du hattest aber gefragt:
Nun soll ich für die mittlere Lufttemperatur, für alle 12 Monate, das Monatsmittel berechnen. ...
Aber wie mache ich das am besten?
Bei dem jetzigen Stand der Aufgabe halte ich

Code: Alles auswählen

MoMi <- tapply(Daten$LUFTTEMPERATUR, Daten$Monat, FUN=mean)
für die beste Variante.

Gruß, Jörg
saskiasal47

Re: Monatsmittel

Beitrag von saskiasal47 »

super, vielen vielen Dank für eure Hilfe :)

ich soll aus den Daten der Lufttemperatur nur boxplots anlegen, also quasi für jeden Monat einen Boxplot und das ganze in einer Grafik

das ganze geht doch am besten auch wieder über eine Schleife oder? und dem Befehl boxplot

boxplot(d,range=1.5), wie klappt das denn am besten für jeden Monat ein Boxplot einzurichten? :?
jogo
Beiträge: 2086
Registriert: Fr Okt 07, 2016 8:25 am

Re: Monatsmittel

Beitrag von jogo »

Code: Alles auswählen

boxplot(iris$Sepal.Length ~ iris$Species)
Bei Dir dann:

Code: Alles auswählen

boxplot(Daten$LUFTTEMPERATUR ~ Daten$Monat)  ## bzw.:
boxplot(LUFTTEMPERATUR ~ Monat, data=Daten) 
Gruß, Jörg
saskiasal47

Re: Monatsmittel

Beitrag von saskiasal47 »

Vielen lieben Dank!

die restlichen Aufgaben bekomme ich jetzt denke ich alleine hin:)

Liebe Grüße
saskiasal47

Re: Monatsmittel

Beitrag von saskiasal47 »

Ich hätte doch nochmal eine Frage :D

ich soll noch das Monatsmittel als Punkt in den Boxpot einfügen, dass habe ich auch soweit hinbekommen.
Als letzte Aufgabe sollen nun alle Monate paarweise miteinander verglichen werden, dazu die Aufgabe:
Vergleichen Sie nun jeweils paarweise alle Monate miteinander. Führen Sie für jeden Paar-Vergleich mit Hilfe des Wilcoxon-Mann_Whitney Rangsummentests einen Signifikanztest auf Gleichheit durch mit dem Siknifikanzniveau 5%. Für welche Paar-Vergleiche lässt sich die Null-Hypothese: „Die beiden Monate unterscheiden sich nicht“ ablehnen. Fassen Sie die Ergebnisse dieser Signifikanztests in einer übersichtlichen Tabelle zusammen, in der Sie kennzeichnen, ob sich die jeweils miteinander verglichenen Monate signifikant voneinander unterscheiden oder nicht.

Das Programm läuft allerdings nicht ganz rund, wenn ich z.B. eine Schleife nehme,dann werde nur noch nie Wahrscheinlichkeiten gespeichert:(
Vllt sieht ja jemand wo das Problem liegt.

Hier mein Programm:
Daten<-read.table(file="Brocken_Gruppe_02.txt",header=T,dec=".",sep="",strip.white=T)

Monat<- Daten[,2]
Tag<-Daten[,3]
Lufttemperatur<-Daten[,4]
d<-data.frame(Monat,Tag,Lufttemperatur)

#Schleife für Mittelwerte der Lufttemperatur
MoMi <- numeric(12)
for (i in 1:12) {
MoMi <- mean(subset(Daten, Monat==i)$LUFTTEMPERATUR)
}
#
#


jpeg(filename="Aufgabe4_Boxplot_.jpg",width=960, height=960, units = "px", pointsize=14, quality=100, bg="white",res=144)
#
boxplot(Lufttemperatur~Monat,range=1.5,main="Luftemperaturen am Brocken im Jahr 2015", xlab="Monat",ylab="Temperatur in °C",ylim=c(-10,25),las=1)
points(c(MoMi),col="blue",pch=19)
#
legend("topleft",legend=c("Mittelwert"),col="blue",pch=19)
#
dev.off()
#
#
# Wenn ich die Schleife dazu nehme speichert er nur noch die Wahrscheinlichkeit :(
#MoVe<-numeric(12)
#for (j in 1:11){
d.x<-subset(Daten, Monat==1)$LUFTTEMPERATUR
d.y<-subset(Daten, Monat==1+1)$LUFTTEMPERATUR
MoVe<-(wilcox.test(d.x,d.y,exact=F,alternative="two.sided",conf.int=T))
#}
#
#
write.table(format(MoVe,digits=3,nsmall=2,width=6),file="Aufgabe4_Signifikanz_t",quote=F,sep=" ",row.name=F)
jogo
Beiträge: 2086
Registriert: Fr Okt 07, 2016 8:25 am

Re: Monatsmittel

Beitrag von jogo »

Hallo Saskia,

bitte lies: viewtopic.php?f=20&t=29
und formatiere den Code in Deiner Nachricht!
Unabhängig davon möchte ich Dir die Funktion pairwise.wilcox.test() empfehlen:

Code: Alles auswählen

pairwise.wilcox.test(iris$Sepal.Length, iris$Species)
Gruß, Jörg
saskiasal47

Re: Monatsmittel

Beitrag von saskiasal47 »

Code: Alles auswählen

Daten<-read.table(file="Brocken_Gruppe_02.txt",header=T,dec=".",sep="",strip.white=T)

Monat<- Daten[,2]
Tag<-Daten[,3]
Lufttemperatur<-Daten[,4]
d<-data.frame(Monat,Tag,Lufttemperatur)

#Schleife für Mittelwerte der Lufttemperatur
MoMi <- numeric(12)
for (i in 1:12) {
MoMi <- mean(subset(Daten, Monat==i)$LUFTTEMPERATUR)
}
#
#


jpeg(filename="Aufgabe4_Boxplot_.jpg",width=960, height=960, units = "px", pointsize=14, quality=100, bg="white",res=144)
#
boxplot(Lufttemperatur~Monat,range=1.5,main="Luftemperaturen am Brocken im Jahr 2015", xlab="Monat",ylab="Temperatur in °C",ylim=c(-10,25),las=1)
points(c(MoMi),col="blue",pch=19)
#
legend("topleft",legend=c("Mittelwert"),col="blue",pch=19)
#
dev.off()
#
#
# Wenn ich die Schleife dazu nehme speichert er nur noch die Wahrscheinlichkeit :(
#MoVe<-numeric(12)
#for (j in 1:11){
d.x<-subset(Daten, Monat==1)$LUFTTEMPERATUR
d.y<-subset(Daten, Monat==1+1)$LUFTTEMPERATUR
MoVe<-(wilcox.test(d.x,d.y,exact=F,alternative="two.sided",conf.int=T))
#}
#
#
write.table(format(MoVe,digits=3,nsmall=2,width=6),file="Aufgabe4_Signifikanz_t",quote=F,sep=" ",row.name=F)]
gibt es noch eine andere Möglichkeit ausser pairwise?
Antworten