劉鑒毅 ，劉建良，林國(guó)源，謝才豐

（1.中國(guó)水產(chǎn)科學(xué)研究院東海水產(chǎn)研究所長(zhǎng)江口漁業(yè)生態(tài)重點(diǎn)實(shí)驗(yàn)室，上海 200090；2.興鯢農(nóng)業(yè)生物科技發(fā)展有限公司，上海 200135；3大龍園農(nóng)業(yè)科技發(fā)展有限公司，廣東 深圳 440300）

性成熟子二代大鯢親體全長(zhǎng)與體質(zhì)量的回歸分析

劉鑒毅1，劉建良2，林國(guó)源2，謝才豐3

（1.中國(guó)水產(chǎn)科學(xué)研究院東海水產(chǎn)研究所長(zhǎng)江口漁業(yè)生態(tài)重點(diǎn)實(shí)驗(yàn)室，上海 200090；2.興鯢農(nóng)業(yè)生物科技發(fā)展有限公司，上海 200135；3大龍園農(nóng)業(yè)科技發(fā)展有限公司，廣東 深圳 440300）

為性成熟子二代大鯢親體的人工繁育選苗提供科學(xué)指導(dǎo)，本文通過對(duì)陜西省漢中市和廣東大有大鯢馴養(yǎng)繁殖有限公司自繁123尾性成熟子二代親體的體質(zhì)量和全長(zhǎng)進(jìn)行實(shí)際測(cè)量，利用matlab分析軟件對(duì)數(shù)據(jù)進(jìn)行回歸建模分析。結(jié)構(gòu)表明，性成熟子二代體質(zhì)量和全長(zhǎng)直接存在極顯著的相關(guān)關(guān)系，子二代養(yǎng)殖大鯢親體的體質(zhì)量和全長(zhǎng)直接的主要關(guān)系可以用冪函數(shù)：y=0.00001129x3.312和三次曲線：y=39400-136.2x＋0.1521x2-0.00004544x3表達(dá)。

大鯢；子二代親體；體質(zhì)量；體長(zhǎng)；性成熟

大鯢（Andrias davidianus）是我國(guó)特有的珍稀兩棲動(dòng)物，是國(guó)家二級(jí)保護(hù)動(dòng)物。在野生大鯢資源量急劇下降的情況下，為保護(hù)野生資源，實(shí)現(xiàn)大鯢資源的增殖，開展大鯢人工繁殖非常重要。國(guó)內(nèi)學(xué)者先后對(duì)大鯢的人工繁育、保護(hù)遺傳和孵化率等[1-4]方面進(jìn)行了研究工作。

在水產(chǎn)養(yǎng)殖中，通常根據(jù)親本的體長(zhǎng)和體質(zhì)量來進(jìn)行質(zhì)量?jī)?yōu)劣的判別，達(dá)到一定體長(zhǎng)規(guī)格的苗種也應(yīng)達(dá)到相應(yīng)的體質(zhì)量，才被認(rèn)同為優(yōu)質(zhì)苗種[5]。作為兩棲類動(dòng)物的大鯢也是這樣。葛蔭榕等[6]對(duì)國(guó)內(nèi)野生大鯢年齡與生長(zhǎng)發(fā)育關(guān)系進(jìn)行了分析，認(rèn)為野外大鯢IV齡時(shí)已經(jīng)基本性成熟，并且可以產(chǎn)卵，最小體質(zhì)量約0.35 kg。但羅慶華等[7]人研究認(rèn)為野外大鯢在V齡左右時(shí)達(dá)到性成熟，而人工養(yǎng)殖大鯢為IV齡左右。王啟軍等[8]和朱忠勝等[9]先后對(duì)養(yǎng)殖大鯢年齡與生長(zhǎng)發(fā)育關(guān)系進(jìn)行了分析，亦認(rèn)為IV齡時(shí)體質(zhì)量和體長(zhǎng)可以作為人工養(yǎng)殖大鯢的繁殖年齡。但是對(duì)于性成熟的養(yǎng)殖個(gè)體，體長(zhǎng)和體質(zhì)量的相關(guān)關(guān)系的報(bào)道尚未見到。因此本文采用123尾子二代性成熟產(chǎn)卵大鯢親體的全長(zhǎng)和體長(zhǎng)的數(shù)據(jù)，采用線性回歸、二次曲線和冪函數(shù)等多種模型構(gòu)建回歸方程，為養(yǎng)殖效果評(píng)定提供安全、快速、準(zhǔn)確的衡量方法，為人工繁育提供科學(xué)指導(dǎo)。

