周益林, 聶 曉, 劉 偉, 范潔茹, 梁革梅, 陸宴輝

(中國農(nóng)業(yè)科學院植物保護研究所,植物病蟲害生物學國家重點實驗室, 北京 100193)

數(shù)理統(tǒng)計學是隨著人類社會發(fā)展和管理的需要而從應用統(tǒng)計學發(fā)展起來的一個分支學科,它是一門通過運用多種模型與技術分析、社會調(diào)查與統(tǒng)計分析來收集數(shù)據(jù)與處理數(shù)據(jù)的方法,對科技前沿、國民經(jīng)濟中復雜或者重大的問題以及社會和政府中的大量問題進行定量分析的學科。此學科從卡爾·皮爾遜(Karl Pearson)[1]于20世紀初創(chuàng)立以來,作為一個實用性極強的學科已廣泛應用于社會各個領域,其在植物保護學科領域中的大量應用始于20世紀60年代,特別是70年代計算機及其相關工程技術的快速發(fā)展,大大促進了數(shù)理統(tǒng)計在植保學科特別是昆蟲生態(tài)學、植物病害流行學及病蟲害監(jiān)測預警等領域中的廣泛應用。郭予元院士作為數(shù)理統(tǒng)計學在植保學科的應用研究先驅(qū)者之一,從20世紀50-80年代就開始在這方面做了大量工作,為數(shù)理統(tǒng)計在植物保護學科中應用研究以及教育普及工作作出了重大貢獻。

1 郭予元院士在植保數(shù)理統(tǒng)計理論和方法中的重要貢獻

20世紀60年代郭予元院士在寧夏永寧農(nóng)校工作期間,計算器/計算機還不普及,大量的數(shù)理統(tǒng)計的計算工作還需要依賴手工計算,其中相關分析是在病蟲害調(diào)查數(shù)據(jù)分析中經(jīng)常使用的重要方法之一。簡單的相關系數(shù)比較容易計算,但多重相關分析由于公式復雜,其手工計算過程比較麻煩,為了簡化這一過程,郭予元院士根據(jù)分解線性回歸方程的原理,創(chuàng)造性地提出了一種緊湊計算法。這種方法不論有多少要分析的因素,都無需查找特殊的公式就可簡捷地把各種相關系數(shù)算出來[2]。在此基礎上郭予元院士于1985年又進一步提出了緊湊消去法解多元回歸聯(lián)立方程的計算方法,這兩種計算相關系數(shù)和回歸系數(shù)的緊湊消去法的特別之處在于把幾階矩陣的中間計算過程串聯(lián)在一起,減少了手工抄寫,大大提高了計算速度,節(jié)省了計算時間,而且涉及的因子越多,其優(yōu)點更加突出,同時計算過程中便于檢驗,能避免出錯[3]。郭予元院士(1981年)通過對科技期刊上已發(fā)表的求致死中量的若干實例數(shù)據(jù)研究分析發(fā)現(xiàn),在統(tǒng)計分析殺蟲劑的致死劑量中采用多項式配線精確度一般優(yōu)于常用的機值分析法,特別是在數(shù)據(jù)點偏離劑量-死亡線的情況下,機值分析法誤差大且結(jié)果很不可靠。但是采用多項式配線的最大問題是運算過程太麻煩,遠遠比不上機值分析法簡便快捷,盡管英國著名統(tǒng)計學家羅納德·費希爾(Ronald Fisher)[4]在1938年提供的正交多項式表可加快系數(shù)變換后的配線速度,但是這種方法解決不了系數(shù)復原的簡化問題,進行高次方程的系數(shù)復原非常費事。郭予元院士根據(jù)正交多項式編制了系數(shù)復原表,使整個配線過程全部表格化,成功解決了這一問題[5-6]。同時他還通過進一步數(shù)據(jù)分析提出了以上兩種方法適用的數(shù)據(jù)特征和各自的局限性。

