Matma_AI_cyber/Projekt_1/T-student-multiple-paired.R

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')

wyniki <- data.frame()
for (i in 1:10) {
  bootstrappedData_Cisnienie <- ourBoot(Cisnienie$Po-Cisnienie$Przed, size, size)
  ttest <- t.test_multiple_paried(bootstrappedData_Cisnienie,m0=0,alternative='less')
  wyniki <- rbind(wyniki, ttest)
}

mean(wyniki$mean)
mean(wyniki$statistic)
mean(wyniki$p.value)

# funkcja wbudowana
#t.test(bootstrappedData_Cisnienie,mu=0,alternative='less')