Tasas de mortalidad

Dinámica de un cohorte

Glosario de términos
Z Tasa instantanea de mortalidad total, coeficiente de mortalidad total o tasa de mortalidad total
M Tasa de mortalidad natural
F Tasa de mortalidad por pesca
q Coeficiente de capturabilidad (medida de la capacidad del arte para capturar peces)
f esfuerzo de pesca
Tr Edad de reclutamiento
Ntr Número de reclutas "R"

Z=F+M

Modelo de extinción exponencial 1

ΔNt/Δt = −Z∗Nt Nt= a *exp -Zt

Modelo de extinción exponencial 2

Nt = Ntr∗exp(−Z∗(t−Tr)) Nt2 = Nt1∗exp(−Z*(t2-t1)) Nt = N0∗exp(−Z)

Z = 1/(t2-t1) * ln (t1/t2)

Ecuación de Baranov o de captura

Ct1,t2 = F/Z∗(Nt1−Nt2)
mortalidad total
Numero de sobrevivientes (Nt) y edad (t) de una cohorte (ejm. estanque)
Curva de extinción exponencial

Nt= N0 *exp -Zt Resumen

si el tiemp si el tiempo çfuese dias se divide entre 365
t (años) Nt Ln Nt Estadísticas de la regresión
0 1000 6.91 Coeficiente 0.999627
1 607 6.41 Coeficient 0.999255
2 368 5.91 R^2 ajust 0.999187
3 223 5.41 Error típic 0.056342
4 135 4.91 Observacio 13
5 82 4.41
6 50 3.91 ANÁLISIS DE VARIANZA
7 30 3.40 Grados de Suma
libertad
de
Promedio
cuadrados
de los cuadrados
F
8 18 2.89 Regresión 1 46.82011 46.82011 14749.41
9 11 2.40 Residuos 11 0.034918 0.003174
10 7 1.95 Total 12 46.85503
11 4 1.39
12 2 0.69 CoeficientesError típicoEstadísticoProbabilidad
t
Intercepció 6.933237 0.029531 234.7783 1.05E-21
Variable X -0.507201 0.004176 -121.4471 1.476E-18

Curva de extinción exponencial

1200 Curva de extinción exponencia
8.00
1000
f(x) = 1025.8096546767 exp( − 0.5072014796668 x ) 7.00
R² = 0.999254763231712 f(x) = − 0.5072014796668 x + 6.93323748676
800 6.00 R² = 0.999254763231712
600 5.00
Nt

Ln Nt

4.00
400
3.00
200 2.00
0 1.00
0 2 4 6 8 10 12 14 0.00
0 2 4 6 8 1
t
t
Entonces: Z = 0.507

Ejemplo. Creación de CCE

N0 1000
Z 1

t Nt Ln Nt %
0 1000 6.91 100
1 368 5.91 36.78794
2 135 4.91 13.53353
3 50 3.91
4 18 2.91
57 1.91
 62 0.91
71 -0.09

8.00

7.00
f(x) = − x + 6.90775527898214
6.00 R² = 1

5.00

4.00
ln Nt

3.00

2.00

1.00

0.00
0 1 2 3 4 5 6 7
t (edad)
stanque)

Valor crítico de F
1.476E-18

Inferior 95%
Superior 95%
Inferior 95.0%
Superior 95.0%
6.86824 6.998235 6.86824 6.998235
-0.516393 -0.498009 -0.516393 -0.498009

e extinción exponencial linealizada

5072014796668 x + 6.93323748676183
9254763231712

4 6 8 10 12 14
t
Uso de estructura de edades del stock (ejem. Crucero de evaluación)

Datos de entrada:
N° o biomasa del stock
Muestra de frecuencia de tallas
Relación peso-talla
Clave edad-talla

Ejemplo en clase
Biomasa: 1000 Kg
Relación peso talla: W = 0.01*L^3

Muestra de frecuencia de tallas:

L (cm) f
10 5
11 10
12 25
13 50
14 90
15 60
16 30
17 15
18 10
19 5
Total 300

Clave edad talla:

talla/edad 1 2 3 4 5 Total
10 3 3
11 25 25
12 15 3 18
13 5 12 17
14 1 26 2 29
15 15 10 25
16 10 22 4 36
17 1 14 24 10 49
18 5 1 20 26
19 5 5
233
Solución:

1. Calculamos la captura por cada talla, necesitamos la muestra de frecuencia de tallas y la relación talla peso
Relación peso talla: a 0.01
 b 3

Peso de c/talla (g) (g)