1 材料和方法

1.1 大鯢

實(shí)驗(yàn)測(cè)量樣本來自陜西漢中和廣東大有大鯢馴養(yǎng)繁殖有限公司自繁個(gè)體，總數(shù)為123尾，測(cè)量對(duì)象都是子二代養(yǎng)殖親體。

1.2 測(cè)量方法

本文體長(zhǎng)為全長(zhǎng)，指自大鯢嘴部最前端至尾末端長(zhǎng)度。用數(shù)顯式電子秤（精確度為0.1 g）稱量大鯢體質(zhì)量，直尺（精確度0.1 cm）測(cè)量大鯢的體長(zhǎng)。

1.3 數(shù)據(jù)處理

本文采用Pearson相關(guān)系數(shù)模型檢驗(yàn)子二代性成熟產(chǎn)卵大鯢的全長(zhǎng)和體長(zhǎng)相關(guān)性。采用線性回歸、二次曲線回歸、冪函數(shù)、三次曲線回歸和指數(shù)回歸等5個(gè)回歸模型建立二代性成熟產(chǎn)卵大鯢全長(zhǎng)和體質(zhì)量模型。采用均方根誤差（Root mean square error，RMSE）

和回歸模型的擬合決定系數(shù)R2作為檢驗(yàn)回歸模型擬合度的指標(biāo)。（1）式中y為實(shí)驗(yàn)測(cè)量值，y?為回歸模型擬合值，n為樣本個(gè)數(shù)。

2 結(jié)果

通過對(duì)123尾子二代性成熟產(chǎn)卵大鯢親體的測(cè)量記錄數(shù)據(jù)進(jìn)行軟件分析，得到全長(zhǎng)與體質(zhì)量的Pearson相關(guān)系數(shù)為0.93；sigf值均為0（即P=0），差異極具有統(tǒng)計(jì)學(xué)意義（P＜0.01）。因此從本次樣本數(shù)據(jù)的分析結(jié)果看，子二代親的全長(zhǎng)與體質(zhì)量二者之間具有極強(qiáng)的正相關(guān)性；將體質(zhì)量和全長(zhǎng)分別設(shè)定為因變量和自變量，經(jīng)過回歸分析后，得到如下的回歸分析結(jié)果（表1）。

表1 回歸模型的參數(shù)和檢驗(yàn)結(jié)果

由表1可以看出，5個(gè)回歸模型的p值均為0，因此所有構(gòu)建的回歸模型方程具有極高的顯著性，不能拒絕模型原假設(shè)。5個(gè)模型方程如下：

線性：y=-14210+23.9x

二次曲線：y=11100－35.25x+0.03392x2

冪函數(shù)：y=0.000001129x3.312

三 次 曲 線 ：y=39440 －136.2x+0.152x2－0.00004544x3

指數(shù)函數(shù)：y=268.3e0.00358x通過模型擬合度指標(biāo)決定系數(shù)和均方根誤差結(jié)果（表1）表明，5種回歸模型大致相當(dāng)，線性模型的擬合優(yōu)劣程度相對(duì)要弱一點(diǎn)，其他4個(gè)模型的R2都在0.890～0.895，但指數(shù)模型的RMSE要明顯高于二次曲線、冪函數(shù)和三次曲線模型。綜合2個(gè)指標(biāo)考慮，冪函數(shù)和三次曲線模型都可以作為子二代性成熟產(chǎn)卵大鯢親體的全長(zhǎng)和體質(zhì)量的回歸模型。圖1給出了冪函數(shù)回歸模型擬合的二區(qū)曲線圖，三次曲線的回歸曲線圖類似。

圖1 二代性成熟產(chǎn)卵大鯢的全長(zhǎng)和體質(zhì)量曲線關(guān)系圖

3 討論

