2022-04-27 12:37:59 +02:00
|
|
|
library('MKinfer')
|
|
|
|
|
library(boot)
|
|
|
|
|
|
2022-04-25 00:05:17 +02:00
|
|
|
#test t-studenta dla jednej próby
|
|
|
|
|
load("Hamulce.RData")
|
|
|
|
|
attach(Hamulce)
|
|
|
|
|
mean(Wynik)
|
|
|
|
|
t.test(Wynik,mu=18.6,alternative='less')
|
2022-05-09 17:33:51 +02:00
|
|
|
boot.t.test(Wynik,mu=18.6,alternative='less', R = 200)
|
2022-04-27 12:37:59 +02:00
|
|
|
|
|
|
|
|
par(mfrow=c(1,2))
|
|
|
|
|
|
|
|
|
|
hist(Wynik,prob=T,main='')
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
meanboot <- function(x,i)mean(x[i])
|
|
|
|
|
bmean=boot(Wynik,meanboot,100)
|
|
|
|
|
|
|
|
|
|
hist(bmean$t,prob=T,main='')
|
|
|
|
|
|
|
|
|
|
t.test(bmean$t,mu=18.6,alternative='less')
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2022-04-25 00:05:17 +02:00
|
|
|
|
|
|
|
|
# dwie próby niezależne
|
|
|
|
|
load("Bledy.RData")
|
|
|
|
|
attach(Bledy)
|
|
|
|
|
tapply(Blad,Kod,mean)
|
|
|
|
|
tapply(Blad,Kod,var)
|
|
|
|
|
var.test(Blad~Kod)
|
|
|
|
|
t.test(Blad~Kod,var.equal=T,alternative='less')
|
|
|
|
|
|
2022-05-08 23:55:40 +02:00
|
|
|
|
|
|
|
|
|
2022-04-25 00:05:17 +02:00
|
|
|
#dwie próby zależne
|
|
|
|
|
load("Cisnienie.RData")
|
|
|
|
|
attach(Cisnienie)
|
|
|
|
|
mean(Po-Przed)
|
|
|
|
|
t.test(Po-Przed,mu=0,alternative='less')
|
|
|
|
|
|
|
|
|
|
#TODO
|
|
|
|
|
# trzeba narysować rozkład statystyki testowej
|
