【正文】 表54 棉花品種與密度試驗(yàn)的皮棉產(chǎn)量結(jié)果(斤/畝)品種 密度(株/畝) 　 3000 4000 5000 A1 102 140 180 130 162 172 A2 130 120 100 140 110 105 A3 105 120 125 92 115 135 A4 180 150 160 160 170 140 　　MINITAB程序:　　 兩向分組組合內(nèi)有重復(fù)觀察值資料　　SET C1　　102 140 180　　130 162 172　　130 120 100　　140 110 105　　105 120 125　　92 115 135　　。五、兩向分組組合內(nèi)有重復(fù)觀察值資料的統(tǒng)計(jì)分析　　〖例4〗4個(gè)棉花品種、3種種植密度,每處理組合有2個(gè)觀察值,其皮棉產(chǎn)量(斤/畝)列于表54所示?！　INITAB計(jì)算結(jié)果:　　FACTOR TYPE LEVELS VALUES　　C3 FIXED 4 1 2 3 4　　C2(C3) RANDOM 3 1 2 3　　ANALYSIS OF VARIANCE FOR C1　　SOURCE DF SS MS F P　　C3 3 　　C2(C3) 8 　　ERROR 12 　　TOTAL 23 　　ANOV命令不再自動(dòng)給出每個(gè)因子的各個(gè)處理的平均值和區(qū)間估計(jì)，若希望獲得這些數(shù)據(jù)，可以使用ANOV命令的子命令MEAN指定輸出哪個(gè)因子的各處理的平均值。不同的結(jié)果輸入順序，編碼的格式也不同。在后面的許多程序中都會(huì)出現(xiàn)這種編碼格式，通過多次的練習(xí)和細(xì)心體會(huì)，相信你能夠掌握這一技巧?！　AND C2.　　 End　　由于本試驗(yàn)中只有處理的區(qū)別，每個(gè)處理中小區(qū)是隨機(jī)設(shè)置的，且完全包含在每個(gè)處理內(nèi)部，所以選用了“處理 小區(qū)（處理）”的巢式模型，并設(shè)小區(qū)為固定因子。試進(jìn)行統(tǒng)計(jì)分析。　　結(jié)果表明：不同棉花品種單株結(jié)鈴數(shù)之間的F＝，,，所以品種間存在著顯著差異。試進(jìn)行統(tǒng)計(jì)分析。由方差分析之后的各處理產(chǎn)量的平均值及其區(qū)間估計(jì)結(jié)果，初步可以看出處理4的產(chǎn)量顯著地高于處理處理5的產(chǎn)量。表51 5種不同的氮肥用量對(duì)小麥籽粒產(chǎn)量的影響N1 N2 N3 N4 N5 22 22 24 24 22 23 24 25 27 23 20 22 23 24 21 20 20 22 23 21 　　MINITAB程序　?。?單向分組組內(nèi)樣本數(shù)相等　　READ C1C5　　22 22 24 24 22　　23 24 25 27 23　　20 22 23 24 21　　20 20 22 23 21　　END　　AOVO C1C5　　LET K1 = 　　LET K2=SQRT(K1/N(C1))　　 End　　MINITAB計(jì)算結(jié)果:　　ANALYSIS OF VARIANCE　　SOURCE DF SS MS F P　　FACTOR 4 　　ERROR 15 　　TOTAL 19 　　INDIVIDUAL 95 PCT CI39。二、單向分組組內(nèi)樣本數(shù)相等資料的方差分析　　〖例1〗用盆栽試驗(yàn)研究不同氮肥用量對(duì)小麥的效應(yīng),每種氮肥用量種4盆,完全隨機(jī)試驗(yàn)設(shè)計(jì),成熟時(shí)測(cè)定每盆小麥籽粒重(g),列于表51。　　在說明模型時(shí)，可以用豎線“|”簡化交互式說明,如A|B代表A和B的主效應(yīng)和其交互效應(yīng)三項(xiàng)；用減號(hào)“”表示刪除某項(xiàng)?！　≡贏NOVA命令中可以使用如下的模型：　　三因子交互式: ANOVA Y=A B C A*B A*C B*C A*B*C　　三因子巢式: ANOVA Y=A B(A) C(A B)　　交互和巢式: ANOVA Y=A B(A) C A*C B*C(A)　　在上面的命令中，等號(hào)后面為選用的模型，因素的交互作用以“*”號(hào)說明,如A*B為A和B兩因素的交互；以括號(hào)說明巢式，如B(A)為B包括在A內(nèi)，C(A B)為C包括在A和B中?！　。ㄈ┒嘁蛩胤讲罘治觥　∮梅? ANOVA 模型　　常用子命令: RANDOM, MEANS　　功能: 進(jìn)行多因素平衡設(shè)計(jì)或單因素的不平衡設(shè)計(jì)的方差分析。該命令的輸出結(jié)果只包括方差分析表。　?。ǘ﹥上蚍纸M資料的方差分析　　用法: TWOWAYAOV C, C(行), C(列)　　TWOWAYAOV C, SUBS C, C [RESIDUALS C [FITS C]]　　子命令: ADDITIVE, MEAN　　功能: 進(jìn)行平衡設(shè)計(jì)的兩向分組資料的方差分析。　　說明: ONEWAYAOV和AOVONEWAY的功能是一樣的,它要求所有的試驗(yàn)數(shù)據(jù)放在一列中,在另一列中存放試驗(yàn)的處理水平。AOVONEWAY 可以作樣本不相等資料的方差分析?！　≌f明：每個(gè)處理水平的資料占一列，至少應(yīng)當(dāng)有兩列數(shù)據(jù)。〖答案：t=，P=〗表4－4 兩種病毒在煙葉上產(chǎn)生的枯斑數(shù)株號(hào) 1 2 3 4 5 6 7 8 病毒A 9 17 31 18 7 8 20 10 病毒B 10 11 18 14 6 7 17 5 第五講 Minitab在方差分析中的應(yīng)用　　這一講我們首先學(xué)習(xí)幾個(gè)MINITAB的方差分析命令，然后舉例說明各種試驗(yàn)結(jié)果的模型確定方法、方差分析的MINITAB程序及計(jì)算結(jié)果解釋?！紅=，一尾P=，達(dá)顯著水平〗　　實(shí)習(xí)5 為測(cè)定A、B兩種病毒對(duì)煙草的致病力，取8株煙草，每一株皆半葉接種A病毒，另半葉接種B病毒，以葉面出現(xiàn)枯斑數(shù)的多少作為致病力強(qiáng)弱的指標(biāo)，得結(jié)果于表4?！即鸢福?，s=，Z=，一尾P=〗　　實(shí)習(xí)4 測(cè)定前作噴灑過某種有機(jī)砷殺雄劑的麥田植株樣本4次，(mg)；測(cè)定對(duì)照（前作未用過有機(jī)砷殺雄劑）的植株樣本3次?，F(xiàn)從中選出一個(gè)株系，取10個(gè)纖維樣品。試測(cè)驗(yàn)施氮肥在產(chǎn)量上是否和未施氮肥差異顯著，并以95%的可靠度估計(jì)施氮肥情況下該品種的畝產(chǎn)范圍?！即鸢福?，s=，CV=〗　　實(shí)習(xí)2 某玉米品種在不施氮肥的情況下，一般產(chǎn)量為250kg/畝。五、實(shí)習(xí)　　實(shí)習(xí)1 測(cè)得1960～1972年間越冬代棉紅鈴蟲在江蘇東臺(tái)的羽化高峰期依次為（以6月30日為0）8,6,10,5,6,6,10,1,12,11,9,1,8。試測(cè)驗(yàn)新品種與當(dāng)?shù)仄贩N之間是否有顯著差異。下面是一個(gè)實(shí)例。 ＃ 計(jì)算F值的累積概率　　F 7 6.　　＃ 不使用聯(lián)合方差的t測(cè)驗(yàn)　　TWOS C1 C2　　 end　　2 計(jì)算結(jié)果與分析　　 （，）　　TWOSAMPLE T FOR C1 VS C2　　N MEAN STDEV SE MEAN　　C1 7 　　C2 8 12　　95 PCT CI FOR MU C1 MU C2： (， 65)　　TTEST MU C1 = MU C2 (VS NE)： T= P= DF= 9　　結(jié)果表明：由于兩個(gè)樣本方差不同質(zhì)，所以使用兩個(gè)樣本各自的方差進(jìn)行t測(cè)驗(yàn)，所以兩個(gè)樣本平均數(shù)差異顯著，也就是說密點(diǎn)播產(chǎn)量顯著不同于撒播產(chǎn)量。試測(cè)驗(yàn)兩種播種方式的小麥產(chǎn)量是否有顯著差異。　　POOLED.　　 end　　2 MINITAB計(jì)算結(jié)果　　 （）　　TWOSAMPLE T FOR C1 VS C2　　N MEAN STDEV SE MEAN　　C1 10 　　C2 10 　　95 PCT CI FOR MU C1 MU C2： (， )　　TTEST MU C1 = MU C2 (VS NE)： T= P= DF= 18　　POOLED STDEV = 　　結(jié)果表明：由于兩個(gè)樣本方差同質(zhì)，所以使用兩個(gè)樣本的聯(lián)合方差進(jìn)行t測(cè)驗(yàn)，所以兩個(gè)樣本平均數(shù)具有顯著差異，也就是說兩種飼料對(duì)豬增重有顯著差異?！　?Example 44　　SET C1　　30, 35, 40, 32, 42, 31, 41, 38, 36, 34　　END　　SET C2　　25, 27, 33, 35, 37, 33, 33, 34, 31, 29　　END　?。?方差同質(zhì)性檢驗(yàn)　　LET K1=STDEV(C1)**2/STDEV(C2)**2　　CDF K1。試測(cè)驗(yàn)兩種飼料對(duì)豬增重有無顯著差異。如果置信水平?jīng)]有指明，假定置信水平為95%；如果沒有ALTERNATIVE子命令，做兩尾測(cè)驗(yàn)?！　±樱篢WOSAMPLE 99 C1 C2 （進(jìn)行CC2兩組數(shù)據(jù)的99%的t測(cè)驗(yàn)）　　　　用法：TWOT [置信水平K] C，C　　子命令：ALTERNATIVE，POOLED　　功能：作兩個(gè)平均數(shù)差數(shù)的t測(cè)驗(yàn)和置信區(qū)間的估計(jì)。當(dāng)P，達(dá)顯著水平；當(dāng)P，達(dá)極顯著水平?！　? MINITAB程序：　　 Example 43　　SET C1　　 　　END　　TTEST MU= C1　　 end　　2 MINITAB輸出結(jié)果：　　TEST OF MU = VS MU . 　　N MEAN STDEV SE MEAN T P VALUE　　C1 4 　　結(jié)果表明：當(dāng)t=，故推斷實(shí)測(cè)結(jié)果與理論值相符合。　　說明：Z代表正態(tài)分布值，P VALUE為其概率值，用以說明Z的顯著性。某日抽樣調(diào)查10包，、(=101)公斤?！　≌f明：其中K為置信概率，如果置信概率沒有指明，那么缺省值為95；如果此值小于1，那么將乘以100?！　≌f明：如果沒有指明K，那么使用μ=0；如果沒有使用子命令，執(zhí)行兩尾測(cè)驗(yàn)。如果置信概率沒有指明，那么缺省值為95；如果此值小于1，那么將乘以100?！　±樱篫TEST 6 2 C1 (檢驗(yàn)C1列中的數(shù)據(jù)是否來自N(6,2)的總體)　　　　用法：ZINTERVAL [置信水平=K%] σ=K， C，...，C　　功能：估計(jì)已知總體方差的置信區(qū)間?！　∽用睿篈LTERNATIVE=K　　如果ALTERNATIVE=1，那么被擇假設(shè)為μK；如果ALTERNATIVE=1，那么被擇假設(shè)為μK。σ為總體的標(biāo)準(zhǔn)差。二、單個(gè)平均數(shù)的假設(shè)測(cè)驗(yàn)　?。ㄒ唬┟睢　≡贛initab中，有4個(gè)命令用來完成單個(gè)平均數(shù)的假設(shè)測(cè)驗(yàn)和置信區(qū)間的估計(jì)，簡介如下：　　　　用法：ZTEST [μ=K] σ=K C，...，C　　功能：在σ已知時(shí)，用正態(tài)分布推斷樣本平均數(shù)是否屬于已知的總體。試計(jì)算它的平均數(shù)、標(biāo)準(zhǔn)差、變異系數(shù)等統(tǒng)計(jì)量。　　說明：將一列的數(shù)據(jù)由小到大排列，Q3為3(N+1)/4位置的值，Q1為(N+1)/4位置的值。第四講 Minitab的基本統(tǒng)計(jì)命令與假設(shè)測(cè)驗(yàn)　　本講介紹了MINITAB的基本統(tǒng)計(jì)命令，包括變數(shù)的統(tǒng)計(jì)描述、單個(gè)平均數(shù)的統(tǒng)計(jì)推斷、兩個(gè)平均數(shù)的假設(shè)測(cè)驗(yàn)等，并給出了具體實(shí)例的MINITAB程序和計(jì)算結(jié)果的分析?！　?、C2兩列中的第2～4行，然后用PRINT顯示CC2兩列，觀察結(jié)果，再用ERASE命令將這兩列刪除，最后用COPY命令將CC4的內(nèi)容復(fù)制到CC2中。表34 DATA EDITOR的功能鍵功能鍵 說明 F1 獲取幫助信息 F2 切換編輯模式 F3 切換輸入方向 F4 切還壓縮與非壓縮顯示模式 F6 相當(dāng)于CTRL+ENTER鍵 F7 插入一個(gè)單元 SHIFT+F7 插入一行 F8 刪除一個(gè)單元 SHIFT+F8 刪除一行 F10 進(jìn)入DATA EDITOR菜單 實(shí)習(xí)　　、READ命令輸入下面兩列數(shù)據(jù)：C1 C2 1 32 2 35 3 29 4 30 5 25 　　，然后用READ命令將其讀入到CC4兩列中。DATA EDITOR的菜單命令列于表33。欲調(diào)整列的顯示格式可以利用DATA EDITOR的菜單命令REFORMAT重新格式化這列?！　?duì)字符型列，其寬度是由輸入的最大字符串長度決定的；對(duì)于數(shù)字列，DATA EDITOR根據(jù)您輸入的數(shù)據(jù)動(dòng)態(tài)地調(diào)整數(shù)據(jù)格式，如果輸入的數(shù)據(jù)沒有小數(shù)點(diǎn)，則為整數(shù)型列，如果輸入的數(shù)值有小數(shù)點(diǎn)，則為浮點(diǎn)型列，如果輸入的數(shù)據(jù)大于1e+9，則為指數(shù)型列?！　∫话銇碚f，DATA EDITOR可以自動(dòng)格式化您輸入的數(shù)據(jù)。按SHIFT+F7鍵，在當(dāng)前行處插入一行；按SHIFT+F8鍵刪除當(dāng)前行?！　h除單個(gè)數(shù)據(jù)單元：移動(dòng)光標(biāo)至指定的單元，按F8鍵，將刪除當(dāng)前單元，這列中的其余單元將上移一個(gè)單元。　　插入缺失值的方法：移動(dòng)光標(biāo)至待插入數(shù)據(jù)的單元，按F7鍵，在當(dāng)前位置插入一個(gè)*，代表缺失值，此列中后續(xù)數(shù)據(jù)順序向下移動(dòng)一個(gè)單元；　　插入新值的方法：移動(dòng)光標(biāo)至待插入數(shù)據(jù)的單元，輸入新值，按F7鍵即輸入新值?！　?)刪除字符　　在活動(dòng)單元中刪除單個(gè)字符的步驟：　　①按F2進(jìn)入EDIT模式；　?、谝苿?dòng)光標(biāo)到正確的位置；　?、郯碊EL鍵，當(dāng)前光標(biāo)處的字符被刪除；　?、馨碋NTER鍵，完成編輯?！　? 編輯已存在的數(shù)據(jù)　　1)插入字符　　在活動(dòng)單元中插入字符的步驟為:　?、侔碏2鍵，進(jìn)入EDIT模式；　?、谝苿?dòng)光標(biāo)到正確的位置；　?、鄞蛉胝_的字符，新字符即插