library('MKinfer')
library(boot)

#test t-studenta dla jednej próby
load("Hamulce.RData")
attach(Hamulce)
mean(Wynik)
t.test(Wynik,mu=18.6,alternative='less')
boot.t.test(Wynik,mu=18.6,alternative='less', R = 100)

par(mfrow=c(1,2))

hist(Wynik,prob=T,main='')


meanboot <- function(x,i)mean(x[i])
bmean=boot(Wynik,meanboot,100)

hist(bmean$t,prob=T,main='')

t.test(bmean$t,mu=18.6,alternative='less')




# dwie próby niezależne
load("Bledy.RData")
attach(Bledy)
tapply(Blad,Kod,mean)
tapply(Blad,Kod,var)
var.test(Blad~Kod)
t.test(Blad~Kod,var.equal=T,alternative='less')

#dwie próby zależne
load("Cisnienie.RData")
attach(Cisnienie)
mean(Po-Przed)
t.test(Po-Przed,mu=0,alternative='less')

#TODO
# trzeba narysować rozkład statystyki testowej