余 純，許 冬，黃 維，雷明洪，萬(wàn) 鵬，*

(1.江西財(cái)經(jīng)大學(xué) 統(tǒng)計(jì)學(xué)院，江西 南昌 330013； 2.農(nóng)業(yè)農(nóng)村部華中作物有害生物綜合治理重點(diǎn)實(shí)驗(yàn)室/農(nóng)作物病蟲(chóng)草害防控湖北省重點(diǎn)實(shí)驗(yàn)室/湖北省農(nóng)業(yè)科學(xué)院植保土肥研究所，湖北 武漢 430064； 3.湖南省安鄉(xiāng)縣植保植檢站，湖南 安鄉(xiāng) 415600)

昆蟲(chóng)的發(fā)生與生理行為受環(huán)境影響，不同種類(lèi)的昆蟲(chóng)，對(duì)外界環(huán)境條件中的溫度、濕度、光照甚至月相等因子的適應(yīng)性各不相同。如棗鐮翅小卷蛾(Ancylissativa)和二化螟(Chilosuppressalis)在高溫時(shí)求偶持續(xù)時(shí)間縮短[1]，異色瓢蟲(chóng)在低溫下的飛行能力較強(qiáng)但搜索能力較弱[2]，棉鈴蟲(chóng)(Helicoverpazea)雌成蟲(chóng)在高溫或低溫條件下性信息素釋放量均會(huì)減少[3]。甜菜夜蛾成蟲(chóng)壽命隨濕度的提高而增加[4]，濕地松粉蚧若蟲(chóng)的存活率隨降雨強(qiáng)度的增加顯著下降[5]。長(zhǎng)光照能促進(jìn)異色瓢蟲(chóng)的求偶交配欲[6]，短日照能減少稻縱卷葉螟的南遷[7]，三化螟在滿月期的撲燈數(shù)顯著少于月虧期，而黑尾葉蟬的誘集量在滿月期顯著高于其他月相期[8]。因此，研究環(huán)境因子對(duì)各種昆蟲(chóng)生理行為的影響，有助于深入認(rèn)識(shí)該昆蟲(chóng)與環(huán)境之間的互作規(guī)律，為其在生產(chǎn)上的治理提供理論依據(jù)。

紅鈴蟲(chóng)為長(zhǎng)江流域棉花上的重要害蟲(chóng)，能?chē)?yán)重危害棉花的產(chǎn)量與品質(zhì)。有研究表明，月相和溫度均能顯著影響紅鈴蟲(chóng)的發(fā)生[9]。為此，本研究采用性信息素誘集的方式，監(jiān)測(cè)了湖南安鄉(xiāng)棉區(qū)的紅鈴蟲(chóng)種群動(dòng)態(tài)，并運(yùn)用非凸懲罰的穩(wěn)健線性回歸模型，分析了月相與溫度對(duì)紅鈴蟲(chóng)發(fā)生量的影響，以期能加深對(duì)紅鈴蟲(chóng)發(fā)生規(guī)律的認(rèn)識(shí)，為該地區(qū)紅鈴蟲(chóng)的綜合治理提供理論依據(jù)。

1 材料與方法

1.1 基于性誘劑的紅鈴蟲(chóng)種群動(dòng)態(tài)監(jiān)測(cè)

2012—2015年，在紅鈴蟲(chóng)的盛發(fā)期，采用紅鈴蟲(chóng)信息素誘集器(北京中捷四方生物科技股份有限公司)在湖南省安鄉(xiāng)縣的棉田內(nèi)進(jìn)行紅鈴蟲(chóng)成蟲(chóng)誘集。按照說(shuō)明書(shū)組裝好誘集器各組件，將其放置于棉田內(nèi)，使誘集器高于棉花植株0.5 m，每天收集所誘成蟲(chóng)并記數(shù)。為保證誘集效果，每月更換1次誘芯。誘集時(shí)間從每年的5月下旬至9月底。

1.2 紅鈴蟲(chóng)監(jiān)測(cè)點(diǎn)氣象數(shù)據(jù)的獲取

紅鈴蟲(chóng)監(jiān)測(cè)點(diǎn)的氣象數(shù)據(jù)從安鄉(xiāng)縣氣象局獲取。

1.3 數(shù)據(jù)處理與分析

