Fisheri testi vastus
Märkus! R Fisheri test leiab sansside suhte 0 vs 1! Kuigi loogiline oleks vastupidi. Kuna aga jagatis on nagunii pöördvõrdeline, siis on võimalik tulemus lihtsalt pöördväärtust võttes ümber arvutada, alljärgneva näite puhul nii ka teeme.
|
cop=as.numeric(I(Cop2>0))
t1=table(cop,Osa);t1 Osa cop Lämmijärv Suurjärv 0 6 15 1 13 4 fisher.test(t1) Fisher's Exact Test for Count Data data: t1 p-value = 0.008138 alternative hypothesis: true odds ratio is not equal to 1 95 percent confidence interval: 0.02146905 0.64249844 sample estimates: odds ratio 0.1313954 1/0.1313954 [1] 7.610617 |
Teisendame arvukused 0-1 tunnuseks-kui esineb, siis on väärtus 1.
Sagedustabel, millest silma järgi võib ütelda, et uuritavat liiki tuleb proovianumasse sagedamini Lämmijärves. Fisheri test . Seos on oluline, šansside suhe ei ole 1. Cop1 mitteleidumise (vaata märkus ülal) šanss Lämmijärves on 0.1313954 korda suurem kui Suurjärves, mis on üpris kohmakas väide. Aga kui võtta sellest pöördväärtus (so 7,61) ja ütelda nii, et liigi Cop2 leidumise šanss Lämmijärves on umber 8 korda suurem, on asi arusaadavam. Ka 95% usalduspiirid peab ümber arvutama! |
t-testiga
t.test(Cop2~Osa)
Welch Two Sample t-test
data: Cop2 by Osa
t = 3.4677, df = 20.542, p-value = 0.002358
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
1.850184 7.412974
sample estimates:
mean in group Lämmijärv mean in group Suurjärv
5.2631579 0.6315789
Lämmijärves proovides on keskmiselt 5,26-0.63=4.63 liigi Cop2 isendit rohkem.
Welch Two Sample t-test
data: Cop2 by Osa
t = 3.4677, df = 20.542, p-value = 0.002358
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
1.850184 7.412974
sample estimates:
mean in group Lämmijärv mean in group Suurjärv
5.2631579 0.6315789
Lämmijärves proovides on keskmiselt 5,26-0.63=4.63 liigi Cop2 isendit rohkem.