成劍波，肖 琳，李云平，楊志敏*，陳玉成

（1.三峽庫區(qū)生態(tài)環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室，重慶 400715；2.西南大學(xué)資源環(huán)境學(xué)院，重慶 400716；3.重慶市煙草公司巫山分公司，重慶 404700）

數(shù)據(jù)挖掘能識(shí)別隱藏于數(shù)據(jù)間和數(shù)據(jù)中的寶貴信息[1]。Meta analysis（或稱整合分析、二次分析等）作為一種有效的數(shù)據(jù)挖掘方法[2]，現(xiàn)已廣泛用于社會(huì)科學(xué)、醫(yī)學(xué)、經(jīng)濟(jì)學(xué)和生態(tài)學(xué)等領(lǐng)域，如對(duì)各種生態(tài)效應(yīng)的研究[3]。它將不同的統(tǒng)計(jì)學(xué)定量法組織到同一架構(gòu)中，可完成某一領(lǐng)域不同研究成果的定量概化[4]。目前，Meta分析往往用于多個(gè)同類研究的結(jié)果的定量總結(jié)（宏觀），如某領(lǐng)域的文獻(xiàn)綜述，其用于單個(gè)研究的結(jié)果的定量表達(dá)鮮有報(bào)道（微觀）。筆者通過Meta分析重金屬對(duì)作物種子萌發(fā)的影響，以實(shí)現(xiàn)Pb2+對(duì)煙草種子萌發(fā)的影響的定量化和概化，探究 Meta分析在研究重金屬對(duì)作物生長(zhǎng)發(fā)育影響中的微觀應(yīng)用。

1 材料與方法

1.1 供試材料

試驗(yàn)所用煙草種子為中國(guó)煙草總公司重慶市公司提供的K326，所用的外源重金屬Pb為實(shí)驗(yàn)室配制的Pb(NO3)2溶液。

1.2 試驗(yàn)方法

將約1 cm厚的脫脂棉平鋪于18個(gè)直徑9 cm的培養(yǎng)皿內(nèi)，再將濾紙鋪在脫脂棉上，按每個(gè)培養(yǎng)皿40粒種子的標(biāo)準(zhǔn)，將種子均勻分布在各個(gè)培養(yǎng)皿內(nèi)。將培養(yǎng)皿均分為C0、C1、C2、C3、C4、C56組，其對(duì)應(yīng)的 Pb2+濃度分別為0、0.2、0.5、1.0、3.0、5.0 mg/L；將培養(yǎng)皿置于恒溫箱中，在23 ℃、80%的相對(duì)濕度和13 h/d的光照條件下，定期加入純水（C0組）和Pb2+溶液，以保持培養(yǎng)床的濕潤(rùn)狀態(tài)；從種子萌芽時(shí)（第 11天）開始，每天記錄每個(gè)濃度的3個(gè)平行樣本CX1、CX2和CX3的已發(fā)芽總數(shù)；于第 20天量取已發(fā)芽種子的根長(zhǎng)和莖長(zhǎng)并結(jié)束試驗(yàn)。

1.3 Meta分析方法