本研究采用穩(wěn)健的回歸分析——基于非凸懲罰的線性回歸模型，研究紅鈴蟲(chóng)種群與月相、溫度之間的關(guān)系，以避免觀察數(shù)據(jù)中極端值或離群值對(duì)模型參數(shù)估計(jì)的影響。在回歸分析前，從3個(gè)方面進(jìn)行考慮：一是完全用陰歷的日期來(lái)代替月相，分析不同月相時(shí)紅鈴蟲(chóng)的發(fā)生量；二是根據(jù)月相，將陰歷的日期拆分成新月、弦月、滿月等時(shí)間段，分別賦值，分析其與紅鈴蟲(chóng)發(fā)生量之間的關(guān)系；第三是除月相外，還需考慮光線的影響，即夜間光線的強(qiáng)弱并結(jié)合月相來(lái)對(duì)陰歷的每天進(jìn)行賦值，再由此分析其與紅鈴蟲(chóng)發(fā)生量之間的關(guān)系。具體賦值情況如下：雨天，0；陰天，1；多云，2；晴天則依月相的不同而賦值不同，其中新月、殘?jiān)?農(nóng)歷初1、初2、29、30)，3；蛾眉月(農(nóng)歷初3、初4、27、28)，4；蛾眉月(農(nóng)歷初5、初6、25、26)，5；弦月(農(nóng)歷初7、初8、23、24)，6；弦月(農(nóng)歷初9、初10、21、22)，7；凸月(農(nóng)歷11、12、19、20)，8；凸月(農(nóng)歷13、14、17、18)，9；滿月(農(nóng)歷15、16)，10。

1.3.1 線性回歸模型

對(duì)于給定的相互獨(dú)立的、具同樣分布特征(independent and identically distributed， iid)的觀測(cè)值(xi，yi)，i=1，2，3，…，n。為了研究因變量(y)和自變量(x)之間的關(guān)系，通常用以下的線性回歸模型：

y=Xβ+，～N(0，σ2I)。

(1)

模型中y=(y1，…，yn)T，是n個(gè)觀測(cè)值yi組成的向量，如果有p(p≥ 1)個(gè)自變量，X就是維度為n×(p+1)的矩陣。β=(β0，β1，…，βp)T是維度為(p+1)×1的未知系數(shù)向量。=(1，…，n)T是隨機(jī)的誤差向量，假定其服從多元正態(tài)分布時(shí)，對(duì)于第i個(gè)觀測(cè)值(xi，yi)，yi和xi滿足如下關(guān)系：i。最常用的估計(jì)系數(shù)向量β的方法是最小二乘法，其工作原理就是使得殘差平方和最小。

Yu等[10]指出，當(dāng)數(shù)據(jù)中出現(xiàn)極端值或離群值時(shí)，最小二乘法不能準(zhǔn)確估計(jì)系數(shù)向量β。在此情況下，需要運(yùn)用穩(wěn)健的線性回歸模型來(lái)處理數(shù)據(jù)并檢測(cè)離群值。

1.3.2 基于非凸懲罰的穩(wěn)健線性回歸模型

She等[11]提出了基于非凸懲罰的穩(wěn)健線性回歸模型：

y=Xβ+γ+，～N(0，σ2I)。

(2)

這個(gè)穩(wěn)健的線性回歸模型(2)是在傳統(tǒng)的線性回歸模型(1)中引入一個(gè)均值漂移參數(shù)γ=(γ1，…，γn)T，n是樣本容量。當(dāng)?shù)趇個(gè)觀測(cè)值為離群值時(shí)，γi值為非零；否則，γi為0。模型(2)通過(guò)對(duì)均值漂移參數(shù)γ的準(zhǔn)確估計(jì)，可達(dá)到同時(shí)估計(jì)參數(shù)β和檢測(cè)離群值的目的?？赏ㄟ^(guò)使如下的目標(biāo)函數(shù)達(dá)到最小來(lái)估計(jì)參數(shù)β和γ。

(3)

式(3)中，pλ(|γi|)是關(guān)于γi的懲罰函數(shù)，λ是調(diào)諧參數(shù)。

1.3.3 懲罰函數(shù)與閾值法則

上文所述目標(biāo)函數(shù)(3)中的懲罰函數(shù)可以有多種選擇，比較普遍的有l(wèi)0范數(shù)和l1范數(shù)，以及SCAD等懲罰函數(shù)[12-15]。每個(gè)懲罰函數(shù)都有其對(duì)應(yīng)的閾值法則(thresholding rule)，比如l0范數(shù)的懲罰函數(shù)對(duì)應(yīng)有hard閾值法則，l1范數(shù)的懲罰函數(shù)對(duì)應(yīng)于soft閾值法則，SCAD的懲罰函數(shù)對(duì)應(yīng)于SCAD閾值法則。本研究中我們采用l0范數(shù)的懲罰函數(shù)：

