用Excel中的VBA編寫卡方測算相關(guān)程序

楊振宇 ，楊海智 ，楊信東

（1.吉林農(nóng)業(yè)大學(xué) 農(nóng)學(xué)院，吉林 長春 130118；2.吉林省農(nóng)業(yè)科學(xué)院大豆所，吉林 公主嶺 136100；3.健康元藥業(yè)集團股份有限公司，廣東 深圳 518057；4.吉林農(nóng)業(yè)大學(xué) 發(fā)展學(xué)院，吉林 長春 130118）

在生命科學(xué)及社會科學(xué)研究中，卡方（χ2）檢驗是最常用的統(tǒng)計分析方法之一[1-8]。 卡方（χ2）檢驗是次數(shù)資料的顯著性檢驗方法，包括適合性檢驗和獨立性檢驗兩類[9]。適合性檢驗用于檢驗?zāi)承誀钣^察次數(shù)與該性狀的理論比率是否符合，如在遺傳分析中研究兩對性狀雜種后代的分離現(xiàn)象是否符合某一特定比率；而獨立性檢驗是用于判斷兩類因子是彼此相關(guān)還是相互獨立的，如：采用某種措施與預(yù)防某種災(zāi)害的關(guān)系。

卡平方的計算很復(fù)雜，但隨著計算機與統(tǒng)計軟件的普及，應(yīng)用計算機計算卡平方成為最準確可靠的方法，各種大型統(tǒng)計軟件如 SAS、SPSS、DPS等均應(yīng)用于卡平方的測算[10-12]，一些研究者也研發(fā)了在Excel上編制運算表來實現(xiàn)卡平方測算的方法[13]。

使用如SAS、SPSS、DPS等大型統(tǒng)計軟件需要有較高花費；采用R統(tǒng)計軟件來編程進行卡平方測算要求研究者有較強學(xué)術(shù)底蘊和計算機編程能力。這在客觀上限制了 SAS、SPSS、DPS、R等統(tǒng)計軟件的普及應(yīng)用。Excel作為 Microsoft Office家族成員，由于多數(shù)學(xué)生和科研人員都具有Excel基礎(chǔ)知識，采用Excel軟件計算卡平方無需考慮運行環(huán)境而受到歡迎。然而，采用編制運算表的方法還是讓人感到有些繁鎖，不易學(xué)習使用。鑒于此，本文采用一種基于 Excel的 VBA編程方法[14-15]，編寫了“卡平方測算”相關(guān)程序，可在所有配備Office 2000以上版本的計算機上使用。“卡平方測算”在VBA程序運算時只需輸入最原始的數(shù)據(jù)，應(yīng)用步驟實行最直觀的人機對話，任何初學(xué)者都可以即學(xué)即用輕松地掌握程序的使用方法；編寫完的程序成為工作模板后，可以任意復(fù)制或通過e-mail郵寄等方式進行傳播、拷貝，因此這種方法受到使用者的一致好評。

本文介紹了該程序的源代碼及使用方法，讓所有不方便使用大型統(tǒng)計軟件的同行都可分享這些程序帶來的所有便利。

1 VBA程序源代碼

1.1 用于適合性檢驗的卡平方計算程序

1.2 用于獨立性檢驗的卡平方計算程序

1.2.1 2×2表的獨立性測驗

1.2.2 2×c表的獨立性測驗

1.2.3 r×c表的獨立性測驗

2 “卡平方測算”VBA程序的應(yīng)用步驟：

（1）運行環(huán)境：Win2003、WinXP、Win2000 等。

（2）VBA 程序的應(yīng)用步驟：打開 Excel“EC50、EC90測算”程序工作簿（在此過程中，若計算機屏幕顯示對話框，則應(yīng)選擇并點擊“啟用宏”），單擊（程序運行命令）按鈕（本程序中該按鈕上面寫著“計算”），則程序開始運行計算機，屏幕依次出現(xiàn)對話框，按對話框所提問題逐個輸入數(shù)據(jù)，輸完后瞬間即在程序工作簿頁面上顯示輸出結(jié)果（包括卡平方值）。

