Статистика
Простые статистические тесты
В этом тьюториале приведены самые простые тесты, данные считыванием так:
chick.w <- read.table("ChickWeight.tab", header=T) sleep.paired <- read.table("sleep.paired.tab“,header=T) rt <- read.table("anaesthetic.reaction.time.tab", sep="\t", header=T)
Файлы тут: http://makarich.fbb.msu.ru/artemov/R/ChickWeight.tab
http://makarich.fbb.msu.ru/artemov/R/sleep.paired.tab
http://makarich.fbb.msu.ru/artemov/R/anaesthetic.reaction.time.tab
Сравнение эмпирического распределения с теоретическим:
- тест на нормальность
ks.test(w.diet.3,"pnorm",mean(w.diet.3),sd(w.diet.3))
Сравнение распределений двух выборок:
ks.test(w.diet.2,w.diet.3)
Тест Shapiro-Wilk (на нормальность):
shapiro.test(w.diet.3)
t-test для независимых выборок
chick.test <- t.test(w.diet.2, w.diet.3, alternative="less")
И для парных:
sleep.test <- t.test(sleep.paired$Drug1, + sleep.paired$Drug2, paired=T, alternative=“less")
U-критерий Манна-Уитни
wilcox.test(rt$Mean.RT~rt$Group, alternative="greater")
Парный:
wilcox.test(test$Q1, test$Q2, paired=T, alternative="greater", exact=F)