在20世紀80年代初,在比較高端、有統(tǒng)計功能的計算器如Casio fx-180p、Casio fx-3600p等還不普及的情況下,郭予元院士(1983年)就開發(fā)出基于單變量函數(shù)式普通電子計算器直線回歸分析程序,由于這類早期型號的計算器(如Casio fx-120)只帶有一個存儲器,還不能實現(xiàn)真正意義的統(tǒng)計學分析如直線回歸分析,他巧妙地利用這單一存儲器創(chuàng)造性地通過對該類型號計算器編程實現(xiàn)快速連續(xù)運算,約1 min可直接求出直線回歸的截距、斜率和回歸誤差[7]。這一絕妙的思路和方法成為他以后利用Casio fx-180p、Casio fx-3600p等型號計算器進行統(tǒng)計學分析方法編程的開端,由此可以看出郭予元院士深厚的數(shù)理統(tǒng)計功底和功夫。隨著Casio fx-180p系列型號計算器的推出和普及,郭予元院士基于此系列型號計算器開發(fā)出幾乎涵蓋植保學科領域常用的數(shù)理統(tǒng)計學方法或模型的編碼程序,如分組資料的平均數(shù)和方差計算、方差同質(zhì)性分析、自相關和偏相關系數(shù)計算、機值法求LD50及卡方檢驗、Leslie矩陣最大特征值和特征向量、邏輯斯蒂曲線擬合及K值求法、3～5次方程求根、病蟲害各種空間分布型擬合(正態(tài)分布、二項分布、泊松分布、負二項分布、奈曼A分布、泊松-二項分布等)、諧波分析、生態(tài)位和病蟲危害指數(shù)等參數(shù)計算。從郭予元院士這些編碼程序可以看出,他把此型號系列計算器的2程序存儲區(qū)共38步可用的功能發(fā)揮到極致,只可惜這些程序編碼當時都沒有在公開期刊或雜志發(fā)表,只保留在郭先生的工作筆記和課程教案中,或者參加過他的數(shù)理統(tǒng)計培訓班或上過他的昆蟲生態(tài)學課程的學生或?qū)W員的筆記中有部分保存,這也是非常遺憾的事情,本文在此展示兩個例子(正態(tài)分布和諧波分析),以體會下這些編碼程序的精妙思路及它的美學之處。

例子1的正態(tài)分布由法國數(shù)學家亞伯拉罕·棣莫弗(Abraham De Moivre)于1733年首次提出,后由德國數(shù)學家Gauss率先將其應用于天文學研究,故正態(tài)分布又叫高斯分布,該分布是許多統(tǒng)計學方法的基礎,應用范圍廣泛。其公式如下:

上式中π為常數(shù),μ為數(shù)學期望,σ2為方差,隨機變量x服從參數(shù)為μ、σ2的正態(tài)分布或高斯(Gauss)分布。

例子2的諧波分析在18世紀和19世紀已經(jīng)奠定了良好的基礎,傅里葉等提出的諧波分析方法仍被廣泛應用。盡管此分析方法起源于電磁學領域,但也可用于對病蟲害時間序列的分析和模擬預測,滿足一定條件(Dirichlet條件)的、以T為周期的時間周期函數(shù)f(x),在連續(xù)點處,可用下述三角函數(shù)的線性組合(傅里葉級數(shù))來表示:

上式稱為f(x)的傅里葉級數(shù),a0為直流分量,r為諧波次數(shù),ar和br為諧波幅度,L為最大諧波次數(shù),x為等距值。

Casio fx-180p求a0,ar,br: 1 Kin 5 MODE 0 INV PCL P1Kout 6 M+ MR sin*ENT Kin+3 Kin 4=Kin+1 MR cos*Kout 4=Kin+2 INV RTN INV P2Kout 5÷÷Kout 2=INV HLT Kout 1=MODE·INV KAC INV Min 2πm/n Kin 6 INV M-n Kin 5 P1y1RUNy2RUN····AC INV P2→arRUN→brKout 3÷N→a0。