L (cm) f wLi (g) peso total proporció biomasa N p/talla
p/talla n p/talla
10 5 10 50 0.0054944 5494.4 549
11 10 13.31 133.1 0.0146262 14626.2 1099
12 25 17.28 432 0.047472 47472.0 2747
13 50 21.97 1098.5 0.120713 120713.0 5494
14 90 27.44 2469.6 0.2713816 271381.6 9890
15 60 33.75 2025 0.222525 222525.0 6593
16 30 40.96 1228.8 0.1350315 135031.5 3297
17 15 49.13 736.95 0.0809826 80982.6 1648
18 10 58.32 583.2 0.0640872 64087.2 1099
19 5 68.59 342.95 0.0376864 37686.4 549
Total 300 340.75 9100.1 1 1000000.0 32967

2. Calculamos la captura por edad

Datos de entrasa: Clave edad talla y numero de individuos por talla
L (cm) 1 2 3 4 5 Total N
10 3 3 549
11 25 25 1098
12 15 3 18 2747
13 5 12 17 5494
14 1 26 2 29 9890
15 15 10 25 6593
16 10 22 4 36 3296
17 1 14 24 10 49 1648
18 5 1 20 26 1098
19 5 5 549
Total 49 67 53 29 35 233 32962

L (cm) 1 2 3 4 5 Total
10 549 0 0 0 0 549
11 1098 0 0 0 0 1098
12 2289 458 0 0 0 2747
13 1616 3878 0 0 0 5494
14 341 8867 682 0 0 9890
15 0 3956 2637 0 0 6593
16 0 916 2014 366 0 3296
17 0 34 471 807 336 1648
18 0 0 211 42 845 1098
Curva
19 0 0 0 0 549 549
12

10
f(x)
8 R² =

6
t
 Curva
12
Total 5893 18108 6016 1216 1730 32962
10
f(x)
Tabla lista para ser usada para el calculo de Z!!! Resumen 8 R² =
Edad Ct ln Ct 6

Ln Ct
1 5893 8.681534654 Estadísticas de la regresión
4
2 18108 9.804100037 Coeficiente de co 0.994403111
2
3 6016 8.702095112 Coeficiente de d 0.988837547
4 1216 7.103023223 R^2 ajustado 0.977675095 0
1.5 2
5 1730 7.455843112 Error típico 0.202926739
Observaciones 3

ANÁLISIS DE VARIANZA
Grados de libertad
Suma de cuadrados
Regresión 1 3.64790797984
Residuos 1 0.041179261225
Total 2 3.689087241065

Coeficientes Error típico

Intercepción 12.58802135 0.446131254138
Variable X 1 -1.350538407 0.143490872923

Entonces: Z: 1.348

RECORDEMOS
¿Como trabajaríamos estos datos para hallar los parámetros de crecimiento de un pez subtropical?

Necesitamos saber la talla media de cada edad

Datos de entrada:
Clave edad-talla
Operación: Calcular el promedio de talla de la edad

En el ejemplo anterior:
L (cm) 1 2 3 4 5
10 3
11 25
12 15 3
13 5 12
14 1 26 2
15 15 10
16 10 22 4
17 1 14 24 10
18 5 1 20 VB
19 5
20.00
Total 49 67 53 29 35 18.00
Talla media 11.51 14.30 16.19 16.90 17.86 16.00
14.00
12.00
Lo que resulta en:
10.00

talla
edad (años) talla media
8.00
(cm)
6.00
1 11.51 4.00
2.00
2 14.30
0.00
3 16.19 0.5 1 1.5 2 2.5 3 3
4 16.90 edad (años)
5 17.86

Usando FW, dt = 1 año: FW

Lt Lt+1
20.00
11.51 14.30 18.00
14.30 16.19 16.00 f(x) = 0.617935756060357
R² = 0.97587065323742
16.19 16.90 14.00
16.90 17.86 12.00
10.00

Lt+1
8.00
6.00
a 7.21 4.00
b 0.62 2.00
0.00
11.0012.0013.0014.0015.0016.00
K 0.48
Lt
Linf 18.88

Usando VB, Linf= 18.88

edad (años) (-ln (1- VB
Lt/Linf)) 3.50
1 0.94
3.00
2 1.42 f(x) = 0.479359064325185
2.50
3 1.95 R² = 0.990745784102252
(-ln (1-Lt/Linf))

4 2.25 2.00

5 2.92 1.50

1.00
a 0.46 0.50
b 0.48
0.00
0.5 1 1.5 2 2.5 3 3.5 4 4.5
K 0.48 edad
t0 -0.96
 0.00
0.5 1 1.5 2 2.5 3 3.5 4 4.5
edad

Entonces:
Linf 18.88
K 0.48
t0 -0.96

¿A que edad hay mayor biomasa?

ación talla peso
 Curva linealizada de captura
12

10
f(x) = − 1.3505384074213 x + 12.5880213463775
8 R² = 0.988837547465365

6
t
 Curva linealizada de captura
12

10
f(x) = − 1.3505384074213 x + 12.5880213463775
8 R² = 0.988837547465365