式中，I(·)為指示函數(shù)(indicator function)。l0范數(shù)的懲罰函數(shù)對(duì)應(yīng)于以下的hard閾值法則：

(4)

式(4)中，Θ為閾值法則函數(shù)，ξ為變量，λ為調(diào)諧參數(shù)。

1.3.4 參數(shù)估計(jì)的具體算法

1.3.5 調(diào)諧參數(shù)λ最優(yōu)取值的選擇

因?yàn)閼土P函數(shù)和閾值法則中包含調(diào)諧參數(shù)λ，很顯然λ控制γi的估計(jì)值，所以調(diào)諧參數(shù)λ的取值對(duì)準(zhǔn)確估計(jì)γi起到很重要的作用。經(jīng)驗(yàn)的做法是在(λmin，λmax)范圍內(nèi)取100個(gè)λ值，其中，λmin是λ的最小取值，使得向量γ=(γ1，…，γn)T中50%的γi值非零，λmax對(duì)應(yīng)的是λ的最大取值，使得向量γ中所有的γi為0。顯然100個(gè)不同的λ值對(duì)應(yīng)100組不同的向量γ=(γ1，…，γn)T的估計(jì)值，我們用BIC(bayesian information criterion)的標(biāo)準(zhǔn)來(lái)選擇最優(yōu)的調(diào)諧參數(shù)λ取值。選定最優(yōu)的調(diào)諧參數(shù)λ取值后，其對(duì)應(yīng)的一組γ估計(jì)值即為最終的參數(shù)γ估計(jì)值。

通過(guò)對(duì)參數(shù)γ和β的準(zhǔn)確估計(jì)，穩(wěn)健線性回歸模型(2)可以實(shí)現(xiàn)同時(shí)進(jìn)行離群值檢測(cè)和模型參數(shù)估計(jì)。

2 結(jié)果與分析

2.1 月相對(duì)誘蛾量的影響

以月相作為自變量(X)，誘蛾量作為因變量(Y)，分別對(duì)2012年、2013年、2014年和2015年的數(shù)據(jù)做誘蛾量與月相之間的線性回歸分析。同時(shí)，為了增加數(shù)據(jù)的代表性，增大樣本容量，2012—2015年的數(shù)據(jù)也被合并進(jìn)行了線性回歸分析，以期擬合出更準(zhǔn)確的回歸模型。考慮到誘蛾量的觀測(cè)值在(0，1 100)之間，變量取值的范圍和方差太大會(huì)給線性回歸的參數(shù)估計(jì)帶來(lái)較大的誤差，故在數(shù)據(jù)分析時(shí)以誘蛾量除以10作為因變量Y的取值。如上文所述，我們建立穩(wěn)健的線性回歸模型為：Y=β0+β1X+γ+。模型中β0、β1和γ為3個(gè)未知的參數(shù)，其中β0為截距，即當(dāng)月相取值為0(下雨天氣)時(shí)的平均燈誘蛾量；β1為斜率，即當(dāng)月相每增加1個(gè)單位時(shí)，平均誘蛾量的變化。均值漂移參數(shù)γ=(γ1，…，γn)T是離群值指示器。如果估計(jì)的γi為非零，則第i個(gè)觀測(cè)值是離群值。除了參數(shù)估計(jì)，本研究還進(jìn)行了模型顯著性t檢驗(yàn)，用于檢驗(yàn)的2個(gè)假設(shè)分別為零假設(shè)H0:β1=0和備擇假設(shè)Ha:β1≠0。

圖1展示了2012年、2013年、2014年和2015年和2012—2015年匯總數(shù)據(jù)中誘蛾量與月相之間關(guān)系的散點(diǎn)圖。為了使散點(diǎn)圖更直觀地呈現(xiàn)模型的特征，各散點(diǎn)圖中添加了線性回歸趨勢(shì)線。圖1顯示誘蛾量與月相存在負(fù)相關(guān)關(guān)系，這一共性在各年份和匯總數(shù)據(jù)中都有明顯的體現(xiàn)。各年份的數(shù)據(jù)和2012—2015年匯總數(shù)據(jù)中都存在離群值，檢測(cè)出的離群值用星號(hào)標(biāo)記。而且，模型擬合出的回歸趨勢(shì)線并沒(méi)受離群值的影響，體現(xiàn)出所用線性回歸模型的穩(wěn)健性。

