Matma_AI_cyber/Projekt_1/T-student-multiple-paired.R
2022-05-10 17:26:35 +02:00

47 lines
1.2 KiB
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')
#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')