diff --git a/Projekt_1/T-student-multiple-paired.R b/Projekt_1/T-student-multiple-paired.R index faea955..b7a0c2e 100644 --- a/Projekt_1/T-student-multiple-paired.R +++ b/Projekt_1/T-student-multiple-paired.R @@ -32,15 +32,23 @@ size <- length(Cisnienie) #dane zwykłe # funkcja wbudowana -t.test(Cisnienie$Po-Cisnienie$Przed,mu=0,alternative='less') +#t.test(Cisnienie$Po-Cisnienie$Przed,mu=0,alternative='less') #nasza funkcja -t.test_multiple_paried(Cisnienie$Po-Cisnienie$Przed,m0=0,alternative='less') +#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) -#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') +#t.test(bootstrappedData_Cisnienie,mu=0,alternative='less') + diff --git a/Projekt_1/T-studnet-multiple.R b/Projekt_1/T-student-multiple.R similarity index 59% rename from Projekt_1/T-studnet-multiple.R rename to Projekt_1/T-student-multiple.R index 4c72095..b08ea0c 100644 --- a/Projekt_1/T-studnet-multiple.R +++ b/Projekt_1/T-student-multiple.R @@ -26,21 +26,24 @@ attach(Bledy) size_Bledy <- length(Bledy$Blad) -bootstrappedData_Bledy_Blad <- ourBoot(Bledy$Blad, 200, size_Bledy) -bootstrappedData_Bledy_Kod <- ourBoot(Bledy$Kod, 200, size_Bledy) +wyniki <- data.frame() +for (i in 1:10){ + bootstrappedData_Bledy_Blad <- ourBoot(Bledy$Blad, size_Bledy, size_Bledy) + bootstrappedData_Bledy_Kod <- ourBoot(Bledy$Kod, size_Bledy, size_Bledy) + ttest <- t.test_multiple(x=bootstrappedData_Bledy_Blad, y=bootstrappedData_Bledy_Kod, alternative='less') + wyniki <- rbind(wyniki, ttest) +} - -#dane zbootstrapowane -#funkcja wbudowana -t.test(x=bootstrappedData_Bledy_Blad, y=bootstrappedData_Bledy_Kod,var.equal=T,alternative='less') -# nasza funkcja -t.test_multiple(x=bootstrappedData_Bledy_Blad, y=bootstrappedData_Bledy_Kod, alternative='less') +mean(wyniki$mean1) +mean(wyniki$mean2) +mean(wyniki$statistic) +mean(wyniki$p.value) #dane normalne #funkcja wbudowana -t.test(x=Blad, y=Kod,var.equal=T,alternative='less') +#t.test(x=Blad, y=Kod,var.equal=T,alternative='less') # nasza funkcja -t.test_multiple(x=Blad, y=Kod, alternative='less') +#t.test_multiple(x=Blad, y=Kod, alternative='less')