1.12 PRIMJER 4
Je li novi lijek značajno utjecao na smanjenje težine ispitanika?
Za 4 ispitanika 3 je puta izmjerena težina u kg (jednom na početku eksperimenta, jednom tijekom eksperimenta i jednom na kraju)
Subjekti | Prvo mjerenje | Drugo mjerenje | Treće mjerenje | ¯yj |
---|---|---|---|---|
Ispitanik 1 | 75 | 74 | 64 | 71 |
Ispitanik 2 | 86 | 85 | 78 | 83 |
Ispitanik 3 | 100 | 99 | 95 | 98 |
Ispitanik 4 | 95 | 86 | 83 | 88 |
¯yi | 89 | 86 | 80 | 85 |
SSBS=3[(71−85)2+(83−85)2+(98−85)2+(88−85)2]=1134SSBG=4[(89−85)2+(86−85)2+(80−85)2]=168SSWS=(75−71)2+(74−71)2+⋯+(100−98)2+⋯+(83−88)2=204SSE=204−168=36SST=1134+168+36=1338
- U konkretnom slučaju primjene testira se značajnost efekta ponavljanja, tj.
F2=SSBG(k−1)SSE(n−1)(k−1)=1682366=846=14
- Podatke iz Tablice 1.6 unesite ručno u R skriptu pomoću generičke funkcije
c()
te ih objedinite u jednom spremniku podataka (data frame) pod nazivomponovljena
kombinirajući naredbedata.frame()
icbind()
=c(1,1,1,2,2,2,3,3,3,4,4,4)
ispitanik=c(1,2,3,1,2,3,1,2,3,1,2,3)
mjerenje=c(75,74,64,86,85,78,100,99,95,95,86,83)
tezina=data.frame(cbind(ispitanik,mjerenje,tezina))
ponovljena ponovljena
## ispitanik mjerenje tezina
## 1 1 1 75
## 2 1 2 74
## 3 1 3 64
## 4 2 1 86
## 5 2 2 85
## 6 2 3 78
## 7 3 1 100
## 8 3 2 99
## 9 3 3 95
## 10 4 1 95
## 11 4 2 86
## 12 4 3 83
- Koristeći naredbu
aov()
provedite jednofaktorsku analizu varijance s ponovljenim mjerenjima, pri čemu je od interesa testirati efekt ponavljanja
=aov(tezina~factor(mjerenje)+Error(factor(ispitanik)),data=ponovljena)
tablica4summary(tablica4)
##
## Error: factor(ispitanik)
## Df Sum Sq Mean Sq F value Pr(>F)
## Residuals 3 1134 378
##
## Error: Within
## Df Sum Sq Mean Sq F value Pr(>F)
## factor(mjerenje) 2 168 84 14 0.0055 **
## Residuals 6 36 6
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1