表1顯示了對(duì)于誘蛾量(Y)與月相(X)的穩(wěn)健線性回歸模型中非零的γ估計(jì)值。從中可以看出：2012年的數(shù)據(jù)中，第100、101、104、105和114個(gè)觀測(cè)值為離群值。2013年的數(shù)據(jù)中，第43和45個(gè)觀測(cè)值為離群值。2014年的數(shù)據(jù)中，第1、100和101個(gè)觀測(cè)值為離群值。2015年的數(shù)據(jù)中，第6個(gè)觀測(cè)值為離群值。2012—2015年的匯總數(shù)據(jù)中，第100和320個(gè)觀測(cè)值為離群值。其他的γ估計(jì)值為0，說(shuō)明其對(duì)應(yīng)的觀測(cè)值為正常值。

表2為誘蛾量(Y)與月相(X)的穩(wěn)健線性回歸分析的參數(shù)估計(jì)和假設(shè)檢驗(yàn)結(jié)果。從表中可以看出，在模型顯著性檢驗(yàn)中，各年份誘蛾量與月相的線性回歸模型都在5%顯著性水平下表現(xiàn)顯著。在各年份的數(shù)據(jù)和匯總數(shù)據(jù)中誘蛾量(Y)與月相(X)之間的線性回歸方程非常接近。而且，由各年份的數(shù)據(jù)和匯總數(shù)據(jù)得出的95%置信區(qū)間相互重疊，說(shuō)明誘蛾量(Y)與月相(X)之間的線性關(guān)系在不同的年份表現(xiàn)出相同的模式，且月相對(duì)誘蛾量的抑制程度在不同的年份沒(méi)有表現(xiàn)出顯著區(qū)別。以2012—2015年的匯總數(shù)據(jù)為例，晚上沒(méi)有月亮?xí)r，平均誘蛾量為169頭。月相在0～10的取值范圍內(nèi)，月相每增加1個(gè)單位，平均誘蛾量下降11頭。

圖中星號(hào)代表離群值。下同。The asterisks (*) represent the outliers. The same as bellow.圖1 2012、2013、2014、2015年和2012—2015年匯總數(shù)據(jù)中誘蛾量與月相之間關(guān)系的散點(diǎn)圖Fig.1 Scatter plots of moth yield and lunar phase in year 2012， 2013， 2014， 2015 and 2012-2015 combined data

表1誘蛾量與月相的穩(wěn)健線性回歸模型參數(shù)γ的估計(jì)

Table1Estimation of parameterγof robust linear regression model with moth yield and lunar phase

年份Year樣本容量Sample size離群值數(shù)Number of outliers非零的γ估計(jì)值Estimated nonzero γ20121165γ100=90.47,γ101=68.63,γ104=58.73,γ105=58.93,γ114=75.932013832γ43=51.31,γ45=47.5720141153γ1=40.32,γ100=28.39,γ101=30.0220151221γ6=112.022012—20154362γ100=91.46,γ320=111.69

表2誘蛾量與月相的穩(wěn)健線性回歸分析結(jié)果

Table2Results of robust linear regression analysis for the model with moth yield and lunar phase

年份Year數(shù)據(jù)量Sample size估計(jì)的回歸方程Estimated regression equationβ的95%置信區(qū)間95% confidence intervalt檢驗(yàn)P值P value of t test20122013201420152012-201511683115120436Y=19.107-1.888XY=23.364-1.500XY=13.280-1.083XY=15.639-1.079XY=16.900-1.178X(-3.294,-0.481)(-2.581,-0.418)(-1.844,-0.323)(-2.077,-0.081)(-1.722,-0.634)0.0090.0070.0060.034<0.001

2.2 溫度對(duì)誘蛾量的影響

把日平均溫度作為自變量(X)，誘蛾量作為因變量(Y)，分別對(duì)2012年、2013年、2014年和2015年的數(shù)據(jù)以及2012—2015年匯總數(shù)據(jù)做誘蛾量與日平均溫度之間的穩(wěn)健線性回歸分析。這部分所建立的線性回歸模型、參數(shù)含義和模型顯著性檢驗(yàn)和2.1節(jié)所述一樣。

