学以致用——使用莱斯利矩阵模型预测蠵龟种群数量的演变(Demographics of the Loggerhead Sea Turtle using Leslie population model)

学了线性代数,感觉像是打开了一个通往新世界的大门,非常神奇。

本文讲述使用莱斯利矩阵模型预测蠵龟种群数量的演变。

首先,蠵龟(Loggerhead Sea Turtle)的种群基本信息如下表所示:

Table 1   Four-Stage Model for Loggerhead Sea Turtle Demographics  
Stage Number Description  (age in years)                        Annual  survivorshipEggs laid per year
1 Eggs, hatchlings (<1) 0.670
2 Juveniles and subadults (1–21)0.740
3Novice breeders (22)0.81127
4Mature breeders (23–54)0.8179

然后,使用以下公式可计算各阶段的留存率(proportion remaining in stage i),

 

再使用以下公式计算各阶段生存下来并成功进入下一阶段的概率,

 

这样,可以计算出Leslie模型中用到一些参数的值如下:

piqidi
00.671
0.739368080.0006319220
00.811
0.8097230370.00027696331

 

各参数的说明如下:

Note:    
1. di represents the duration of the ith stage  
2. si is the annual survivorship rate for that stage  
3. pi is the proportion remaining in stage i the following year  
4. e2, e3, e4 denote the average number of eggs laid by a member of stage i (i = 2, 3, 4) in 1 year
5. qi proportion of the population that will survive and move into stage i + 1 

 莱斯利模型如下,

Leslie Matrix Model:   
p1e2e3e4
q1p200
0q2p30
00q3p4

 将上述已知及计算出的参数填入上述莱斯利模型,得到如下莱斯利矩阵:

Leslie Matrix:   
0012779
0.670.739400
00.00063200
000.810.8077

然后,使用矩阵乘法及VBA宏代码,计算出以下结果:

Stage Number1234
year0 (Initial Population)2000003000005001500
year11820003558201901617
year21518733850332251460
year31439153864482431361
year41383803821632441296
year51333723752862421244
year61290103668462371201
year71249783576832321162
year81212623482062261126
year91176563387092201093
year101142873292712141061
year111109973200352081030
year121077863110022021000
year13104654302171197971
year14101728293543191944
year1598833285203186917
year1696065277097180891
year1793249269249175865
year1890560261560170840
year1987950254073165816
year2085419246788161793
year2183094239706156771
year2280721232912151749
year2378348226298147727
year2476102219818143706
year2573935213522139686
year2671847207415135667
year2769838201500131648
year2867829195781127629
year2965820190206124611
year3064017184738120594
year3162166179487117577
year3260442174364113561
year3358670169421110545
year3457025164579107529
year3555380159896104514
year3653814155332101499
year375224815090898485
year385076114658895471
year394927414239793457
year404791413830290444
year414650613436387432
year424517713050785419
year434389612676582407
year444256712314080395
year454136511957078384
year464024211612576373
year473911911282573363
year483794810963371352
year493682510648869342
year503578110341067332

 这和教材中的结果很接近(具体差异应该是参数的细微差异导致的)。

结论:

50年以后,这种动物所剩数量将大幅减少,看来这种动物属于濒危动物了。 

 

参考:

教材中的完整案例如下,

 

 

评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值