Meta分析實(shí)為一類數(shù)學(xué)和統(tǒng)計(jì)學(xué)方法的多層次組織方式，其作用因方法和方法組合方式不同而不同，且實(shí)用范圍廣。因此，它本身也是一項(xiàng)研究,需要認(rèn)真設(shè)計(jì)[5]。一般地，Meta分析通過效應(yīng)指標(biāo)量化研究結(jié)果（數(shù)據(jù)）所含的規(guī)律，通過加權(quán)整合概化同類研究的不同結(jié)果，進(jìn)而較直觀簡(jiǎn)略地表達(dá)客觀規(guī)律。上述過程存在著效應(yīng)指標(biāo)的選擇、各待整合效應(yīng)間異質(zhì)性的分析、效應(yīng)整合方式的選擇等主要環(huán)節(jié)，同時(shí)也隱含著偏性等問題?，F(xiàn)有 lnR、Hedges's d、risk ratio等指標(biāo)和 RevMan、Stata、CMA等軟件可幫助完成Meta分析?；诎l(fā)芽數(shù)、根長(zhǎng)、莖長(zhǎng)和體長(zhǎng)（葉除外）數(shù)據(jù)，為表征Pb2+對(duì)種子萌發(fā)過程（萌芽+發(fā)育）的影響，本研究在計(jì)算 Pb2+組對(duì)比純水組而得到的分效應(yīng)e（lnOR和lnR）后，借助CMA[6]整合了分效應(yīng)e得到了綜合效應(yīng)E。以本試驗(yàn)為例，視加Pb2+液為處理，加純水為對(duì)照[7]，簡(jiǎn)介本文涉及的效應(yīng)指標(biāo)機(jī)會(huì)odds、機(jī)會(huì)比OR和響應(yīng)比 R[6]：odds = n已萌芽/n未萌芽，可表現(xiàn)出某環(huán)境下種子萌芽相對(duì)不萌芽?jī)?yōu)勢(shì)或萌芽的機(jī)會(huì)；OR =oddsT/oddsC，可表現(xiàn)出Pb2+液對(duì)萌芽機(jī)會(huì)的影響；可表現(xiàn)出 Pb2+對(duì)根發(fā)育（或莖發(fā)育、體長(zhǎng)）的影響。在表征萌發(fā)過程時(shí)，之所以納入體長(zhǎng)，是因?yàn)椋海?）雖單粒種子的體長(zhǎng) = 根長(zhǎng)+莖長(zhǎng)，但數(shù)據(jù)預(yù)處理會(huì)丟掉各數(shù)據(jù)組的離群值和異常離群值，則單粒種子的根長(zhǎng)、莖長(zhǎng)和體長(zhǎng)不一定同時(shí)被采納，故大量種子的體長(zhǎng) ≠ 根長(zhǎng)+莖長(zhǎng)；（2）相對(duì)于根、莖長(zhǎng)，體長(zhǎng)更能反映種子對(duì)包衣營(yíng)養(yǎng)物質(zhì)的吸收和轉(zhuǎn)化能力；（3）從更多角度去分析結(jié)果，可更客觀地還原試驗(yàn)過程。在量化Pb2+液的影響時(shí)，之所以未直接使用OR和R而采用lnOR和lnR，是因?yàn)閇7-9]：（a）對(duì)數(shù)化指標(biāo)能平等反映正效應(yīng)和負(fù)效應(yīng)，如發(fā)芽機(jī)會(huì)增至100倍或減為1/100其效應(yīng)大小均為4.6；（b）對(duì)數(shù)化指標(biāo)能等效反映T和C而不偏重任一方，對(duì)數(shù)化前，若分子/分母變?yōu)椋ǚ肿印?00）/分母則該變化效應(yīng)值為﹢99（由 100-1算得），若變?yōu)榉肿?（分母×100）則該變化效應(yīng)值僅為-0.99而非-99，對(duì)數(shù)化后，兩種變化的效應(yīng)大小均為 4.6；（c）由比值型指標(biāo)得到的效應(yīng)值構(gòu)成的小樣本為偏分布而非正態(tài)分布，對(duì)數(shù)化后的樣本會(huì)更接近正態(tài)分布，對(duì)客觀規(guī)律更具表征力。

Pb2+在每個(gè)濃度處對(duì)種子萌芽影響的計(jì)算過程分 3步：（1）在 CMA中，視加 Pb2+液為 T處理;加純水為 C處理，選用 lnOR從 10個(gè)時(shí)間尺度（10～19 d）計(jì)算 3個(gè)培養(yǎng)皿的發(fā)芽總數(shù)數(shù)據(jù)（圖1a），得到分效應(yīng)e10～e19；（2）在CMA中判定分效應(yīng)間的異質(zhì)性[10]；（3）結(jié)合異質(zhì)性判定結(jié)果，用CMA的隨機(jī)模型將同質(zhì)分效應(yīng)整合為綜合效應(yīng)E。Pb2+在每個(gè)濃度處對(duì)種子萌發(fā)（根發(fā)育、莖發(fā)育或體長(zhǎng)）影響的計(jì)算過程分5步：（1）用SPSS分別丟掉3個(gè)培養(yǎng)皿的長(zhǎng)度數(shù)據(jù)中的離群值和異常離群值；（2）視加Pb2+液為T處理，加純水為C處理，為每個(gè)培養(yǎng)皿設(shè)計(jì)3個(gè)對(duì)比（圖1b）；（3）根據(jù)式1和2[7]，用Excel算出3×3個(gè)分效應(yīng)e1～e9（lnR）及其誤差（νlnR）；（4）在CMA中判定分效應(yīng)間異質(zhì)性[10]；（5）用CMA的隨機(jī)模型將同質(zhì)分效應(yīng)整合為綜合效應(yīng)E。

