Verfügbare Informationen zu "zettel 6 aufg 27"
Qualität des Beitrags: Beteiligte Poster: Lucie - daniel - Tommy Gun - rage Forum: das Problem R aus dem Unterforum: Vergleichen Antworten: 4 Forum gestartet am: Freitag 18.05.2007 Sprache: deutsch Link zum Originaltopic: zettel 6 aufg 27 Letzte Antwort: vor 16 Jahren, 10 Monaten, 17 Tagen, 17 Stunden, 47 Minuten
Alle Beiträge und Antworten zu "zettel 6 aufg 27"
Re: zettel 6 aufg 27
Lucie - 30.05.2007, 17:53zettel 6 aufg 27
so, stell mal meine lösung von 27 rein, hab aber keine ahnung ob das so richtig ist, deshalb bitte nochmal nachprüfen oder vergleichen...und sonst halt nen bissl abändern wär nett :D
#Aufgabe 27
#(a)
#alles löschen
rm(list=ls(all=TRUE))
#alle Graphiken speichern (zum Blättern mit PgUp und PgDn)
windows(record=T)
n=1
m=500
n1 = rnorm(m)
x1 = n1^2
n=5
m=500
n1 = rnorm(m)
n2 = rnorm(m)
n3 = rnorm(m)
n4 = rnorm(m)
n5 = rnorm(m)
x5 = n1^2 + n2^2 + n3^2 + n4^2 + n5^2
#Histogramm der X^2-Werte
#mit prob=T ist die Gesamtfläche der Säulen gleich Eins
hist(x1,prob=T)
hist(x5,prob=T)
#Dichte der chi^2-Verteilung mit n Freiheitsgraden
chi1=dchisq(x1,1) #n=1
#Kontrolle
chi1
chi2=dchisq(x5,5) #n=5
#Kontrolle
chi2
#erstellen eines data.frame mit Varaiblen chi1 und chi2
Daten = data.frame(chi1, chi2)
#Kontrolle
Daten
#speichern der Daten als csv-Datei mit dem Namen "aufg27"
write.csv2(Daten, file="aufg27.csv")
#Kontrolle des speicherns mit wiederholtem einlesen der Datei1
Daten = read.csv2("aufg27",header=TRUE) #funkts nicht!!!!
#Daten im R-eigenen Format als binäre Datei speichern
save(chi1, chi2, file="aufg27.RData")
#Datei kann mit load("") wieder eingelesen werden
load("aufg27.RData")
###(b)
##Exeldatensatz "klausur" einlesen und vorher im csv-Format verfügbar machen
##Umbenennen der eingelesenen Datei in Data und speichern
Daten = read.csv2("klausur0607.csv",header=T)
#Kontrolle
Daten
#Variablen verfügbar machen
attach (Daten)
Data = Daten$Klausur
#Kontrolle
Data
#(1) normaler QQ-Plot mit vertauschten Achsen durch datax=FALSE
qqnorm(Data,main="Normaler QQ-Plot")
#Hilfslinie mit qqline, so ist die Normalenverteilung besser erkennbar
qqline(Data,datax=FALSE,col="green")
##(2) per Hand
#Quantile der Standardnormalverteilung für
#die in der Übung gegebenen alpha_i
n=86
alpha = 1:n/(n+1)
#Kontrolle
alpha
y = qnorm(alpha)
#Kontrolle
y
x = sort(Data)
#Kontrolle
x
plot(y,x) #nicht ganz genau das errechnete Bild von R, deshalb:
alpha2 = (1:n-1/2)/n
#Kontrolle
alpha2
y2 = qnorm(alpha2)
#Kontrolle
y2
plot(y2,x)
#oder auch
alpha3 = ppoints(n)
#Kontrolle
alpha3
y3 = qnorm(alpha3)
#Kontrolle
y3
plot(y3,x) ##plot(y2,x) und plot(y3,x) gleich mit dem von errechneten von R
##Teil 2
#(1)
#Datei ist im R-eigenen Format als binäre Datei gespeichert und wird wieder eingelesen
load("aufg27.RData")
#abrufen der Werte für x^2-Werte zu n=5
x5
#normaler QQ-Plot zu den x^2-Werten zu n=5,
#mit vertauschten Achsen durch datax=FALSE
qqnorm(x5,main="Normaler QQ-Plot")
#Hilfslinie mit qqline, so ist die Normalenverteilung besser erkennbar
qqline(x5,datax=FALSE, col="green")
#abrufen von chi2
chi2
#theoretischen Quantile der chi^2-Verteilung
z = qchisq(chi2, df=5)
#QQ-Plot mit zwei Datensätzen (chi2 und z)
qqplot(chi2,z)
so, das wars, viel Spaß dabei :P
Re: zettel 6 aufg 27
daniel - 31.05.2007, 12:08
huhu
Da claudia noch nicht im forum regestriert ist poste ich hier ihre genialen verbesserungsforschläge!
hey! bei dir ging das einlesen bei 27a nicht, weil du die dateiendung vergessen hattest:n readcsv2("aufg27.csv")
und ich glaube, dass : save(chi1, chi2, file="aufg27.RData")
muss so heißen: save(x1, x5, file="aufg27.RData")
load("aufg27.RData")
Re: zettel 6 aufg 27
Tommy Gun - 31.05.2007, 20:38
Hi,
sollte die Dichte der X^2-Verteilung nicht ins Histogramm eingezeichnet werden? Unser Tutor meinte, das gehe mit "curve(...)".
Hab damit rumprobiert, habs aber nicht hinbekommen, also hab ich "lines(dchisq(x1,1))" bzw. "lines(dchisq(x5,1))" genommen. Dann zeichnet R nen Graphen ins Histogramm, ob das allerdings richtig ist, weiß ich auch nicht. Sieht jedenfall etwas merkwürdig aus.
Außerdem bin ich auch der Meinung, dass x1 und x5 gespeichert und im b)-Teil verwendet werden sollen.
Re: zettel 6 aufg 27
rage - 31.05.2007, 20:40
Hey Daniel! Hast du auch einen tollen Vorschlag für die 27b? ich bekomme den qqplot nicht hin!
lg, Rahel
Mit folgendem Code, können Sie den Beitrag ganz bequem auf ihrer Homepage verlinken
Weitere Beiträge aus dem Forum das Problem R
36c - gepostet von ClaudiaB am Donnerstag 14.06.2007
36 a) - gepostet von Tomcio am Dienstag 12.06.2007
Ähnliche Beiträge wie "zettel 6 aufg 27"
netter zettel am auto - djder1 (Montag 05.11.2007)
e-fensterheber - Richta (Donnerstag 31.07.2008)
deutsch-zettel-dasparfum - Anonymous (Dienstag 03.05.2005)
Einige Zettel an der Tarvene - linguist (Dienstag 10.05.2005)
Mathe-Zettel - abi2010stg (Samstag 17.09.2005)
HA: Zettel vom 5.9.07 Aufg. 6 - kowalski (Dienstag 18.09.2007)
Warum die Zettel nicht immer gelb trotz neuem Beitrag - Eldariel (Samstag 29.09.2007)
Zettel - Bernhard (Samstag 03.12.2005)
Physik-Zettel - july (Mittwoch 17.10.2007)
Ein Zettel am Schwarzen Brett - Nyachelle (Samstag 09.02.2008)