s444455/Matma_AI_cyber
1
1
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

student dwie

Browse Source
This commit is contained in:
Jakub Adamski 2022-05-09 12:34:50 +02:00
parent 11f948c114
commit fa03e4d6fe
1 changed files with 24 additions and 16 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

38
Projekt_1/T-studnet-multiple.R
View File

@ -1,26 +1,34 @@
t.test_knownvar <- function(x, y, V1, V2, m0 = 0, alpha = 0.05, alternative = "two.sided") { t.test_multiple <- function(x, y, m0=0, alternative = "two.sided") {
M1 <- mean(x) M1 <- mean(x)
M2 <- mean(y) M2 <- mean(y)
n1 <- length(x) n1 <- length(x)
n2 <- length(y) n2 <- length(y)
sigma1 <- sqrt(V1) S <- sqrt((var(x) / n1) + (var(y) / n2))
sigma2 <- sqrt(V2)
S <- sqrt((V1 / n1) + (V2 / n2))
statistic <- (M1 - M2 - m0) / S statistic <- (M1 - M2 - m0) / S
p <- if (alternative == "two.sided") { p <- if (alternative == "two.sided") {
2 * pnorm(abs(statistic), lower.tail = FALSE) 2 * pt(q=abs(statistic), df=n1+n2-2, lower.tail=FALSE)
} else if (alternative == "less") { } else if (alternative == "less") {
pnorm(statistic, lower.tail = TRUE) pt(q=statistic, df=n1+n2-2, lower.tail = TRUE)
} else { } else {
pnorm(statistic, lower.tail = FALSE) pt(q=statistic, df=n1+n2-2, lower.tail = FALSE)
} }
LCL <- (M1 - M2 - S * qnorm(1 - alpha / 2))
UCL <- (M1 - M2 + S * qnorm(1 - alpha / 2)) value <- list(mean1 = M1, mean2 = M2, df=n1+n2-2, m0 = m0, statistic = statistic, p.value = p, alternative = alternative)
value <- list(mean1 = M1, mean2 = M2, m0 = m0, sigma1 = sigma1, sigma2 = sigma2, S = S, statistic = statistic, p.value = p, LCL = LCL, UCL = UCL, alternative = alternative)
# print(sprintf("P-value = %g",p))
# print(sprintf("Lower %.2f%% Confidence Limit = %g",
# alpha, LCL))
# print(sprintf("Upper %.2f%% Confidence Limit = %g",
# alpha, UCL))
return(value) return(value)
} }
#dwie próby zależne
load("Cisnienie.RData")
attach(Cisnienie)
t.test(x=Po, y=Przed, alternative='less')
t.test_multiple(x=Po, y=Przed, alternative='less')