目前在水產(chǎn)養(yǎng)殖中描述體長(zhǎng)與體質(zhì)量關(guān)系時(shí)應(yīng)用最多的是冪函數(shù)關(guān)系（W=bxa）[5-9]，本文的研究結(jié)果與葛蔭榕等前人研究相似，均成冪函數(shù)關(guān)系。但本文的冪函數(shù)表達(dá)式的參數(shù)與前人的研究結(jié)果不一致，主要是因?yàn)椴蓸訑?shù)據(jù)的不同。本文的采樣數(shù)據(jù)測(cè)量對(duì)象是子二代養(yǎng)殖大鯢親體，測(cè)量的全長(zhǎng)，王文林等和葛蔭榕等的報(bào)道結(jié)果顯示，他們均是研究的大鯢體長(zhǎng)和體質(zhì)量關(guān)系[9-10]。在規(guī)模化養(yǎng)殖中，養(yǎng)殖人員往往對(duì)全長(zhǎng)和體長(zhǎng)混為一談，因此不同的人員測(cè)量，測(cè)量的結(jié)果往往不一樣，給模型統(tǒng)計(jì)帶來極大的影響。在實(shí)際測(cè)量中，包含尾部的測(cè)量往往更容易更準(zhǔn)確，因此本文采用全長(zhǎng)數(shù)據(jù)建立和體質(zhì)量的回歸關(guān)系。

此外本文的采樣數(shù)據(jù)是子二代性成熟產(chǎn)卵大鯢，是經(jīng)過長(zhǎng)期養(yǎng)殖的個(gè)體。王文林等和葛蔭榕等的采樣的數(shù)據(jù)是野生個(gè)體，包含不同年齡。人工養(yǎng)殖的大鯢性成熟年齡和野外的研究結(jié)果并不一致，因?yàn)樵谌斯ゐB(yǎng)殖條件下，水溫相對(duì)恒定，餌料充足，基本克服冬眠的現(xiàn)象，相同年齡的人工養(yǎng)殖大鯢生長(zhǎng)比野外的更快，體質(zhì)量和體長(zhǎng)的等生長(zhǎng)參數(shù)比野外更大。王啟軍、朱忠勝等[7-8]人雖然也是養(yǎng)殖個(gè)體，但他們測(cè)量的對(duì)象包含為成年和未性成熟的小個(gè)體。大鯢在不同年齡段的體質(zhì)量和體長(zhǎng)的生長(zhǎng)速度是不一樣的。朱忠勝和于振海[8，11]研究表明I-II齡段的養(yǎng)殖大鯢體質(zhì)量和體長(zhǎng)的增速最快，之后增速逐漸下降。因此其體質(zhì)量和體長(zhǎng)關(guān)系和成熟的個(gè)體不一樣。

王啟軍等[7]人采用不同模型建立養(yǎng)殖大鯢全長(zhǎng)和體質(zhì)量關(guān)系，認(rèn)為冪函數(shù)模型擬合最優(yōu)，其次是二次曲線模型和三次曲線模型，冪函數(shù)曲線應(yīng)該作為養(yǎng)殖大鯢全長(zhǎng)和體質(zhì)量關(guān)系的最佳曲線。本文冪函數(shù)建立的模型和二次曲線、三次曲線模型相似，這種差別主要還是本文采用數(shù)據(jù)是子二代性成熟產(chǎn)卵大鯢，這期間的個(gè)體，尤其是大齡個(gè)體的體質(zhì)量和體長(zhǎng)接近生長(zhǎng)極限，其體質(zhì)量和體長(zhǎng)的增長(zhǎng)絕對(duì)量和增長(zhǎng)倍數(shù)都下降，階段性的回歸關(guān)系也符合拋物線函數(shù)。而王啟軍等人采用數(shù)據(jù)包含未成熟的小個(gè)體，在生長(zhǎng)拐點(diǎn)前，體質(zhì)量和體長(zhǎng)的回歸關(guān)系更符合冪函數(shù)。

無論是在野生狀態(tài)下還是在人工養(yǎng)殖狀態(tài)下，大鯢的生長(zhǎng)趨勢(shì)都是符合冪函數(shù)的生長(zhǎng)規(guī)律。根據(jù)本文建立的回歸模型，對(duì)親本子二代親體選育過程中，可以在水池中測(cè)量個(gè)體全長(zhǎng)，然后估計(jì)個(gè)體體質(zhì)量，避免了以往的傳統(tǒng)捕捉、稱量過程所造成的應(yīng)激反應(yīng)和應(yīng)激性疾病的發(fā)生，為今后大鯢養(yǎng)殖效果評(píng)定提供了一個(gè)安全、快速、準(zhǔn)確的方法。

