留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于线性混合模型的不同群体柔鱼胴长体质量关系的异质性研究

韩霈武 李楠 方舟 陈新军

引用本文:
Citation:

基于线性混合模型的不同群体柔鱼胴长体质量关系的异质性研究

    作者简介: 韩霈武 (1996—),男,硕士研究生,研究方向为大洋性头足类生物学。E-mail:2718344838@qq.com;
    通讯作者: 方舟, zfang@shou.edu.cn
  • 中图分类号: S 931.1

Heterogeneity of mantle length-body mass relationship in different Ommastrephes bartramii populations based on linear mixed model

    Corresponding author: Zhou FANG, zfang@shou.edu.cn ;
  • CLC number: S 931.1

  • 摘要: 为研究不同群体柔鱼 (Ommastrephes bartramii) 生长的异质性,根据2016—2018年在北太平洋采集的柔鱼样本,对该物种胴长和体质量的性别和月间差异进行了分析。构建了广义线性模型 (Generalized linear model, GLM) 和9个线性混合模型 (Linear mixed model, LMM),用于研究不同群体柔鱼的胴长和体质量关系 (W=aLb) 性别及时间上的差异。结果显示,东、西部群体不同性别柔鱼胴长、体质量均存在极显著性差异 (P<0.01);各月份间也均存在极显著性差异 (P<0.01)。东部群体b的固定值 (2.763 4) 小于3,表示其为负异速增长,肥满度随胴长增加而减少,胴体趋于细长;西部群体b的固定值 (3.117 19) 大于3,表示其为正异速增长,肥满度随胴长增加而增大,胴体趋于宽厚。研究表明,性别和月份对不同群体柔鱼胴长和体质量关系均有显著影响,线性混合模型能将性别和月份的异质性通过随机效应在单个模型中更准确、方便地体现出来,进一步证实了此模型在个体生长异质性研究中的优势。
  • 图 1  北太平洋柔鱼不同群体采样站点

    Figure 1.  Sampling station of different stocks O. bartramii in North Pacific Ocean

    图 2  东部群体柔鱼不同性别胴长体质量频率分布图

    Figure 2.  Frequency distribution of mantle length and body mass of O. bartramii of different sex for eastern population

    图 3  西部群体柔鱼不同性别胴长体质量频率分布图

    Figure 3.  Frequency distribution of mantle length and body mass of O. bartramii of different sex for western population

    图 4  柔鱼不同月份胴长体质量分布图

    Figure 4.  Distribution of mantle length and body mass of O. bartramii in different months

    图 5  不同群体柔鱼体质量-胴长关系

    Figure 5.  Mantle length-body mass relationship for different O. bartramii populations

    表 1  不同月份柔鱼样本基本信息

    Table 1.  Basic information of O. bartramii samples in different months

    群体
    Population
    年份
    Year
    月份
    Month
    样本数量
    Number of sample
    胴长范围
    Range of mantle length/mm
    东部群体 Eastern population2016—20175394163~500
    6456212~516
    西部群体 Western population2016—20187124187~290
    8283195~357
    9331199~386
    10 165190~377
    11 105215~421
    下载: 导出CSV

    表 2  不同群体柔鱼胴长-体质量关系模型的拟合效果指标

    Table 2.  Fitting effect index (AIC) of mantle length-body mass relationship models of different O. bartramii populations

    种群
    Population
    模型缩写
    Model abbreviation
    模型
    Model
    对数转化
    Log-transformed
    赤池信息量准则
    AIC
    均方根误差
    RMSE
    东部群体
    Eastern population
    LWRm = a×Lbln(m)=ln(a)+b×ln(L)−466.4020.183 265
    RSSm=a×L(b+RSS)ln(m)=ln(a)+(b+RSS)×ln(L)−472.3380.181 617
    RSIm=a×exp(RSILbln(m)=[ln(a)+RSI]+b×ln(L)−472.270.181 641
    RSI&Sm=a×exp(RSIL(b+RSS)ln(m)=[ln(a)+RSI]+(b+RSS)×ln(L)−469.0570.181 428
    RMSm=a×L(b+RMS)ln(m)=ln(a)+(b+RMS)×ln(L)−489.1070.179 784
    RMIm=a×exp(RMILbln(m)=[ln(a)+RMI]+b×ln(L)−492.7720.179367
    RMI&Sm=a×exp(RMIL(b+ RMS)ln(m)=[ln(a)+RMI]+(b+RMS)×ln(L)−595.1370.166 468
    RS&MSm=a×L(b+RSS+RMS)ln(m)=ln(a)+(b+RSS + RMS)×ln(L)−494.4520.178 558
    RS&MIm=a×exp(RSI)×exp(RMILbln(m)=[ln(a)+RSI + RMI]+b×ln(L)−497.6940.178 211
    RS&MI&Sm=a×exp(RSI)×exp(RMSL(b+RSI+ RMS)ln(m)=[ln(a) + RSI + RMS]+(b+RSI + RMI)×ln(L)−596.1380.165 394
    西部群体
    Western population
    LWRm = a×Lbln(m)=ln(a)+b×ln(L)−1 795.630.099 003
    RSSm=a×L(b+RSS)ln(m)=ln(a)+(b+RSS)×ln(L)−1 797.150.098 614
    RSIm=a×exp(RSILbln(m)=[ln(a)+RSI]+b×ln(L)−1 797.310.098 603
    RSI&Sm=a×exp(RSIL(b+RSS)ln(m)=[ln(a)+RSI]+(b+RSS)×ln(L)−1 796.720.098 512
    RMSm=a×L(b+RMS)ln(m)=ln(a)+(b+RMS)×ln(L)−1 807.040.097 626
    RMIm=a×exp(RMILbln(m)=[ln(a)+RMI]+b×ln(L)−1 807.410.097 615
    RMI&Sm=a×exp(RMIL(b+ RMS)ln(m)=[ln(a)+RMI]+(b+RMS)×ln(L)−1 807.590.097 215
    RS&MSm=a×L(b+RSS+RMS)ln(m)=ln(a)+(b+RSS + RMS)×ln(L)−1 808.140.097 347
    RS&MIm=a×exp(RSI)×exp(RMILbln(m)=[ln(a)+RSI + RMI]+b×ln(L)−1 808.560.097 333
    RS&MI&Sm=a×exp(RSI)×exp(RMSL(b+RSI+ RMS)ln(m)=[ln(a) + RSI + RMS]+(b+RSI + RMI)×ln(L)−1 808.090.096 976
    注:第一列为柔鱼的不同群体,第二列为第三、四列模型的缩写 RSS、RSI、RSI&S、RMS、RMI、RMI&S、RS&MS、RS&MI 和 RS&MI&S 分别为斜率 b 的性别随机效应、截距 ln(a) 的性别随机效应、截距 ln(a) 与斜率 b的性别随机效应、斜率 b 的月份随机效应、截距 ln(a) 的月份随机效应、截距 ln(a) 与斜率 b 的月份随机效应、斜率 b的性别与月份随机效应、截距 ln(a) 的性别与月份随机效应,以及截距 ln(a) 与斜率 b 的性别与月份随机效应Note: The first column shows different O. bartramii populations. The second column shows the abbreviations of models detailed in the third and fourth columns. RSS, RSI, RSI&S, RMS, RMI, RMI&S, RS&MS, RS&MI和RS&MI&S indicate random effects on slope(b) in sexual difference, random effects on intercept [ln(a)] in sexual difference, random effects on intercept [ln(a)] and slope(b) in sexual difference, random effects on slope(b) in months, random effects on intercept [ln(a)] in months, random effects on intercept [ln(a)] and slope(b) in months, random effects on slope(b) in sexual difference and months, random effects on intercept [ln(a)] fin sexual difference and months, and random effects on intercept [ln(a)] and slope(b) in sexual difference and months, respectively.
    下载: 导出CSV

    表 4  不同月份东西部群体柔鱼胴长与体质量关系

    Table 4.  Mantle length-body mass relationship for different O. bartramii populations in different months

    群体
    Population
    月份
    Month
    胴长与体质量关系
    Mantle length-mass relationship
    样本数量
    Number
    判定系数R2
    Coefficient of determination
    P
    东部群体
    Eastern population
    5m =1.9×10−3L2.295 73940.670 60.001
    6m =1×10−5L3.160 14560.973 80.001
    西部群体
    Western population
    7m = 9×10−6L3.205 31240.889 50.001
    8m = 2×10−5L3.090 92830.934 50.001
    9m = 2×10−5L3.059 73310.952 00.001
    10 m = 1×10−5L3.178 31650.972 80.001
    11 m = 5×10−5L2.902 81050.947 70.001
    下载: 导出CSV

    表 3  不同性别东西部群体柔鱼胴长与体质量关系

    Table 3.  Mantle length-body mass relationship of O. bartramii of different sex for different populations

    群体
    Population
    性别
    Sex
    胴长与体质量关系
    Mantle length-mass relationship
    样本数量
    Number
    判定系数R2
    Coefficient of determination
    P
    东部群体
    Eastern population
    雌性m=3×10−5L2.981 87440.892 90.001
    雄性m=2×10−4L2.621 01060.794 30.001
    西部群体
    Western population
    雌性m=1×10−5L3.159 95120.968 10.001
    雄性m =2×10−5L3.103 54960.941 50.001
    下载: 导出CSV
  • [1] 王尧耕, 陈新军. 世界大洋性经济柔鱼类资源及其渔业[M]. 北京: 海洋出版社, 2005: 124-137.
    [2] MURATA M. Oceanic resources of squids[J]. Mar Freshw Behav Phy, 1990, 18(1): 19-71. doi:  10.1080/10236249009378779
    [3] BOWER J R, ICHⅡ T. The red flying squid (Ommastrephes bartramii): a review of recent research and the fishery in Japan[J]. Fish Res, 2005, 76(1): 39-55. doi:  10.1016/j.fishres.2005.05.009
    [4] 陈新军. 渔业资源与渔场学[M]. 北京: 海洋出版社, 2004: 29-212.
    [5] YATSU A, TANAKA H, MORI J. Population structure of the neon flying squid, Ommastrephes bartramii, in the North Pacific Ocean[J]. Jpn Mar Fish Resour Res Cent, 1998: 31-48.
    [6] FANG Z, LIU B, LI J, et al. Stock identification of neon flying squid (Ommastrephes bartramii) in the North Pacific Ocean on the basis of beak and statolith morphology[J]. Sci Mar, 2014, 78(2): 239-248. doi:  10.3989/scimar.03991.06A
    [7] 方舟, 陈新军, 瞿俊跃, 等. 北太平洋柔鱼角质颚形态及生长年间差异[J]. 上海海洋大学学报, 2020, 29(1): 109-120. doi:  10.12024/jsou.20181202485
    [8] WATANABE H, KUBODERA T, ICHⅡ T, et al. Feeding habits of neon flying squid Ommastrephes bartramii in the transitional region of the central North Pacific[J]. Mar Ecol Prog Ser, 2004, 266(1): 173-184.
    [9] CHEN C S, CHIU T S. Variations of life history parameters in two geographical groups of the neon flying squid, Ommastrephes bartramii, from the North Pacific[J]. Fish Res, 2003, 63(3): 349-366. doi:  10.1016/S0165-7836(03)00101-2
    [10] YATSU A. Distribution of flying squid, Ommastrephes bartramii, in the north Pacific based on a jigging survey, 1976−1983[J]. Bull Natl Res Inst Far Seas Fish, 1992(29): 13-37.
    [11] YATSU A, MORI J. Early growth of the autumn cohort of neon flying squid, Ommastrephes bartramii, in the North Pacific Ocean[J]. Fish Res, 2000, 45(2): 189-194. doi:  10.1016/S0165-7836(99)00112-5
    [12] SUN M C, ZHANG S H, QIAN W G. Approach on biological characteristics of Ommastrephes bartramii on the central and eastern squid fishing grounds in the north pacific[J]. Mar Tieher, 2001, 23(1): 21-24.
    [13] YATSU A, MIDORIKAWA S, SHIMADA T, et al. Age and growth of the neon flying squid, Ommastrephes bartramii, in the North Pacific Ocean[J]. Fish Res, 1997, 29(3): 257-270. doi:  10.1016/S0165-7836(96)00541-3
    [14] MCINTYRE R B A D. Cephalopod Life Cycles. Vol. 1. Species Accountsby P. P. Boyle[J]. J Anim Ecol, 1986, 55(1): 395-396.
    [15] GRANADOS-AMORES J, SALINAS-ZAVALA C A, FLORES-ORTEGA J R, et al. Length-weight relationship and condition factor for 7 loliginid squid species in Mexican waters[J]. Cienc Mar, 2019, 45(4): 175-180. doi:  10.7773/cm.v45i4.3015
    [16] VASCONCELOS P, PEREIRA F, CARVALHO A N, et al. Weight-length relationships and relative growth of the cuttlefish (Sepia officinalis): causes and effects of hypoallometry[J]. Thalassas: Int J Mar Sci, 2018, 32(4): 323-331.
    [17] 衷思剑, 麻秋云, 刘淑德, 等. 基于线性混合效应模型的黄鮟鱇体长体重关系的时空差异[J]. 中国水产科学, 2018, 25(6): 152-160.
    [18] XU H, SUN Y J, WANG X J, et al. Linear mixed-effects models to describe individual tree crown width for China-fir in Fujian province, southeast China[J]. PLoS One, 2015, 10(4): e0122257. doi:  10.1371/journal.pone.0122257
    [19] MA Q Y, JIAO Y, REN Y P. Linear mixed-effects models to describe length-weight relationships for yellow croaker (Larimichthy spolyactis) along the north coast of China[J]. PLoS One, 2017, 12(2): e0171811. doi:  10.1371/journal.pone.0171811
    [20] 瞿俊跃, 方舟, 陈新军. 基于线性混合模型的蓝点马鲛叉长体重生长关系的月间及性别差异[J]. 中国水产科学, 2020: 1-10.
    [21] 廖为耕. 浅谈北太平洋之海洋环境及赤鱿之洄游[J]. 中国水产(台刊), 1980, 457: 39-47.
    [22] LIPINSKI M R, UNDERHILL L G. Sexual maturation in squid: quantum or continuum[J]. S Afr J Mar Sci, 1995, 15: 207-223. doi:  10.2989/02577619509504844
    [23] 李春喜, 邵云, 姜丽娜. 生物统计学[M]. 北京: 科学出版社, 2008: 85-110.
    [24] GUEORGUIEVA R, KRYSTAL J H. Move over anova: progress in analyzing repeated-measures data and its reflection in papers published in the archives of general psychiatry[J]. Arch Gen Psychiat, 2004, 61(3): 310-317. doi:  10.1001/archpsyc.61.3.310
    [25] 王斌会. 多元统计分析及R语言建模[M]. 广州: 暨南大学出版社, 2010, 157-158.
    [26] KEYS A B. The weight-length relationship in fishes[J]. P Natl A Sci USA, 1928, 14(12): 922-925. doi:  10.1073/pnas.14.12.922
    [27] BURNHAM K P, ANDERSON D R. Model selection and multi-model inference: a practical information theoretic approach[M]. 2nd ed. New York: Springer Press, 2002: 484-488.
    [28] BERNAL-RUSIEL J L, GREVE D N, REUTER M, et al. Statistical analysis of longitudinal neuroimage data with Linear Mixed Effects models[J]. NeuroImage, 2013, 66(2): 249-260.
    [29] 黄洪亮, 郑元甲, 程家骅. 北太平洋海区柔鱼生物学特征研究[J]. 海洋渔业, 2003, 25(3): 21-24.
    [30] ICHⅡ T, MAHAPATRA K, OKAMURA H, et al. Stock assessment of the autumn cohort of neon flying squid (Ommastrephes bartramii) in the North Pacific based on past large-scale high seas driftnet fishery data[J]. Fish Res, 2006, 78(2): 286-297.
    [31] KATUGIN O N. Patterns of genetic variability and population structure in the north pacific squids Ommastrephes bartramii, Todarodes pacificus, and berryteuthis magister[J]. Bull Mar Sci, 2002, 71(1): 383-420.
    [32] HILE R. Age and growth of the cisco, Leucichthys artedi (Le Sueur), in the lakes of the northeastern highlands, Wisconsin[J]. Bull Bur Fish, 1936, 48(11): 211-317.
    [33] KEITH A, BIGELOW, KEVIN C, et al. Hatch dates and growth of Ommastrephes bartramii paralarvae from Hawaiian waters as determined from statolith analysis[J]. Rec Adv Fish Biol, 1993: 15-24.
    [34] 唐峰华, 伍玉梅, 樊伟. 2010年北太平洋海区巴特柔鱼(Ommastrephe bartrami)的生物学特征[J]. 海洋科学进展, 2012, 30(2): 257-265. doi:  10.3969/j.issn.1671-6647.2012.02.013
    [35] FORSYTHE J W. Accounting for the effect of temperature on squid growth in nature: from hypothesis to practice[J]. Mar Freshw Res, 2004, 55(4): 331-339. doi:  10.1071/MF03146
    [36] ICHⅡ T, MAHAPATRA K, SAKAI M, et al. Life history of the neon flying squid: effect of the oceanographic regime in the North Pacific Ocean[J]. Mar Ecol Prog Ser, 2009, 378: 1-11. doi:  10.3354/meps07873
    [37] 章寒, 郑基, 虞聪达, 等. 东太平洋赤道海域茎柔鱼主要生物学特性比较研究[J]. 中国水产科学, 2019, 26(4): 745-755.
    [38] DOMINGUES P M, SYKES A, ANDRADE J P. The effects of temperature in the life cycle of two consecutive generations of the cuttlefish Sepia officinalis (Linnaeus, 1758), cultured in the Algarve (South Portugal)[J]. Aquacult Int, 2002, 10(3): 207-220. doi:  10.1023/A:1022148802078
    [39] DOMINGUES P M, POIRIER R, DICKEL L, et al. Effects of culture density and live prey on growth and survival of juvenile cuttlefish, Sepia officinalis[J]. Aquacult Int, 2003, 11(3): 225-242. doi:  10.1023/A:1024803802486
    [40] 陈锦辉, 王学昉, 田思泉, 等. 基于混合效应模型的印度洋大眼金枪鱼生长特征异质性分析[J]. 中国水产科学, 2020, 27(5): 570-578.
  • [1] 崔雪森唐峰华周为峰吴祖立杨胜龙化成君 . 基于支持向量机的西北太平洋柔鱼渔场预报模型构建. 南方水产科学, doi: 10.3969/j.issn.2095-0780.2016.05.001
    [2] 党莹超戴小杰吴峰 . 北太平洋金枪鱼延绳钓钓钩垂直分布及浸泡时间对渔获物的影响. 南方水产科学, doi: 10.12131/20190252
    [3] 杨铭霞陈新军冯永玖官文江 . 利用多种空间插值法分析西北太平洋柔鱼渔场. 南方水产科学, doi: 10.3969/j.issn.2095-0780.2015.02.011
    [4] 李佳佳汪金涛陈新军雷林关长涛 . 不同气候模态下西北太平洋柔鱼冬春生群资源时空分布变化研究. 南方水产科学, doi: 10.12131/20190175
    [5] 王旭徐衡邹莉杜肖王凤丽郑伟韩志强水柏年 . 3个单齿螺地理群体形态差异分析. 南方水产科学,
    [6] 孙立元郭华阳朱彩艳马振华江世贵张殿昌 . 卵形鲳鲹育种群体遗传多样性分析. 南方水产科学, doi: doi:10.3969/j.issn.2095-0780.2014.02.010
    [7] 韩真肖永双高天翔 . 中国近海9个小黄鱼群体的形态学比较研究. 南方水产科学, doi: 10.3969/j.issn.2095-0780.2012.03.004
    [8] 董在杰刘念傅建军朱文彬王兰梅苏胜彦 . 6个野生与选育鲤群体的微卫星遗传分析. 南方水产科学, doi: 10.3969/j.issn.2095-0780.2018.04.006
    [9] 侯娟周为峰樊伟 . 基于集成学习的南太平洋长鳍金枪鱼渔场预报模型研究. 南方水产科学, doi: 10.12131/20200022
    [10] 薛嘉伦樊伟唐峰华郭刚刚汤伟张胜茂 . 基于最大熵模型预测西北太平洋公海鲐潜在栖息地分布. 南方水产科学, doi: 10.3969/j.issn.2095-0780.2018.01.012
    [11] 张嘉容杨晓明戴小杰邹莉瑾 . 南太平洋长鳍金枪鱼延绳钓渔获率与环境因子的关系研究. 南方水产科学, doi: 10.12131/20190178
    [12] 党莹超陈屹林戴小杰李云凯吴峰刘攀 . 北太平洋亚热带海域大眼金枪鱼秋季摄食习性的初步研究. 南方水产科学, doi: 10.12131/20190184
    [13] 郭刚刚张胜茂樊伟陈新军杨胜龙 . 南太平洋长鳍金枪鱼垂直活动水层空间分析. 南方水产科学, doi: 10.3969/j.issn.2095-0780.2016.05.016
    [14] 徐博张衡唐峰华隋芯张瑛瑛侯刚 . 基于GAM的北太平洋公海围网主要渔获种类渔场重心与环境因子的关系分析. 南方水产科学, doi: 10.12131/20200006
    [15] 彭露江艳娥徐姗楠陈作志 . 南海北部陆架区金线鱼矢耳石形态及其与体长、体质量关系. 南方水产科学, doi: 10.12131/20180049
    [16] 李鹏许柳雄周成王学昉唐浩刘伟 . 中西太平洋金枪鱼围网自由鱼群渔场重心变动及其与南方涛动指数的关系. 南方水产科学, doi: 10.12131/20190158
    [17] 翟天晨戴小杰朱江峰陈丽雯 . 南太平洋雌性大眼金枪鱼性腺成熟等级研究. 南方水产科学, doi: 10.3969/j.issn.2016.01.014
    [18] 刘攀戴小杰王杰高春霞吴峰王学昉 . 热带西太平洋帆蜥鱼胃含物组成分析. 南方水产科学, doi: 10.12131/20180141
    [19] 陈春光张敏邹晓荣陆奇巍许啸梁严威 . 东南太平洋智利竹筴鱼中心渔场的月间变动研究. 南方水产科学, doi: 10.3969/j.issn.2095-0780.2014.05.009
    [20] 郑晓春戴小杰朱江峰陈彦高春霞吴峰张艳波 . 太平洋中东部海域大眼金枪鱼胃含物分析. 南方水产科学, doi: 10.3969/j.issn.2095-0780.2015.01.011
  • 加载中
计量
  • 文章访问数:  135
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-05-29
  • 录用日期:  2020-06-20

基于线性混合模型的不同群体柔鱼胴长体质量关系的异质性研究

    作者简介:韩霈武 (1996—),男,硕士研究生,研究方向为大洋性头足类生物学。E-mail:2718344838@qq.com
    通讯作者: 方舟, zfang@shou.edu.cn
  • 1. 上海海洋大学海洋科学学院,上海 201306
  • 2. 大洋渔业资源可持续开发教育部重点实验室/国家远洋渔业工程技术研究中心/农业农村部大洋渔业开发重点实验室/农业农村部大洋渔业资源环境科学观测实验站,上海 201306

摘要: 为研究不同群体柔鱼 (Ommastrephes bartramii) 生长的异质性,根据2016—2018年在北太平洋采集的柔鱼样本,对该物种胴长和体质量的性别和月间差异进行了分析。构建了广义线性模型 (Generalized linear model, GLM) 和9个线性混合模型 (Linear mixed model, LMM),用于研究不同群体柔鱼的胴长和体质量关系 (W=aLb) 性别及时间上的差异。结果显示,东、西部群体不同性别柔鱼胴长、体质量均存在极显著性差异 (P<0.01);各月份间也均存在极显著性差异 (P<0.01)。东部群体b的固定值 (2.763 4) 小于3,表示其为负异速增长,肥满度随胴长增加而减少,胴体趋于细长;西部群体b的固定值 (3.117 19) 大于3,表示其为正异速增长,肥满度随胴长增加而增大,胴体趋于宽厚。研究表明,性别和月份对不同群体柔鱼胴长和体质量关系均有显著影响,线性混合模型能将性别和月份的异质性通过随机效应在单个模型中更准确、方便地体现出来,进一步证实了此模型在个体生长异质性研究中的优势。

English Abstract

参考文献 (40)

返回顶部

目录

    /

    返回文章
    返回