圖2報(bào)告了各年份和2012—2015年匯總數(shù)據(jù)中誘蛾量與日平均溫度的散點(diǎn)圖。和圖1展示的一樣，在各圖中加入線性回歸趨勢(shì)線。圖2顯示誘蛾量與溫度存在負(fù)相關(guān)的關(guān)系，這一共性在各年份和匯總數(shù)據(jù)中都有明顯的體現(xiàn)。各年份的數(shù)據(jù)和2012—2015年匯總數(shù)據(jù)中都存在離群值，檢測(cè)出的離群值用星號(hào)標(biāo)記。而且，和前文所述一樣，模型擬合出的回歸直線并沒(méi)有受離群值的影響，體現(xiàn)出所用的線性回歸模型的穩(wěn)健性。

圖2 2012、2013、2014、2015年和2012—2015年匯總數(shù)據(jù)誘蛾量與日平均溫度的散點(diǎn)圖Fig.2 Scatter plots of moth yield and daily mean temperature in year 2012， 2013， 2014， 2015 and 2012-2015 combined data

表3報(bào)告了對(duì)于誘蛾量(Y)與溫度(X)的穩(wěn)健線性回歸模型中非零的γ估計(jì)值。從表3可以看出：2012年的數(shù)據(jù)中，第100、101、104、105和114個(gè)觀測(cè)值為離群值；2013年的數(shù)據(jù)中，第43、45、53、54和55個(gè)觀測(cè)值為離群值；2014年的數(shù)據(jù)中，第1、99和100個(gè)觀測(cè)值為離群值；2015年的數(shù)據(jù)中，第6個(gè)觀測(cè)值為離群值；2012—2015年的匯總數(shù)據(jù)中，第100和320個(gè)觀測(cè)值為離群值。

表4報(bào)告了對(duì)于誘蛾量(Y)與溫度(X)的穩(wěn)健線性回歸分析結(jié)果。從表4可以看出，各年份和匯總數(shù)據(jù)擬合的模型在5%顯著性水平下都表現(xiàn)為顯著。在各年份數(shù)據(jù)中誘蛾量(Y)與溫度(X)之間的線性回歸方程非常接近。而且，由各年份數(shù)據(jù)和匯總數(shù)據(jù)得出的β1的95%置信區(qū)間相互重疊，說(shuō)明誘蛾量(Y)與溫度(X)之間的線性關(guān)系在不同年份表現(xiàn)出相同的模式，且溫度對(duì)誘蛾量的抑制程度在不同年份沒(méi)有表現(xiàn)出顯著區(qū)別。以2012—2015年的匯總數(shù)據(jù)為例，當(dāng)日平均溫度為15～35 ℃時(shí)，溫度每增加1 ℃，平均誘蛾量減少10頭。

2.3 月相和溫度對(duì)誘蛾量的影響

在分析誘蛾量與單個(gè)因素(月相或溫度)之間關(guān)系的基礎(chǔ)上，我們還考慮在誘蛾量(Y)與月相(X1)之間的線性回歸模型中加入第2個(gè)變量——溫度(X2)。得到如下穩(wěn)健二元線性回歸模型：Y=β0+β1X1+β2X2+γ+。在對(duì)二元線性回歸模型進(jìn)行顯著性假設(shè)檢驗(yàn)中，月相(X1)在2012、2013、2014和2015年均表現(xiàn)為不顯著；月相和溫度在2012—2015年匯總數(shù)據(jù)中均表現(xiàn)為顯著。所以，我們只在2012—2015年匯總數(shù)據(jù)中同時(shí)考慮誘蛾量(Y)與月相(X1)和溫度(X2)以建立線性回歸模型。

