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

##########################
#                        #
#  Jednovyberovy t-test  #
#                        #
##########################

ubytek<-c(2.00,1.20,-0.12,2.35,0.40,0.56,4.23,2.35,1.11,-0.02,0.00,3.10,0.05,2.14,2.04,3.10,4.05,1.88,3.02,0.57,-1.78,2.00,2.54,1.21)
referencni<-2.12

# prumerny ubytek hmotnosti
mean(ubytek)

# smerodatna odchylka hmotnosti
sd(ubytek)

# stupne volnosti
length(ubytek)-1

# t na hladině 95 % spolehlivosti
qt(1-0.05/2,length(ubytek)-1)

# t na hladině 90 % spolehlivosti
qt(1-0.10/2,length(ubytek)-1)

# vypocet p-hodnoty
t.test(ubytek,mu=2.12)$p.value


#########################
#                       #
#  Dvouvyberovy t-test  #
#                       #
#########################

podzim2010<-c(75,44,62,57,54,59,65,62,68,70,68,59,48,42,28,90,53,76,74,68)
jaro2011<-c(63,78,83,66,60,63,98,82,79,78,79,71,93,52,38,85,41,26,89,95,83,88,68,62,59)

# p-hodnota
t.test(podzim2010,jaro2011,var.equal=TRUE)

# prumer podzim 2010
mean(podzim2010)

# prumer jaro 2011
mean(jaro2011)


###################
#                 #
#  Parovy t-test  #
#                 #
###################

pbr<-c(2.51,2.52,1.93,1.42,2.03,3.10,1.01,1.09,2.35,2.31,2.06,1.92,1.75,2.19,2.35,2.50,2.37,2.49,2.03,1.75,1.03,1.70,2.82,2.94,1.18,2.84,1.60)
psr<-c(2.46,2.32,1.78,1.44,2.32,3.12,1.19,1.27,2.29,2.14,2.15,2.18,1.54,2.38,1.94,2.67,2.00,2.55,2.03,1.65,1.35,1.63,2.73,2.89,1.07,2.63,1.62)

# prumerny prospech pred bitim
mean(pbr)

# prumerny prospech po biti
mean(psr)

# smerodatna odchylka pred bitim
sd(pbr)

# smerodatna odchylka po biti
sd(psr)

# rozdil prospechu
mean(psr)-mean(pbr)

# s (upravene)
sqrt(sd(pbr)^2+sd(psr)^2-2*cov(pbr,psr))

# stupne volnosti
length(pbr)-1

# 1-α/2
1-0.05/2

# t na hladině 95 % spolehlivosti
qt(1-0.05/2,length(pbr)-1)

# vypoctene t
(mean(pbr)-mean(psr))/sqrt(sd(pbr)^2+sd(psr)^2-2*cov(pbr,psr))*sqrt(length(pbr)-1)

# n potrebne k vyznamnosti (pozor spoustu „nekonecny“ cyklus)
t<-(mean(pbr)-mean(psr))/sqrt(sd(pbr)^2+sd(psr)^2-2*cov(pbr,psr))*sqrt(length(pbr)-1)
stupne_volnosti<-1
while (qt(1-0.05/2,stupne_volnosti)>t) {
  stupne_volnosti<-stupne_volnosti+1
}
stupne_volnosti

# p-hodnota
t.test(pbr,psr,paired=TRUE)