# Nastavení pracovního adresáře
setwd("C:/Jirka/Environment/Výuka/Bi7541 Analýza dat na PC/2020 podzim")

# Načtení dat ve formátu csv oddělených středníkem
data<-read.csv2(file="05_statistica.csv",header=TRUE,dec=",",sep=";",skip=0)

# Počet radku a sloupcu souboru
nrow(data)
ncol(data)

# Kontrola pojmenovani sloupcu
colnames(data)

# T-test vlivu pohlavi na delku hospitalizace
class(data$los)
data$los<-as.numeric(as.character(data$los))

muzi<-data[which(data$gender=="Male"),]
zeny<-data[which(data$gender=="Female"),]

t.test(muzi$los,zeny$los,paired=FALSE,alternative=c("two.sided"),na.rm=TRUE)$p.value

# T-test vlivu koureni na delku hospitalizace
  kuraci<-data[which(data$smoker=="Yes"),]
nekuraci<-data[which(data$smoker=="No"),]

t.test(kuraci$los,nekuraci$los,paired=FALSE,alternative=c("two.sided"),na.rm=TRUE)[3]

# Sloupce r1, r2, r3 a rt
ncol(data)

data$r1<-runif(nrow(data),0,1)
data$r2<-runif(nrow(data),0,1)
data$r3<-runif(nrow(data),0,1)

data$rt<-data$r1+data$r2+data$r3

data$rt<-(data$rt-mean(data$rt))/sd(data$rt)

# Histogram
deleni<-c(0:20)*(max(data$rt)-min(data$rt))/20+min(data$rt)
hist(data$rt,breaks=deleni)

install.packages("genasis")
library(genasis)
genhistogram (data$rt,breaks=19,emboss=3,col="orange")