Porovnávanie priemerov v programovom prostredí R
 
 
Step: Ak si vyberieme nasledujúci druh testu dvojvzorkový Z test...
  Ak si vyberieme nasledujúci druh testu dvojvzorkový Z test
Potom príkaz bude vyzerať nasledovne:
library(TeachingDemos)
z.test(y,mu,sd)

Vysvetlivky k príkazom

y - kontinuálna premenná
mu - populačný priemer
sd - smerodajná odchýlka
Príklad
Napríklad ak chceme zistiť, či priemerná pôrodná hmotnosť novorodencov vo vybranej pôrodnici za uplynulý kalendárny rok bola 2,5 kg, pričom vieme, že smerodajná odchýlka údajov o pôrodnej hmotnosti bola 0,75 kg, môžeme využiť jednovzorkový z-test:
library(TeachingDemos)
z.test(porodna_hmotnost, mu=2.5, sd=0.75)


H0: priemerná pôrodná hmotnosť novorodencov = 2,5
HA: priemerná pôrodná hmotnosť novorodencov ≠ 2,5

Priemerná pôrodná hmotnosť novorodencov vo vybranej pôrodnici za uplynulý kalendárny rok bola štatisticky významne odlišná od stanovenej hodnoty pôrodnej hmotnosti (2,5 kg) na hladine významnosti 0,05 (prijímame alternatívnu hypotézu, pretože p< 0,05). Priemerná pôrodná hmotnosť novorodencov bola 2,8 kg.
 
 
Step: Ak si vyberieme nasledujúci druh testu jednovzorkový T test...
  Ak si vyberieme nasledujúci druh testu jednovzorkový T test
Potom príkaz bude vyzerať nasledovne:
t.test(y,mu)

Vysvetlivky k príkazom

y - kontinuálna premenná
mu - populačný priemer
Príklad
Napríklad ak chceme zistiť, či priemerné trvanie expozície ionizujúcim žiarením u vybraných respondentov bolo 32 minút, pričom nepoznáme smerodajnú odchýlku údajov o expozícii, môžeme využiť jednovzorkový t-test:
t.test(trvanie_expozicie_minuty, mu=32)



H0: priemerné trvanie expozície ionizujúcim žiarením = 32
HA: priemerné trvanie expozície ionizujúcim žiarením ≠ 32

Priemerné trvanie expozície ionizujúcim žiarením u vybraných respondentov nebolo štatisticky významne odlišné od stanovenej hodnoty (32 minút) na hladine významnosti 0,05. Prijímame teda nulovú hypotézu (p>0,05).
 
 
Step: Ak si vyberieme nasledujúci druh testu dvojvzorkový T test...
  Ak si vyberieme nasledujúci druh testu dvojvzorkový T test
Potom príkaz bude vyzerať nasledovne:
t.test(y~x)

alebo
t.test(y1,y2)


Vysvetlivky k príkazom

y - kontinuálna premenná
x - kategorická premenná s 2 kategóriami
Príklad
Je rozdiel v pôrodnej hmotnosti medzi chlapcami a dievčatami vo vybranej pôrodnici za uplynulý kalendárny rok?
t.test(porodna_hmotnost~pohlavie)



H0: nie je štatisticky významný rozdiel v priemernej pôrodnej hmotnosti medzi chlapcami a dievčatami.
HA: je štatisticky významný rozdiel v priemernej pôrodnej hmotnosti medzi chlapcami a dievčatami

Hypotézu sme testovali na hladine významnosti 0,05, pričom sme nezistili štatisticky významný rozdiel v priemernej pôrodnej hmotnosti medzi chlapcami a dievčatami (prijímame nulovú hypotézu, pretože p>0,05). Priemerná pôrodná hmotnosť chlapcov bola 2,86 kg a dievčat 2,92 kg.
 
 
Step: Ak si vyberieme nasledujúci druh testu párový T test...
  Ak si vyberieme nasledujúci druh testu párový T test
Potom príkaz bude vyzerať nasledovne:
t.test(y1,y2,paired=TRUE)

Vysvetlivky k príkazom

y - kontinuálna premenná
Príklad
Zníži sa hladina glukózy v krvi diabetikov po diétnych opatreniach?
t.test(glukozaPred, glukozaPo, paired=T)
boxplot(glukozaPred, glukozaPo)



H0: neexistuje štatisticky významný rozdiel v priemernej hladine glukózy v krvi respondentov pred a po zavedení diétnych opatrení
HA: existuje štatisticky významný rozdiel v priemernej hladine glukózy v krvi respondentov pred a po zavedení diétnych opatrení

Priemerná hladina glukózy v krvi respondentov sa štatisticky významne znížila po zavedení diétnych opatrení. Kým priemerná hladina glukózy pred diétou bola 8,5 mmol/l, po zavedení diétnych opatrení sa znížila na 6,6 mmol/l. Prijímame alternatívnu hypotézu (p<0,05). Graficky môžeme tento rozdiel vyjadriť s využitím boxplotu.
 
 
Step: Ak si vyberieme nasledujúci druh testu ANOVA...
  Ak si vyberieme nasledujúci druh testu ANOVA
Zníži sa hladina glukózy v krvi diabetikov po diétnych opatreniach?
anova(lm(y~x))

alebo
summary(aov(y~x))
pairwise.t.test(y,x)


Vysvetlivky k príkazom

y - kontinuálna premenná
x - kategorická premenná s 3 a viac kategóriami
Príklad
Existuje štatisticky významný rozdiel v priemernej pôrodnej hmotnosti novorodencov vo vybranej pôrodnici za uplynulý kalendárny rok podľa vzdelania ich matiek?
summary(aov(porodna_hmotnost~vzdelanie_matky))


H0: neexistuje štatisticky významný rozdiel v priemernej pôrodnej hmotnosti novorodencov podľa vzdelania ich matiek
HA: existuje štatisticky významný rozdiel v priemernej pôrodnej hmotnosti novorodencov aspoň medzi 2 skupinami matiek podľa ich vzdelania

Pri testovaní hypotézy na hladine významnosti 0,05 sme zistili, že neexistuje štatisticky významný rozdiel v priemernej pôrodnej hmotnosti novorodencov vo vybranej pôrodnici za uplynulý kalendárny rok podľa vzdelania matiek. Prijali sme nulovú hypotézu (p<0,05).