6
Ln Ct

4 Curva de captura
2 20000
18000
0 16000
1.5 2 2.5 3 3.5 4 4.5 14000
t (edad) 12000
10000

N
8000
6000
4000
Promedio de los cuadradosF Valor crítico de F 2000
0
3.64790798 88.58604723 0.067386241
0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5
0.041179261
edad

Estadístico t Probabilidad Inferior 95% Superior 95%Inferior 95.0%Superior 95.0%

28.21595938 0.022552961 6.919386292 18.2566564 6.919386292 18.2566564
-9.41201611 0.067386241 -3.17376282 0.472686002 -3.17376282 0.472686002

ubtropical?
 VB

Column B

1 1.5 2 2.5 3 3.5 4 4.5 5 5.5

edad (años)

FW

f(x) = 0.617935756060357 x + 7.21205209866178

R² = 0.97587065323742

Lt+1
Linear (Lt+1)

0012.0013.0014.0015.0016.0017.0018.00
Lt

VB

f(x) = 0.479359064325185 x + 0.457994009629751

R² = 0.990745784102252

Column B
Linear (Column B)

5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5

edad
5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5
edad
captura

3 3.5 4 4.5 5 5.5

edad
Curva de captura linealizada con intervalos de tiempo constante (ejm. Desembarqu
Usando estructura de edades de la captura
C(t,t+∆t) = g - Z*t

La captura por edades se calcula por cada año

Edad 1974 1975 1976 1977 1978 1979 1980 promedio Observaciones
0 599 239 424 664 685 478 330 488.43 no se considera
1 679 860 431 1004 418 607 288 612.43 no se considera
2 1097 390 1071 532 335 464 323 601.71 no se considera
3 275 298 159 269 203 211 243 236.86 si
4 40 54 75 32 69 86 80 62.29 si
5 6 9 13 18 8 25 31 15.71 si
6 1 8 3 5 5 3 8 4.71 si
7 6 1 1 1 1 2.00 (no se considera)

C(a,t,t+1) numero de peces capturados en el año (a), de edad entre t y t+1 (en millones). Ejemplo C(1976,3,4) = 159

1974 1975 1976 1977 1978 1979 1980 promedio

Z 1.875 1.264 1.366 1.253 1.327 1.400 1.119 1.372

Edad 1974 1975 1976 1977 1978 1979 1980 promedio

0 6.395262 5.476464 6.049733 6.498282 6.529419 6.169611 5.799093 6.191193
1 6.520621 6.756932 6.066108 6.911747 6.035481 6.408529 5.66296 6.417432
2 7.000334 5.966147 6.976348 6.276643 5.814131 6.139885 5.777652 6.399783
3 5.616771 5.697093 5.068904 5.594711 5.313206 5.351858 5.493061 5.467457
4 3.688879 3.988984 4.317488 3.465736 4.234107 4.454347 4.382027 4.131732
5 1.791759 2.197225 2.564949 2.890372 2.079442 3.218876 3.433987 2.75457
6 0 2.079442 1.098612 1.609438 1.609438 1.098612 2.079442 1.550597
7 0 0 0.693147

Curva de captura linealizada

8
7
6
5 f(x) = − 1.87474332778856 x + 11.2106974803007
4
Ln C

3
2
1
0
0 1 2 3 4 5 6 7 8
edad
 2
1
0
0 1 2 3 4 5 6 7 8
edad

b= -Z -1.8747
Z 1.8747
te (ejm. Desembarque en un añ)

599
679
1097
275
40
6
1
6
1974 Linear (1974)
1975 Linear (1975)
6

5 f(x)
f(x) == −− 1.87474332778856
1.26447153037048 xx ++ 11.2106974803007
9.18080779968857
R² = 0.999720610732292
R² = 0.910405399683352
4
Ln (Ct)

0
2.5 3 3.5 4 4.5 5 5.5 6 6.5
t (edad en años)
La curva de captura linealizada basada en datos de talla

ln(C(L1,L2)/Δt) = c−Z∗t(L1+L2)/2

Para hallar Z se necesitaFrecuencia de tallas de las capturas y los parámetros de crecimiento

Ecuación inversa de VBt= t0- ln(1-L/Linf)/k

Ejemplo: Composición de tallas para el salmonete dato

Linf 23.1
K 0.59

L1-L2 L1 L2 C(L1-L2) t(L1) t(L2) ∆t

6-7 6 7 3 0.510 0.612 0.102

7-8 7 8 143 0.612 0.721 0.109
8-9 8 9 271 0.721 0.837 0.116
9-10 9 10 318 0.837 0.961 0.125
10-11 10 11 416 0.961 1.096 0.135
11 12 488 1.096 1.242 0.146
12 13 614 1.242 1.402 0.160
13 14 613 1.402 1.579 0.177
14 15 493 1.579 1.776 0.197
15 16 278 1.776 2.000 0.223
16 17 93 2.000 2.257 0.257
17 18 73 2.257 2.560 0.303
18 19 7 2.560 2.930 0.370
19 20 2 2.930 3.404 0.474
20 21 2 3.404 4.064 0.660
21 22 0 4.064 5.160 1.096
22 23 1 5.160 9.224 4.064
23 23.09 1 9.224 13.127 3.903

Por lo que tenemos:

t[(L1+L2)/2] Ln C((L1,L2)/∆t) Curva de captura linealizada basada en los datos de
=X =Y composición por tallas
0.56 3.38 10.00
0.67 7.18
0.78 7.76 8.00 f(x) = − 4.16163496283054 x + 14.7689220679855
0.90 7.84 R² = 0.975201722672616
ln C(L1,L2)/Δt

6.00
1.03 8.04
1.17 8.11 4.00
1.32 8.25
2.00 Resumen
1.49 8.15
1.68 7.82 Parte utilizada
0.00de la curva
1.89 7.13 0.00 2.00 Estadísticas
4.00 de la regresión8.00
6.00 10.00 12.00
2.13 5.89 -2.00 Coeficiente 0.988922
2.41 5.48 Coeficient 0.977968
Edad relativa t[(L1+L2)/2]
2.75 2.94 R^2 ajust 0.97246
3.17 1.44 Error típic 0.408891
3.73 1.11 Observacio 6
4.61 Err:502
7.19 -1.40 ANÁLISIS DE VARIANZA
11.18 -1.36 Grados de Suma
libertad
de
Promedio
cuadrados
de los cuadrados
F
Regresión 1 29.68519 29.68519 177.5517
Residuos 4 0.668767 0.167192
Total 5 30.35396

CoeficientesError típicoEstadísticoProbabilidad
t
Intercepció 15.36853 0.787236 19.52214 4.06E-05
Z 4.2392 1.490555 -4.388788 0.329369 -13.32485 0.000183

Resultados de datos de probabilidad

Percentil 8.15159
8.333333 1.439964
25 2.940375
41.66667 5.482927
58.33333 5.890119
75 7.126691
91.66667 7.82351
 Gráfico de probabilidad normal
8.15159000311126

10

0
0 10 20 30 40 50 60 70 80 90 100
Muestra percentil

Valor crítico de F
0.000183

Inferior 95%
Superior 95%
Inferior 95.0%
Superior 95.0%
13.18282 17.55425 13.18282 17.55425
-5.303262 -3.474314 -5.303262 -3.474314
Gráfico de Z respecto al esfuerzo para la estimación separada de F y M

Lcm talla media de la captura

Lc talla de retención del 50%

Datos para la estimación de M y q en Selanoides leptolepis del golfo de Tailandia.

K 1.16
Lc 10
Linf 20

Z obtenido por la ecuación de Beverton y Holt (1956), basada en tallas

Za= K* (Linf-Lcm)/(Lcm-Lc)

Método de Paloheimo (1961)

Za=M + q*fa Regresión ec 4.6.1

En millones de h de arrastre
año esfuerzo (ftalla media Za= K* (Linf-Lcm)/(Lcm-Lc)
1966 2.08 13.25 2.41
1967 2.08 13.01 2.69
1968 3.50 12.99 2.72
1969 3.60 13.07 2.62
1970 3.80 12.37 3.73
1971
1972 7.19 12.30 3.88
1973 9.94 12.01 4.61
1974 6.06 12.60 3.30

Mortalidad total vs
5.00
4.50
f(x) = 0.25335063201037
4.00 R² = 0.806891491289058
3.50
3.00
2.50
Z

2.00
1.50
1.00
0.50
0.00
1.00 2.00 3.00 4.00 5.00 6.00
esfuerzo (millones de ho
Entonce tenemos:
M 2.04
q 0.25

Con estos datos calculamos el F

año esfuerzo (ftalla media Za= K* (LinM F
1966 2.08 13.25 2.41 2.0351 0.37
1967 2.08 13.01 2.69 2.0351 0.66
1968 3.50 12.99 2.72 2.0351 0.68
1969 3.60 13.07 2.62 2.0351 0.58
1970 3.80 12.37 3.73 2.0351 1.70
1971 2.0351
1972 7.19 12.30 3.88 2.0351 1.85
1973 9.94 12.01 4.61 2.0351 2.58
1974 6.06 12.60 3.30 2.0351 1.27
de F y M

Mortalidad total vs esfuerzo

0
0
f(x) = 0.25335063201037 x + 2.03514571908251
0 R² = 0.806891491289058
0
0
0
0
0
0
0
0
1.00 2.00 3.00 4.00 5.00 6.00 7.00 8.00 9.00 10.00 11.00
esfuerzo (millones de horas de arrastre)