load("diagnozaOS.Rdata")
#attach(dane)
#head(dane)

#dane[, "bp124"]
#dane <- dane[which(! is.na(dane[, "bp124"])), ]

# Wczytanie danych
readData <- function() {
  dane <<- dane[which(! is.na(dane[, "bp124"])), ]
  dane <- na.omit(dane)
}

#readData()

lB = 0.05
uB = 0.95

# Obcięcie skrajnych wartosci
trimData <- function(lowerBound = lB, upperBound = uB) {
  dane <<- dane[dane[, "bp124"] > quantile(dane[, "bp124"], probs = lowerBound), ]
  dane <<- dane[dane[, "bp124"] < quantile(dane[, "bp124"], probs = upperBound), ]
}

readData()
trimData(lowerBound = lB, upperBound = uB)

# Nowe nazwy kolumn
dane_wartosciowe <- dane[, c("bp124", "cp124", "bp129_11", "bp125", "bp129_5", "bp126_1", "bp130", "bp131")]
colnames(dane_wartosciowe) <- c("godzin_tygodniowo", "ocena_umiejetnosci", "sciaganie_softu_muzyki", "czy_korzystales_z_internetu",
                    "przegladanie_www", "w_pracy", "ile_godzin_tygodniowo_w_internecie", "od_kiedy_uzywa_internetu")

fmla <- as.formula("log(godzin_tygodniowo) ~ ocena_umiejetnosci + sciaganie_softu_muzyki + czy_korzystales_z_internetu + przegladanie_www + w_pracy + ile_godzin_tygodniowo_w_internecie")
model_liniowy <- lm(formula = fmla, data = dane_wartosciowe)

# Interesujaca zmienna: 658 bp131 "Kiedy zaczšł Pan korzystać z Internetu"