圖3展示了2012—2015年的匯總數(shù)據(jù)中誘蛾量(Y)與月相(X1)和日平均溫度(X2)三維散點(diǎn)圖，在圖中2個(gè)離群值用星號(hào)標(biāo)注出來(lái)。從圖中的長(zhǎng)方形趨勢(shì)平面來(lái)看，月相和溫度都與誘蛾量呈負(fù)相關(guān)，說(shuō)明月相和溫度都對(duì)誘蛾量有抑制作用。對(duì)2012—2015年的匯總數(shù)據(jù)中誘蛾量(Y)與月相(X1)和日平均溫度(X2)建立穩(wěn)健的線性回歸模型進(jìn)行參數(shù)估計(jì)，非零的估計(jì)為：γ100=92.90，γ320=113.49，顯示第100和320個(gè)觀測(cè)值為離群值。在對(duì)模型的顯著性檢驗(yàn)中，二元線性回歸模型在5%顯著性水平下表現(xiàn)顯著(F檢驗(yàn)的P值為2.957e-07<0.05)。而且月相(X1)和溫度(X2)2個(gè)變量在單獨(dú)顯著性檢驗(yàn)中都表現(xiàn)為顯著(對(duì)X1和X2的t檢驗(yàn)P值分別為0.020 2和0.000 1)。估計(jì)出的線性回歸模型方程為Y=39.066 - 0.680X1-0.842X2。當(dāng)溫度保持恒定時(shí)，月相每增加1個(gè)單位，平均誘蛾量減少7頭；當(dāng)月相保持恒定時(shí)，日平均溫度每增加1 ℃，平均誘蛾量減少8頭。

表3誘蛾量與溫度的穩(wěn)健線性回歸模型參數(shù)γ的估計(jì)

Table3Estimation of parameterγof robust linear regression model with moth yield and temperature

年份Year數(shù)據(jù)量Sample size離群值數(shù)Number of outliers非零的γ估計(jì)值Estimated nonzero γ20122013201420152012—20151168311512243655312γ100=94.65, γ101=74.14, γ104=59.22, γ105=55.24, γ114=71γ43=57.37, γ45=53.12, γ53=38.48, γ54=41.09, γ55=48.84 γ1=42.83, γ99=29.82, γ100=29.99γ6=116.78γ100=93.28, γ320=115.05

表4誘蛾量與溫度的穩(wěn)健線性回歸分析結(jié)果

Table4Results of robust linear regression analysis for the model with moth yield and temperature

年份Year樣本容量Sample size估計(jì)的回歸方程Estimated regression equationβ的95%置信區(qū)間95% confidence intervalt檢驗(yàn)P值P value of t test20122013201420152012—201511683115122436Y=59.645-1.674XY=36.644-0.6973XY=29.627-0.725XY=40.234-1.004XY=42.484-1.023X(-2.403,-0.944)(-1.208,-0.186)(-1.188,-0.263)(-1.879,-0.129)(-1.422,-0.625)0.00010.00810.00240.0249<0.0001

圖3 2012—2015年匯總數(shù)據(jù)誘蛾量與月相和日平均溫度的三維散點(diǎn)圖Fig.3 Three dimensional scatter plot of moth yield and lunar phase and daily mean temperature for 2012-2015 combined data

3 結(jié)論與討論

盛承發(fā)等[9]通過(guò)圓形分析表明：從全月來(lái)看，紅鈴蟲(chóng)性誘劑誘蛾高峰集中角無(wú)統(tǒng)計(jì)意義；從分上、下半月來(lái)看，紅鈴蟲(chóng)的性誘蛾高峰數(shù)與月相有一定關(guān)系，即望、朔蛾峰均少，上弦(初7、初8)、下弦(22～23日)蛾峰均較多。但他僅分析了紅鈴蟲(chóng)的蛾峰與月相的關(guān)系。本研究將紅鈴蟲(chóng)全年的發(fā)生動(dòng)態(tài)與月相之間的關(guān)系進(jìn)行了分析，研究表明，紅鈴蟲(chóng)發(fā)生量在全年內(nèi)與月相無(wú)顯著相關(guān)性，但進(jìn)一步將天氣因素考慮進(jìn)去后發(fā)現(xiàn)，夜間光線對(duì)紅鈴蟲(chóng)發(fā)生量的影響更大，即夜間光線越強(qiáng)，紅鈴蟲(chóng)的誘蛾量就越小，反之則大。該結(jié)果表明，月相對(duì)紅鈴蟲(chóng)成蟲(chóng)活動(dòng)的影響主要由夜間光線引起，月相本身(月球引力)并不影響紅鈴蟲(chóng)的發(fā)生。