郭予元院士早期在數(shù)理統(tǒng)計學理論和方法方面的貢獻中其“正交多項式配線求殺蟲劑的致死中量”和“多元回歸分析的因子相關選擇法”兩項成果分別獲得了1983年寧夏技術改進獎二等獎和寧夏優(yōu)秀科技成果獎三等獎。

2 郭予元院士在植保數(shù)理統(tǒng)計學應用研究中的重要貢獻

20世紀50年代中期郭予元院士在寧夏回族自治區(qū)試驗農(nóng)場做稻瘟病(病原菌Pyriculariaoryzae)發(fā)生規(guī)律、預測預報和防治研究的時候,就開始在試驗研究中應用相關的數(shù)理統(tǒng)計學方法。他從調(diào)查稻瘟病發(fā)生規(guī)律入手,在田間調(diào)查和病菌孢子捕捉技術收集數(shù)據(jù)基礎上,采用多元回歸分析方法對多年稻瘟病病情和氣象資料進行分析,通過對分析結(jié)果的實地驗證,證實關鍵時期雨量和雨日的乘積與稻瘟病發(fā)生程度極顯著相關,由此制定了預測稻瘟病病情和確定最佳防治時間的中長期預測模型,此模型可準確測報稻瘟病的發(fā)生,有效地指導對該病的防治。當時,采用數(shù)理統(tǒng)計方法定量化研究病蟲害規(guī)律并建立病蟲害預測模型的研究在國內(nèi)是非常少見的,只可惜限于當時客觀條件,這些有重要意義的結(jié)果大多都沒能在相關學術期刊上正式發(fā)表[8]。

1982年底郭予元院士調(diào)入中國農(nóng)業(yè)科學院植物保護研究所,此后在植物保護研究所工作的35年間,在他的引領下他和團隊的老師及學生使用數(shù)理統(tǒng)計的方法解決了不少植保學科數(shù)據(jù)分析中的難題、障礙和科學問題,取得許多豐碩的研究成果,為數(shù)理統(tǒng)計在植保學科中的應用作出了重大貢獻。此文鑒于篇幅所限,在此只介紹郭予元先生在害蟲生命表和種群動態(tài)及多病蟲危害產(chǎn)量損失和防治指標兩個方面的研究工作。

2.1 害蟲生命表和種群動態(tài)的研究

20世紀80年代開始郭予元院士和團隊的研究人員就開展了麥田黏蟲Mythimnaseparata生命表的研究,在對黏蟲發(fā)育階段致死因子K值和種群趨勢指數(shù)分析的基礎上,首次采用了歐式距離相似程度的方法分析其關鍵因子,篩選出了影響?zhàn)はx種群的關鍵因子為蛹的被寄生和1～3齡幼蟲的被捕食,這2個因子是造成一代黏蟲年度間變動的主要原因。此研究采用的歐式距離分析法比前人常用的圖解法、b值分析法、r2值分析方法等所得結(jié)果更精確,還可避免采用逐步回歸分析存在的缺陷,即不能排除選中高負值r的因子可能性(這種情況是不符合生物學邏輯,即由此因子造成的死亡率越高,最終殘蟲率也越高,是不合理的)[11]。1989年郭予元院士和團隊的研究人員通過田間試驗首次建立了麥長管蚜Sitobionavenae(目前已更名為荻草谷網(wǎng)蚜Sitobionmiscanthi)田間自然種群的生命表,此研究在對麥長管蚜各代種群趨勢指數(shù)分析的基礎上,采用聚類分析中的距離系數(shù)方法對影響麥長管蚜種群的影響因子K值進行了分析,明確了影響種群波動的主要因素是天氣條件和天敵,低齡若蚜主要受風雨影響,高齡若蚜和成蚜則受天敵和風雨共同作用[12]。隨后還建立了棉鈴蟲Helicoverpaarmigera自然種群生命表,通過關鍵因子分析明確了影響全年棉鈴蟲種群數(shù)量變動的關鍵因子是2齡期被寄生,年度間關鍵因子是5齡期被捕食,明確了化學防治對種群I值的影響,給出了相應的化學防治策略[13]。1994年通過對田間麥無網(wǎng)長管蚜Metopolophiumdirhodum系統(tǒng)調(diào)查,分析了該蚜蟲的時空分布格局和動態(tài)變化,通過對生態(tài)位的定量分析,明確了麥無網(wǎng)長管蚜、禾谷縊管蚜Rhopalosiphumpadi和麥長管蚜的種間競爭關系,在用相關分析和通徑分析明確麥無網(wǎng)長管蚜種群消長關鍵因子的基礎上,通過多元回歸分析法建立了基于關鍵因子的預測模型[14]。

