library('MKinfer')
library(boot)

t.test_multiple_paried <- function(x, m0=0, alternative = "two.sided") {
  M <- mean(x)
  n <- length(x)
  #warianncja
  V <- var(x)
  sigma <- sqrt(V)
  S <- sqrt(V / n)
  statistic <- (M - m0) / S
  p <- if (alternative == "two.sided") {
    2 * pt(q=abs(statistic), df=length(x)-1, lower.tail=FALSE)
  } else if (alternative == "less") {
    pt(q=statistic, df=length(x)-1, lower.tail = TRUE)
  } else {
    pt(q=statistic, df=length(x)-1, lower.tail = FALSE)
  }
  
  value <- list(mean = M, m0 = m0, df=length(x)-1, statistic = statistic, p.value = p,  alternative = alternative)
  return(value)
} 


#dwie próby zależne
load("Cisnienie.RData")
attach(Cisnienie)
mean(Po-Przed)

#t.test(Po-Przed,mu=0,alternative='less')
size <- length(Cisnienie)

#dane zwykłe 
# funkcja wbudowana
t.test(Cisnienie$Po-Cisnienie$Przed,mu=0,alternative='less')
#nasza funkcja
t.test_multiple_paried(Cisnienie$Po-Cisnienie$Przed,m0=0,alternative='less')

#dane zbootstrpowane
bootstrappedData_Cisnienie <- ourBoot(Cisnienie$Po-Cisnienie$Przed, 200, size)
# funkcja wbudowana
t.test(bootstrappedData_Cisnienie,mu=0,alternative='less')
#nasza funkcja
t.test_multiple_paried(bootstrappedData_Cisnienie,m0=0,alternative='less')