R Project [Regessionsgerade]

[Minimax83]

Hallo,

Ich brauche mal eure Hilfe. Na dem ganzen Suchen im Internet habe ich schon Kopfschmerzen bekommen...

Also:

Nach langem Suchen habe ich es bis hier geschafft. Ich hab x und y werte und lasse diese zeichnen:

cx <- c(1,2,3,4)
cy <- c(0.1,0.5,0.6,0.9)

plot(cx, cy, type = "p", main="Lineare Regression")
abline(lm(cy ~ cx))

variable = lm(cy ~ cx)
variable$coefficients[1]
variable$coefficients[2]

Wie kann ich:

1. die komplette Geradengleichung und das Bestimmtheitsmaß in den Plot zeichnen?
mit summary kann ich zwar R^2 etc. sehen, aber habe keine ahnung, wie ich weiter vorgehen muss.
2. wie sagen, dass die gerade mal durch den Ursprung gehen soll?
3. Wie kann ich einstellen, dass die gerade vor dem 1. x-wert und nach dem 2. x-wert nicht mehr weitergezeichnet wird?

4. ps: Wie zeichne ich eine quadratische trendlinie in den plot?

danka,
minimax

 

[FrankR]

Hallo,
1. die Funktion um Text in einen Plot zu bringen lautet >text()<. Es besteht hiermit zusammen mit der Konkatenation (>paste()<) der Koeffizienten die Möglichkeit die Gleichung in den Plot zu bekommen.
2. Zunächst die Koeffizienten anpassen und erst danach plotten.

Frank

 

[Minimax83]

Danke. Das mit der Funktionsgleichung zeichnen habe ich jetzt hinbekommen.

Am wichtigsten wäre noch für mich, wie ich bei der Quadratischen Funktion R² ausrechnen kann.

Beispiel:

y = 0.1108x² + 22.94x -0.2765 ist das ergebnis einer quad ratischen regression. wenn ich mir das ergebnis mit summary anschaue, bekomme ich:

---------------------
Call:
lm(formula = cy ~ cx + I(cx^2))

Residuals:
1 2 3 4
-0.475 2.224 -2.738 0.989

Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.2765 5.2094 -0.05 0.966
cx 22.9381 0.7549 30.39 0.021 *
I(cx^2) 0.1108 0.0227 4.88 0.129
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 3.69 on 1 degrees of freedom
Multiple R-squared: 1, Adjusted R-squared: 1
F-statistic: 1.01e+04 on 2 and 1 DF, p-value: 0.00702


---------------------

Dort steht: Multiple R-squared: 1
Doch wie komme ich an diesen Wert, so dass ich ihn in eine Variable schreiben kann?

 

[FrankR]

Hallo,
Auf das Bestimmtheitsmaß kannst du via summary(lm(...))$r.squared zugreifen.

lm(cy~0+cx) liefert Regressionsgerade, die durch den Ursprung läuft.

Frank