2.2 多病蟲害危害產(chǎn)量損失及復合防治指標

20世紀80年代開始郭予元院士以及團隊的科研人員還開展了棉鈴蟲對棉花的危害和防治指標的研究,組建了不同年份、不同肥料水平及不同代棉鈴蟲為害棉花的產(chǎn)量損失率模型,在這項研究過程中采用協(xié)方差分析成功地解決了模型之間的斜率和截距差異性比較,明確了不同模型的適用性范圍和穩(wěn)定性差異,最終分別建立了適用于三種情況,即高產(chǎn)二代或中產(chǎn)二代或高產(chǎn)、中產(chǎn)、低產(chǎn)三代及低產(chǎn)二代的棉鈴蟲損失率估計模型[15]。1988年通過田間自然蟲源輔助人工控制培養(yǎng)研究了4種麥蚜混合種群對小麥產(chǎn)量的影響,此研究采用偏相關分析方法,成功地解決了穗數(shù)對小麥產(chǎn)量的影響及不同蚜種對小麥產(chǎn)量的影響問題,從而明確了4種蚜蟲混合種群中麥長管蚜和禾谷縊管蚜對產(chǎn)量有明顯影響,并通過多元回歸分析方法獲得了復合蚜量的動態(tài)防治指標[16]。郭予元院士等于1994年通過田間開放式小區(qū)試驗,巧妙地采用了回歸最優(yōu)設計的“206”田間試驗設計,對河南省二代棉鈴蟲和葉螨(朱砂葉螨Tetranychuscinnabarinus和截型葉螨T.truncatus的復合種群)為害棉花的復合防治指標進行了研究,在數(shù)據(jù)分析中用偏相關分析發(fā)現(xiàn)黃萎病對葉螨有一定干擾作用,因此通過協(xié)方差方法排除了黃萎病的干擾,并用主成分分析明確了棉鈴蟲和葉螨的一次和二次型數(shù)值,尤其是二者的交互作用,建立了棉鈴蟲和葉螨復合為害的二次曲線方程,并組建了一系列復合防治指標[17]。另外在棉鈴蟲預測預報研究方面,在華北棉區(qū)用掃網(wǎng)法和目測查蟲法調(diào)查麥田一代棉鈴蟲幼蟲發(fā)生,用比較法和生命表分析法預測二代棉鈴蟲在棉田的發(fā)生程度和蟲口密度,經(jīng)1991年-1994年應用檢驗,其預測結(jié)果均與實際發(fā)生情況吻合[18]。

3 郭予元院士在植保數(shù)理統(tǒng)計學教學和普及中的重要貢獻

郭予元院士值得一書的還有他在植保數(shù)理統(tǒng)計學的教學和普及方面的重要貢獻。多年來他通過在全國各地舉辦植保數(shù)理統(tǒng)計培訓班和研究生授課,大大提升了我國從事植保學科的一代又一代科技人員在試驗設計、科學地統(tǒng)計和分析試驗數(shù)據(jù)的水平,解決了一些科研和生產(chǎn)及相關的數(shù)理統(tǒng)計中的重大問題和復雜問題,為科學地預測和防治病蟲害以及提高我國植保的科學水平都起到重要的作用。