[1]陽(yáng)愛生，卞偉，劉運(yùn)清，等.大鯢胚胎發(fā)育的初步研究[J].動(dòng)物學(xué)報(bào)，1983（1）：42-47.

[2]章克家，王小明，吳巍，等.大鯢保護(hù)生物學(xué)及其研究進(jìn)展[J].生物多樣性，2002，10（3）：291-297.

[3]劉鑒毅，譚永安，盧興孫，等.中國(guó)大鯢子二代規(guī)?；斯し趸夹g(shù)的研究[J].經(jīng)濟(jì)動(dòng)物學(xué)報(bào)，2005，（3）：152-155，159.

[4]鄧捷，王啟軍，趙虎，等.池養(yǎng)成年大鯢體質(zhì)量和性別發(fā)育與水溫的相關(guān)性研究 [J].基因組學(xué)與應(yīng)用生物學(xué)，2016（9）：2342-2346.

[5]陳友明，李瀟軒，王明華，等.斑點(diǎn)叉尾鮰苗種全長(zhǎng)與體質(zhì)量的關(guān)系[J].水產(chǎn)養(yǎng)殖，2010，31（5）：13-15.

[6]葛蔭榕，鄭合勛，李繼海.大鯢年齡與生長(zhǎng)的初步研究[J].河南師范大學(xué)學(xué)報(bào)（自然科學(xué)版），1995，23（1）：59-63.

[7]羅慶華，盧成英，劉清波.中國(guó)大鯢繁殖生物學(xué)研究進(jìn)展[J].四川動(dòng)物，2005，（3）：417-420.

[8]王啟軍，趙虎，張紅星，等.人工養(yǎng)殖大鯢全長(zhǎng)與體質(zhì)量關(guān)系的回歸分析[J].基因組學(xué)與應(yīng)用生物學(xué)，2012，31（4）：381-384

[9]朱忠勝，李媛，姚俊杰，等.貴州人工養(yǎng)殖大鯢年齡與體質(zhì)量和體長(zhǎng)的關(guān)系[J].貴州農(nóng)業(yè)科學(xué)，2015，（9）：142-144.

[10]王文林，詹克慧，陳平，等．池養(yǎng)大鯢生長(zhǎng)發(fā)育的初步研究[J].淡水漁業(yè)，1994（4）：20-22.

[11]于振海，陳有光，靖瑩，等．中國(guó)大鯢地下室微流水人工養(yǎng)殖技術(shù)初步研究[J].上海海洋大學(xué)學(xué)報(bào)，2013（1）：60-65.

江蘇宿遷泗陽(yáng)縣工業(yè)化養(yǎng)殖每年畝產(chǎn)成魚30萬千克

通常，每667 m2水面年產(chǎn)成魚1 500 kg，就算高產(chǎn)魚塘了?？墒?，令人驚訝的是，江蘇宿遷市泗陽(yáng)縣綠興水產(chǎn)養(yǎng)殖專業(yè)合作社每667 m2水面年產(chǎn)成魚高達(dá)30萬kg!

泗陽(yáng)縣綠興水產(chǎn)養(yǎng)殖專業(yè)合作社養(yǎng)魚高產(chǎn)，源于工廠化的生態(tài)養(yǎng)殖模式。2017年9月，筆者采訪時(shí)，正趕上合作社起魚。工人們站在魚池兩邊，從西頭拉著漁網(wǎng)慢慢后退，將魚趕進(jìn)東側(cè)的凈化區(qū)。工人按動(dòng)電鈕后，凈化區(qū)的自動(dòng)裝置由南向北推進(jìn)，再將魚趕進(jìn)北頭的網(wǎng)箱里。最后，工人撈起網(wǎng)箱里的魚過磅，整個(gè)起魚過程不到半個(gè)小時(shí)。

