rm(list=ls(all=T))

# Daten einlesen
matphys <- read.table("D://Uebungen//2010_Sommersemester//Wirtschaftsstatistik_SS_10//Blatt 5//matphys.txt", header = TRUE)

# Daten als Streudiagramm darstellen
plot(matphys)

# lineare Regression mit "Punkte in Mathematik" als Ausgangsvariable
x <- matphys$Mathematik
y <- matphys$Physik
beta <- cov(x, y) / cov(x, x)
alpha <- mean(y) - beta * mean(x)
abline(alpha, beta, lwd = 2)

# lineare Regression mit "Punkte in Physik" als Ausgangsvariable
x <- matphys$Physik
y <- matphys$Mathematik
beta <- cov(x, y) / cov(x, x)
alpha <- mean(y) - beta * mean(x)
abline(alpha, beta, lty = 3, lwd = 2)