與光線相比較，溫度對(duì)昆蟲(chóng)發(fā)生的影響更大。一般而言，影響昆蟲(chóng)生長(zhǎng)發(fā)育的是有效積溫。在其發(fā)育起點(diǎn)溫度之上時(shí)，隨溫度的上升，昆蟲(chóng)生長(zhǎng)發(fā)育速率加快，導(dǎo)致昆蟲(chóng)的種群數(shù)量也隨之上升，但溫度過(guò)高也會(huì)增加其死亡率。在田間條件下，紅鈴蟲(chóng)的盛發(fā)期為每年的6月中下旬至10月上旬，在此期間，如果天氣呈現(xiàn)出高溫高濕條件，就會(huì)促進(jìn)紅鈴蟲(chóng)的發(fā)生，從而導(dǎo)致其發(fā)生量大，為害加重。但從本研究結(jié)果來(lái)看，溫度為15～35 ℃時(shí)，隨溫度的上升，紅鈴蟲(chóng)成蟲(chóng)的活動(dòng)顯著下降，即較高的溫度對(duì)紅鈴蟲(chóng)成蟲(chóng)活動(dòng)有抑制作用。由于溫度能影響棉鈴蟲(chóng)的性信息素分泌，如在22～30 ℃條件下，溫度越低，棉鈴蟲(chóng)性信息素的分泌量越多[16]。因此，在本研究中，高溫對(duì)紅鈴蟲(chóng)誘蛾量的抑制可能與紅鈴蟲(chóng)的性信息素分泌有關(guān)。

線性回歸模型是研究因變量與2個(gè)甚至幾個(gè)變量(因子)之間關(guān)系的主要統(tǒng)計(jì)分析方法，通常采用最小二乘法來(lái)估計(jì)回歸模型中的參數(shù)。但是最小二乘法對(duì)離群值很敏感，從而影響數(shù)據(jù)分析的結(jié)果。本研究采用了基于非凸懲罰的線性回歸模型來(lái)進(jìn)行分析，該方法具備同時(shí)進(jìn)行離群值檢測(cè)和穩(wěn)健參數(shù)估計(jì)的功能。分析結(jié)果表明，2012—2015年紅鈴蟲(chóng)的誘蛾量均存在離群值，顯示其不適合用最小二乘法分析。但在處理離群值時(shí)，本研究并沒(méi)有全部舍棄，而是根據(jù)實(shí)際情況進(jìn)行了分析，比如在2012年的數(shù)據(jù)中，分析出有離群值5個(gè)，分別為第100、101、104、105和114。本研究?jī)H舍棄了第100個(gè)數(shù)據(jù)，即9月7日的誘蛾量。原因是：第一，該數(shù)值在本年中最高，顯示其為蛾峰值；第二，其實(shí)際值與預(yù)測(cè)值之間的殘差最大。表明這一天的誘蛾量主要受紅鈴蟲(chóng)成蟲(chóng)發(fā)生高峰的影響，故本研究舍棄了該值。其他幾個(gè)離群值雖然實(shí)際值與預(yù)測(cè)值的殘差較大，但它們低于蛾峰值，且其與前后幾天的誘蛾量差異比較小，表明其主要受環(huán)境影響所致，故加以保留。同理，對(duì)2013年的數(shù)據(jù)僅舍棄了第43個(gè)，2014年的數(shù)據(jù)沒(méi)有舍棄，2015年的數(shù)據(jù)舍棄了第6個(gè)。這些值都為當(dāng)年的蛾峰值。經(jīng)此處理后，各年份的回歸模型t測(cè)驗(yàn)的P值均有顯著提升，進(jìn)一步證明了文中對(duì)離群值處理的合理性。同理，在處理溫度與誘蛾量回歸的離群值時(shí)，也僅舍棄了蛾高峰期的數(shù)據(jù)，其他均予以保留。

紅鈴蟲(chóng)發(fā)生除受月相影響外，受降雨的影響也較大。在一定條件下，紅鈴蟲(chóng)成蟲(chóng)的活動(dòng)與雨量呈正相關(guān)，即在較小雨量下，弱光和高濕有利于紅鈴蟲(chóng)成蟲(chóng)的活動(dòng)與產(chǎn)卵；但當(dāng)雨量過(guò)大時(shí)，紅鈴蟲(chóng)成蟲(chóng)的活動(dòng)受到抑制[17]。說(shuō)明紅鈴蟲(chóng)的發(fā)生非常復(fù)雜，是多種因素共同作用的結(jié)果。如能更系統(tǒng)地研究這些因子對(duì)紅鈴蟲(chóng)發(fā)生與行為的影響，將有助于深化對(duì)紅鈴蟲(chóng)發(fā)生規(guī)律的認(rèn)識(shí)，從而在生產(chǎn)上制定合適的防治對(duì)策。