Autor |
Nachricht |
cubicglobe
|
Betreff des Beitrags: EDV3 Streudiagramm Verfasst: 18.06.2012, 12:21 |
|
Eprouvettenschüttler/in |
|
|
Registriert: 12.09.2009, 11:46 Beiträge: 113
|
ich habe ein problem mit dem streudiagramm. das hab ich bisher eingegen: Zitat: > edv3 <- read.xls(file.choose()) > modell <- lm(leit ~ temp, data=edv3) > summary (modell)
> hilfsw <- data.frame(temp=seq(min(edv3$temp), max(edv3$temp), lenght=100)) > konfidenzband <- predict(modell, hilfsw, intervall="conf") > vorhersageband <- predict(modell, hilfsw, intervall="pred")
> plot(edv3$temp, edv3$leit) > abline (modell) > matlines (hilfsw$temp, konfidenzband, lty=c("blank", "dotted", "dotted"), col=c("black", "red", "red")) > matlines (hilfsw$temp, vorhersageband, lty=c("blank", "dashed", "dashed"), col=c("black", "green", green")) er generiert zwar das diagramm mit plot und zeigt auch die abline an, aber die linien für konfidenzband und vorhersageband erscheinen nicht in der grafik. wo ist der fehler?
|
|
|
|
|
newbie18
|
Betreff des Beitrags: Re: EDV3 Streudiagramm Verfasst: 18.06.2012, 12:35 |
|
Mikroskopierer/in |
|
|
Registriert: 18.06.2010, 09:39 Beiträge: 186
|
Hey versuchs mal mit dem hier du musst dann nur noch deine Daten eingeben Lg newbie scatterplot( ertrag~ grasanteil, reg.line=lm, smooth=FALSE, spread=FALSE, boxplots=FALSE, span=0.5, data= gras) hilfswerte <- data.frame( grasanteil= seq(min( gras$ grasanteil), max( gras$ grasanteil), length=100)) konfidenzband <- predict(modell, hilfswerte, interval="conf") vorhersageband <- predict(modell, hilfswerte, interval="pred") matlines(hilfswerte$ grasanteil, konfidenzband, lty=c("blank", "dotted", "dotted"), col=c("black", "red", "red")) matlines(hilfswerte$ grasanteil, vorhersageband, lty=c("blank", "dashed", "dashed"), col=c("black", "green", "green")) P.S. Nicht vergessen alles zu markieren, was aufgeführt werden soll! Und wenn du dich gar nicht auskennst einfach mal hier gucken https://moodle.boku.ac.at/mod/resource/ ... p?id=76881
|
|
|
|
|
cubicglobe
|
Betreff des Beitrags: Re: EDV3 Streudiagramm Verfasst: 18.06.2012, 12:42 |
|
Eprouvettenschüttler/in |
|
|
Registriert: 12.09.2009, 11:46 Beiträge: 113
|
danke, das sieht schon sehr gut aus so. weisst du vielleicht auch noch wie man das zentriert? das vorhersageband ist nicht komplett am plot. und was meinst du mit markieren? die achsenbeschriftung oder wie?
|
|
|
|
|
newbie18
|
Betreff des Beitrags: Re: EDV3 Streudiagramm Verfasst: 18.06.2012, 12:49 |
|
Mikroskopierer/in |
|
|
Registriert: 18.06.2010, 09:39 Beiträge: 186
|
cubicglobe hat geschrieben: danke, das sieht schon sehr gut aus so. weisst du vielleicht auch noch wie man das zentriert? das vorhersageband ist nicht komplett am plot. und was meinst du mit markieren? die achsenbeschriftung oder wie? Hey ich bin mir nicht sicher mit dem zentrieren aber vielleicht eventuell bei length eine kleinere zahl eingeben Mit markieren meinte ich, den Input also das was du hineinkopiert hast zu markieren, weil nur dann der Befehl ausgeführt wird. Lg
|
|
|
|
|
cubicglobe
|
Betreff des Beitrags: Re: EDV3 Streudiagramm Verfasst: 18.06.2012, 13:13 |
|
Eprouvettenschüttler/in |
|
|
Registriert: 12.09.2009, 11:46 Beiträge: 113
|
lenght ist es nicht. kann mann das eventuell mit ylim machen?
|
|
|
|
|
newbie18
|
Betreff des Beitrags: Re: EDV3 Streudiagramm Verfasst: 18.06.2012, 13:45 |
|
Mikroskopierer/in |
|
|
Registriert: 18.06.2010, 09:39 Beiträge: 186
|
Hm weiß ich leider nicht, bei mir war das Vorhersageintervall und Kofidenzintervall ziemlich nahe an der Regressionslinie Lg
|
|
|
|
|
cubicglobe
|
Betreff des Beitrags: Re: EDV3 Streudiagramm Verfasst: 18.06.2012, 13:48 |
|
Eprouvettenschüttler/in |
|
|
Registriert: 12.09.2009, 11:46 Beiträge: 113
|
ok. meines sieht dzt. so aus:
Dateianhänge: |
scatter.png [6.21 KiB]
653-mal heruntergeladen
|
|
|
|
|
|
newbie18
|
Betreff des Beitrags: Re: EDV3 Streudiagramm Verfasst: 18.06.2012, 13:50 |
|
Mikroskopierer/in |
|
|
Registriert: 18.06.2010, 09:39 Beiträge: 186
|
Naja immerhin sieht man die eine Seite Ich denke das ist nicht so schlimm wenn man das nicht sieht
|
|
|
|
|
cubicglobe
|
Betreff des Beitrags: Re: EDV3 Streudiagramm Verfasst: 18.06.2012, 13:59 |
|
Eprouvettenschüttler/in |
|
|
Registriert: 12.09.2009, 11:46 Beiträge: 113
|
das ist btw meine angabe: Um die Prozesssteuerung eines groen Lebensmittelherstellers zu optimieren, wird unter anderem auch dieWarmeleitfahigkeit [mW/(cmK)] von Staubzucker in Abhangigkeit von der Verarbeitungstemperatur [o C] untersucht (Code ZU-3, Aufgabencode 32391): Temperatur 30 30 30 34 34 34 38 38 38 60 Leitfaehigkeit 9.03 9.22 9.12 9.22 9.26 8.92 9.24 9.00 9.03 9.27 Temperatur 60 60 64 64 64 68 68 68 Leitfaehigkeit 9.15 8.96 9.39 8.83 9.20 9.72 9.06 9.49 Geben Sie ein lineares Modell der Form y = a + bx fur den Zusammenhang zwischen Temperatur und Leitfahigkeit an. Schatzen Sie die zugehorigen Koezienten a und b sowie die Fehlervarianz. Testen Sie auerdem die Nullhypothesen H0 : a = 0 und H0 : b = 0 bei einem Risiko 1. Art von = 0:05. (4) Fertigen Sie ein Streudiagramm der Daten an, das die geschatzte Regressionsgerade, das Kondenzband und das Vorhersageband enthalt. (1) Beurteilen Sie an Hand numerischer und grascher Kriterien, wie gut die Modellgerade den Zusammenhang zwischen den beiden Variablen beschreibt. (2) Uberprufen Sie an Hand einer geeigneten graschen Darstellung, ob die fur die Regression notwendige Voraussetzung erfullt ist. (1) Geben Sie ein 95 %-iges Kondenzintervall fur die erwartete Warmeleitfahigkeit bei 50o C an. (1) Geben Sie ein 99 %-iges Vorhersageintervall fur einen gemessenen Wert der Warmeleitf ahigkeit bei 45o C an. [schei? encoding ]
|
|
|
|
|
cubicglobe
|
Betreff des Beitrags: Re: EDV3 Streudiagramm Verfasst: 18.06.2012, 14:19 |
|
Eprouvettenschüttler/in |
|
|
Registriert: 12.09.2009, 11:46 Beiträge: 113
|
so hat es jetzt geklappt:
> plot (edv3$temp, edv3$leit, xlab="Temp", ylab="Leitfähigkeit", main="Leitfähigkeit von Staubzucker", ylim=c(min(vorhersageband[,2]), max(vorhersageband[,3]))) > abline (modell) > matlines(hilfswerte$temp, konfidenzband, lty=c("blank","dotted","dotted"), col=c("black", "red", "red")) > matlines(hilfswerte$temp, vorhersageband, lty=c("blank", "dashed", "dashed"), col=c("black", "green", "green"))
|
|
|
|
|
newbie18
|
Betreff des Beitrags: Re: EDV3 Streudiagramm Verfasst: 18.06.2012, 14:33 |
|
Mikroskopierer/in |
|
|
Registriert: 18.06.2010, 09:39 Beiträge: 186
|
Hast du da nur die Zahlen neben dem Vorhersageintervall hingeschrieben? Muss man sich merken
|
|
|
|
|
sminim
|
Betreff des Beitrags: Re: EDV3 Streudiagramm Verfasst: 19.06.2012, 12:41 |
|
Eprouvettenschüttler/in |
|
|
Registriert: 07.10.2010, 19:54 Beiträge: 176
|
bei mir sieht das ganze sehr komisch aus. kann das stimmen? habs mal angehängt
Dateianhänge: |
bild10.jpeg [ 24.42 KiB | 617-mal betrachtet ]
|
|
|
|
|
|
newbie18
|
Betreff des Beitrags: Re: EDV3 Streudiagramm Verfasst: 19.06.2012, 12:59 |
|
Mikroskopierer/in |
|
|
Registriert: 18.06.2010, 09:39 Beiträge: 186
|
|
|
|
|
weiland
|
Betreff des Beitrags: Re: EDV3 Streudiagramm Verfasst: 19.06.2012, 13:01 |
|
Mikroskopierer/in |
|
|
Registriert: 09.10.2011, 08:02 Beiträge: 189
|
sminim hat geschrieben: bei mir sieht das ganze sehr komisch aus. kann das stimmen? habs mal angehängt Deine Werte korrelieren sehr gut, daher sind die Bänder recht schmal! Ist das eine Wiederfindungsfunktion die du da geplottet hast?!
|
|
|
|
|
sminim
|
Betreff des Beitrags: Re: EDV3 Streudiagramm Verfasst: 19.06.2012, 13:09 |
|
Eprouvettenschüttler/in |
|
|
Registriert: 07.10.2010, 19:54 Beiträge: 176
|
ok danke. dachte ich mir eh. also die aufgabe war: Fertigen Sie ein Streudiagramm der Daten an, das die geschatzte Regressionsgerade, das Kondenzband und das Vorhersageband enthalt.
und ich habe es so gelöst: > Hilfswerte <- data.frame(Standard=seq(min(Kalibrierung$Standard),max(Kalibrierung$Standard), length=100)) > Konfidenzband <- predict(Modell, Hilfswerte, interval ="conf") > Vorhersageband <- predict(Modell, Hilfswerte, interval="pred") > plot(Kalibrierung$Standard, Kalibrierung$Messung, xlab="Standard [%]", ylab="Messung [%]", main="Standard und Messung von Kalibrierung", ylim=c(min(Vorhersageband[,2]),max(Vorhersageband[,3]))) > abline (Modell) > matlines(Hilfswerte$Standard, Konfidenzband, lty=c("blank", "dotted", "dotted"), col=c("black", "red", "red")) > matlines(Hilfswerte$Standard, Vorhersageband, lty=c("blank", "dashed", "dashed"), col=c("black", "blue", "blue"))
beim auswerten tu ich mir etwas schwer. ich würde es so machen: Anhand des Wertes "Multiple R-squared" aus r^2=0.9982 kann man ablesen, dass rund 99.8% der Variabilität der y-Werte durch das Modell erklärt werden. Das ist sehr hoch und deshalb sind die Bänder recht schmal. Es sind 11 Werte auf der Grafik erkennbar sind obwohl 20 eingegeben wurden. Das rührt daher, dass manche Wertepaare sich überlagern, also gleich sind.
|
|
|
|
|
|