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

# 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[,34]<-runif(nrow(data),0,1)
data[,35]<-runif(nrow(data),0,1)
data[,36]<-runif(nrow(data),0,1)

data[,37]<-data[,34]+data[,35]+data[,36]

data[,37]<-(data[,37]-mean(data[,37]))/sd(data[,37])

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

install.packages("genasis")
library(genasis)
genhistogram (data[,37],breaks=19,emboss=3,col="orange")