圖1 種子萌發(fā)中Pb2+所產(chǎn)生效應(yīng)的構(gòu)造Fig.1 Calculation design for quantifying Pb2+ effect on seeds bourgeon

2 結(jié)果與討論

2.1 Pb2+對(duì)煙草種子萌芽的影響

0.2 和0.5 mg/L的各分效應(yīng)組均無異質(zhì)性，1、3和5 mg/L的各分效應(yīng)組的I2均在25%～50%之間屬中度異質(zhì)（表1），說明對(duì)各組e10-e19的整合均有效，所得的E均能表征Pb2+對(duì)種子萌芽的影響。每個(gè)E的95%CI的上限都小于0，表明各濃度 Pb2+液均能降低種子萌芽的機(jī)會(huì)，能明顯抑制并推遲種子萌芽（表1和圖2）。在（0，1.5 mg/L）內(nèi)，抑制作用隨濃度增加而快速增大。在1.5 mg/L處，抑制作用升至最大，效應(yīng)值達(dá)到-1.58。在（1.5，5.0 mg/L）內(nèi)，抑制作用緩慢回落至較穩(wěn)定水平，效應(yīng)值定在-0.8左右，但仍比0.2 mg/L大80%。

表1 Pb2+對(duì)種子萌芽的影響Table1 Pb2+ effect on seed germination

圖2 Pb2+對(duì)種子萌芽的影響Fig.2 Pb2+ effect on seed germination

2.2 Pb2+對(duì)煙草種子發(fā)育的影響

各濃度3個(gè)效應(yīng)組I2均低于50%，說明各效應(yīng)組均無顯著異質(zhì)性，對(duì)各組e1-e9的整合均有效，所得的E均能表征Pb2+對(duì)種子發(fā)育的影響（表2）。

根的響應(yīng)均為負(fù)，表明各濃度Pb2+溶液均會(huì)抑制根的發(fā)育（表2和圖3a）。所有濃度（除0.2 mg/L外）的響應(yīng)上限均小于 0，說明它們對(duì)根的發(fā)育有明顯的抑制作用。從0.2 mg/L開始，抑制作用增強(qiáng)。到0.7 mg/L處，響應(yīng)變?yōu)?0.8，抑制作用達(dá)到較大水平。在隨后的濃度翻倍后，響應(yīng)減為-0.3。在隨后的濃度再翻倍后，抑制作用再次增大。3 mg/L后，抑制作用繼續(xù)增大但增速變緩。故認(rèn)為0.7 mg/L和1.5 mg/L為抑制作用的拐點(diǎn)。

各濃度Pb2+溶液對(duì)莖的發(fā)育均有明顯的抑制作用，抑制作用的變化在不同濃度區(qū)間差異較大（表2和圖 3b）。經(jīng)過 0.2～0.7 mg/L后抑制作用增加109.68%，而經(jīng)過0.7～5.0 mg/L后抑制作用僅僅增加17.85%。到0.7 mg/L處，響應(yīng)為-0.33，抑制作用較大。到5 mg/L響應(yīng)為-0.38，抑制作用最大。

體長(zhǎng)的響應(yīng)趨勢(shì)雖與根、莖大體相同，但仍有差別（表2和圖3）。0.5、1.0、3.0和5 mg/L對(duì)體長(zhǎng)的抑制作用大小處于根、莖之間，但0.2 mg/L對(duì)體長(zhǎng)的抑制作用大小卻大于根和莖，說明0.2 mg/L處對(duì)體長(zhǎng)的影響較乏規(guī)律性，或說明在0.2 mg/L處引入了較多隨機(jī)因素。