3.1 植保數(shù)理統(tǒng)計學培訓班

20世紀80年代初郭先生在實際研究和工作中發(fā)現(xiàn),由于我國教育體系受文化大革命的影響,植保領域的研究人員特別是從事病蟲害預測預報等方面的年輕工作者普遍存在數(shù)理統(tǒng)計知識薄弱的問題,嚴重制約植保學科領域科學研究的發(fā)展和研究水平的提升。他把自己多年來堅持自學的數(shù)理統(tǒng)計知識和經(jīng)驗總結(jié)整理出適用于植?？茖W試驗研究的一套統(tǒng)計分析方法—《植保數(shù)理統(tǒng)計學方法》,并自己刻板油印,先后在我國7個省20多個農(nóng)業(yè)科研和推廣單位舉辦多次數(shù)理統(tǒng)計講習班,不但使當時許多植?？蒲泻屯茝V人員的研究論文數(shù)據(jù)分析的科學性得到明顯提高,更重要的是培養(yǎng)了一大批懂數(shù)理統(tǒng)計的植保領域的專家和學者,這些人員至今還活躍在各自的研究領域并發(fā)揮著重要的作用。

郭予元院士油印本的《植保數(shù)理統(tǒng)計學方法》內(nèi)容廣泛,不僅涵蓋了當時植保學科領域所有能用到的基本統(tǒng)計學方法,如回歸分析、相關分析和通徑分析、曲線方程擬合、病情發(fā)生趨勢預測方法、列聯(lián)表分析法、聚類分析和判別分析、t檢驗和方差分析、正交試驗、卡方分析、協(xié)方差分析、空間分布型和抽樣技術、病蟲危害特性和防治指標等,還包括了當時一些剛剛出現(xiàn)的數(shù)學或統(tǒng)計新興領域的方法如模糊數(shù)學、灰色理論等,還有一些不太常見的統(tǒng)計學方法或當時比較新的田間試驗設計方法如題總相關分析、嶺回歸方法、灰色突變長期預測模型、回歸正交旋轉(zhuǎn)設計等也在培訓班中多有涉及。一些學員或?qū)W生在學習后就很快在自己研究中應用,如有學員1989年在小麥白粉病(病原菌:Blumeriagraminisf.sp.tritici)田間藥劑防治試驗中就采用了回歸正交旋轉(zhuǎn)設計,這種試驗設計的優(yōu)點是試驗處理少,獲得的信息多,精確度高,而且計算簡便,從藥劑試驗結(jié)果看,它不但能獲得每種藥劑的防效大小和藥劑之間交互作用的有無,更重要的是可以獲得每種藥劑用量與病情之間的關系模型,從而在防治時可根據(jù)實際情況,計算出最適用藥量[19]。

