推 dan0818:多謝老師 03/30 21:37
如果你還沒有R
請到http://cran.r-project.org/下載
選Windows (95 and later)->base->R-2.2.1-win32.exe即可!
------------------------------------------------------------------
你要輸入的資料請先用"記事本"打好
格式如下:
(株號) (冠區) (葉長)
number region length
1 1 15.7
... ... ...
葉長葉寬要分兩個檔案儲存
-------------------------------------------------------------------
使用R:
打開R後，先叫出檔案，檔案名稱請自己令(我先用我的說明)
> alnuslen=read.table("C:/Documents and Settings/睿涵/My Documents/涵報告/赤楊
葉長.txt",header=T)
(我令葉長的檔名為alnuslen，後面的路徑請依據你檔案位置填寫)
將檔案作anova分析
你可以令檔名，也可以不令
> len_aov=aov(length~number*region,data=alnuslen)
把anova結果叫出來
> anova(len_aov)
會出現下列結果
Analysis of Variance Table
Response: length
Df Sum Sq Mean Sq F value Pr(>F)
number 4 50399 12600 33.9888 < 2.2e-16 ***
region 4 159344 39836 107.4615 < 2.2e-16 ***
number:region 16 46551 2909 7.8484 < 2.2e-16 ***
Residuals 1475 546784 371
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
這裡要注意一件事
就是number及region的編號若是用阿拉伯數字
請先把它改成"character"
>fix(alnuslen)
在number及region的框框上按一下
將numeric改成character
再將視窗關掉即可
比較方便的方法就是使用英文字編號
這樣就不用再做更改
若number:region的p value太小即推翻
表二者有交感
接下來要找出誰跟誰有顯著差異(即p value<0.05者，若a=0.05)
這時候要用TukeyHSD做檢測
他會幫你每一項都做比較
> TukeyHSD(len_aov)
此時會出現一大串資料如下:
Tukey multiple comparisons of means
95% family-wise confidence level
Fit: aov(formula = length ~ number * region, data = alnuslen)
$"number:region"
diff lwr upr p adj
2:1-1:1 -8.9500000 -21.8331744 3.93317443 0.6604584
... ... ... ... ...
我們要看的是最後一欄的值
而上式所列為:第二株的第一區與第一株的第一區之差異，其p value為0.6604584
接下來我的作法是
利用excell幫我排序
找出同一棵樹不同冠區 及 同一冠區不同顆樹 的p value
如此可以得知誰和誰有顯著差異
葉寬的計算以此類推
然後
你就可以寫報告了^^
---------------------------------------------------------------------
得出的結果...
是有常識就可以知道的...orz
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 218.166.235.158