表2 Pb2+對(duì)種子幼芽的影響Table2 Pb2+ effect on seed growth

圖3 Pb2+對(duì)種子幼芽的影響Fig.3 Pb2+ effect on seed growth

2.3 Pb2+對(duì)煙草種子萌發(fā)影響的波動(dòng)性分析

從各效應(yīng)圖（圖2和圖3）可看出，各效應(yīng)值的95%CI上下限間構(gòu)成了較為明顯的作用帶，作用帶中的任何一點(diǎn)均可能代表真實(shí)作用。作用帶越寬，該作用越不穩(wěn)定。隨著Pb2+液對(duì)種子萌發(fā)抑制作用的增大，作用帶寬也趨向于增大，Pb2+作用的波動(dòng)性隨之增大（圖4）。大部分散點(diǎn)在趨勢(shì)線以下，說明在0.2～5 mg/L內(nèi)上述波動(dòng)性總體較為保守。0.2 mg/L的散點(diǎn)多在趨勢(shì)線以上，說明在0.2 mg/L處上述波動(dòng)性較為激烈，這與研究體長(zhǎng)的響應(yīng)時(shí)的發(fā)現(xiàn)相符。

圖4 作用帶寬隨作用大小的變化Fig.4 Volatility values of Pb2+ effects

3 結(jié) 論

初步認(rèn)為 Meta分析能量化并概化重金屬對(duì)作物生長(zhǎng)發(fā)育的影響，可用于重金屬對(duì)作物生長(zhǎng)發(fā)育影響的研究。0.2～5.0 mg/L的Pb2+液會(huì)抑制煙草種子的萌芽和發(fā)育，對(duì)萌芽和發(fā)育的抑制作用均在0.2 mg/L處最小。隨著濃度的增大，對(duì)萌芽的抑制作用呈“增→減→穩(wěn)”，對(duì)根、莖發(fā)育的抑制作用均呈“增→減→增”且對(duì)根的抑制約為對(duì)莖的 1.7～3.4倍，種子對(duì)營(yíng)養(yǎng)物的轉(zhuǎn)化受到抑制且末端抑制作用為起始處的 5.7倍。此外，Pb2+作用的波動(dòng)性隨作用增大而增大，總體較為保守，但在0.2 mg/L處較為激烈。

[1]McCue C.Data Mining and Predictive Analysis[M].Burlington: Butterworth-Heinemann, 2007: 25.

[2]E.Bayam, J Liebowitz, W Agresti.Older drivers and accidents: A meta analysis and data mining application on traffic accident data[J].Expert Systems with Applications, 2005, 29(3): 598-629.

[3]Gavin Stewart.Meta-analysis in applied ecology[J].Biology Letters, 2010(6): 78-81.

[4]Arnqvist G, D Wooster.Meta-analysis: synthesizing research findings in ecology and evolution[J].Trends in Ecology & Evolution, 1995, 10(6): 236-240.

[5]彭少麟，鄭鳳英.Meta分析及MetaWin軟件[J].土壤與環(huán)境，1999，8（4）：295-299.

[6]Borenstein M, Hedges L, Higgins J, et al.Comprehensive Meta Analysis, version 2.0: introduction to the program[CP/OL].Comprehensive Meta Analysis, 2005:[2011-7-20].http://www.meta- analysis.com.

[7]H Hillebrand.Meta-analysis in Ecology [J/OL].Encyclopedia of life science, 2008: [2011-7-20].http://onlinelibrary.wiley.com/doi/10.1002/9780470015902.a0003272.

[8]Larry V.Hedges, Jessica Gurevitch, Peter S.Curtis.The meta-analysis of response ratios in experimental ecology[J].Ecology, 1999, 80(4): 1150-1156.

[9]Chen Xuezhong, Yin Xiangchu.Application of the Load/Unload Response Ratio Theory to Predictions of Moderate Earthquakes (7.0＞M≥6.0)[J].Earthquake Research in China, 1995, 9(3): 225-229.

[10]Julian P T Higgins, Simon G Thompson, Jonathan J Deeks,et al.Measuring inconsistency in meta-analyses[J].BMJ,2003, 327(7414): 557-560.