Commit 1d12514c authored by Iago Mosqueira's avatar Iago Mosqueira
Browse files

WKFlatNSCS 2020

parent f88d2f8a
---
title:
author: Iago MOSQUEIRA (WMR) <iago.mosqueira@wur.nl>
date: 2019-10-07 10:33
documentclass: article
output:
pdf_document:
fig_caption: yes
number_sections: yes
toc: yes
tags:
abstract:
license: Creative Commons Attribution-ShareAlike 4.0 International Public License
---
# TASKS
- REWRITE data.R
- LOAD data from Datras + InterCatch
- CREATE AAP package
# STRUCTURE
- data/
- data.R
- input/
- input.R
- model/
- model.R
- output/
- output.R
- report/
- report.R
- report/report.Rmd
- upload.R
- Makefile
# SA checklist
1. Exploratory data analysis, trends and likely signals and signs for later
2. Model structures to be explored
3. Setting up the model to run (initial guess etc)
4. Model run and convergence
5. Full suite of model fits and diagnostics
6. Potential modifications after initial fitting
7. Stock status summary
8. Projections etc.
# ICES Suggested details to deal with in benchmark assessments.
- What have been the key areas of focus for the benchmark?
- Have comments in previous technical minutes been addressed?
- Have environmental issues and species interactions been taken into account?
- Are the catch statistics representative of the removals from the stocks (discards, underreporting, misreporting)?
- Has the quality of the survey data been evaluated and presented in the report?
- Has the quality of fishery dependent data (catch, effort, cpue) been evaluated and presented in the report?
- Has the sampling level (commercial catch, discards, surveys) been presented and evaluated?
## Stock assessment modelling:
- Has the sensitivity to model assumptions (different model formulations) been explored?
- Has the sensitivity of parameter estimates to data been evaluated and presented (e.g. bootstrap analysis)?
- Has retrospective analysis been carried out and presented?
- Is the final assessment acceptable to: 1) show the historical stock development 2) show the present status, and 3) be used for a short term forecast?
- Are the recruitment estimates presented and credible?
- If a short term forecast is presented:
- Is the short term forecast consistent with the stock assessment output?
- Has the sensitivity of the short term forecast been discussed (including the current year assumption)?
- Are the (previous) reference points (still) appropriate?
# TMP.R - DESC
# /TMP.R
# Copyright Iago MOSQUEIRA (WMR), 2020
# Author: Iago MOSQUEIRA (WMR) <iago.mosqueira@wur.nl>
#
# Distributed under the terms of the EUPL-1.2
library(FLa4a)
load('data/data.RData')
# CHECK residuals()
control <- AAP.control(pGrp=TRUE, qplat.surveys=7, qplat.Fmatrix=8,
Fage.knots=6, Ftime.knots=22, Wtime.knots=5, mcmc=FALSE)
run <- aap(stock, indices[c("BTS-ISIS", "SNS")], control=control,
wkdir=tempfile())
x <- run@index.res[[1]][-10,ac(1985:2018)] -
(log(index(indices[[1]])) - log(run@index.hat[[1]][-10,ac(1985:2018)]))
# TEST effect of maturity ogive on fit
control <- AAP.control(pGrp=TRUE, qplat.surveys=7, qplat.Fmatrix=8,
Fage.knots=6, Ftime.knots=22, Wtime.knots=5, mcmc=FALSE)
stk <- stock
mat(stk) <- c(0,0.5,0.9,rep(1,7))
run <- aap(stock, indices[c("BTS-ISIS", "SNS")], control=control,
wkdir=tempfile())
residuals <- FLQuants(c(run@index.res, list(
catch.n=residuals(catch.n(run) + 0.5, catch.n(stock)),
catch=residuals(catch(run) + 0.5, catch(stock)))))
ggplot(residuals, aes(x=year, y=as.factor(age))) +
geom_point(aes(size=abs(data), fill=as.factor(sign(data))),
shape=21, alpha=0.4, na.rm=TRUE) + facet_wrap(~qname) +
scale_fill_manual(values=c("black", NA, "white")) +
scale_size(range = c(0.1, 5)) +
theme(legend.position="none") + xlab("") + ylab("")
# TEST effect of control changes on residual patterns
control <- AAP.control(pGrp=TRUE, qplat.surveys=7, qplat.Fmatrix=8,
Fage.knots=6, Ftime.knots=22, Wtime.knots=5, mcmc=FALSE)
run1 <- aap(stock, indices[c("BTS-ISIS", "SNS")], control=control,
wkdir=tempfile())
control <- AAP.control(pGrp=TRUE, qplat.surveys=7, qplat.Fmatrix=8,
Fage.knots=6, Ftime.knots=25, Wtime.knots=5, mcmc=FALSE)
run2 <- aap(stock, indices[c("BTS-ISIS", "SNS")], control=control,
wkdir=tempfile())
# NPIN
# NPIN
pin <- btsisis@stdfile
idx <- which(pin$name == 'log_sel_coff1')
npin <- rbind(pin[1:idx[length(idx)],], pin[seq(idx[1], length=12),], pin[165:441,])
npin$index <- seq(1, nrow(npin))
npin[seq(165, length=12),]
npin[seq(160, length=22),]
control <- AAP.control(pGrp=TRUE, qplat.surveys=7, qplat.Fmatrix=8,
Fage.knots=6, Ftime.knots=26, Wtime.knots=5, mcmc=FALSE)
pbtsisis <- aap(stock, indices[c("BTS-ISIS", "SNS")], control=control,
wkdir="model/aap/pin", stdfile=npin)
......@@ -64,4 +64,4 @@ sol-nsea in ICES area 27.4: IC discard numbers (thousands)
9454 9166 3913 1991 1528 415 15 50 8 1 0 0 0 0 0
5145 5338 5048 1393 291 536 226 4 1 0 1 0 0 0 0
6083 4171 3633 2712 469 89 342 138 0 0 0 0 0 0 0
1970 7801 1720 1464 1195 100 15 126 36 0 0 0 0 0 0
\ No newline at end of file
2928 7760 1704 1448 1186 98 15 125 36 0 0 0 0 0 0
\ No newline at end of file
sol-nsea - ICES area 27.4: 01BTS-ISIS 02SNS
102
103
BTS-ISIS
1985 2018
1 1 0.66 0.75
......@@ -41,53 +41,86 @@ BTS-ISIS
SNS
1970 2018
1 1 0.66 0.75
1 6
1 5410 734 238 35 4.0 0.0
1 903 1831 113 3 28.9 0.0
1 1455 272 149 -1 28.3 0.0
1 5587 935 84 37 13.0 0.0
1 2348 361 65 -1 0.0 4.4
1 525 865 177 18 0.0 17.1
1 1399 74 229 27 5.7 0.0
1 3743 776 104 43 31.7 3.9
1 1548 1355 294 28 99.4 13.3
1 94 408 301 78 0.0 16.7
1 4313 89 109 61 3.3 0.0
1 3737 1413 50 20 0.0 0.0
1 5857 1146 228 7 10.0 0.0
1 2621 1123 121 40 0.0 19.7
1 2493 1100 318 74 8.0 0.0
1 3619 716 167 49 4.4 0.0
1 3705 458 69 31 16.7 0.0
1 1948 944 65 21 0.0 0.0
1 11227 594 282 82 10.2 15.5
1 2831 5005 208 53 18.2 18.6
1 2856 1120 914 100 49.6 12.5
1 1254 2529 514 624 27.2 35.8
1 11114 144 360 195 284.8 20.0
1 1291 3420 154 213 0.0 191.7
1 652 498 934 10 59.3 0.0
1 1362 224 143 411 7.1 31.1
1 218 349 30 36 90.0 10.0
1 10279 154 190 27 58.1 230.0
1 4095 3126 142 99 0.0 10.0
1 1649 972 456 10 20.7 0.0
1 1639 126 166 118 0.0 2.0
1 970 655 107 36 56.2 0.0
1 7548 379 195 -1 30.8 19.2
1 -1 -1 -1 -1 -1 -1
1 1370 624 393 69 53.1 7.5
1 568 163 124 -1 21.3 6.7
1 2726 117 25 30 0.0 0.0
1 849 911 33 40 14.4 0.0
1 1259 259 325 -1 10.0 0.0
1 1932 344 62 103 0.0 0.0
1 2637 237 67 42 23.2 0.0
1 1248 884 211 112 0.0 38.0
1 -1 -1 -1 -1 -1 -1
1 967 427 491 179 50.8 7.6
1 2849 448 45 60 34 0
1 3192 2334 138 160 162 151
1 734 623 495 110 17 43
1 956.7 204.3 209.6 209.7 41.6 5.2
1 1002.3 482.4 163.1 94.1 82.4 5.7
\ No newline at end of file
0 6
1 623.1 5410.3 734.4 237.7 35.4 4.0 0.0
1 10685.1 902.7 1831.1 113.4 2.9 28.9 0.0
1 16.0 1454.7 272.3 148.6 0.0 28.3 0.0
1 895.7 5587.2 935.3 83.8 37.3 13.0 0.0
1 174.4 2347.9 361.4 65.2 0.0 0.0 4.4
1 577.5 525.4 864.5 177.0 17.5 0.0 17.1
1 464.6 1399.4 73.6 229.1 26.7 5.7 0.0
1 1585.0 3742.9 776.1 103.8 43.1 31.7 3.9
1 10370.5 1547.7 1354.7 294.1 28.0 99.4 13.3
1 3922.7 93.8 408.3 300.8 76.9 0.0 16.7
1 5145.8 4312.9 88.9 109.3 61.3 3.3 0.0
1 3240.7 3737.2 1413.1 50.0 20.0 0.0 0.0
1 2147.0 5856.5 1146.2 227.8 6.7 10.0 0.0
1 769.1 2621.1 1123.3 120.6 39.9 0.0 19.7
1 3334.0 2493.1 1099.9 318.3 74.4 8.0 0.0
1 2713.4 3619.4 715.6 167.1 49.3 4.4 0.0
1 742.0 3705.1 457.6 69.2 31.4 16.7 0.0
1 13610.1 1947.9 943.7 64.8 21.3 0.0 0.0
1 522.7 11226.7 593.8 281.6 81.5 10.2 15.5
1 1743.4 2830.7 5005.0 207.6 53.1 18.2 18.6
1 50.8 2856.2 1119.5 914.3 100.4 49.6 12.5
1 3639.7 1253.6 2529.1 513.8 623.9 27.2 35.8
1 302.9 11114.0 144.4 360.4 194.9 284.8 20.0
1 231.3 1290.8 3419.6 153.8 212.8 0.0 191.7
1 4692.7 651.8 498.3 934.1 10.2 59.3 0.0
1 1374.9 1362.1 223.7 142.8 411.1 7.1 31.1
1 2322.3 218.4 349.1 29.6 35.5 90.0 10.0
1 803.0 10279.3 153.6 189.8 26.5 58.1 230.0
1 327.9 4094.6 3126.4 141.7 98.7 0.0 10.0
1 2187.9 1648.9 971.8 455.6 10.0 20.7 0.0
1 70.0 1639.2 125.9 166.3 118.0 0.0 2.0
1 8340.0 970.3 655.4 106.7 35.5 56.2 0.0
1 1127.7 7547.5 379.0 195.3 0.0 30.8 19.2
1 -1 -1 -1 -1 -1 -1 -1
1 162.0 1369.5 624.4 393.0 68.9 53.1 7.5
1 305.0 568.1 162.9 124.0 0.0 21.3 6.7
1 16.0 2726.4 117.1 25.0 30.0 0.0 0.0
1 466.9 848.6 911.0 33.3 39.5 14.4 0.0
1 754.7 1259.1 258.5 325.3 0.0 10.0 0.0
1 2291.0 1931.6 344.4 61.7 102.7 0.0 0.0
1 333.9 2636.9 237.1 67.1 42.2 23.2 0.0
1 136.3 1248.0 883.9 211.3 111.8 0.0 38.0
1 144.7 226.6 159.5 54.0 18.0 0.0 0.0
1 237.3 967.4 426.6 490.5 179.3 50.8 7.6
1 126.0 2849.0 448.2 44.8 60.0 33.6 0.0
1 109.7 3192.0 2333.9 137.8 159.9 162.4 150.6
1 373.2 733.8 623.3 494.6 109.8 16.7 42.9
1 205.9 956.7 204.3 209.6 209.7 41.6 5.2
1 6574.9 1002.3 482.4 163.1 94.1 82.4 5.7
DFS
1990 2018
1 1 0.66 0.75
0 0
1 6.38099979420552
1 167.562790818552
1 9.26602772844516
1 15.3239778292221
1 22.0632417416225
1 7.0647780276121
1 40.2717398630735
1 26.9395747980212
1 -1
1 -1
1 9.5041330402662
1 51.4241905001998
1 58.582992521043
1 10.6093395105163
1 31.2517760468482
1 40.9870116446901
1 12.5667007889326
1 13.7274846340426
1 11.7676192904652
1 27.3315119316548
1 42.8619685862924
1 12.1299798388105
1 11.2261423473025
1 44.8188382604357
1 23.616080619542
1 7.44835180402519
1 12.2755392345859
1 20.9656062744883
1 56.7482834387156
......@@ -63,6 +63,5 @@ sol-nsea in ICES area 27.4: IC landings numbers (thousands)
371 8111 6916 22942 11440 2591 1808 620 840 459 0 0 0 0 0
201 10512 16589 4738 14756 6157 1470 562 393 398 15 45 12 42 33
119 6151 24249 11489 4475 8994 4495 774 278 140 524 102 14 38 36
297 4839 17561 16756 5900 2117 3738 2002 443 866 0 0 0 0 0
416 4928 17641 16818 5909 2118 3745 2005 443 204 44 140 18 5 87
343 11207 9192 11818 10120 3846 1043 1945 797 404 0 0 0 0 0
\ No newline at end of file
416 4928 17641 16818 5909 2118 3745 2005 443 204 44 140 18 5 87
331 11141 9184 11994 10095 3918 1096 1942 804 158 88 118 26 36 10
\ No newline at end of file
......@@ -49,7 +49,7 @@ sol-nsea in ICES area 27.4: mean weight discards (kilograms)
0.059 0.082 0.096 0.104 0.112 0.106 0.120 0.131 0.137 0.163 0.164 0.164 0.164 0.164 0.164
0.059 0.082 0.096 0.104 0.112 0.106 0.120 0.131 0.137 0.163 0.164 0.164 0.164 0.164 0.164
0.046 0.068 0.084 0.091 0.096 0.11 0.124 0.137 0.137 0.169 0.169 0.169 0.169 0.169 0.169
0.054 0.087 0.1 0.107 0.114 0.11 0.124 0.137 0.137 0.169 0.169 0.169 0.169 0.169 0.169
0.054 0.087 0.1 0.107 0.114 0.11 0.124 0.137 0.137 0.169 0.169 0.169 0.169 0.169 0.169
0.065 0.089 0.103 0.111 0.118 0.095 0.124 0.137 0.137 0.169 0.169 0.169 0.169 0.169 0.169
0.068 0.089 0.104 0.109 0.114 0.103 0.107 0.137 0.137 0.169 0.169 0.169 0.169 0.169 0.169
0.066 0.082 0.099 0.109 0.108 0.115 0.113 0.121 0.137 0.169 0.169 0.169 0.169 0.169 0.169
......@@ -64,4 +64,4 @@ sol-nsea in ICES area 27.4: mean weight discards (kilograms)
0.032 0.076 0.095 0.087 0.105 0.117 0.132 0.124 0.159 0.199 0.000 0.000 0.000 0.141 0.180
0.024 0.073 0.087 0.095 0.114 0.108 0.124 0.221 0.214 0.197 0.209 0.235 0.741 0.487 0.507
0.047 0.073 0.086 0.086 0.097 0.124 0.111 0.113 0.287 0.223 0.488 0.340 0.376 0.565 0.379
0.052 0.069 0.086 0.091 0.097 0.102 0.101 0.105 0.127 0 0 0 0 0 0
\ No newline at end of file
0.035 0.069 0.086 0.091 0.097 0.103 0.102 0.105 0.0127 0 0 0 0 0 0
\ No newline at end of file
......@@ -64,4 +64,4 @@ sol-nsea in ICES area 27.4: mean weight of landings
0.145 0.169 0.205 0.240 0.263 0.274 0.304 0.293 0.330 0.330 0.141 0.141 0.141 0.346 0.546
0.143 0.175 0.200 0.236 0.265 0.275 0.273 0.294 0.325 0.397 0.252 0.251 0.627 0.520 0.493
0.109 0.168 0.190 0.226 0.276 0.274 0.313 0.309 0.280 0.311 0.471 0.334 0.456 0.736 0.362
0.123 0.166 0.199 0.235 0.257 0.265 0.243 0.258 0.268 0.286 0.321 0.172 0.522 0.245 0.486
\ No newline at end of file
0.123 0.165 0.198 0.233 0.256 0.263 0.242 0.258 0.268 0.285 0.321 0.180 0.488 0.250 0.399
\ No newline at end of file
......@@ -64,4 +64,4 @@ sol-nsea in ICES area 27.4: mean weight in stock (kilograms)
0.070 0.075 0.142 0.148 0.227 0.244 0.263 0.288 0.370 0.317 0.642 0.682 0.380 0.398 0.741
0.010 0.067 0.151 0.186 0.232 0.248 0.236 0.261 0.221 0.359 0.227 0.221 0.179 0.522 0.716
0.021 0.074 0.131 0.174 0.231 0.242 0.249 0.217 0.233 0.223 0.338 0.608 0.347 0.660 0.321
0.025 0.085 0.146 0.180 0.205 0.237 0.227 0.219 0.259 0.237 0.627 0.566 0.542 0.434 0.406
\ No newline at end of file
0.026 0.084 0.146 0.180 0.205 0.237 0.228 0.219 0.260 0.242 0.591 0.526 0.525 0.423 0.406
\ No newline at end of file
This diff is collapsed.
"","year","age","data","lo","up"
"1",1985,1,885.245267245989,401.417340840198,1995.31575321807
"2",1986,1,3484.01573358945,1977.83509711318,5879.30488612688
"3",1987,1,697.445572809923,382.651176274574,1298.94447929789
"4",1988,1,7154.22177222469,4401.44385372812,11386.5857019332
"5",1989,1,2160.26655674738,1471.88742261744,3301.99479139265
"6",1990,1,2878.2991519854,2076.59184027215,3962.37258619541
"7",1991,1,1253.75621150213,894.14588890932,1749.81767294643
"8",1992,1,17427.5482139418,13071.8665792417,23300.3360666063
"9",1993,1,5237.21849357082,3909.24740201558,7162.62131814864
"10",1994,1,3787.71566951735,2777.01020985792,5206.97507685888
"11",1995,1,7208.45104077995,5474.22610574098,9981.60971739648
"12",1996,1,1588.96870961093,1076.86120987878,2276.2194476237
"13",1997,1,16126.9421962396,11784.2339452113,22195.2020616325
"14",1998,1,2807.05761504071,1890.02794997336,4110.62863864486
"15",1999,1,2621.02319241893,1839.23167374387,3770.2339960633
"16",2000,1,2851.19908739495,2048.51591863913,3988.14814351972
"17",2001,1,2086.82864919004,1423.48857796841,3086.94354954885
"18",2002,1,3235.55759512912,2271.35990240477,4757.45726507869
"19",2003,1,2878.76475133545,2104.15046154157,4022.27228986262
"20",2004,1,994.880704041486,747.200337675367,1339.26209597156
"21",2005,1,1624.11738327746,1182.95269128316,2244.9102560112
"22",2006,1,4261.32686715748,3161.98886843776,5811.43802648199
"23",2007,1,2099.34665483228,1603.45451077164,2727.58533070216
"24",2008,1,2920.97187736084,2225.11321051406,3890.55196904906
"25",2009,1,3193.88552074751,2443.18402522334,4264.03568988006
"26",2010,1,3583.86186016809,2636.52437905925,4873.1807578358
"27",2011,1,2960.32268752066,2147.87581193049,4004.84832592836
"28",2012,1,1356.83573682516,1014.60986106737,1875.93785047227
"29",2013,1,1708.74381170432,1338.83485699953,2308.40427658764
"30",2014,1,4024.32843439171,2970.23158410783,5508.52729600901
"31",2015,1,3163.04285744059,2342.81621654212,4138.83609577776
"32",2016,1,1665.83428059843,1271.90874792618,2165.51617512977
"33",2017,1,6513.22820467539,4807.73906474979,8738.73892710758
"34",2018,1,3506.40747133835,2695.10385026748,4696.35247218188
"35",2019,1,15366.9312458669,11875.1441241968,20215.4979537047
"36",1985,2,817.438903128104,540.648526631306,1220.75852741301
"37",1986,2,1730.66935597855,1186.65455119492,2513.32858742292
"38",1987,2,2199.16663123192,1570.32592076231,3040.9453257453
"39",1988,2,834.690589024656,601.979249384768,1185.29290140677
"40",1989,2,6155.9890414012,4629.25999444573,8154.26065583893
"41",1990,2,2586.0184972998,2084.89705794353,3308.08611226265
"42",1991,2,3404.9552747872,2713.95689554414,4380.88186710778
"43",1992,2,3269.32321704898,2591.38110761935,4152.84006193633
"44",1993,2,11274.6247966547,9065.21104197336,14428.5487392064
"45",1994,2,2738.65867974224,2225.05226118994,3456.8177890392
"46",1995,2,2405.218191486,1914.08348736204,3058.89852583012
"47",1996,2,2326.07428413028,1841.02416398526,2982.6999757574
"48",1997,2,1051.6253964729,828.224945415454,1400.79758327596
"49",1998,2,5093.13247027982,4007.94971987439,6711.65858531083
"50",1999,2,1603.54472844094,1289.01740846595,2053.59918107687
"51",2000,2,874.504257723054,697.948311241222,1126.69423787227
"52",2001,2,1372.19887572455,1060.49120110937,1813.98314805903
"53",2002,2,647.985426586012,515.061338137283,844.909289866143
"54",2003,2,1459.80893528502,1146.76219344992,1861.68058861387
"55",2004,2,1120.36084131959,915.91369639623,1399.26795796916
"56",2005,2,878.222108528244,711.3877309528,1098.80315860818
"57",2006,2,734.286058544723,594.394843205329,918.437562100286
"58",2007,2,3010.36059479447,2541.92194976206,3688.04926685822
"59",2008,2,1449.89406562256,1194.92290046661,1803.23389465374
"60",2009,2,1256.28386381268,1041.11462875969,1548.22868694721
"61",2010,2,1510.62435047336,1218.02692760373,1961.32707261567
"62",2011,2,2706.19079982552,2211.07761922572,3367.89481641217
"63",2012,2,3884.65923762828,3277.51418723183,4727.37621633698
"64",2013,2,877.69936800471,727.655086516262,1065.59672031764
"65",2014,2,2101.01837916173,1693.45402157152,2637.13665696994
"66",2015,2,2587.32107843426,2151.80805741465,3226.28116130571
"67",2016,2,2059.0688176355,1724.97804339985,2536.12887640748
"68",2017,2,1388.9583215711,1142.16368234646,1711.0357131139
"69",2018,2,2161.89624173943,1769.36229088666,2644.83117812926
"70",2019,2,1897.03168548775,1562.28536056808,2335.92714367681
"71",1985,3,782.771743805392,589.776650202583,1082.89230509795
"72",1986,3,642.634105379691,468.114392591008,884.707225306
"73",1987,3,501.503963413311,386.830195997114,651.194827938094
"74",1988,3,760.858350015694,564.163177308481,1014.95482763952
"75",1989,3,651.144320894508,507.660639458296,843.053376043267
"76",1990,3,4910.92202897093,4053.07785648405,6139.73708797783
"77",1991,3,1400.0427068075,1146.8548251285,1751.67445245908
"78",1992,3,3671.38326022095,2981.7674079212,4513.68299091158
"79",1993,3,588.036486444184,475.952808495043,723.048571170687
"80",1994,3,6551.91567126433,5468.28847566805,8044.26690981588
"81",1995,3,1872.13126448901,1529.21768091873,2321.01605523148
"82",1996,3,490.789675945628,398.176330159525,611.272973321633
"83",1997,3,809.044966865757,657.056483562706,1036.61332957673
"84",1998,3,199.773934151247,149.624158418936,263.874836304752
"85",1999,3,2004.16475539777,1628.08619530297,2447.93938607226
"86",2000,3,692.127427036746,566.930846697121,846.793379064722
"87",2001,3,464.158930853127,373.305724885494,594.708836590657
"88",2002,3,574.18762100698,458.804077839069,724.76796090927
"89",2003,3,394.322082896816,319.539266564911,490.796845128407
"90",2004,3,863.074428502919,724.338452948847,1046.23941411343
"91",2005,3,532.87863721785,440.096651027568,645.138904820799
"92",2006,3,241.809933701052,193.554613234562,308.161625488152
"93",2007,3,387.277384379109,318.942847092617,467.904332154186
"94",2008,3,1401.6737850886,1184.86498562276,1686.84521532026
"95",2009,3,694.089548896337,577.125510953634,845.018028803138
"96",2010,3,482.939123736065,394.452818361278,603.493293370734
"97",2011,3,686.858189290339,579.530365494126,841.199941332626
"98",2012,3,1471.50554994168,1262.37669260305,1737.0621901725
"99",2013,3,2052.35827753475,1758.09023392841,2418.79427153548
"100",2014,3,447.567466340266,364.057527667391,548.754744375886
"101",2015,3,1391.63547543337,1200.48240967096,1664.69654212137
"102",2016,3,1386.66245327325,1191.39436016763,1666.1116568041
"103",2017,3,1260.25246086381,1079.76326412246,1491.30677654083
"104",2018,3,614.977392515483,520.410018912056,744.927604797235
"105",2019,3,1262.66153291423,1058.99372026513,1506.33349262537
"106",1985,4,420.417242580565,312.28953518104,562.447606027209
"107",1986,4,399.686840636214,301.399806025177,546.018678796798
"108",1987,4,168.334516889007,126.533041362242,226.907962039662
"109",1988,4,224.501512857991,165.288793038757,294.022294806903
"110",1989,4,554.173684049579,435.91534818489,698.204527609017
"111",1990,4,378.597696105176,313.254422317975,474.143163852849
"112",1991,4,1899.09003065404,1580.13366413738,2360.53445832733
"113",1992,4,721.28865484504,597.922920319081,889.125422054875
"114",1993,4,1869.56064278673,1548.94276448936,2294.63524246467
"115",1994,4,77.7415359552511,57.9302310146899,104.271142900224
"116",1995,4,1924.6973772929,1591.64630147349,2335.59697047747
"117",1996,4,482.902978562259,397.066598126649,601.852451218943
"118",1997,4,191.960870462513,151.88111426108,247.559151871896
"119",1998,4,245.383117503245,192.54811365995,314.490729869866
"120",1999,4,48.3226699998385,35.2771393920622,66.7317568077787
"121",2000,4,366.669834209028,301.660541880956,448.7728099043
"122",2001,4,473.802246488073,377.172679644102,591.237390030885
"123",2002,4,211.781312902552,172.121390418403,267.33549533206
"124",2003,4,257.620979632062,210.166299934167,319.98583935783
"125",2004,4,167.266352523474,140.228398652096,203.501907801577
"126",2005,4,371.787056823297,309.033692705441,451.949535544506
"127",2006,4,381.629170792465,316.72766874831,463.113474346264
"128",2007,4,115.971985374537,94.0929057449056,143.521223076109
"129",2008,4,211.614649662865,176.708294816439,254.873448954296
"130",2009,4,802.970936684493,685.101754608795,951.562147342033
"131",2010,4,266.443011883227,219.693380953603,329.135099693595
"132",2011,4,233.248386205913,193.240800543413,284.618032161289
"133",2012,4,344.183445660832,297.038570500247,407.035812504095
"134",2013,4,592.84997914782,508.948893756541,698.65112102811
"135",2014,4,825.492658858397,696.510303008863,1017.96492532587
"136",2015,4,374.888732088268,316.271060452153,448.579215725869
"137",2016,4,689.404056534758,594.643340210926,809.428283615581
"138",2017,4,627.072206061204,539.997194131454,737.85256494288
"139",2018,4,597.701553262438,515.950689549303,705.422938902458
"140",2019,4,336.881249748993,284.763116375845,412.190941338215
"141",1985,5,189.689168068214,144.587627174244,254.034225406716
"142",1986,5,271.621700468519,210.968395340927,369.48920464668
"143",1987,5,155.621255684513,117.762152489681,211.184076551429
"144",1988,5,84.2456170335031,60.1829950144419,111.554457520308
"145",1989,5,105.896653268125,81.8867201041758,139.375536381797
"146",1990,5,229.783566785458,191.276295307883,286.396200140548
"147",1991,5,102.013630900384,80.5454117570586,132.182014187491
"148",1992,5,1094.6915718993,908.40752715701,1346.28872745855
"149",1993,5,703.712965989135,593.915937518459,863.572837445675
"150",1994,5,378.164428843057,318.13641388173,464.621656548629
"151",1995,5,194.599100994728,158.954554436343,244.744286288332
"152",1996,5,720.17984731712,600.276957318717,894.374116177783
"153",1997,5,221.213967665484,182.887493089514,276.463600045162
"154",1998,5,59.9876860889166,38.9615793596352,93.6679110183275
"155",1999,5,131.664843636464,107.880477801757,168.287700404878
"156",2000,5,73.0206240585843,52.5498344928392,101.794710194417
"157",2001,5,239.62329768346,194.81816707179,303.026158982786
"158",2002,5,97.4998658597499,77.6880763196279,125.145380326495
"159",2003,5,76.9306879913088,60.1655965064006,99.6997464385869
"160",2004,5,130.169866262261,109.288467586528,158.040665347964
"161",2005,5,84.4914376027517,68.5115558468604,105.873880157424
"162",2006,5,201.197050519759,170.219585364412,248.091974703527
"163",2007,5,152.737162369493,127.404053592375,191.577192870635
"164",2008,5,66.6538961357514,53.5008209897404,84.2186583026469
"165",2009,5,116.057501439432,95.6167614632916,144.104977434487
"166",2010,5,248.583807020562,210.405977239539,305.670242094165
"167",2011,5,189.363653098402,159.351440683695,233.404628971101
"168",2012,5,162.981809238246,139.550969878072,197.177177771542
"169",2013,5,189.505975965697,160.233965370773,228.485554027361
"170",2014,5,346.23875020661,296.167419642282,424.166497139657
"171",2015,5,683.571196925805,596.668992843595,814.993113905771
"172",2016,5,202.860986085943,173.466320952725,243.243394195638
"173",2017,5,264.833048285563,228.908024862236,317.20397052952
"174",2018,5,195.787433775128,168.18323947144,237.259001905901
"175",2019,5,235.14155120013,201.390159632031,289.534571434398
"176",1985,6,82.9432737850167,58.6112621917546,116.7975478387
"177",1986,6,107.908989539864,74.2684452608543,160.218120633441
"178",1987,6,115.686264099727,89.5961269134677,154.489147752355
"179",1988,6,80.7690174267657,60.6529501598271,110.66360578076
"180",1989,6,62.0042547051705,47.415300495431,87.3352113889435
"181",1990,6,133.109005252938,108.604646508505,173.640276855682
"182",1991,6,53.4984429195477,36.6704739630874,77.8408733935214
"183",1992,6,17.2930016065626,9.46593858314941,31.1945135946346
"184",1993,6,1250.0416991901,1053.92002668995,1531.16860676653
"185",1994,6,60.8117263526566,47.1496902636166,81.0906045175716
"186",1995,6,266.658473427614,224.553599946951,333.174283027705
"187",1996,6,89.934715737354,71.5658621722292,116.43689873339
"188",1997,6,158.565063978857,130.127539379415,201.671465486241
"189",1998,6,51.466233212728,36.9778254095791,72.7726366681311
"190",1999,6,20.0125083967438,12.6153697069228,34.1723300600795
"191",2000,6,24.4286955983514,15.6937747912768,37.0832252108037
"192",2001,6,34.8246339464476,21.5725671394018,56.006605105203
"193",2002,6,113.98778765299,93.404624846646,144.304830149335
"194",2003,6,58.3948015387814,43.767932508052,80.0412055512566
"195",2004,6,37.4289200654364,29.2159172495064,50.1998449034011
"196",2005,6,73.8382205168454,62.1451821935404,93.6777184123103
"197",2006,6,50.9915685970536,41.0418837431329,67.9315161449616
"198",2007,6,143.129778753423,121.743274859014,175.655778392814
"199",2008,6,72.7955731170764,59.3985800307077,92.2466900533596
"200",2009,6,40.7649583100512,31.1388969392982,56.5525163025519
"201",2010,6,79.2094624143805,64.0386295541864,104.127759933662
"202",2011,6,186.229521989385,158.234072979396,225.889150957922
"203",2012,6,99.448512438135,84.3100662999029,121.233372581765
"204",2013,6,51.1913843254037,41.6610031536543,66.2397505412743
"205",2014,6,79.7290661519387,64.8822740512077,100.625689727086
"206",2015,6,225.153967176972,193.937215589032,273.83360513849
"207",2016,6,357.702968513379,313.080873595206,423.942194942794
"208",2017,6,88.9491673408867,75.9311438711128,109.184773884409
"209",2018,6,134.231577022942,116.398971388706,163.738088379328
"210",2019,6,82.2851628236639,70.4329705955083,102.639002293071
"211",1985,7,0,0,0
"212",1986,7,60.9812791903388,41.6901415864002,88.0102934591576
"213",1987,7,95.7373158732541,70.9239225715573,125.662540951065
"214",1988,7,72.4629744884188,52.2515026871996,98.6486071078635
"215",1989,7,56.5443009891271,40.3175217217491,78.4483985348214
"216",1990,7,35.1401725851834,22.1588162240482,56.6049438791096
"217",1991,7,31.078819698857,17.7044812830979,55.2616685349055
"218",1992,7,60.0184699348616,44.2138832444863,80.1925160934518
"219",1993,7,41.8694005894146,29.452091833952,57.0202393401081
"220",1994,7,335.141618299323,274.44177257596,404.792411009699
"221",1995,7,83.7064026010709,64.3556018820615,107.667013037009
"222",1996,7,118.426354683274,93.6425051856571,146.434281174111
"223",1997,7,28.0947478135694,18.3660519825776,43.6069708038055
"224",1998,7,115.823565583034,89.2435320785745,152.250544754967
"225",1999,7,14.9143731578283,6.85203659908655,32.1419597959186
"226",2000,7,16.9416029008832,9.28921467312485,31.9562731684502
"227",2001,7,5.86225673610374,1.46241918055286,19.5368741284354
"228",2002,7,22.1943307572673,12.1948264831613,42.1436909587224
"229",2003,7,66.929458879729,50.4271063191133,90.8532431890891
"230",2004,7,22.8149843499021,14.8297202607428,35.035984505381
"231",2005,7,34.9185084252709,23.3764182462961,53.6639669226369
"232",2006,7,62.1997698548161,45.173523019824,84.3730986828511
"233",2007,7,30.8974496207097,21.0102271824885,45.9733299008719
"234",2008,7,95.9692454093312,77.2661313718598,119.741091488006
"235",2009,7,89.2430460950006,70.4233291600916,110.938491227714
"236",2010,7,22.4319723964378,14.8497149324567,35.6365270008875
"237",2011,7,30.1714961143009,19.8885946226712,47.7916197558601
"238",2012,7,59.8977376439963,47.1013352737865,75.9752026575584
"239",2013,7,55.7687292998964,40.3966765693263,76.2314819846394
"240",2014,7,25.6335592611312,16.9542701044364,37.858695799586
"241",2015,7,104.370552059383,85.4775773698537,126.931261843305
"242",2016,7,104.633513025617,88.1698185158077,124.267022112991
"243",2017,7,122.618191722756,101.955460377298,146.058658102741
"244",2018,7,46.8772600405679,34.9696910206496,61.4319896116942
"245",2019,7,69.8578879699497,56.6202236272072,86.8363280888972
"246",1985,8,0,0,0
"247",1986,8,0,0,0
"248",1987,8,57.3703020649246,40.011215636992,87.3521022892934
"249",1988,8,42.3626388722229,28.9728841471531,68.7199884436062
"250",1989,8,1.87017828797231,0.340207820944419,13.4452891552401
"251",1990,8,28.3622492394677,17.2571186897828,54.3019241734278
"252",1991,8,26.7694869278082,13.2702673499415,50.0234353302902
"253",1992,8,13.0458953136381,5.6109300576018,29.3391524883311
"254",1993,8,60.8508901641745,45.1349100526201,86.7335535542976
"255",1994,8,50.8640313043414,36.739365893158,73.9458607159241
"256",1995,8,156.343442548305,123.935696991333,208.734234506122
"257",1996,8,28.0037969547093,16.9838325173624,50.8397082834334
"258",1997,8,16.4014622069341,8.25568534061435,32.6230051464355
"259",1998,8,13.789544957763,5.18481622398387,34.9060713852702
"260",1999,8,56.8083760571431,39.1165828844477,85.95260137547
"261",2000,8,2.54221707367671,0.905085714236905,8.0185338043606
"262",2001,8,13.8786548577972,5.45871369922171,36.7534070136326
"263",2002,8,18.2635108584771,10.0370336466421,33.7359843826447
"264",2003,8,6.5771563553275,2.54691763701511,19.9614860772291
"265",2004,8,10.6772367723609,4.72504804796756,24.3434838907356
"266",2005,8,11.1304905249256,5.25884301748409,24.6863494151375
"267",2006,8,21.0095488435824,12.6024453884834,36.8800508214461
"268",2007,8,27.4399242915766,16.6610484446145,46.8632279131564
"269",2008,8,11.7936105415766,6.19277635648824,23.3540919368291
"270",2009,8,63.2358827563054,45.6450396831728,93.3901212937533
"271",2010,8,19.8129018594757,12.0288917566789,34.9118220730455
"272",2011,8,15.0104714127924,8.44656342732611,27.7706289630407
"273",2012,8,20.8300078945041,13.078017516063,35.0539747091289
"274",2013,8,47.951518366234,31.7960183416512,76.9766844421693
"275",2014,8,28.0282302278271,17.8526530894749,48.6740466004068
"276",2015,8,25.8269030049775,16.8638155058191,44.1583879351743
"277",2016,8,23.2930816930189,13.3307269791261,40.3336668931417
"278",2017,8,60.3369127443743,43.1525331894087,87.0237588504865
"279",2018,8,71.2790132571986,56.493718098496,99.8737213147402
"280",2019,8,16.7197926594053,9.90850212154539,29.5675260292362