郭予元院士授課通俗實用,所用的實例都是植保領域田間試驗或病蟲害發(fā)生流行的相關數(shù)據(jù)實例,每次上課都根據(jù)學員的實際情況舉例,而且在課堂上實時采用Casio fx-3600p、Casio fx-180p計算器編程,進行實際操作計算,使學員能學到數(shù)理統(tǒng)計的真本領,尤其重要的是他在授課時不但教各種統(tǒng)計學方法的使用,還特別強調(diào)使用時需要注意的問題,為此他還專門整理油印了相關的材料如《病蟲測報中的幾個問題》《棉蟲試驗的幾種統(tǒng)計分析方法》等供學員學習。比如在《病蟲測報中的幾個問題》中就強調(diào)了預報因子在選擇時不能只憑生物學邏輯或簡單相關系數(shù)挑選因子,一定還要注意因子之間的自相關性和主次關系;在對有重復使用的數(shù)據(jù)進行回歸分析時,不但要檢驗回歸的顯著性,還要檢驗重復之間失擬的顯著性;在《棉蟲試驗的幾種統(tǒng)計分析方法》中強調(diào)了數(shù)據(jù)轉(zhuǎn)換的問題以及數(shù)據(jù)轉(zhuǎn)換與數(shù)據(jù)分布型的關系,例如在對棉蚜Aphisgossypii田間種群量數(shù)據(jù)進行方差分析時,應先把數(shù)據(jù)轉(zhuǎn)換成對數(shù)值,這是因為棉蚜在田間是聚集分布的。統(tǒng)計分析數(shù)據(jù)轉(zhuǎn)換和數(shù)據(jù)分布型在當前依然是植保數(shù)理統(tǒng)計中的薄弱環(huán)節(jié)或被忽視的問題。

郭予元院士不但通過舉辦培訓班傳授數(shù)理統(tǒng)計知識,國內(nèi)植保同行或?qū)W生也經(jīng)常通過信件向他請教相關問題,郭先生都是有問必答,甚至親自幫助不少老師處理和分析他們的試驗數(shù)據(jù)。如1986年前后新疆植保站的榮麗君老師就數(shù)次給郭先生寫信請教多年度害蟲為害產(chǎn)量損失數(shù)據(jù)的分析和模型的組建問題,郭先生都一一回信做了解答,而且還親自對他的試驗數(shù)據(jù)進行了統(tǒng)計分析,通過采用協(xié)方差的分析方法,解決了黃地老虎Agrotissegetum、冬麥地老虎Caradrinaauguroides、麥長管蚜、棉鈴蟲和牧草盲蝽Lyguspratensis等5種害蟲不同年份模型差異性的檢驗問題,并建立了適用不同情況或條件下相關蟲害為害產(chǎn)量損失的損失估計模型,相關的文章也在學術期刊上及時發(fā)表[20-22]。

3.2 昆蟲生態(tài)學研究生課程

郭予元院士10多年來一直擔任中國農(nóng)業(yè)科學院研究生院《昆蟲生態(tài)學》課程的主講老師,他的課程授課內(nèi)容豐富、涵蓋范圍廣,包括昆蟲種群的時間動態(tài)、空間模式、種間關系以及群落生態(tài)系統(tǒng)等諸多方面,其課程內(nèi)容的最大特點是每個講授內(nèi)容或章節(jié)都有實際的例子,而且如果涉及數(shù)理統(tǒng)計內(nèi)容,都有課堂實時分析計算,使同學們不但學到了昆蟲生態(tài)學的理論知識,而且還學到數(shù)理統(tǒng)計的實戰(zhàn)方法,并在自己的科研中很快就能得到應用。特別是在講授種群的空間圖式(正態(tài)分布、泊松分布、均勻分布、二項分布和負二項分布、奈曼A型分布、泊松-二項分布等P-E分布型)、種群增長模型(邏輯斯蒂曲線、指數(shù)增長曲線、Gompertz曲線、冪函數(shù)曲線、差數(shù)指數(shù)曲線、差數(shù)冪函數(shù)曲線、雙曲線、麥肯齊模型等)、種間關系(競爭模型、捕食模型等)以及這些模型的擬合和檢驗時,郭先生都有實際的植保例子供學習和分析,最有特色的是在課堂上用他開發(fā)的基于計算器編程的程序進行實時計算,這種教學方法和模式,大大提升了同學們昆蟲生態(tài)學知識和相應的數(shù)理統(tǒng)計學知識水平。

4 郭予元院士團隊對數(shù)理統(tǒng)計在植保領域應用后續(xù)工作的影響

