108 年 國立臺灣大學農藝系暨研究所生物統計甲組《試驗設計學》
第 1 題
某研究員測試四種不同製劑(A、B、C、D)的生物活性並獲得部分變方分析表如下。
| 變異來源 | 自由度 | 平方和 |
|---|---|---|
| 製劑 | 3 | 64.42 |
| 誤差 | 26 | 62.12 |
| 總和 | 29 | 126.54 |
(1) 說明該研究員採用的試驗設計為何?此試驗共採用多少試驗單位?以 、、、 代表四種製劑的重複數,說明這些試驗單位應如何進行適當的隨機配置。(10分)
(2) 根據上列部分變方分析表提供的資訊,在顯著水準 下進行 檢定以比較四種製劑的生物活性是否相同 ()。(5分)
(3) 四種製劑中 A、B 屬於廠商提供的製劑,C、D 屬於實驗室自行開發的製劑。若分別以 、、、 代表四種製劑的平均生物活性,以 、、、 代表四種製劑的重複數,列式並提供適當說明應如何比較「廠商提供的製劑」與「實驗室自行開發的製劑」在顯著水準 下是否有顯著差異。(10分)
登入後即可作答並保存紀錄。
核心觀念
本題考查:
- 由變方分析表判斷試驗設計與試驗單位數。
- 單因子完全隨機設計的變方分析與 檢定。
- 以事前規劃的線性對比,比較兩組製劑的平均生物活性。
單因子完全隨機設計中,若有 個處理、總試驗單位數為 ,則:
(1) 試驗設計、試驗單位與隨機配置
試驗設計
製劑是唯一的試驗因子,共有四個處理水準:
由於變方分析表只有「製劑」與「誤差」兩個來源,沒有區集、批次或其他阻隔因素,因此採用:
單因子完全隨機設計(Completely Randomized Design, CRD)
其統計模式為:
其中 通常假設彼此獨立,且服從:
試驗單位數
由總和自由度:
因此:
也可由誤差自由度驗算:
所以本試驗共採用 30 個試驗單位。
四種製劑的重複數須滿足:
且每個 均為正整數。
適當的隨機配置
先將 30 個性質相近的試驗單位編號為 ,再將其中:
- 個配置給製劑 A;
- 個配置給製劑 B;
- 個配置給製劑 C;
- 個配置給製劑 D。
配置方式應採完全隨機,例如將四種製劑標籤依各自重複數寫入 30 張籤中,充分混合後隨機抽取,或使用亂數表、電腦亂數進行隨機分派。
每一個試驗單位接受 A、B、C、D 的機率分別為:
隨機配置的目的,是避免研究員主觀選擇試驗單位所造成的系統性偏差。
(2) 四種製劑平均生物活性是否相同
假設
欲檢定四種製劑的平均生物活性是否相同:
計算均方
製劑均方為:
誤差均方為:
因此檢定統計量:
判定
題目給定:
由於:
故拒絕虛無假設 。
結論
在顯著水準 下,有充分證據顯示四種製劑的平均生物活性並不完全相同。此檢定只能說明至少有一組平均數不同,不能直接指出是哪幾組不同;若要進一步判定,須進行適當的多重比較。
(3) 比較廠商製劑與實驗室製劑
核心觀念
A、B 屬於廠商製劑,C、D 屬於實驗室製劑。欲比較的是兩個群組的平均生物活性,因此應使用線性對比,而不是直接比較四個處理的個別平均數。
廠商製劑的加權平均為:
實驗室製劑的加權平均為:
加權是必要的,因為四種製劑的重複數不一定相同。
建立對比
定義欲檢定的母體對比為:
其估計量為:
檢定假設為:
第 2 題
某研究室受託調查 15 處森林的土壤化學成分,這 15 處森林依據樹木的物種分為 (A) 20-30% (B) 30-50% (C) 50-80% 三種森林覆蓋程度,每種程度各 5 重複。由於研究時間緊迫,研究室派出 5 位研究生各自負責 3 處森林採樣取回土壤樣本一份,並自行對採樣進行化學分析。5 位研究生包含經驗豐富的博士班學生及剛進入研究室的碩士班一年級學生,對於化學分析的熟練程度不同。
(1) 針對上述情境,說明應如何指派 5 位研究生採樣的森林覆蓋程度種類及順序,以適當控制森林覆蓋程度及研究生造成的試驗誤差。(10分)
(2) 根據題(1)的分派方式,列出適當的變方分析表中的「變異來源」(SOV) 和「自由度」(DF),以分析比較三種不同森林覆蓋程度的土壤化學組成是否有顯著差異。(10分)
(3) 若每個森林採樣取回「兩份」土壤樣本時,題(2)的變方分析表應該如何改寫?(5分)
登入後即可作答並保存紀錄。
核心觀念
本題同時考查:
- 隨機完全區集設計(Randomized Complete Block Design, RCBD)
- 區集效應控制
- 採樣順序的隨機化與平衡
- 重複樣本與實驗單位的區分
- 變方分析表的自由度配置
三種森林覆蓋程度是研究處理:
- A:20–30%
- B:30–50%
- C:50–80%
5 位研究生的採樣及化學分析熟練程度不同,因此「研究生」會造成額外試驗誤差,應將研究生作為區集。每位研究生負責 3 處森林,且每種覆蓋程度各採樣 1 處,正好形成一個完整區集。
(1)森林覆蓋程度及採樣順序的指派
1. 以研究生作為區集
令第 位研究生為第 個區集,每位研究生均採集:
- 1 處 A 類森林
- 1 處 B 類森林
- 1 處 C 類森林
因此,每一位研究生都接觸三種覆蓋程度,研究生的經驗差異便不會集中在某一種覆蓋程度。
例如可安排如下:
| 研究生 | 第 1 處 | 第 2 處 | 第 3 處 |
|---|---|---|---|
| 1 | A | B | C |
| 2 | B | C | A |
| 3 | C | A | B |
| 4 | A | C | B |
| 5 | C | B | A |
實際執行時,應將各覆蓋程度內的實際森林位置以隨機方式分派給研究生,避免森林地理位置、交通距離或其他環境因素與研究生混淆。
2. 採樣順序應隨機化或予以平衡
不可讓所有研究生都依照相同順序採樣,例如全部採用 A → B → C。因為研究生可能在採樣過程中產生:
- 疲勞效應
- 熟練度逐漸提高的學習效應
- 時間延遲造成的環境差異
- 儀器或分析流程的時間效應
因此,三種覆蓋程度的採樣順序應在研究生之間隨機化,並盡量使 A、B、C 分散於第 1、2、3 個採樣位置。
上述安排中:
- 第 1 次採樣:A、B、C、A、C
- 第 2 次採樣:B、C、A、C、B
- 第 3 次採樣:C、A、B、B、A
各覆蓋程度在不同採樣位置皆有分布,可降低採樣順序效應與覆蓋程度效應的混淆。
3. 設計表示
此設計可表示為:
其中:
- :第 位研究生採集第 種覆蓋程度森林所得的土壤化學分析值
- :總平均
- :第 位研究生的區集效應
- :第 種森林覆蓋程度效應
- :隨機誤差
欲檢定的假設為:
亦即三種森林覆蓋程度的平均土壤化學組成相同。
對立假設為:
(2)變方分析表中的變異來源與自由度
共有:
筆觀測值,因此總自由度為:
在隨機完全區集設計中:
- 研究生區集自由度:
- 覆蓋程度自由度:
- 誤差自由度:
因此變方分析表如下:
| 變異來源(SOV) | 自由度(DF) |
|---|---|
| 研究生(區集) | |
| 森林覆蓋程度 | |
| 誤差 | |
| 總和 |
自由度驗算:
覆蓋程度的檢定統計量為:
其中:
第 3 題
某農藝學家擬探討不同培養介質對洋菇產量及品質的影響,試驗將於 3 間菇室中進行;每間菇室內有 4 座床架,每座床架由上至下可擺置 4 席面積為 6 平方公尺的菇床,而擬參試的培養介質共有 4 種配方。根據所給的條件:
(1) 請幫忙設計一個重複拉丁方來進行此試驗,詳細說明試驗進行的步驟。(10分)
(2) 請寫出其ANOVA 分析的統計線型模式(statistical linear model),及ANOVA表中的SOV和DF。(10分)
登入後即可作答並保存紀錄。
核心觀念
本題考查「重複拉丁方設計(replicated Latin square design)」。
有 4 種培養介質,因此每一個拉丁方為 :
- 列因子:菇室內的 4 座床架;
- 行因子:床架由上至下的 4 個位置;
- 處理因子:4 種培養介質;
- 重複因子:3 間菇室。
拉丁方的基本要求是:每一種培養介質在每一列及每一行中均恰好出現一次。如此可同時控制床架差異與上下位置差異,提高處理比較的精確度。
每間菇室有 個菇床,三間共 個實驗單位;每種培養介質共使用 個菇床。
(1) 重複拉丁方設計與試驗步驟
第一步:界定實驗單位與因子
每一席面積 平方公尺的菇床為一個實驗單位。
配置如下:
- 重複:菇室,共 個;
- 列:每間菇室的床架,共 列;
- 行:每座床架由上至下的位置,共 行;
- 處理:培養介質配方 。
第二步:建立三個 拉丁方
可採用下列配置。表中每個字母代表一種培養介質:
菇室一
| 床架\位置 | 1 | 2 | 3 | 4 |
|---|---|---|---|---|
| 1 | A | B | C | D |
| 2 | B | C | D | A |
| 3 | C | D | A | B |
| 4 | D | A | B | C |
菇室二
| 床架\位置 | 1 | 2 | 3 | 4 |
|---|---|---|---|---|
| 1 | A | C | D | B |
| 2 | C | D | B | A |
| 3 | D | B | A | C |
| 4 | B | A | C | D |
菇室三
| 床架\位置 | 1 | 2 | 3 | 4 |
|---|---|---|---|---|
| 1 | B | D | A | C |
| 2 | D | A | C | B |
| 3 | A | C | B | D |
| 4 | C | B | D | A |
每一個菇室內,每種介質在每座床架及每個上下位置均出現一次。
實際試驗時, 應以隨機方式對應到四種實際培養介質配方;三個菇室也應各自獨立進行拉丁方隨機化。
第三步:執行隨機化
每間菇室的隨機化可依下列順序進行:
- 將四種培養介質以隨機方式標記為 。
- 隨機排列四座床架,作為拉丁方的四列。
- 隨機排列床架由上至下的四個位置,作為拉丁方的四行。
- 將隨機化後的拉丁方配置於該菇室的 席菇床。
- 三間菇室各使用一個拉丁方,並以菇室作為三個重複。
第四步:管理試驗條件
除培養介質外,菇種、播種量、含水量、溫度、濕度、通風、光照、培養時間及採收標準應盡量一致。
各菇床應記錄:
- 產量,例如每席菇床總產量;
- 品質指標,例如菇體大小、商品率、硬度或其他預先訂定的品質測定值。
產量與品質屬於不同反應變數,通常分別進行 ANOVA 分析。
(2) 統計線型模式、SOV 與 DF
統計線型模式
令:
- :第 間菇室、第 座床架、第 個上下位置之觀測值;
- :總平均;
- :第 間菇室效果;
- :第 間菇室內第 座床架效果;
- :第 個上下位置效果;
- :配置於該菇床的第 種培養介質效果;
第 4 題
某農藝學家擬在 RCBD 下,執行兩重複 r=2 的 A,B 兩因子試驗,其試級(level)數分別為 a=2 及 b=3。下列有三個田間配置圖:
配置圖 A:
Block I
Block II
A1B1
A1B2 A1B3
A2B2
A2B3 A2B1
A2B2
A2B1 A2B3
A1B1
A1B2 A1B3
配置圖 B:
Block I
Block II
A1B1
A1B2 A1B3
A2B1
A2B2 A2B3
A1B3
A1B1 A1B2
A2B3
A2B1 A2B2
配置圖 C:
Block I
Block II
A2B1
A1B1 A1B3
A2B3 A2B2 A1B2
A1B2 A1B2 A2B1
A1B1
A2B3 A2B3
(1) 請選擇複因子設計(factorial design)的田間配置圖,並寫出其ANOVA 分析的統計線型模式(statistical linear model);及ANOVA表中的SOV和DF。(10分)
(2) 請選擇裂區設計(split-plot design)的田間配置圖,並寫出其ANOVA 分析的統計線型模式(statistical linear model);及ANOVA表中的SOV和DF。(10分)
(3) 請選擇條區設計(strip-plot design)的田間配置圖,並寫出其ANOVA 分析的統計線型模式(statistical linear model);及ANOVA表中的SOV和DF。(10分)
登入後即可作答並保存紀錄。
核心觀念
本題比較三種二因子試驗的田間配置:
- 複因子設計(factorial design):每一個區集內,、 的所有處理組合各出現一次,且處理組合可完全隨機配置。
- 裂區設計(split-plot design):先將主因子配置於較大主區,再將次因子配置於主區內的次區。因此同一主因子水準的各個 水準必須形成同一主區內的次區。
- 條區設計(strip-plot design): 與 分別以兩組互相垂直的條帶配置,兩因子的條帶交會處形成一個試驗小區。
本題有:
總自由度為:
解題方法與配置圖判斷
配置圖 B:複因子設計
配置圖 B 中,每個區集都包含完整的六個處理組合:
而且各處理組合在區集內各出現一次,符合 RCBD 下複因子試驗的配置原則。
因此:
配置圖 A:裂區設計
配置圖 A 中,在每一個區集內, 的同一水準之下配置 三個次區。例如:
因此先配置 主區,再於每一個 主區內配置 次區,符合裂區設計。
配置圖 C:條區設計
配置圖 C 的配置方式是將 的不同水準形成一組條帶, 的不同水準形成另一組與其交叉的條帶;每一個 條帶與 條帶的交會處形成一個觀測小區。
因此, 與 都不是「主區—次區」的階層關係,而是兩種條帶相交的配置,符合條區設計。
(1)複因子設計:配置圖 B
統計線型模式
令:
- :第 個區集、 的第 水準、 的第 水準之觀測值;
- ;
- ;
- 。
RCBD 下複因子設計的模型為:
其中:
- :總平均;
- :區集效果;
- :因子 效果;
- :因子 效果;
- : 交互作用;
- :隨機誤差,通常假設獨立且服從 。
ANOVA 表
| 變異來源 SOV | 自由度 DF |
|---|---|
| 區集 Block | |
| 誤差 Error | |
| 總和 Total |
因此複因子設計中,、 及 均以共同的實驗誤差項檢定。
(2)裂區設計:配置圖 A
統計線型模式
裂區設計中, 配置於主區, 配置於主區內的次區。因此必須有兩層誤差:
- 主區誤差:區集與 的交互作用;
- 次區誤差:次區層級的隨機誤差。
模型為:
其中:
- :主區誤差,代表區集內 主區的隨機變異;
- :次區誤差。
主區誤差可表示為:
因此也可寫成:
第 4 題
某農藝學家擬在 RCBD 下,執行兩重複 r=2 的 A,B 兩因子試驗,其試級(level)數分別為 a=2 及 b=3。下列有三個田間配置圖:
配置圖 A:
Block I
Block II
A1B1
A1B2 A1B3
A2B2
A2B3 A2B1
A2B2
A2B1 A2B3
A1B1
A1B2 A1B3
配置圖 B:
Block I
Block II
A1B1
A1B2 A1B3
A2B1
A2B2 A2B3
A1B3
A1B1 A1B2
A2B3
A2B1 A2B2
配置圖 C:
Block I
Block II
A2B1
A1B1 A1B3
A2B3 A2B2 A1B2
A1B2 A1B2 A2B1
A1B1
A2B3 A2B3
(1) 請選擇複因子設計(factorial design)的田間配置圖,並寫出其ANOVA 分析的統計線型模式(statistical linear model);及ANOVA表中的SOV和DF。(10分)
(2) 請選擇裂區設計(split-plot design)的田間配置圖,並寫出其ANOVA 分析的統計線型模式(statistical linear model);及ANOVA表中的SOV和DF。(10分)
(3) 請選擇條區設計(strip-plot design)的田間配置圖,並寫出其ANOVA 分析的統計線型模式(statistical linear model);及ANOVA表中的SOV和DF。(10分)
登入後即可作答並保存紀錄。
核心觀念
本題比較三種二因子試驗的田間配置:
- 複因子設計(factorial design):每一個區集內,、 的所有處理組合各出現一次,且處理組合可完全隨機配置。
- 裂區設計(split-plot design):先將主因子配置於較大主區,再將次因子配置於主區內的次區。因此同一主因子水準的各個 水準必須形成同一主區內的次區。
- 條區設計(strip-plot design): 與 分別以兩組互相垂直的條帶配置,兩因子的條帶交會處形成一個試驗小區。
本題有:
總自由度為:
解題方法與配置圖判斷
配置圖 B:複因子設計
配置圖 B 中,每個區集都包含完整的六個處理組合:
而且各處理組合在區集內各出現一次,符合 RCBD 下複因子試驗的配置原則。
因此:
配置圖 A:裂區設計
配置圖 A 中,在每一個區集內, 的同一水準之下配置 三個次區。例如:
因此先配置 主區,再於每一個 主區內配置 次區,符合裂區設計。
配置圖 C:條區設計
配置圖 C 的配置方式是將 的不同水準形成一組條帶, 的不同水準形成另一組與其交叉的條帶;每一個 條帶與 條帶的交會處形成一個觀測小區。
因此, 與 都不是「主區—次區」的階層關係,而是兩種條帶相交的配置,符合條區設計。
(1)複因子設計:配置圖 B
統計線型模式
令:
- :第 個區集、 的第 水準、 的第 水準之觀測值;
- ;
- ;
- 。
RCBD 下複因子設計的模型為:
其中:
- :總平均;
- :區集效果;
- :因子 效果;
- :因子 效果;
- : 交互作用;
- :隨機誤差,通常假設獨立且服從 。
ANOVA 表
| 變異來源 SOV | 自由度 DF |
|---|---|
| 區集 Block | |
| 誤差 Error | |
| 總和 Total |
因此複因子設計中,、 及 均以共同的實驗誤差項檢定。
(2)裂區設計:配置圖 A
統計線型模式
裂區設計中, 配置於主區, 配置於主區內的次區。因此必須有兩層誤差:
- 主區誤差:區集與 的交互作用;
- 次區誤差:次區層級的隨機誤差。
模型為:
其中:
- :主區誤差,代表區集內 主區的隨機變異;
- :次區誤差。
主區誤差可表示為:
因此也可寫成: