daten <- read.csv("C:/Dokumente und Einstellungen/w6reto2/Desktop/desktopJena/teaching/empirie/data/cars.csv",header=TRUE, sep=";")
  # replace pathname with C:/data for instance
  
   result <- lm(S ~ T + E + P + H,data=daten)
   summary(result)
   
  
#   Coefficients:
#              Estimate Std. Error t value Pr(>|t|)    
# (Intercept)  9.2239226  1.9354666   4.766 3.67e-05 ***
# T           -0.7865682  0.5902559  -1.333   0.1918    
# E            7.9060736  3.6601435   2.160   0.0381 *  
# P            0.0004080  0.0004954   0.824   0.4161    
# H           -0.0189712  0.0026741  -7.094 4.01e-08 ***

# Adjusted R-squared: 0.6742

# Variable H in inverser Form

daten$h = 1/daten$H

   result <- lm(S ~ T + E + P + h,data=daten)
   summary(result)

    # Ergebnis: Bessere T-Werte, höhere Schätzgüte

# Alternative Spezifikation mit M/H, anstelle von M + 1/H

daten$p <- daten$P/daten$H

   result <- lm(S ~ T + E + p,data=daten)
   summary(result)
   
    # Ergebnis: Sehr hoher T-Wert für p, noch höhere Schätzgüte, Effekt von M und H nicht trennbar   

# Alternative Spezifikation mit M, 1/H und M/H, anstelle von M + 1/H

   result <- lm(S ~ T + E + P + h + p,data=daten)
   summary(result)