（3）“卡平方測算”VBA程序的下載及拷貝方法：①下載源程序的方法參見文獻[6-7]；②“卡平方測算”文件的再拷貝：上述工作完成后，此Excel文件即成為一個工作模板，可以任意復(fù)制、粘貼或通過e-mail郵寄等方式進行傳播和拷貝。

本Excel文件所占內(nèi)存約為60 KB，與其他有類似功能的程序相比要小得多。上述4個程序也可以組合成一個較大的程序，但從應(yīng)用角度考慮，這樣做會浪費計算機的內(nèi)部資源（因為事先聲明了較多在當次運算中并不使用的變量），故本文未這樣處理。

[1]武曉玲，周斌，孫石，等.大豆對大豆疫霉菌株 Pm14抗性的遺傳分析及基因定位 [J].中國農(nóng)業(yè)科學(xué)，2011，44（3）：456-460.

[2]王保通，李強，胡茂林，等.小麥品種Libellula和N.strampelli抗條銹病主效、微效基因遺傳分析[J].植物病理學(xué)報，2010，40（3）：300-306.

[3]張宏，任志龍，胡銀崗，等.陜麥 139抗條銹病基因遺傳分析[J].作物學(xué)報，2010，36（1）：109-114.

[4]何麗華，牛寶龍，齊曉朋，等.棉鈴蟲成蟲體色突變體的發(fā)現(xiàn)及其遺傳分析[J].核農(nóng)學(xué)報，2007（4）：397-400.

[5]楊振宇，王曉麗，張曉波，等.部分抗 SMV大豆品種成株抗性基因?qū)?shù)分析 [J].吉林農(nóng)業(yè)大學(xué)學(xué)報，2011，33（6）：591-594.

[6]龔瑞，楊炬，黎唏，等.2007-2010年度寧夏流感監(jiān)測結(jié)果分析[J].寧夏醫(yī)學(xué)雜志，2011，33（3）：222-224.

[7]錢峰.基于卡方檢驗的國內(nèi)外知識管理研究熱點比較[J].情報雜志，2008（9）：56-58.

[8]徐向陽.卡方檢驗在學(xué)生成績差異性分析中的應(yīng)用[J].常州技術(shù)師范學(xué)院學(xué)報，2001，7（4）：13-16.

[9]蓋鈞鎰.試驗統(tǒng)計方法[M].北京，中國農(nóng)業(yè)出版社，2000.

[10]詹秋文.Excel和SAS在生物統(tǒng)計學(xué)的應(yīng)用比較 [J].生物學(xué)雜志，2009，26（1）：74-75，83.

[11]向窮，施樹良，李鈺.常用統(tǒng)計軟件在生物統(tǒng)計中的應(yīng)用比較[J].現(xiàn)代生物醫(yī)學(xué)進展，2009，9（9）：1775-1777，1789.

[12]唐啟義，馮明光.實用統(tǒng)計分析及其DPS數(shù)據(jù)處理系統(tǒng)[M].北京：科學(xué)出版社，2002：188-95.

[13]譚永強，余華強，陳橋生，等.利用 Excel軟件建立卡方檢驗分析模板在農(nóng)業(yè)統(tǒng)計中的應(yīng)用 [J].湖北農(nóng)業(yè)科學(xué)，2010，49（12）：3192-3195.

[14]龔沛曾，陸慰民.Visual Basic程序設(shè)計教程 （6.0版）[M].北京：高等教育出版社，2001.

[15]李曉玫，楊小平.Excel中的 VBA程序設(shè)計[J].四川師范大學(xué)學(xué)報（自然科學(xué)版），2004（4）：423-426.

[16]馬海霞，劉 影，王艷紅，等.用 EXCEL中的 VBA編寫“多項式的三角函數(shù)擬合單峰曲線”程序 [J].菌物研究，2009，7（3-4）：195-200.

[17]段顯德，王艷紅，楊信東.用 EXCEL中的 VBA編寫“試卷分析”程序[J].通化師范學(xué)院學(xué)報，2010，31（8）：52-53.