setwd("C:/Jirka/Environment/Výuka/E7541 Analýza dat na PC/2022 podzim")

library("xlsx")

#######################
#                     # 
#  Mann-Whitney test  #
#                     # 
#######################

data0<-read.xlsx("10_neparametricke_testy_zadani.xlsx",2)
pozitivni<-as.numeric(data0[6:16,2])
negativni<-as.numeric(data0[6:16,4])

wilcox.test(pozitivni,negativni,paired=FALSE)


##############################
#                            # 
#  Wilcoxonuv (parovy) test  #
#                            # 
##############################

data1<-read.xlsx("10_neparametricke_testy_zadani.xlsx",3)
data<-data1[7:32,3:5]
colnames(data)<-c("datum","digitel","ps1")
data$datum<-as.Date(as.numeric(data$datum),origin="1899-12-30")
data$digitel<-as.numeric(data$digitel)
data$ps1    <-as.numeric(data$ps1)

shapiro.test(data$ps1)
shapiro.test(data$digitel)

hist(data$ps1)
hist(data$digitel)

wilcox.test(data$ps1,data$digitel,paired=TRUE)
boxplot(data$ps1,data$digitel,xaxt="n")
axis(1,labels=c("PS-1","Digitel"),at=c(1,2))

# data lze ovsem i zlogaritmovat
hist(log10(data$ps1))
hist(log10(data$digitel))

shapiro.test(log10(data$ps1))
shapiro.test(log10(data$digitel))

t.test(data$ps1,data$digitel,paired=TRUE)


#####################
#                   # 
#  Znamenkovy test  #
#                   # 
#####################

binom.test(13,31)
binom.test(12,31)