郭予元院士在植物保護研究領域特別是發(fā)展棉鈴蟲綜合防治策略和技術體系、有效解決生產(chǎn)抗性棉鈴蟲的種群治理問題和促進農(nóng)作物病蟲害綜合防治理論的創(chuàng)新等方面做出重要貢獻,在他的帶動和引領下,年輕的一代植?？蒲泄ぷ髡?尤其是研究團隊繼承和發(fā)揚他們這一輩老科學家思想和精神,青出于藍而勝于藍,在棉花害蟲研究領域做出了世界水平研究成果,在這些科學研究中也大量地使用了數(shù)理數(shù)據(jù)統(tǒng)計的方法,其研究結(jié)果先后在《科學》《自然》等國際著名期刊發(fā)表。例如棉鈴蟲研究團隊通過對我國北方6省過去10年100多個樣點棉鈴蟲在多種作物上卵(多代)和幼蟲發(fā)生數(shù)據(jù)的回歸分析發(fā)現(xiàn),我國華北地區(qū)大規(guī)模種植Bt棉,不僅降低了棉花上棉鈴蟲的數(shù)量,而且減輕了周邊其他非轉(zhuǎn)基因農(nóng)作物(玉米、大豆和花生)上的棉鈴蟲為害程度。這項研究結(jié)果是世界首個針對轉(zhuǎn)基因作物抗蟲性進行的大規(guī)模長期跟蹤研究,該論文曾作為2008年《科學》雜志的封面文章[23];進一步的研究工作又發(fā)現(xiàn)雖然大范圍種植Bt棉大大降低了棉鈴蟲發(fā)生和危害,但另一類害蟲盲蝽種群在棉花中快速上升,通過對長期多點收集的數(shù)據(jù)進行模型擬合,結(jié)果顯示,Bt棉種植比例與該害蟲的種群密度呈線性正相關,與盲蝽殺蟲劑的使用次數(shù)呈非線性指數(shù)增長關系,通過對種群密度與多因子包括殺蟲劑的使用次數(shù)、主要氣象因子溫度、降雨等的回歸和相關分析表明,該害蟲的種群密度增長主要是種植Bt棉帶來的棉鈴蟲殺蟲劑打藥次數(shù)減少造成的。另外研究團隊對我國華北大范圍種植Bt棉后該害蟲在棉花、蘋果、葡萄、桃、梨、棗等作物上的發(fā)生危害嚴重度多年多點發(fā)生數(shù)據(jù)的模型模擬發(fā)現(xiàn),其危害嚴重度與Bt棉的種植比例也呈指數(shù)增長關系,該研究論文2010年再次登上《科學》雜志[24]。2012年研究團隊又闡明了Bt棉種植后由于打藥次數(shù)減少導致的天敵昆蟲種群和棉蚜種群等的變化趨勢,論文發(fā)表于《自然》雜志[25],同年該研究入選“中國科學十大進展”。這篇論文在數(shù)據(jù)分析中也采用了數(shù)理統(tǒng)計的多種方法如多因子方差分析、逐步回歸、相關分析、直線或曲線模型擬合等方法。

郭予元院士不僅在植?？蒲猩献龀隽司薮蟮呢暙I,而且在數(shù)理統(tǒng)計學的理論和方法研究、數(shù)理統(tǒng)計在植保研究領域的應用以及數(shù)理統(tǒng)計學知識的教學普及和推廣方面做出了重要的貢獻。近年來隨著數(shù)理統(tǒng)計學本身的發(fā)展,如廣義線性模型、加性模型、混合效應模型、隨機系數(shù)模型等發(fā)展,同時一些重要數(shù)理支持軟件的發(fā)展和改進如SAS、SPSS等,特別是開源軟件R軟件的快速發(fā)展,隨著大數(shù)據(jù)、機器學習、AI人工智能技術時代的來臨,相信郭予元院士奠基和引領的植保數(shù)理統(tǒng)計學必將在植保領域發(fā)揮越來越重要的作用,并產(chǎn)出更多更大的研究成果。