##########################################################################################
install.packages("nortest")       # nacteni balicku pro testovani normality
library("nortest")
install.packages("DescTools")     # nacteni balicku pro znamenkovy test
library("DescTools")
install.packages("car")           # nacteni balicku pro Levenuv test
library(car)
# Nastaveni pracovniho adresare - napr.:
setwd("c:/Users/littnerovas/Documents/VYUKA_2019/2020/vicerozmerky/02/")
getwd()
kardio <- read.table(file="Data_pro_R.txt", sep="\t", dec=".", header=T) # nacteni dat
dim(kardio)                       # pocet radku a sloupcu datove matice
head(kardio)                      # prvnich 10 radku tabulky
str(kardio)                       # sumarni informace o jednotlivych promennych
is(kardio)				                # typ objektu
names(kardio)
hist(kardio$sys_tlak)
# Histogram s krivkou hustoty normalniho rozdeleni
x <- kardio$sys_tlak
h<-hist(x, col="gray", xlab="vyska",main="Histogram s krivkou hustoty normalniho rozdeleni")
xfit<-seq(min(x),max(x),length=length(unique(x)))
yfit<-dnorm(xfit,mean=mean(x),sd=sd(x))
yfit <- yfit*diff(h$mids[1:2])*length(x)
lines(xfit, yfit, col="red", lwd=2)
boxplot(kardio$sys_tlak)
qqnorm(kardio$sys_tlak) # Q-Q graf: teoreticky (za predpokladu X~N) vs. pozorovany kvantil
# Kolmogorovuv - Smirnovuv test: Lilieforsovu variantu testu pocitame vzdy, kdyz nezname stredni hodnotu a rozptyl dat
lillie.test(kardio$sys_tlak)
dim(kardio)