泗陽(yáng)縣綠興水產(chǎn)養(yǎng)殖專業(yè)合作社擁有水面80×667 m2，其中魚池面積不足3×667 m2?！拔覀兊墓S化生態(tài)養(yǎng)魚技術(shù)是從國(guó)外引進(jìn)的，共有6項(xiàng)專利技術(shù)，能夠?qū)崿F(xiàn)高產(chǎn)，而且具有生態(tài)水循環(huán)功能。”合作社董事長(zhǎng)周秀紅介紹。在漁場(chǎng)西側(cè)，可以看到，由南到北排列著寬5 m、長(zhǎng)22 m、大小一致的16個(gè)池口。魚池的西首，安裝著水循環(huán)系統(tǒng)，7臺(tái)不停轉(zhuǎn)動(dòng)的泵機(jī)產(chǎn)生的高強(qiáng)度氣體推動(dòng)著水流，源源不斷流進(jìn)魚池內(nèi)，緩緩流向東出口，池內(nèi)眾多魚則逆流聚攏。

循環(huán)水養(yǎng)魚，水質(zhì)如何保證?合作社飼養(yǎng)員周家明介紹，在7臺(tái)泵機(jī)“以氣推流”作用下，水體以每秒0.3～0.7 m的速度流動(dòng)，池水流動(dòng)，池魚不會(huì)缺氧，就不會(huì)造成魚浮頭。從魚池流出的水體先后進(jìn)入一級(jí)凈化區(qū)、二級(jí)凈化區(qū)、生物凈化區(qū)等6道凈化區(qū)，通過一圈循環(huán)凈化后，水質(zhì)恢復(fù)清澈，又成為魚池的新水源。

生態(tài)環(huán)保是工廠化高產(chǎn)養(yǎng)魚的新概念。周家明說，池魚排放的糞便很多，如果不清理就會(huì)影響水質(zhì)，也保證不了魚的品質(zhì)。因此，合作社在每個(gè)魚池都裝了吸污系統(tǒng)，喂食后一個(gè)半小時(shí)到兩個(gè)小時(shí)，工作人員開動(dòng)吸污系統(tǒng)，魚池底下的魚糞通過8根管道被吸到100 m外的沉淀池。沉淀后的魚糞被集中到島上再次沉淀，然后作為農(nóng)業(yè)生產(chǎn)的有機(jī)肥，實(shí)現(xiàn)了工廠化養(yǎng)魚的零排放。

工廠化生態(tài)養(yǎng)魚每年畝產(chǎn)30萬kg，究竟是怎么實(shí)現(xiàn)的?周家明說，工廠化養(yǎng)魚的最大特點(diǎn)就是用地少，高密度。合作社16個(gè)魚池，每個(gè)不到123 m2，分別放養(yǎng)鯽魚、花鰱、白鰱、鯉魚和草魚，每池每批可起成魚2萬kg左右，每年可出欄成魚3批。

“我們引進(jìn)工廠化生態(tài)養(yǎng)魚技術(shù)不是光顧自己賺錢，而是想給村民做個(gè)示范，讓大家都來學(xué)習(xí)，一起賺錢?！敝苄慵t說。

（www.bbwfish.com）

Regression analysis of whole length and body weight of sexual maturity giant salamanders(Andrias davidianus)F2

Liu Jianyi1，Liu Jianliang2，Lin Guoyuan2，Xie Caifeng3

（1.Key Laboratory of Fisheries Ecology of the Yangtze Estuary，China Academy of Fishery Sciences，College of Fisheries and Life Science，Shanghai 200090，China；2.Shanghai Xing salamander agricultural biotechnology development Limited，Shanghai 200135，China；3.Da Long Park Agricultural Science and Technology Development Limited，Shenzhen 440300，China）

In order to support the artificial breeding of sexual maturity giant samanders(Andrias davidianus)F2，this paper analyzed the whole length as well as the weight of 123 sexual maturity giant samanders F2using regression model by matlab software.The results indicated that there was a significant correlative relation between the weight and the whole length of sexual maturity giant samanders F2.The relationship between the weight and the whole length of sexual maturity giant samanders F2could be modeled by Power，Quadratic and cubic polynomial curve.

giant samanders；F2parent；weight；whole length；sexual maturity

Q174

A

1004-2091（2017）10-0034-04

10.3969/j.issn.1004-2091.2017.10.007

上?？萍寂d農(nóng)項(xiàng)目[滬農(nóng)科字（2015）第18號(hào)]；深圳市科技計(jì)劃項(xiàng)目（CXZZ20120828143309692）

劉鑒毅（1965-），男，研究員，主要從事淡水珍稀水生動(dòng)物繁殖生理研究.E-mail:liujy@ecsf.ac.cn

2016-11-26）