Comparing assessment results with truth: first pass
Sarah Gaichas and Christine Stawitz
01 October, 2019
So your SS model finally ran!
Now we get to plot outputs using r4ss
library(here)## here() starts at /Users/sgaichas/Documents/0_Data/Atlantis/Poseidon/poseidon-devlibrary(r4ss)
library(atlantisom)
library(dplyr)##
## Attaching package: 'dplyr'## The following objects are masked from 'package:stats':
##
## filter, lag## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, unionlibrary(ggplot2)
library(ggthemes)#Use the right Atlantis output
source(here("config/CC3Config.R"))
#Set species name
species_ss <- "Pacific_sardine" #change to species_ss if we want both
#Directory with SS files
model_dir <- "Sardine_SS_files/"
#Name of SS data file
datfile_name <- "sardEM_3_3.dat"
#ssinput.data <- r4ss::SS_readdat_3.30(paste0("./inst/extdata/", model_dir, datfile_name))
#Name of SS control file
ctlfile_name <- "sardEM_3_3.ctl"
# this does not appear to read correctly; has 20 ages for maturity
# but the model reads it correcty, so don't use for now
#ssinput.ctl <-r4ss::SS_readctl_3.30(paste0("./inst/extdata/", model_dir, ctlfile_name))Basic SS Plots
replist <- r4ss::SS_output(dir = here(paste0("./inst/extdata/", model_dir)), verbose=TRUE, printstats=TRUE,covar=FALSE)## Getting header info from:
## /Users/sgaichas/Documents/0_Data/Atlantis/Poseidon/poseidon-dev/./inst/extdata/Sardine_SS_files//Report.sso## Note: this function tested on SS versions 3.24 and 3.30.
## You are using 3.30.14.0 which SHOULD work with this R code.## Report file time: Tue Oct 1 15:48:44 2019
## Reading full report file
## Got all columns. To speed code, use ncols=61 in the future.
## Got Report file
## Got Forecast-report file
## Setting minimum biomass threshhold to 0.25
## based on US west coast assumption associated with biomass target of 0.4.
## (can replace or override in SS_plots by setting 'minbthresh')
## Got log file. There were NO temporary files were written in this run.
## Got warning file.
## There were 5 warnings in /Users/sgaichas/Documents/0_Data/Atlantis/Poseidon/poseidon-dev/./inst/extdata/Sardine_SS_files//warning.sso
## Finished reading files
## CompReport file separated by this code as follows (rows = Ncomps*Nbins):
## 4800 rows of length comp data,
## 0 rows of generalized size comp data,
## 500 rows of age comp data,
## 19990 rows of conditional age-at-length data,
## 0 rows of ghost fleet age comp data,
## 0 rows of ghost fleet conditional age-at-length data,
## 0 rows of ghost fleet length comp data,
## 0 rows of mean length at age data,
## 0 rows of mean weight at age data,
## 0 rows of 'TAG1' comp data, and
## 0 rows of 'TAG2' comp data.
## Finished dimensioning
## You skipped the covar file## Finished primary run statistics list## Statistics shown below (to turn off, change input to printstats=FALSE)
## $SS_version
## [1] "3.30.14.04-safe;_2019_08_14;_Stock_Synthesis_by_Richard_Methot_(NOAA)_using_ADMB_12.0"
##
## $SS_versionshort
## [1] "3.30"
##
## $SS_versionNumeric
## [1] 3.3
##
## $StartTime
## [1] "StartTime: Tue Oct 1 15:48:44 2019"
##
## $RunTime
## [1] "0 hours, 1 minutes, 10 seconds."
##
## $Files_used
## [1] "Data_File: sardEM_3_3.dat Control_File: sardEM_3_3.ctl"
##
## $Nwarnings
## [1] 5
##
## $warnings
## [1] "#V3.30.14.04-safe;_2019_08_14;_Stock_Synthesis_by_Richard_Methot_(NOAA)_using_ADMB_12.0"
## [2] "#Stock Synthesis (SS) is a work of the U.S. Government and is not subject to copyright protection in the United States."
## [3] "#Foreign copyrights may apply. See copyright.txt for more information."
## [4] "Tue Oct 1 15:48:44 2019"
## [5] ""
## [6] " catch is 0.0 in endyr; this can cause problem in the benchmark and forecast calculations"
## [7] "note: Std_yrmax=-1 in starter, so no variance output for forecast quantities after endyr+1 "
## [8] "You have just one settlement event; recommend changing to recr_dist_method 4 which takes no recr_dist parameters"
## [9] " Minimum size bin is:_12; which is >10cm, which is large for use as size-at-age 0.0 recruitment"
## [10] "Final gradient: 3.8274 is larger than final_conv: 0.0001"
## [11] " N warnings: 5"
## [12] "Number_of_active_parameters_on_or_near_bounds: 0"
##
## $likelihoods_used
## values lambdas
## TOTAL 20059.4000000000014551915228 NA
## Catch 0.0000000000000000477525 NA
## Equil_catch 0.0000000000000000000000 1
## Survey 716.8880000000000336513040 NA
## Length_comp 19203.0999999999985448084772 NA
## Age_comp 55.7323999999999983856469 NA
## Recruitment 83.7015999999999991132427 1
## InitEQ_Regime 0.0000000000000000000000 1
## Forecast_Recruitment 0.0000000000000000000000 1
## Parm_priors 0.0000000000000000000000 1
## Parm_softbounds 0.0001623969999999999971 NA
## Parm_devs 0.0000000000000000000000 1
## Crash_Pen 0.0000000000000000000000 1
##
## $likelihoods_laplace
## values lambdas
## NA NA NA
## Crash_Pen 0 1
##
## $likelihoods_by_fleet
## Label ALL Fishery Survey
## 141 Catch_lambda NA 1.00000e+00 1.0000
## 142 Catch_like 4.77525e-17 4.77525e-17 0.0000
## 143 Surv_lambda NA 0.00000e+00 1.0000
## 144 Surv_like 7.16888e+02 0.00000e+00 716.8880
## 145 Length_lambda NA 1.00000e+00 1.0000
## 146 Length_like 1.92031e+04 7.89269e+03 11310.4000
## 147 Age_lambda NA 1.00000e+00 1.0000
## 148 Age_like 5.57324e+01 3.87906e+01 16.9418
##
## $Parm_devs_detail
## [1] NA
##
## $N_estimated_parameters
## [1] 56
##
## $table_of_phases
##
## -99 -5 -3 -2 -1 1 2 3 4
## 1 1 17 1 1 51 1 2 2
##
## $estimated_non_dev_parameters
## Value Phase Min Max Init Status
## VonBert_K_Fem_GP_1 0.2734270 3 0.05 0.99 0.40 OK
## SR_LN(R0) 15.8324000 1 3.00 25.00 20.45 OK
## Age_inflection_Fishery(1) 1.9807600 2 0.00 5.00 2.00 OK
## Age_95%width_Fishery(1) 0.0123213 3 -20.00 10.00 1.00 OK
## Parm_StDev Gradient Pr_type Prior Pr_SD
## VonBert_K_Fem_GP_1 0.00122423 1.232300 No_prior NA NA
## SR_LN(R0) 0.01894840 0.265441 No_prior NA NA
## Age_inflection_Fishery(1) 13.87890000 0.423146 No_prior NA NA
## Age_95%width_Fishery(1) 8.89175000 3.827400 No_prior NA NA
## Pr_Like Afterbound
## VonBert_K_Fem_GP_1 NA OK
## SR_LN(R0) NA OK
## Age_inflection_Fishery(1) NA OK
## Age_95%width_Fishery(1) NA OK
##
## $maximum_gradient_component
## [1] 3.8274
##
## $parameters_with_highest_gradients
## Value Gradient
## Age_95%width_Fishery(1) 0.0123213 3.827400
## VonBert_K_Fem_GP_1 0.2734270 1.232300
## Age_inflection_Fishery(1) 1.9807600 0.423146
## SR_LN(R0) 15.8324000 0.265441
## Main_RecrDev_54 3.1754100 0.121431
##
## $Length_comp_Eff_N_tuning_check
## Factor Fleet Recommend_var_adj # N Npos min_inputN max_inputN
## 1607 4 1 0.1160580 # 50 50 756 965
## 1608 4 2 0.0474779 # 50 50 983 998
## mean_adj_inputN mean_effN HarMean Curr_Var_Adj Fleet_name
## 1607 939.02 291.5290 108.9810 1 Fishery
## 1608 990.30 79.9557 47.0174 1 Survey
##
## $Age_comp_Eff_N_tuning_check
## Factor Fleet Recommend_var_adj # N Npos min_inputN max_inputN
## 3662 5 1 2.12096 # 50 50 1.000 1.00
## 3663 5 2 122.30200 # 1999 1999 0.001 0.18
## mean_adj_inputN mean_effN HarMean Curr_Var_Adj Fleet_name
## 3662 1.0000000 7.6958 2.12096 1 Fishery
## 3663 0.0249825 14.5296 3.05541 1 Survey
##
## $SBzero
## [1] 335385
##
## $current_depletion
## [1] 1.484573
##
## $last_years_SPR
## [1] 1
##
## $SPRratioLabel
## [1] "1-SPR"
##
## $sigma_R_in
## [1] 0.727
##
## $sigma_R_info
## period N_devs SD_of_devs Var_of_devs mean_SE mean_SEsquared
## 1 Main 50 1.343749 1.805661 0.1055551 0.02072205
## 2 Early+Main 50 1.343749 1.805661 0.1055551 0.02072205
## 3 Early+Main+Late 51 1.330244 1.769548 0.1177403 0.03067905
## sqrt_sum_of_components SD_of_devs_over_sigma_R sqrt_sum_over_sigma_R
## 1 1.351438 1.848348 1.858924
## 2 1.351438 1.848348 1.858924
## 3 1.341725 1.829771 1.845565
## alternative_sigma_R
## 1 1.351438
## 2 1.351438
## 3 1.341725
##
## $rmse_table
## ERA N RMSE RMSE_over_sigmaR mean_BiasAdj
## 1215 main 50 1.33024 3.34806 0
## 1216 early 0 0.00000 0.00000 0
##
## completed SS_outputSS_plots(replist)## Warning in SS_plots(replist): covar information unavailable, changing
## 'uncertainty' to FALSE## Changing 'forecastplot' input to FALSE because all years up to endyr+1 are included by default
## Finished defining objects
## Plots will be written to PNG files in the directory:
## /Users/sgaichas/Documents/0_Data/Atlantis/Poseidon/poseidon-dev/./inst/extdata/Sardine_SS_files//plots
## Starting biology plots (group 1)## Starting selectivity and retention plots (group 2)## Starting timeseries plots (group 3)## finished time series subplot 1: Total biomass (mt)## finished time series subplot 4: Summary biomass (mt)## finished time series subplot 7: Spawning output## finished time series subplot 9: Fraction of unfished## finished time series subplot 11: Age-0 recruits (1,000s)## Plotting Summary F
## Starting recruitment deviation plots (group 4)## Skipping estimation of recruitment bias adjustment (group 5) because uncertainty=FALSE
## Starting spawner-recruit curve plot (group 6)## Starting catch plots (group 7)
## Note: catch is in numbers for some, but not all fleets,
## so be careful interpreting catch plots.## finished catch subplot 1: landings## finished catch subplot 2: landings stacked## finished catch subplot 3: observed and expected landings (if different)## finished catch subplot 9: harvest rate
## Starting SPR plots (group 8)## skipped SPR ratio timeseries: requires both sprtarg>0 and uncertainty=TRUE.## Skipping discard plot (group 9) because no discard data
## Skipping mean weight plot (group 10) because no mean weight data
## Starting index plots (group 11)## Starting numbers at age plots (group 12)## Starting length comp data plots (group 13)## Starting age comp data plots (group 14)## Starting conditional comp data plots (group 15)## Starting fit to length comp plots (group 16)## [1] "Francis Weights - len: Fishery: 0.0095 (0.007-0.0155)"## [1] "Francis Weights - len: Survey: 0.0282 (0.0205-0.0457)"## Starting fit to age comp plots (group 17)## [1] "Francis Weights - age: Fishery: 5.3884 (4.0256-8.4571)"## Starting fit to conditional age-at-length comp plots (group 18)## Starting conditional age-at-length diagnostic plots (group 19)## [1] "Francis CAA Weights: Survey: 23.99891 (15.46668-55.88717)"## Starting mean length-at-age and mean weight-at-age plots (group 20)
## Skipping tag plots (group 21) because no tag data in model
## Starting yield plots (group 22)## Skipping movement plots (group 23) because no movement in model
## Starting data range plots (group 24)## Finished all requested plots in SS_plots function
## Starting diagnostic tables (group 25)
## Wrote table of info on PNG files to:
## /Users/sgaichas/Documents/0_Data/Atlantis/Poseidon/poseidon-dev/./inst/extdata/Sardine_SS_files//plots/plotInfoTable_01-10-2019_22.01.50.csv
## Running 'SS_html':
## By default, this function will look in the directory where PNG files were created
## for CSV files with the name 'plotInfoTable...' written by 'SS_plots.'
## HTML files are written to link to these plots and put in the same directory.
##
## Removing duplicate rows in combined plotInfoTable based on multiple CSV files
## Home HTML file with output will be:
## /Users/sgaichas/Documents/0_Data/Atlantis/Poseidon/poseidon-dev/./inst/extdata/Sardine_SS_files//plots/SS_output.html
## Opening HTML file in your default web-browser.replist$timeseries## Area Yr Era Seas Bio_all Bio_smry SpawnBio Recruit_0 SpawnBio_GP:1
## 948 1 -1 VIRG 1 780801 780801 670770 7515010 670770
## 949 1 0 INIT 1 780801 780801 670770 7515010 670770
## 950 1 1 TIME 1 780801 780801 670770 0 670770
## 951 1 2 TIME 1 780801 780801 670770 7515010 670770
## 952 1 3 TIME 1 780801 780801 670770 7515010 670770
## 953 1 4 TIME 1 780801 780801 670770 7515010 670770
## 954 1 5 TIME 1 780801 780801 670770 7515010 670770
## 955 1 6 TIME 1 780801 780801 670770 7515010 670770
## 956 1 7 TIME 1 780801 780801 670770 7515010 670770
## 957 1 8 TIME 1 780801 780801 670770 7515010 670770
## 958 1 9 TIME 1 780801 780801 670770 7515010 670770
## 959 1 10 TIME 1 780801 780801 670770 7515010 670770
## 960 1 11 TIME 1 780801 780801 670770 7515010 670770
## 961 1 12 TIME 1 780801 780801 670770 7515010 670770
## 962 1 13 TIME 1 780801 780801 670770 7515010 670770
## 963 1 14 TIME 1 780801 780801 670770 7515010 670770
## 964 1 15 TIME 1 780801 780801 670770 7515010 670770
## 965 1 16 TIME 1 780801 780801 670770 7515010 670770
## 966 1 17 TIME 1 780801 780801 670770 7515010 670770
## 967 1 18 TIME 1 780801 780801 670770 7515010 670770
## 968 1 19 TIME 1 780801 780801 670770 7515010 670770
## 969 1 20 TIME 1 780801 780801 670770 7515010 670770
## 970 1 21 TIME 1 780801 780801 670770 7515010 670770
## 971 1 22 TIME 1 780801 780801 670770 7515010 670770
## 972 1 23 TIME 1 780801 780801 670770 7515010 670770
## 973 1 24 TIME 1 780801 780801 670770 7515010 670770
## 974 1 25 TIME 1 780801 780801 670770 7515010 670770
## 975 1 26 TIME 1 780801 780801 670770 7515010 670770
## 976 1 27 TIME 1 780801 780801 670770 7515010 670770
## 977 1 28 TIME 1 780801 780801 670770 7515010 670770
## 978 1 29 TIME 1 780801 780801 670770 7515010 670770
## 979 1 30 TIME 1 780801 780801 670770 7515010 670770
## 980 1 31 TIME 1 756802 756802 665659 524983 665659
## 981 1 32 TIME 1 690362 690362 623378 2827400 623378
## 982 1 33 TIME 1 604603 604603 539935 3988500 539935
## 983 1 34 TIME 1 519518 519518 460652 1469640 460652
## 984 1 35 TIME 1 445503 445503 396419 2466290 396419
## 985 1 36 TIME 1 381139 381139 334165 3049790 334165
## 986 1 37 TIME 1 325942 325942 287224 395109 287224
## 987 1 38 TIME 1 280660 280660 249930 1729700 249930
## 988 1 39 TIME 1 252762 252762 209170 6266950 209170
## 989 1 40 TIME 1 287801 287801 197295 15581500 197295
## 990 1 41 TIME 1 366777 366777 262996 4411850 262996
## 991 1 42 TIME 1 470148 470148 381126 7458430 381126
## 992 1 43 TIME 1 639361 639361 462595 35057900 462595
## 993 1 44 TIME 1 806949 806949 609345 1013420 609345
## 994 1 45 TIME 1 960723 960723 835664 5072260 835664
## 995 1 46 TIME 1 1040790 1040790 897299 18964400 897299
## 996 1 47 TIME 1 1078800 1078800 901135 10013800 901135
## 997 1 48 TIME 1 1105120 1105120 952573 6411250 952573
## 998 1 49 TIME 1 1092950 1092950 963305 7258590 963305
## 999 1 50 TIME 1 1066070 1066070 917430 14845700 917430
## 1000 1 51 TIME 1 1041010 1041010 884436 7938280 884436
## 1001 1 52 TIME 1 1029420 1029420 894723 6379850 894723
## 1002 1 53 TIME 1 1001240 1001240 866814 11899000 866814
## 1003 1 54 TIME 1 1038180 1038180 839932 26633500 839932
## 1004 1 55 TIME 1 1702200 1702200 949069 180363000 949069
## 1005 1 56 TIME 1 2595370 2595370 1646100 22285800 1646100
## 1006 1 57 TIME 1 3734760 3734760 3034160 51408000 3034160
## 1007 1 58 TIME 1 4350270 4350270 3802420 15069600 3802420
## 1008 1 59 TIME 1 4514660 4514660 4038440 21672700 4038440
## 1009 1 60 TIME 1 4377480 4377480 3875360 40500100 3875360
## 1010 1 61 TIME 1 4062990 4062990 3590330 9674330 3590330
## 1011 1 62 TIME 1 3735730 3735730 3329950 25465600 3329950
## 1012 1 63 TIME 1 3360590 3360590 2983620 16571300 2983620
## 1013 1 64 TIME 1 3045410 3045410 2713740 12931200 2713740
## 1014 1 65 TIME 1 2646410 2646410 2385890 5931850 2385890
## 1015 1 66 TIME 1 2300840 2300840 2070260 14242400 2070260
## 1016 1 67 TIME 1 2248440 2248440 1787800 87193300 1787800
## 1017 1 68 TIME 1 2341650 2341650 1832330 3100600 1832330
## 1018 1 69 TIME 1 2508050 2508050 2217860 1860220 2217860
## 1019 1 70 TIME 1 2399640 2399640 2223730 2166430 2223730
## 1020 1 71 TIME 1 2104730 2104730 1941620 3918870 1941620
## 1021 1 72 TIME 1 1760200 1760200 1591560 11511400 1591560
## 1022 1 73 TIME 1 1590550 1590550 1310150 43677100 1310150
## 1023 1 74 TIME 1 1524330 1524330 1242650 622436 1242650
## 1024 1 75 TIME 1 1514370 1514370 1353990 1162800 1353990
## 1025 1 76 TIME 1 1412320 1412320 1301910 3769370 1301910
## 1026 1 77 TIME 1 1215100 1215100 1097240 6691000 1097240
## 1027 1 78 TIME 1 1019400 1019400 910811 3561130 910811
## 1028 1 79 TIME 1 874551 874551 773544 6746820 773544
## 1029 1 80 TIME 1 1107940 1107940 697751 106816000 697751
## 1030 1 81 FORE 1 1529970 1529970 995807 7518830 995807
## SmryBio_SX:1_GP:1 SmryNum_SX:1_GP:1 sel(B):_1 dead(B):_1 retain(B):_1
## 948 780801 22794900 0.00 0.00 0.00
## 949 780801 22794900 0.00 0.00 0.00
## 950 780801 22794900 0.00 0.00 0.00
## 951 780801 22794900 0.00 0.00 0.00
## 952 780801 22794900 0.00 0.00 0.00
## 953 780801 22794900 0.00 0.00 0.00
## 954 780801 22794900 0.00 0.00 0.00
## 955 780801 22794900 0.00 0.00 0.00
## 956 780801 22794900 0.00 0.00 0.00
## 957 780801 22794900 0.00 0.00 0.00
## 958 780801 22794900 0.00 0.00 0.00
## 959 780801 22794900 0.00 0.00 0.00
## 960 780801 22794900 0.00 0.00 0.00
## 961 780801 22794900 0.00 0.00 0.00
## 962 780801 22794900 0.00 0.00 0.00
## 963 780801 22794900 0.00 0.00 0.00
## 964 780801 22794900 0.00 0.00 0.00
## 965 780801 22794900 0.00 0.00 0.00
## 966 780801 22794900 0.00 0.00 0.00
## 967 780801 22794900 0.00 0.00 0.00
## 968 780801 22794900 0.00 0.00 0.00
## 969 780801 22794900 0.00 0.00 0.00
## 970 780801 22794900 0.00 0.00 0.00
## 971 780801 22794900 0.00 0.00 0.00
## 972 780801 22794900 0.00 0.00 0.00
## 973 780801 22794900 0.00 0.00 0.00
## 974 780801 22794900 0.00 0.00 0.00
## 975 780801 22794900 0.00 0.00 0.00
## 976 780801 22794900 0.00 0.00 0.00
## 977 780801 22794900 0.00 0.00 0.00
## 978 780801 22794900 0.00 0.00 0.00
## 979 780801 22794900 0.00 0.00 0.00
## 980 756802 15804800 18726.50 18726.50 18726.50
## 981 690362 13161800 18808.70 18808.70 18808.70
## 982 604603 12608000 16706.50 16706.50 16706.50
## 983 519518 9745830 15619.40 15619.40 15619.40
## 984 445503 8819750 16786.90 16786.90 16786.90
## 985 381139 8785770 14845.90 14845.90 14845.90
## 986 325942 6119520 12645.70 12645.70 12645.70
## 987 280660 5674330 14587.40 14587.40 14587.40
## 988 252762 9919670 12951.50 12951.50 12951.50
## 989 287801 22084600 7868.48 7868.48 7868.48
## 990 366777 19072000 5729.07 5729.07 5729.07
## 991 470148 20100200 8190.15 8190.15 8190.15
## 992 639361 48383500 3143.99 3143.99 3143.99
## 993 806949 33393200 2645.96 2645.96 2645.96
## 994 960723 27392400 4790.62 4790.62 4790.62
## 995 1040790 37251800 6583.58 6583.58 6583.58
## 996 1078800 34899900 5944.72 5944.72 5944.72
## 997 1105120 29710700 7140.06 7140.06 7140.06
## 998 1092950 27070200 8692.41 8692.41 8692.41
## 999 1066070 32879800 9203.40 9203.40 9203.40
## 1000 1041010 29864000 7616.38 7616.38 7616.38
## 1001 1029420 26285900 8152.40 8152.40 8152.40
## 1002 1001240 29408000 8704.40 8704.40 8704.40
## 1003 1038180 46235800 8039.88 8039.88 8039.88
## 1004 1702200 211241000 44707.60 44707.60 44707.60
## 1005 2595370 163057000 11197.40 11197.40 11197.40
## 1006 3734760 160363000 6753.33 6753.33 6753.33
## 1007 4350270 122436000 5095.49 5095.49 5095.49
## 1008 4514660 103661000 4722.43 4722.43 4722.43
## 1009 4377480 109926000 9981.42 9981.42 9981.42
## 1010 4062990 83245000 11827.10 11827.10 11827.10
## 1011 3735730 81115400 23497.30 23497.30 23497.30
## 1012 3360590 70688500 22660.90 22660.90 22660.90
## 1013 3045410 60051000 20145.30 20145.30 20145.30
## 1014 2646410 45958700 18218.10 18218.10 18218.10
## 1015 2300840 44850000 14376.30 14376.30 14376.30
## 1016 2248440 117115000 13144.20 13144.20 13144.20
## 1017 2341650 81459700 5980.30 5980.30 5980.30
## 1018 2508050 56336800 13878.20 13878.20 13878.20
## 1019 2399640 39727400 24180.50 24180.50 24180.50
## 1020 2104730 30281900 32440.20 32440.20 32440.20
## 1021 1760200 31512400 25004.50 25004.50 25004.50
## 1022 1590550 64585300 18469.30 18469.30 18469.30
## 1023 1524330 43720100 11794.40 11794.40 11794.40
## 1024 1514370 30249400 16533.40 16533.40 16533.40
## 1025 1412320 23827500 19437.50 19437.50 19437.50
## 1026 1215100 22463600 26175.80 26175.80 26175.80
## 1027 1019400 18367200 15482.20 15482.20 15482.20
## 1028 874551 18894100 20013.10 20013.10 20013.10
## 1029 1107940 119274000 0.00 0.00 0.00
## 1030 1529970 87470300 0.00 0.00 0.00
## sel(N):_1 dead(N):_1 retain(N):_1 obs_cat:_1 F:_1 SSB_vir_LH
## 948 0 0 0 0 0.00000000 717012
## 949 0 0 0 0 0.00000000 717012
## 950 0 0 0 0 0.00000000 717012
## 951 0 0 0 0 0.00000000 717012
## 952 0 0 0 0 0.00000000 717012
## 953 0 0 0 0 0.00000000 717012
## 954 0 0 0 0 0.00000000 717012
## 955 0 0 0 0 0.00000000 717012
## 956 0 0 0 0 0.00000000 717012
## 957 0 0 0 0 0.00000000 717012
## 958 0 0 0 0 0.00000000 717012
## 959 0 0 0 0 0.00000000 717012
## 960 0 0 0 0 0.00000000 717012
## 961 0 0 0 0 0.00000000 717012
## 962 0 0 0 0 0.00000000 717012
## 963 0 0 0 0 0.00000000 717012
## 964 0 0 0 0 0.00000000 717012
## 965 0 0 0 0 0.00000000 717012
## 966 0 0 0 0 0.00000000 717012
## 967 0 0 0 0 0.00000000 717012
## 968 0 0 0 0 0.00000000 717012
## 969 0 0 0 0 0.00000000 717012
## 970 0 0 0 0 0.00000000 717012
## 971 0 0 0 0 0.00000000 717012
## 972 0 0 0 0 0.00000000 717012
## 973 0 0 0 0 0.00000000 717012
## 974 0 0 0 0 0.00000000 717012
## 975 0 0 0 0 0.00000000 717012
## 976 0 0 0 0 0.00000000 717012
## 977 0 0 0 0 0.00000000 717012
## 978 0 0 0 0 0.00000000 717012
## 979 0 0 0 0 0.00000000 717012
## 980 319865 319865 319865 319865 0.02578960 714612
## 981 249922 249922 249922 249922 0.02976120 669662
## 982 215679 215679 215679 215679 0.03086530 580126
## 983 220781 220781 220781 220781 0.03297180 495119
## 984 216810 216810 216810 216810 0.04228980 426739
## 985 203003 203003 203003 203003 0.04395030 359820
## 986 193817 193817 193817 193817 0.04203750 309181
## 987 185898 185898 185898 185898 0.05879990 269785
## 988 180176 180176 180176 180176 0.06177890 225881
## 989 176713 176713 176713 176713 0.03370460 212053
## 990 175377 175377 175377 175377 0.01471820 281812
## 991 182145 182145 182145 182145 0.01766680 408593
## 992 64694 64694 64694 64694 0.00592886 494962
## 993 78720 78720 78720 78720 0.00297530 651675
## 994 91349 91349 91349 91349 0.00497869 894014
## 995 103954 103954 103954 103954 0.00693209 960260
## 996 116472 116472 116472 116472 0.00572272 964173
## 997 127937 127937 127937 127937 0.00670227 1019380
## 998 137312 137312 137312 137312 0.00846070 1031160
## 999 140571 140571 140571 140571 0.00952496 982226
## 1000 138208 138208 138208 138208 0.00771030 946620
## 1001 136455 136455 136455 136455 0.00837186 957736
## 1002 135849 135849 135849 135849 0.00947842 928029
## 1003 140378 140378 140378 140378 0.00875973 899133
## 1004 1018980 1018980 1018980 1018980 0.04104210 1021030
## 1005 424892 424892 424892 424892 0.00370006 1760660
## 1006 157952 157952 157952 157952 0.00176277 3244280
## 1007 101970 101970 101970 101970 0.00115664 4065330
## 1008 73612 73612 73612 73612 0.00109124 4317630
## 1009 141543 141543 141543 141543 0.00248166 4144230
## 1010 185664 185664 185664 185664 0.00307760 3839810
## 1011 315029 315029 315029 315029 0.00689851 3563600
## 1012 324810 324810 324810 324810 0.00733011 3193200
## 1013 278613 278613 278613 278613 0.00721515 2904300
## 1014 245273 245273 245273 245273 0.00747680 2553540
## 1015 174959 174959 174959 174959 0.00696702 2215540
## 1016 178687 178687 178687 178687 0.00729344 1913340
## 1017 156851 156851 156851 156851 0.00244962 1959450
## 1018 249209 249209 249209 249209 0.00556691 2372960
## 1019 328456 328456 328456 328456 0.01066620 2381260
## 1020 366142 366142 366142 366142 0.01699410 2081360
## 1021 264745 264745 264745 264745 0.01620270 1705870
## 1022 240040 240040 240040 240040 0.01407260 1403720
## 1023 270350 270350 270350 270350 0.00769022 1330020
## 1024 268981 268981 268981 268981 0.01128080 1450060
## 1025 245319 245319 245319 245319 0.01494860 1395130
## 1026 309845 309845 309845 309845 0.02414220 1177600
## 1027 202432 202432 202432 202432 0.01676890 976322
## 1028 254846 254846 254846 254846 0.02581230 830425
## 1029 0 0 0 0 0.00000000 745853
## 1030 0 0 0 NA 0.00000000 1064460
## ABC_buffer
## 948 NA
## 949 NA
## 950 NA
## 951 NA
## 952 NA
## 953 NA
## 954 NA
## 955 NA
## 956 NA
## 957 NA
## 958 NA
## 959 NA
## 960 NA
## 961 NA
## 962 NA
## 963 NA
## 964 NA
## 965 NA
## 966 NA
## 967 NA
## 968 NA
## 969 NA
## 970 NA
## 971 NA
## 972 NA
## 973 NA
## 974 NA
## 975 NA
## 976 NA
## 977 NA
## 978 NA
## 979 NA
## 980 NA
## 981 NA
## 982 NA
## 983 NA
## 984 NA
## 985 NA
## 986 NA
## 987 NA
## 988 NA
## 989 NA
## 990 NA
## 991 NA
## 992 NA
## 993 NA
## 994 NA
## 995 NA
## 996 NA
## 997 NA
## 998 NA
## 999 NA
## 1000 NA
## 1001 NA
## 1002 NA
## 1003 NA
## 1004 NA
## 1005 NA
## 1006 NA
## 1007 NA
## 1008 NA
## 1009 NA
## 1010 NA
## 1011 NA
## 1012 NA
## 1013 NA
## 1014 NA
## 1015 NA
## 1016 NA
## 1017 NA
## 1018 NA
## 1019 NA
## 1020 NA
## 1021 NA
## 1022 NA
## 1023 NA
## 1024 NA
## 1025 NA
## 1026 NA
## 1027 NA
## 1028 NA
## 1029 NA
## 1030 0replist$wtatage## Yr Seas Sex Bio_Pattern BirthSeas Fleet 0 1
## 1 1 1 1 1 1 1 0.00343331 0.006000580
## 2 1 1 1 1 1 2 0.00343331 0.006000580
## 3 1 1 1 1 1 -2 0.00000000 0.000343331
## 4 1 1 1 1 1 0 0.00343331 0.003433310
## 5 1 1 1 1 1 -1 0.00343331 0.006000580
## 6 2 1 1 1 1 1 0.00343331 0.006000580
## 7 2 1 1 1 1 2 0.00343331 0.006000580
## 8 2 1 1 1 1 -2 0.00000000 0.000343331
## 9 2 1 1 1 1 0 0.00343331 0.003433310
## 10 2 1 1 1 1 -1 0.00343331 0.006000580
## 11 3 1 1 1 1 1 0.00343331 0.006000580
## 12 3 1 1 1 1 2 0.00343331 0.006000580
## 13 3 1 1 1 1 -2 0.00000000 0.000343331
## 14 3 1 1 1 1 0 0.00343331 0.003433310
## 15 3 1 1 1 1 -1 0.00343331 0.006000580
## 16 4 1 1 1 1 1 0.00343331 0.006000580
## 17 4 1 1 1 1 2 0.00343331 0.006000580
## 18 4 1 1 1 1 -2 0.00000000 0.000343331
## 19 4 1 1 1 1 0 0.00343331 0.003433310
## 20 4 1 1 1 1 -1 0.00343331 0.006000580
## 21 5 1 1 1 1 1 0.00343331 0.006000580
## 22 5 1 1 1 1 2 0.00343331 0.006000580
## 23 5 1 1 1 1 -2 0.00000000 0.000343331
## 24 5 1 1 1 1 0 0.00343331 0.003433310
## 25 5 1 1 1 1 -1 0.00343331 0.006000580
## 26 6 1 1 1 1 1 0.00343331 0.006000580
## 27 6 1 1 1 1 2 0.00343331 0.006000580
## 28 6 1 1 1 1 -2 0.00000000 0.000343331
## 29 6 1 1 1 1 0 0.00343331 0.003433310
## 30 6 1 1 1 1 -1 0.00343331 0.006000580
## 31 7 1 1 1 1 1 0.00343331 0.006000580
## 32 7 1 1 1 1 2 0.00343331 0.006000580
## 33 7 1 1 1 1 -2 0.00000000 0.000343331
## 34 7 1 1 1 1 0 0.00343331 0.003433310
## 35 7 1 1 1 1 -1 0.00343331 0.006000580
## 36 8 1 1 1 1 1 0.00343331 0.006000580
## 37 8 1 1 1 1 2 0.00343331 0.006000580
## 38 8 1 1 1 1 -2 0.00000000 0.000343331
## 39 8 1 1 1 1 0 0.00343331 0.003433310
## 40 8 1 1 1 1 -1 0.00343331 0.006000580
## 41 9 1 1 1 1 1 0.00343331 0.006000580
## 42 9 1 1 1 1 2 0.00343331 0.006000580
## 43 9 1 1 1 1 -2 0.00000000 0.000343331
## 44 9 1 1 1 1 0 0.00343331 0.003433310
## 45 9 1 1 1 1 -1 0.00343331 0.006000580
## 46 10 1 1 1 1 1 0.00343331 0.006000580
## 47 10 1 1 1 1 2 0.00343331 0.006000580
## 48 10 1 1 1 1 -2 0.00000000 0.000343331
## 49 10 1 1 1 1 0 0.00343331 0.003433310
## 50 10 1 1 1 1 -1 0.00343331 0.006000580
## 51 11 1 1 1 1 1 0.00343331 0.006000580
## 52 11 1 1 1 1 2 0.00343331 0.006000580
## 53 11 1 1 1 1 -2 0.00000000 0.000343331
## 54 11 1 1 1 1 0 0.00343331 0.003433310
## 55 11 1 1 1 1 -1 0.00343331 0.006000580
## 56 12 1 1 1 1 1 0.00343331 0.006000580
## 57 12 1 1 1 1 2 0.00343331 0.006000580
## 58 12 1 1 1 1 -2 0.00000000 0.000343331
## 59 12 1 1 1 1 0 0.00343331 0.003433310
## 60 12 1 1 1 1 -1 0.00343331 0.006000580
## 61 13 1 1 1 1 1 0.00343331 0.006000580
## 62 13 1 1 1 1 2 0.00343331 0.006000580
## 63 13 1 1 1 1 -2 0.00000000 0.000343331
## 64 13 1 1 1 1 0 0.00343331 0.003433310
## 65 13 1 1 1 1 -1 0.00343331 0.006000580
## 66 14 1 1 1 1 1 0.00343331 0.006000580
## 67 14 1 1 1 1 2 0.00343331 0.006000580
## 68 14 1 1 1 1 -2 0.00000000 0.000343331
## 69 14 1 1 1 1 0 0.00343331 0.003433310
## 70 14 1 1 1 1 -1 0.00343331 0.006000580
## 71 15 1 1 1 1 1 0.00343331 0.006000580
## 72 15 1 1 1 1 2 0.00343331 0.006000580
## 73 15 1 1 1 1 -2 0.00000000 0.000343331
## 74 15 1 1 1 1 0 0.00343331 0.003433310
## 75 15 1 1 1 1 -1 0.00343331 0.006000580
## 76 16 1 1 1 1 1 0.00343331 0.006000580
## 77 16 1 1 1 1 2 0.00343331 0.006000580
## 78 16 1 1 1 1 -2 0.00000000 0.000343331
## 79 16 1 1 1 1 0 0.00343331 0.003433310
## 80 16 1 1 1 1 -1 0.00343331 0.006000580
## 81 17 1 1 1 1 1 0.00343331 0.006000580
## 82 17 1 1 1 1 2 0.00343331 0.006000580
## 83 17 1 1 1 1 -2 0.00000000 0.000343331
## 84 17 1 1 1 1 0 0.00343331 0.003433310
## 85 17 1 1 1 1 -1 0.00343331 0.006000580
## 86 18 1 1 1 1 1 0.00343331 0.006000580
## 87 18 1 1 1 1 2 0.00343331 0.006000580
## 88 18 1 1 1 1 -2 0.00000000 0.000343331
## 89 18 1 1 1 1 0 0.00343331 0.003433310
## 90 18 1 1 1 1 -1 0.00343331 0.006000580
## 91 19 1 1 1 1 1 0.00343331 0.006000580
## 92 19 1 1 1 1 2 0.00343331 0.006000580
## 93 19 1 1 1 1 -2 0.00000000 0.000343331
## 94 19 1 1 1 1 0 0.00343331 0.003433310
## 95 19 1 1 1 1 -1 0.00343331 0.006000580
## 96 20 1 1 1 1 1 0.00343331 0.006000580
## 97 20 1 1 1 1 2 0.00343331 0.006000580
## 98 20 1 1 1 1 -2 0.00000000 0.000343331
## 99 20 1 1 1 1 0 0.00343331 0.003433310
## 100 20 1 1 1 1 -1 0.00343331 0.006000580
## 101 21 1 1 1 1 1 0.00343331 0.006000580
## 102 21 1 1 1 1 2 0.00343331 0.006000580
## 103 21 1 1 1 1 -2 0.00000000 0.000343331
## 104 21 1 1 1 1 0 0.00343331 0.003433310
## 105 21 1 1 1 1 -1 0.00343331 0.006000580
## 106 22 1 1 1 1 1 0.00343331 0.006000580
## 107 22 1 1 1 1 2 0.00343331 0.006000580
## 108 22 1 1 1 1 -2 0.00000000 0.000343331
## 109 22 1 1 1 1 0 0.00343331 0.003433310
## 110 22 1 1 1 1 -1 0.00343331 0.006000580
## 111 23 1 1 1 1 1 0.00343331 0.006000580
## 112 23 1 1 1 1 2 0.00343331 0.006000580
## 113 23 1 1 1 1 -2 0.00000000 0.000343331
## 114 23 1 1 1 1 0 0.00343331 0.003433310
## 115 23 1 1 1 1 -1 0.00343331 0.006000580
## 116 24 1 1 1 1 1 0.00343331 0.006000580
## 117 24 1 1 1 1 2 0.00343331 0.006000580
## 118 24 1 1 1 1 -2 0.00000000 0.000343331
## 119 24 1 1 1 1 0 0.00343331 0.003433310
## 120 24 1 1 1 1 -1 0.00343331 0.006000580
## 121 25 1 1 1 1 1 0.00343331 0.006000580
## 122 25 1 1 1 1 2 0.00343331 0.006000580
## 123 25 1 1 1 1 -2 0.00000000 0.000343331
## 124 25 1 1 1 1 0 0.00343331 0.003433310
## 125 25 1 1 1 1 -1 0.00343331 0.006000580
## 126 26 1 1 1 1 1 0.00343331 0.006000580
## 127 26 1 1 1 1 2 0.00343331 0.006000580
## 128 26 1 1 1 1 -2 0.00000000 0.000343331
## 129 26 1 1 1 1 0 0.00343331 0.003433310
## 130 26 1 1 1 1 -1 0.00343331 0.006000580
## 131 27 1 1 1 1 1 0.00343331 0.006000580
## 132 27 1 1 1 1 2 0.00343331 0.006000580
## 133 27 1 1 1 1 -2 0.00000000 0.000343331
## 134 27 1 1 1 1 0 0.00343331 0.003433310
## 135 27 1 1 1 1 -1 0.00343331 0.006000580
## 136 28 1 1 1 1 1 0.00343331 0.006000580
## 137 28 1 1 1 1 2 0.00343331 0.006000580
## 138 28 1 1 1 1 -2 0.00000000 0.000343331
## 139 28 1 1 1 1 0 0.00343331 0.003433310
## 140 28 1 1 1 1 -1 0.00343331 0.006000580
## 141 29 1 1 1 1 1 0.00343331 0.006000580
## 142 29 1 1 1 1 2 0.00343331 0.006000580
## 143 29 1 1 1 1 -2 0.00000000 0.000343331
## 144 29 1 1 1 1 0 0.00343331 0.003433310
## 145 29 1 1 1 1 -1 0.00343331 0.006000580
## 146 30 1 1 1 1 1 0.00343331 0.006000580
## 147 30 1 1 1 1 2 0.00343331 0.006000580
## 148 30 1 1 1 1 -2 0.00000000 0.000343331
## 149 30 1 1 1 1 0 0.00343331 0.003433310
## 150 30 1 1 1 1 -1 0.00343331 0.006000580
## 151 31 1 1 1 1 1 0.00343331 0.006000580
## 152 31 1 1 1 1 2 0.00343331 0.006000580
## 153 31 1 1 1 1 -2 0.00000000 0.000343331
## 154 31 1 1 1 1 0 0.00343331 0.003433310
## 155 31 1 1 1 1 -1 0.00343331 0.006000580
## 156 32 1 1 1 1 1 0.00343331 0.006000580
## 157 32 1 1 1 1 2 0.00343331 0.006000580
## 158 32 1 1 1 1 -2 0.00000000 0.000343331
## 159 32 1 1 1 1 0 0.00343331 0.003433310
## 160 32 1 1 1 1 -1 0.00343331 0.006000580
## 161 33 1 1 1 1 1 0.00343331 0.006000580
## 162 33 1 1 1 1 2 0.00343331 0.006000580
## 163 33 1 1 1 1 -2 0.00000000 0.000343331
## 164 33 1 1 1 1 0 0.00343331 0.003433310
## 165 33 1 1 1 1 -1 0.00343331 0.006000580
## 166 34 1 1 1 1 1 0.00343331 0.006000580
## 167 34 1 1 1 1 2 0.00343331 0.006000580
## 168 34 1 1 1 1 -2 0.00000000 0.000343331
## 169 34 1 1 1 1 0 0.00343331 0.003433310
## 170 34 1 1 1 1 -1 0.00343331 0.006000580
## 171 35 1 1 1 1 1 0.00343331 0.006000580
## 172 35 1 1 1 1 2 0.00343331 0.006000580
## 173 35 1 1 1 1 -2 0.00000000 0.000343331
## 174 35 1 1 1 1 0 0.00343331 0.003433310
## 175 35 1 1 1 1 -1 0.00343331 0.006000580
## 176 36 1 1 1 1 1 0.00343331 0.006000580
## 177 36 1 1 1 1 2 0.00343331 0.006000580
## 178 36 1 1 1 1 -2 0.00000000 0.000343331
## 179 36 1 1 1 1 0 0.00343331 0.003433310
## 180 36 1 1 1 1 -1 0.00343331 0.006000580
## 181 37 1 1 1 1 1 0.00343331 0.006000580
## 182 37 1 1 1 1 2 0.00343331 0.006000580
## 183 37 1 1 1 1 -2 0.00000000 0.000343331
## 184 37 1 1 1 1 0 0.00343331 0.003433310
## 185 37 1 1 1 1 -1 0.00343331 0.006000580
## 186 38 1 1 1 1 1 0.00343331 0.006000580
## 187 38 1 1 1 1 2 0.00343331 0.006000580
## 188 38 1 1 1 1 -2 0.00000000 0.000343331
## 189 38 1 1 1 1 0 0.00343331 0.003433310
## 190 38 1 1 1 1 -1 0.00343331 0.006000580
## 191 39 1 1 1 1 1 0.00343331 0.006000580
## 192 39 1 1 1 1 2 0.00343331 0.006000580
## 193 39 1 1 1 1 -2 0.00000000 0.000343331
## 194 39 1 1 1 1 0 0.00343331 0.003433310
## 195 39 1 1 1 1 -1 0.00343331 0.006000580
## 196 40 1 1 1 1 1 0.00343331 0.006000580
## 197 40 1 1 1 1 2 0.00343331 0.006000580
## 198 40 1 1 1 1 -2 0.00000000 0.000343331
## 199 40 1 1 1 1 0 0.00343331 0.003433310
## 200 40 1 1 1 1 -1 0.00343331 0.006000580
## 201 41 1 1 1 1 1 0.00343331 0.006000580
## 202 41 1 1 1 1 2 0.00343331 0.006000580
## 203 41 1 1 1 1 -2 0.00000000 0.000343331
## 204 41 1 1 1 1 0 0.00343331 0.003433310
## 205 41 1 1 1 1 -1 0.00343331 0.006000580
## 206 42 1 1 1 1 1 0.00343331 0.006000580
## 207 42 1 1 1 1 2 0.00343331 0.006000580
## 208 42 1 1 1 1 -2 0.00000000 0.000343331
## 209 42 1 1 1 1 0 0.00343331 0.003433310
## 210 42 1 1 1 1 -1 0.00343331 0.006000580
## 211 43 1 1 1 1 1 0.00343331 0.006000580
## 212 43 1 1 1 1 2 0.00343331 0.006000580
## 213 43 1 1 1 1 -2 0.00000000 0.000343331
## 214 43 1 1 1 1 0 0.00343331 0.003433310
## 215 43 1 1 1 1 -1 0.00343331 0.006000580
## 216 44 1 1 1 1 1 0.00343331 0.006000580
## 217 44 1 1 1 1 2 0.00343331 0.006000580
## 218 44 1 1 1 1 -2 0.00000000 0.000343331
## 219 44 1 1 1 1 0 0.00343331 0.003433310
## 220 44 1 1 1 1 -1 0.00343331 0.006000580
## 221 45 1 1 1 1 1 0.00343331 0.006000580
## 222 45 1 1 1 1 2 0.00343331 0.006000580
## 223 45 1 1 1 1 -2 0.00000000 0.000343331
## 224 45 1 1 1 1 0 0.00343331 0.003433310
## 225 45 1 1 1 1 -1 0.00343331 0.006000580
## 226 46 1 1 1 1 1 0.00343331 0.006000580
## 227 46 1 1 1 1 2 0.00343331 0.006000580
## 228 46 1 1 1 1 -2 0.00000000 0.000343331
## 229 46 1 1 1 1 0 0.00343331 0.003433310
## 230 46 1 1 1 1 -1 0.00343331 0.006000580
## 231 47 1 1 1 1 1 0.00343331 0.006000580
## 232 47 1 1 1 1 2 0.00343331 0.006000580
## 233 47 1 1 1 1 -2 0.00000000 0.000343331
## 234 47 1 1 1 1 0 0.00343331 0.003433310
## 235 47 1 1 1 1 -1 0.00343331 0.006000580
## 236 48 1 1 1 1 1 0.00343331 0.006000580
## 237 48 1 1 1 1 2 0.00343331 0.006000580
## 238 48 1 1 1 1 -2 0.00000000 0.000343331
## 239 48 1 1 1 1 0 0.00343331 0.003433310
## 240 48 1 1 1 1 -1 0.00343331 0.006000580
## 241 49 1 1 1 1 1 0.00343331 0.006000580
## 242 49 1 1 1 1 2 0.00343331 0.006000580
## 243 49 1 1 1 1 -2 0.00000000 0.000343331
## 244 49 1 1 1 1 0 0.00343331 0.003433310
## 245 49 1 1 1 1 -1 0.00343331 0.006000580
## 246 50 1 1 1 1 1 0.00343331 0.006000580
## 247 50 1 1 1 1 2 0.00343331 0.006000580
## 248 50 1 1 1 1 -2 0.00000000 0.000343331
## 249 50 1 1 1 1 0 0.00343331 0.003433310
## 250 50 1 1 1 1 -1 0.00343331 0.006000580
## 251 51 1 1 1 1 1 0.00343331 0.006000580
## 252 51 1 1 1 1 2 0.00343331 0.006000580
## 253 51 1 1 1 1 -2 0.00000000 0.000343331
## 254 51 1 1 1 1 0 0.00343331 0.003433310
## 255 51 1 1 1 1 -1 0.00343331 0.006000580
## 256 52 1 1 1 1 1 0.00343331 0.006000580
## 257 52 1 1 1 1 2 0.00343331 0.006000580
## 258 52 1 1 1 1 -2 0.00000000 0.000343331
## 259 52 1 1 1 1 0 0.00343331 0.003433310
## 260 52 1 1 1 1 -1 0.00343331 0.006000580
## 261 53 1 1 1 1 1 0.00343331 0.006000580
## 262 53 1 1 1 1 2 0.00343331 0.006000580
## 263 53 1 1 1 1 -2 0.00000000 0.000343331
## 264 53 1 1 1 1 0 0.00343331 0.003433310
## 265 53 1 1 1 1 -1 0.00343331 0.006000580
## 266 54 1 1 1 1 1 0.00343331 0.006000580
## 267 54 1 1 1 1 2 0.00343331 0.006000580
## 268 54 1 1 1 1 -2 0.00000000 0.000343331
## 269 54 1 1 1 1 0 0.00343331 0.003433310
## 270 54 1 1 1 1 -1 0.00343331 0.006000580
## 271 55 1 1 1 1 1 0.00343331 0.006000580
## 272 55 1 1 1 1 2 0.00343331 0.006000580
## 273 55 1 1 1 1 -2 0.00000000 0.000343331
## 274 55 1 1 1 1 0 0.00343331 0.003433310
## 275 55 1 1 1 1 -1 0.00343331 0.006000580
## 276 56 1 1 1 1 1 0.00343331 0.006000580
## 277 56 1 1 1 1 2 0.00343331 0.006000580
## 278 56 1 1 1 1 -2 0.00000000 0.000343331
## 279 56 1 1 1 1 0 0.00343331 0.003433310
## 280 56 1 1 1 1 -1 0.00343331 0.006000580
## 281 57 1 1 1 1 1 0.00343331 0.006000580
## 282 57 1 1 1 1 2 0.00343331 0.006000580
## 283 57 1 1 1 1 -2 0.00000000 0.000343331
## 284 57 1 1 1 1 0 0.00343331 0.003433310
## 285 57 1 1 1 1 -1 0.00343331 0.006000580
## 286 58 1 1 1 1 1 0.00343331 0.006000580
## 287 58 1 1 1 1 2 0.00343331 0.006000580
## 288 58 1 1 1 1 -2 0.00000000 0.000343331
## 289 58 1 1 1 1 0 0.00343331 0.003433310
## 290 58 1 1 1 1 -1 0.00343331 0.006000580
## 291 59 1 1 1 1 1 0.00343331 0.006000580
## 292 59 1 1 1 1 2 0.00343331 0.006000580
## 293 59 1 1 1 1 -2 0.00000000 0.000343331
## 294 59 1 1 1 1 0 0.00343331 0.003433310
## 295 59 1 1 1 1 -1 0.00343331 0.006000580
## 296 60 1 1 1 1 1 0.00343331 0.006000580
## 297 60 1 1 1 1 2 0.00343331 0.006000580
## 298 60 1 1 1 1 -2 0.00000000 0.000343331
## 299 60 1 1 1 1 0 0.00343331 0.003433310
## 300 60 1 1 1 1 -1 0.00343331 0.006000580
## 301 61 1 1 1 1 1 0.00343331 0.006000580
## 302 61 1 1 1 1 2 0.00343331 0.006000580
## 303 61 1 1 1 1 -2 0.00000000 0.000343331
## 304 61 1 1 1 1 0 0.00343331 0.003433310
## 305 61 1 1 1 1 -1 0.00343331 0.006000580
## 306 62 1 1 1 1 1 0.00343331 0.006000580
## 307 62 1 1 1 1 2 0.00343331 0.006000580
## 308 62 1 1 1 1 -2 0.00000000 0.000343331
## 309 62 1 1 1 1 0 0.00343331 0.003433310
## 310 62 1 1 1 1 -1 0.00343331 0.006000580
## 311 63 1 1 1 1 1 0.00343331 0.006000580
## 312 63 1 1 1 1 2 0.00343331 0.006000580
## 313 63 1 1 1 1 -2 0.00000000 0.000343331
## 314 63 1 1 1 1 0 0.00343331 0.003433310
## 315 63 1 1 1 1 -1 0.00343331 0.006000580
## 316 64 1 1 1 1 1 0.00343331 0.006000580
## 317 64 1 1 1 1 2 0.00343331 0.006000580
## 318 64 1 1 1 1 -2 0.00000000 0.000343331
## 319 64 1 1 1 1 0 0.00343331 0.003433310
## 320 64 1 1 1 1 -1 0.00343331 0.006000580
## 321 65 1 1 1 1 1 0.00343331 0.006000580
## 322 65 1 1 1 1 2 0.00343331 0.006000580
## 323 65 1 1 1 1 -2 0.00000000 0.000343331
## 324 65 1 1 1 1 0 0.00343331 0.003433310
## 325 65 1 1 1 1 -1 0.00343331 0.006000580
## 326 66 1 1 1 1 1 0.00343331 0.006000580
## 327 66 1 1 1 1 2 0.00343331 0.006000580
## 328 66 1 1 1 1 -2 0.00000000 0.000343331
## 329 66 1 1 1 1 0 0.00343331 0.003433310
## 330 66 1 1 1 1 -1 0.00343331 0.006000580
## 331 67 1 1 1 1 1 0.00343331 0.006000580
## 332 67 1 1 1 1 2 0.00343331 0.006000580
## 333 67 1 1 1 1 -2 0.00000000 0.000343331
## 334 67 1 1 1 1 0 0.00343331 0.003433310
## 335 67 1 1 1 1 -1 0.00343331 0.006000580
## 336 68 1 1 1 1 1 0.00343331 0.006000580
## 337 68 1 1 1 1 2 0.00343331 0.006000580
## 338 68 1 1 1 1 -2 0.00000000 0.000343331
## 339 68 1 1 1 1 0 0.00343331 0.003433310
## 340 68 1 1 1 1 -1 0.00343331 0.006000580
## 341 69 1 1 1 1 1 0.00343331 0.006000580
## 342 69 1 1 1 1 2 0.00343331 0.006000580
## 343 69 1 1 1 1 -2 0.00000000 0.000343331
## 344 69 1 1 1 1 0 0.00343331 0.003433310
## 345 69 1 1 1 1 -1 0.00343331 0.006000580
## 346 70 1 1 1 1 1 0.00343331 0.006000580
## 347 70 1 1 1 1 2 0.00343331 0.006000580
## 348 70 1 1 1 1 -2 0.00000000 0.000343331
## 349 70 1 1 1 1 0 0.00343331 0.003433310
## 350 70 1 1 1 1 -1 0.00343331 0.006000580
## 351 71 1 1 1 1 1 0.00343331 0.006000580
## 352 71 1 1 1 1 2 0.00343331 0.006000580
## 353 71 1 1 1 1 -2 0.00000000 0.000343331
## 354 71 1 1 1 1 0 0.00343331 0.003433310
## 355 71 1 1 1 1 -1 0.00343331 0.006000580
## 356 72 1 1 1 1 1 0.00343331 0.006000580
## 357 72 1 1 1 1 2 0.00343331 0.006000580
## 358 72 1 1 1 1 -2 0.00000000 0.000343331
## 359 72 1 1 1 1 0 0.00343331 0.003433310
## 360 72 1 1 1 1 -1 0.00343331 0.006000580
## 361 73 1 1 1 1 1 0.00343331 0.006000580
## 362 73 1 1 1 1 2 0.00343331 0.006000580
## 363 73 1 1 1 1 -2 0.00000000 0.000343331
## 364 73 1 1 1 1 0 0.00343331 0.003433310
## 365 73 1 1 1 1 -1 0.00343331 0.006000580
## 366 74 1 1 1 1 1 0.00343331 0.006000580
## 367 74 1 1 1 1 2 0.00343331 0.006000580
## 368 74 1 1 1 1 -2 0.00000000 0.000343331
## 369 74 1 1 1 1 0 0.00343331 0.003433310
## 370 74 1 1 1 1 -1 0.00343331 0.006000580
## 371 75 1 1 1 1 1 0.00343331 0.006000580
## 372 75 1 1 1 1 2 0.00343331 0.006000580
## 373 75 1 1 1 1 -2 0.00000000 0.000343331
## 374 75 1 1 1 1 0 0.00343331 0.003433310
## 375 75 1 1 1 1 -1 0.00343331 0.006000580
## 376 76 1 1 1 1 1 0.00343331 0.006000580
## 377 76 1 1 1 1 2 0.00343331 0.006000580
## 378 76 1 1 1 1 -2 0.00000000 0.000343331
## 379 76 1 1 1 1 0 0.00343331 0.003433310
## 380 76 1 1 1 1 -1 0.00343331 0.006000580
## 381 77 1 1 1 1 1 0.00343331 0.006000580
## 382 77 1 1 1 1 2 0.00343331 0.006000580
## 383 77 1 1 1 1 -2 0.00000000 0.000343331
## 384 77 1 1 1 1 0 0.00343331 0.003433310
## 385 77 1 1 1 1 -1 0.00343331 0.006000580
## 386 78 1 1 1 1 1 0.00343331 0.006000580
## 387 78 1 1 1 1 2 0.00343331 0.006000580
## 388 78 1 1 1 1 -2 0.00000000 0.000343331
## 389 78 1 1 1 1 0 0.00343331 0.003433310
## 390 78 1 1 1 1 -1 0.00343331 0.006000580
## 391 79 1 1 1 1 1 0.00343331 0.006000580
## 392 79 1 1 1 1 2 0.00343331 0.006000580
## 393 79 1 1 1 1 -2 0.00000000 0.000343331
## 394 79 1 1 1 1 0 0.00343331 0.003433310
## 395 79 1 1 1 1 -1 0.00343331 0.006000580
## 396 80 1 1 1 1 1 0.00343331 0.006000580
## 397 80 1 1 1 1 2 0.00343331 0.006000580
## 398 80 1 1 1 1 -2 0.00000000 0.000343331
## 399 80 1 1 1 1 0 0.00343331 0.003433310
## 400 80 1 1 1 1 -1 0.00343331 0.006000580
## 401 81 1 1 1 1 1 0.00343331 0.006000580
## 402 81 1 1 1 1 2 0.00343331 0.006000580
## 403 81 1 1 1 1 -2 0.00000000 0.000343331
## 404 81 1 1 1 1 0 0.00343331 0.003433310
## 405 81 1 1 1 1 -1 0.00343331 0.006000580
## 406 -9999 1 1 1 1 0 0.00343331 0.006000580
## 2 3 4 5 6 7 8
## 1 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 2 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 3 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 4 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 5 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 6 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 7 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 8 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 9 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 10 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 11 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 12 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 13 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 14 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 15 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 16 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 17 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 18 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 19 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 20 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 21 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 22 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 23 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 24 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 25 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 26 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 27 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 28 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 29 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 30 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 31 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 32 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 33 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 34 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 35 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 36 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 37 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 38 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 39 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 40 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 41 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 42 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 43 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 44 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 45 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 46 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 47 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 48 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 49 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 50 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 51 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 52 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 53 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 54 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 55 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 56 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 57 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 58 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 59 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 60 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 61 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 62 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 63 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 64 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 65 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 66 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 67 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 68 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 69 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 70 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 71 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 72 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 73 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 74 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 75 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 76 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 77 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 78 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 79 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 80 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 81 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 82 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 83 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 84 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 85 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 86 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 87 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 88 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 89 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 90 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 91 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 92 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 93 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 94 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 95 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 96 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 97 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 98 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 99 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 100 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 101 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 102 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 103 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 104 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 105 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 106 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 107 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 108 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 109 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 110 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 111 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 112 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 113 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 114 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 115 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 116 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 117 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 118 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 119 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 120 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 121 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 122 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 123 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 124 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 125 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 126 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 127 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 128 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 129 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 130 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 131 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 132 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 133 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 134 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 135 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 136 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 137 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 138 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 139 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 140 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 141 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 142 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 143 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 144 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 145 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 146 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 147 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 148 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 149 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 150 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 151 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 152 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 153 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 154 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 155 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 156 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 157 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 158 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 159 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 160 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 161 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 162 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 163 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 164 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 165 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 166 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 167 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 168 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 169 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 170 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 171 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 172 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 173 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 174 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 175 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 176 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 177 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 178 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 179 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 180 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 181 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 182 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 183 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 184 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 185 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 186 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 187 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 188 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 189 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 190 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 191 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 192 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 193 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 194 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 195 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 196 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 197 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 198 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 199 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 200 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 201 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 202 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 203 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 204 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 205 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 206 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 207 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 208 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 209 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 210 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 211 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 212 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 213 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 214 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 215 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 216 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 217 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 218 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 219 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 220 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 221 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 222 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 223 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 224 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 225 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 226 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 227 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 228 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 229 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 230 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 231 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 232 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 233 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 234 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 235 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 236 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 237 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 238 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 239 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 240 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 241 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 242 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 243 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 244 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 245 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 246 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 247 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 248 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 249 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 250 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 251 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 252 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 253 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 254 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 255 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 256 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 257 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 258 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 259 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 260 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 261 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 262 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 263 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 264 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 265 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 266 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 267 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 268 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 269 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 270 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 271 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 272 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 273 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 274 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 275 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 276 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 277 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 278 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 279 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 280 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 281 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 282 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 283 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 284 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 285 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 286 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 287 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 288 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 289 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 290 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 291 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 292 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 293 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 294 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 295 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 296 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 297 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 298 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 299 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 300 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 301 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 302 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 303 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 304 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 305 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 306 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 307 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 308 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 309 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 310 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 311 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 312 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 313 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 314 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 315 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 316 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 317 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 318 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 319 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 320 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 321 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 322 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 323 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 324 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 325 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 326 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 327 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 328 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 329 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 330 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 331 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 332 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 333 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 334 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 335 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 336 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 337 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 338 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 339 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 340 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 341 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 342 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 343 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 344 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 345 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 346 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 347 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 348 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 349 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 350 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 351 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 352 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 353 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 354 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 355 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 356 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 357 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 358 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 359 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 360 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 361 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 362 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 363 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 364 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 365 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 366 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 367 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 368 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 369 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 370 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 371 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 372 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 373 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 374 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 375 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 376 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 377 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 378 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 379 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 380 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 381 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 382 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 383 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 384 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 385 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 386 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 387 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 388 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 389 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 390 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 391 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 392 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 393 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 394 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 395 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 396 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 397 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 398 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 399 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 400 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 401 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 402 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 403 0.00605058 0.0268895 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 404 0.01210120 0.0298772 0.0512501 0.0729178 0.092935 0.110394 0.125053
## 405 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 406 0.02027710 0.0403544 0.0621976 0.0832091 0.102011 0.118072 0.131363
## 9 10 comment
## 1 0.142120 0.159599 #wt_flt_1
## 2 0.142120 0.159599 #wt_flt_2
## 3 0.137038 0.156718 #fecundity
## 4 0.137038 0.156718 #popwt_beg
## 5 0.142120 0.159599 #popwt_mid
## 6 0.142120 0.159599 #wt_flt_1
## 7 0.142120 0.159599 #wt_flt_2
## 8 0.137038 0.156718 #fecundity
## 9 0.137038 0.156718 #popwt_beg
## 10 0.142120 0.159599 #popwt_mid
## 11 0.142120 0.159599 #wt_flt_1
## 12 0.142120 0.159599 #wt_flt_2
## 13 0.137038 0.156718 #fecundity
## 14 0.137038 0.156718 #popwt_beg
## 15 0.142120 0.159599 #popwt_mid
## 16 0.142120 0.159599 #wt_flt_1
## 17 0.142120 0.159599 #wt_flt_2
## 18 0.137038 0.156718 #fecundity
## 19 0.137038 0.156718 #popwt_beg
## 20 0.142120 0.159599 #popwt_mid
## 21 0.142120 0.159599 #wt_flt_1
## 22 0.142120 0.159599 #wt_flt_2
## 23 0.137038 0.156718 #fecundity
## 24 0.137038 0.156718 #popwt_beg
## 25 0.142120 0.159599 #popwt_mid
## 26 0.142120 0.159599 #wt_flt_1
## 27 0.142120 0.159599 #wt_flt_2
## 28 0.137038 0.156718 #fecundity
## 29 0.137038 0.156718 #popwt_beg
## 30 0.142120 0.159599 #popwt_mid
## 31 0.142120 0.159599 #wt_flt_1
## 32 0.142120 0.159599 #wt_flt_2
## 33 0.137038 0.156718 #fecundity
## 34 0.137038 0.156718 #popwt_beg
## 35 0.142120 0.159599 #popwt_mid
## 36 0.142120 0.159599 #wt_flt_1
## 37 0.142120 0.159599 #wt_flt_2
## 38 0.137038 0.156718 #fecundity
## 39 0.137038 0.156718 #popwt_beg
## 40 0.142120 0.159599 #popwt_mid
## 41 0.142120 0.159599 #wt_flt_1
## 42 0.142120 0.159599 #wt_flt_2
## 43 0.137038 0.156718 #fecundity
## 44 0.137038 0.156718 #popwt_beg
## 45 0.142120 0.159599 #popwt_mid
## 46 0.142120 0.159599 #wt_flt_1
## 47 0.142120 0.159599 #wt_flt_2
## 48 0.137038 0.156718 #fecundity
## 49 0.137038 0.156718 #popwt_beg
## 50 0.142120 0.159599 #popwt_mid
## 51 0.142120 0.159599 #wt_flt_1
## 52 0.142120 0.159599 #wt_flt_2
## 53 0.137038 0.156718 #fecundity
## 54 0.137038 0.156718 #popwt_beg
## 55 0.142120 0.159599 #popwt_mid
## 56 0.142120 0.159599 #wt_flt_1
## 57 0.142120 0.159599 #wt_flt_2
## 58 0.137038 0.156718 #fecundity
## 59 0.137038 0.156718 #popwt_beg
## 60 0.142120 0.159599 #popwt_mid
## 61 0.142120 0.159599 #wt_flt_1
## 62 0.142120 0.159599 #wt_flt_2
## 63 0.137038 0.156718 #fecundity
## 64 0.137038 0.156718 #popwt_beg
## 65 0.142120 0.159599 #popwt_mid
## 66 0.142120 0.159599 #wt_flt_1
## 67 0.142120 0.159599 #wt_flt_2
## 68 0.137038 0.156718 #fecundity
## 69 0.137038 0.156718 #popwt_beg
## 70 0.142120 0.159599 #popwt_mid
## 71 0.142120 0.159599 #wt_flt_1
## 72 0.142120 0.159599 #wt_flt_2
## 73 0.137038 0.156718 #fecundity
## 74 0.137038 0.156718 #popwt_beg
## 75 0.142120 0.159599 #popwt_mid
## 76 0.142120 0.159599 #wt_flt_1
## 77 0.142120 0.159599 #wt_flt_2
## 78 0.137038 0.156718 #fecundity
## 79 0.137038 0.156718 #popwt_beg
## 80 0.142120 0.159599 #popwt_mid
## 81 0.142120 0.159599 #wt_flt_1
## 82 0.142120 0.159599 #wt_flt_2
## 83 0.137038 0.156718 #fecundity
## 84 0.137038 0.156718 #popwt_beg
## 85 0.142120 0.159599 #popwt_mid
## 86 0.142120 0.159599 #wt_flt_1
## 87 0.142120 0.159599 #wt_flt_2
## 88 0.137038 0.156718 #fecundity
## 89 0.137038 0.156718 #popwt_beg
## 90 0.142120 0.159599 #popwt_mid
## 91 0.142120 0.159599 #wt_flt_1
## 92 0.142120 0.159599 #wt_flt_2
## 93 0.137038 0.156718 #fecundity
## 94 0.137038 0.156718 #popwt_beg
## 95 0.142120 0.159599 #popwt_mid
## 96 0.142120 0.159599 #wt_flt_1
## 97 0.142120 0.159599 #wt_flt_2
## 98 0.137038 0.156718 #fecundity
## 99 0.137038 0.156718 #popwt_beg
## 100 0.142120 0.159599 #popwt_mid
## 101 0.142120 0.159599 #wt_flt_1
## 102 0.142120 0.159599 #wt_flt_2
## 103 0.137038 0.156718 #fecundity
## 104 0.137038 0.156718 #popwt_beg
## 105 0.142120 0.159599 #popwt_mid
## 106 0.142120 0.159599 #wt_flt_1
## 107 0.142120 0.159599 #wt_flt_2
## 108 0.137038 0.156718 #fecundity
## 109 0.137038 0.156718 #popwt_beg
## 110 0.142120 0.159599 #popwt_mid
## 111 0.142120 0.159599 #wt_flt_1
## 112 0.142120 0.159599 #wt_flt_2
## 113 0.137038 0.156718 #fecundity
## 114 0.137038 0.156718 #popwt_beg
## 115 0.142120 0.159599 #popwt_mid
## 116 0.142120 0.159599 #wt_flt_1
## 117 0.142120 0.159599 #wt_flt_2
## 118 0.137038 0.156718 #fecundity
## 119 0.137038 0.156718 #popwt_beg
## 120 0.142120 0.159599 #popwt_mid
## 121 0.142120 0.159599 #wt_flt_1
## 122 0.142120 0.159599 #wt_flt_2
## 123 0.137038 0.156718 #fecundity
## 124 0.137038 0.156718 #popwt_beg
## 125 0.142120 0.159599 #popwt_mid
## 126 0.142120 0.159599 #wt_flt_1
## 127 0.142120 0.159599 #wt_flt_2
## 128 0.137038 0.156718 #fecundity
## 129 0.137038 0.156718 #popwt_beg
## 130 0.142120 0.159599 #popwt_mid
## 131 0.142120 0.159599 #wt_flt_1
## 132 0.142120 0.159599 #wt_flt_2
## 133 0.137038 0.156718 #fecundity
## 134 0.137038 0.156718 #popwt_beg
## 135 0.142120 0.159599 #popwt_mid
## 136 0.142120 0.159599 #wt_flt_1
## 137 0.142120 0.159599 #wt_flt_2
## 138 0.137038 0.156718 #fecundity
## 139 0.137038 0.156718 #popwt_beg
## 140 0.142120 0.159599 #popwt_mid
## 141 0.142120 0.159599 #wt_flt_1
## 142 0.142120 0.159599 #wt_flt_2
## 143 0.137038 0.156718 #fecundity
## 144 0.137038 0.156718 #popwt_beg
## 145 0.142120 0.159599 #popwt_mid
## 146 0.142120 0.159599 #wt_flt_1
## 147 0.142120 0.159599 #wt_flt_2
## 148 0.137038 0.156718 #fecundity
## 149 0.137038 0.156718 #popwt_beg
## 150 0.142120 0.159599 #popwt_mid
## 151 0.142120 0.159599 #wt_flt_1
## 152 0.142120 0.159599 #wt_flt_2
## 153 0.137038 0.156718 #fecundity
## 154 0.137038 0.156718 #popwt_beg
## 155 0.142120 0.159599 #popwt_mid
## 156 0.142120 0.159599 #wt_flt_1
## 157 0.142120 0.159599 #wt_flt_2
## 158 0.137038 0.156718 #fecundity
## 159 0.137038 0.156718 #popwt_beg
## 160 0.142120 0.159599 #popwt_mid
## 161 0.142120 0.159599 #wt_flt_1
## 162 0.142120 0.159599 #wt_flt_2
## 163 0.137038 0.156718 #fecundity
## 164 0.137038 0.156718 #popwt_beg
## 165 0.142120 0.159599 #popwt_mid
## 166 0.142120 0.159599 #wt_flt_1
## 167 0.142120 0.159599 #wt_flt_2
## 168 0.137038 0.156718 #fecundity
## 169 0.137038 0.156718 #popwt_beg
## 170 0.142120 0.159599 #popwt_mid
## 171 0.142120 0.159599 #wt_flt_1
## 172 0.142120 0.159599 #wt_flt_2
## 173 0.137038 0.156718 #fecundity
## 174 0.137038 0.156718 #popwt_beg
## 175 0.142120 0.159599 #popwt_mid
## 176 0.142120 0.159599 #wt_flt_1
## 177 0.142120 0.159599 #wt_flt_2
## 178 0.137038 0.156718 #fecundity
## 179 0.137038 0.156718 #popwt_beg
## 180 0.142120 0.159599 #popwt_mid
## 181 0.142120 0.159599 #wt_flt_1
## 182 0.142120 0.159599 #wt_flt_2
## 183 0.137038 0.156718 #fecundity
## 184 0.137038 0.156718 #popwt_beg
## 185 0.142120 0.159599 #popwt_mid
## 186 0.142120 0.159599 #wt_flt_1
## 187 0.142120 0.159599 #wt_flt_2
## 188 0.137038 0.156718 #fecundity
## 189 0.137038 0.156718 #popwt_beg
## 190 0.142120 0.159599 #popwt_mid
## 191 0.142120 0.159599 #wt_flt_1
## 192 0.142120 0.159599 #wt_flt_2
## 193 0.137038 0.156718 #fecundity
## 194 0.137038 0.156718 #popwt_beg
## 195 0.142120 0.159599 #popwt_mid
## 196 0.142120 0.159599 #wt_flt_1
## 197 0.142120 0.159599 #wt_flt_2
## 198 0.137038 0.156718 #fecundity
## 199 0.137038 0.156718 #popwt_beg
## 200 0.142120 0.159599 #popwt_mid
## 201 0.142120 0.159599 #wt_flt_1
## 202 0.142120 0.159599 #wt_flt_2
## 203 0.137038 0.156718 #fecundity
## 204 0.137038 0.156718 #popwt_beg
## 205 0.142120 0.159599 #popwt_mid
## 206 0.142120 0.159599 #wt_flt_1
## 207 0.142120 0.159599 #wt_flt_2
## 208 0.137038 0.156718 #fecundity
## 209 0.137038 0.156718 #popwt_beg
## 210 0.142120 0.159599 #popwt_mid
## 211 0.142120 0.159599 #wt_flt_1
## 212 0.142120 0.159599 #wt_flt_2
## 213 0.137038 0.156718 #fecundity
## 214 0.137038 0.156718 #popwt_beg
## 215 0.142120 0.159599 #popwt_mid
## 216 0.142120 0.159599 #wt_flt_1
## 217 0.142120 0.159599 #wt_flt_2
## 218 0.137038 0.156718 #fecundity
## 219 0.137038 0.156718 #popwt_beg
## 220 0.142120 0.159599 #popwt_mid
## 221 0.142120 0.159599 #wt_flt_1
## 222 0.142120 0.159599 #wt_flt_2
## 223 0.137038 0.156718 #fecundity
## 224 0.137038 0.156718 #popwt_beg
## 225 0.142120 0.159599 #popwt_mid
## 226 0.142120 0.159599 #wt_flt_1
## 227 0.142120 0.159599 #wt_flt_2
## 228 0.137038 0.156718 #fecundity
## 229 0.137038 0.156718 #popwt_beg
## 230 0.142120 0.159599 #popwt_mid
## 231 0.142120 0.159599 #wt_flt_1
## 232 0.142120 0.159599 #wt_flt_2
## 233 0.137038 0.156718 #fecundity
## 234 0.137038 0.156718 #popwt_beg
## 235 0.142120 0.159599 #popwt_mid
## 236 0.142120 0.159599 #wt_flt_1
## 237 0.142120 0.159599 #wt_flt_2
## 238 0.137038 0.156718 #fecundity
## 239 0.137038 0.156718 #popwt_beg
## 240 0.142120 0.159599 #popwt_mid
## 241 0.142120 0.159599 #wt_flt_1
## 242 0.142120 0.159599 #wt_flt_2
## 243 0.137038 0.156718 #fecundity
## 244 0.137038 0.156718 #popwt_beg
## 245 0.142120 0.159599 #popwt_mid
## 246 0.142120 0.159599 #wt_flt_1
## 247 0.142120 0.159599 #wt_flt_2
## 248 0.137038 0.156718 #fecundity
## 249 0.137038 0.156718 #popwt_beg
## 250 0.142120 0.159599 #popwt_mid
## 251 0.142120 0.159599 #wt_flt_1
## 252 0.142120 0.159599 #wt_flt_2
## 253 0.137038 0.156718 #fecundity
## 254 0.137038 0.156718 #popwt_beg
## 255 0.142120 0.159599 #popwt_mid
## 256 0.142120 0.159599 #wt_flt_1
## 257 0.142120 0.159599 #wt_flt_2
## 258 0.137038 0.156718 #fecundity
## 259 0.137038 0.156718 #popwt_beg
## 260 0.142120 0.159599 #popwt_mid
## 261 0.142120 0.159599 #wt_flt_1
## 262 0.142120 0.159599 #wt_flt_2
## 263 0.137038 0.156718 #fecundity
## 264 0.137038 0.156718 #popwt_beg
## 265 0.142120 0.159599 #popwt_mid
## 266 0.142120 0.159599 #wt_flt_1
## 267 0.142120 0.159599 #wt_flt_2
## 268 0.137038 0.156718 #fecundity
## 269 0.137038 0.156718 #popwt_beg
## 270 0.142120 0.159599 #popwt_mid
## 271 0.142120 0.159599 #wt_flt_1
## 272 0.142120 0.159599 #wt_flt_2
## 273 0.137038 0.156718 #fecundity
## 274 0.137038 0.156718 #popwt_beg
## 275 0.142120 0.159599 #popwt_mid
## 276 0.142120 0.159599 #wt_flt_1
## 277 0.142120 0.159599 #wt_flt_2
## 278 0.137038 0.156718 #fecundity
## 279 0.137038 0.156718 #popwt_beg
## 280 0.142120 0.159599 #popwt_mid
## 281 0.142120 0.159599 #wt_flt_1
## 282 0.142120 0.159599 #wt_flt_2
## 283 0.137038 0.156718 #fecundity
## 284 0.137038 0.156718 #popwt_beg
## 285 0.142120 0.159599 #popwt_mid
## 286 0.142120 0.159599 #wt_flt_1
## 287 0.142120 0.159599 #wt_flt_2
## 288 0.137038 0.156718 #fecundity
## 289 0.137038 0.156718 #popwt_beg
## 290 0.142120 0.159599 #popwt_mid
## 291 0.142120 0.159599 #wt_flt_1
## 292 0.142120 0.159599 #wt_flt_2
## 293 0.137038 0.156718 #fecundity
## 294 0.137038 0.156718 #popwt_beg
## 295 0.142120 0.159599 #popwt_mid
## 296 0.142120 0.159599 #wt_flt_1
## 297 0.142120 0.159599 #wt_flt_2
## 298 0.137038 0.156718 #fecundity
## 299 0.137038 0.156718 #popwt_beg
## 300 0.142120 0.159599 #popwt_mid
## 301 0.142120 0.159599 #wt_flt_1
## 302 0.142120 0.159599 #wt_flt_2
## 303 0.137038 0.156718 #fecundity
## 304 0.137038 0.156718 #popwt_beg
## 305 0.142120 0.159599 #popwt_mid
## 306 0.142120 0.159599 #wt_flt_1
## 307 0.142120 0.159599 #wt_flt_2
## 308 0.137038 0.156718 #fecundity
## 309 0.137038 0.156718 #popwt_beg
## 310 0.142120 0.159599 #popwt_mid
## 311 0.142120 0.159599 #wt_flt_1
## 312 0.142120 0.159599 #wt_flt_2
## 313 0.137038 0.156718 #fecundity
## 314 0.137038 0.156718 #popwt_beg
## 315 0.142120 0.159599 #popwt_mid
## 316 0.142120 0.159599 #wt_flt_1
## 317 0.142120 0.159599 #wt_flt_2
## 318 0.137038 0.156718 #fecundity
## 319 0.137038 0.156718 #popwt_beg
## 320 0.142120 0.159599 #popwt_mid
## 321 0.142120 0.159599 #wt_flt_1
## 322 0.142120 0.159599 #wt_flt_2
## 323 0.137038 0.156718 #fecundity
## 324 0.137038 0.156718 #popwt_beg
## 325 0.142120 0.159599 #popwt_mid
## 326 0.142120 0.159599 #wt_flt_1
## 327 0.142120 0.159599 #wt_flt_2
## 328 0.137038 0.156718 #fecundity
## 329 0.137038 0.156718 #popwt_beg
## 330 0.142120 0.159599 #popwt_mid
## 331 0.142120 0.159599 #wt_flt_1
## 332 0.142120 0.159599 #wt_flt_2
## 333 0.137038 0.156718 #fecundity
## 334 0.137038 0.156718 #popwt_beg
## 335 0.142120 0.159599 #popwt_mid
## 336 0.142120 0.159599 #wt_flt_1
## 337 0.142120 0.159599 #wt_flt_2
## 338 0.137038 0.156718 #fecundity
## 339 0.137038 0.156718 #popwt_beg
## 340 0.142120 0.159599 #popwt_mid
## 341 0.142120 0.159599 #wt_flt_1
## 342 0.142120 0.159599 #wt_flt_2
## 343 0.137038 0.156718 #fecundity
## 344 0.137038 0.156718 #popwt_beg
## 345 0.142120 0.159599 #popwt_mid
## 346 0.142120 0.159599 #wt_flt_1
## 347 0.142120 0.159599 #wt_flt_2
## 348 0.137038 0.156718 #fecundity
## 349 0.137038 0.156718 #popwt_beg
## 350 0.142120 0.159599 #popwt_mid
## 351 0.142120 0.159599 #wt_flt_1
## 352 0.142120 0.159599 #wt_flt_2
## 353 0.137038 0.156718 #fecundity
## 354 0.137038 0.156718 #popwt_beg
## 355 0.142120 0.159599 #popwt_mid
## 356 0.142120 0.159599 #wt_flt_1
## 357 0.142120 0.159599 #wt_flt_2
## 358 0.137038 0.156718 #fecundity
## 359 0.137038 0.156718 #popwt_beg
## 360 0.142120 0.159599 #popwt_mid
## 361 0.142120 0.159599 #wt_flt_1
## 362 0.142120 0.159599 #wt_flt_2
## 363 0.137038 0.156718 #fecundity
## 364 0.137038 0.156718 #popwt_beg
## 365 0.142120 0.159599 #popwt_mid
## 366 0.142120 0.159599 #wt_flt_1
## 367 0.142120 0.159599 #wt_flt_2
## 368 0.137038 0.156718 #fecundity
## 369 0.137038 0.156718 #popwt_beg
## 370 0.142120 0.159599 #popwt_mid
## 371 0.142120 0.159599 #wt_flt_1
## 372 0.142120 0.159599 #wt_flt_2
## 373 0.137038 0.156718 #fecundity
## 374 0.137038 0.156718 #popwt_beg
## 375 0.142120 0.159599 #popwt_mid
## 376 0.142120 0.159599 #wt_flt_1
## 377 0.142120 0.159599 #wt_flt_2
## 378 0.137038 0.156718 #fecundity
## 379 0.137038 0.156718 #popwt_beg
## 380 0.142120 0.159599 #popwt_mid
## 381 0.142120 0.159599 #wt_flt_1
## 382 0.142120 0.159599 #wt_flt_2
## 383 0.137038 0.156718 #fecundity
## 384 0.137038 0.156718 #popwt_beg
## 385 0.142120 0.159599 #popwt_mid
## 386 0.142120 0.159599 #wt_flt_1
## 387 0.142120 0.159599 #wt_flt_2
## 388 0.137038 0.156718 #fecundity
## 389 0.137038 0.156718 #popwt_beg
## 390 0.142120 0.159599 #popwt_mid
## 391 0.142120 0.159599 #wt_flt_1
## 392 0.142120 0.159599 #wt_flt_2
## 393 0.137038 0.156718 #fecundity
## 394 0.137038 0.156718 #popwt_beg
## 395 0.142120 0.159599 #popwt_mid
## 396 0.142120 0.159599 #wt_flt_1
## 397 0.142120 0.159599 #wt_flt_2
## 398 0.137038 0.156718 #fecundity
## 399 0.137038 0.156718 #popwt_beg
## 400 0.142120 0.159599 #popwt_mid
## 401 0.142120 0.159599 #wt_flt_1
## 402 0.142120 0.159599 #wt_flt_2
## 403 0.137038 0.156718 #fecundity
## 404 0.137038 0.156718 #popwt_beg
## 405 0.142120 0.159599 #popwt_mid
## 406 0.142120 0.159599 #terminatorInitial skill assessment
This code is from the multispecies model skill assessment work using Hydra as an operating model and presented in 2015. Original file is hydra_sim_wrapper_modtest.R
#skill assessment metrics to evaluate model fit
# following from http://heuristically.wordpress.com/2013/07/12/calculate-rmse-and-mae-in-r-and-sas/
# where error = modeled-true or relative error = modeled/true
# Function that returns Root Mean Squared Error
rmse <- function(error)
{
sqrt(mean(error^2, na.rm=T))
}
# Function that returns Mean Absolute Error (SAME AS AAE FOR US)
aae <- function(error)
{
mean(abs(error), na.rm=T)
}
# Function that returns Mean Absolute Error (SAME AS AE FOR US)
ae <- function(error)
{
mean(error, na.rm=T)
}
# Function that returns Modeling Efficiency MEF
mef <- function(obs, error)
{
obsavg <- mean(obs, na.rm=T)
obserr <- obs - obsavg
(sum(obserr^2)-sum(error^2, na.rm=T))/sum(obserr^2)
}Reading in the truth–these files were generated and saved previously here. Recruitment and mortality are detailed here.
# generalized timesteps all models
runpar <- atlantisom::load_runprm(d.name, run.prm.file)
noutsteps <- runpar$tstop/runpar$outputstep
stepperyr <- if(runpar$outputstepunit=="days") 365/runpar$toutinc
midptyr <- round(median(seq(0,stepperyr)))
#Load functional groups
funct.groups <- load_fgs(dir=d.name,
file_fgs = functional.groups.file)
#Get just the names of active functional groups
funct.group.names <- funct.groups %>%
filter(IsTurnedOn == 1) %>%
select(Name) %>%
.$Name
trueB <- readRDS(file.path(d.name, paste0(scenario.name,"surveyBcensus.rds")))
trueB_ss <- trueB %>%
filter(species==species_ss)
trueN <- readRDS(file.path(d.name, paste0(scenario.name,"surveyNcensus.rds")))
trueN_ss <- trueN %>%
filter(species==species_ss)
Natage <- readRDS(file.path(d.name, paste0(scenario.name, "natage_census_sard_hake.rds")))
Natage_ss <- Natage %>%
filter(species==species_ss)
lengthwt_census_ss <- readRDS(file.path(d.name, paste0(scenario.name, "lengthwt_census_sard_hake.rds")))
len_ss <- lengthwt_census_ss$natlength %>%
filter(species==species_ss)
YOY <- atlantisom::load_yoy(d.name, paste0("output", scenario.name, "YOY.txt"))
# load biolprm in some initialize file?
biol <- atlantisom::load_biolprm(d.name, biol.prm.file)## Warning in `[<-.data.frame`(`*tmp*`, , 3:12, value = structure(list(`1` =
## structure(c(469L, : provided 13 variables to replace 10 variables# get code matching species name to split YOY file
code_ss <- funct.groups$Code[which(funct.groups$Name == species_ss)]
# cut to a single species in YOY file
YOY_ss <- YOY %>%
select(Time, paste0(code_ss, ".0"))
# mg carbon converted to wet weight in tonnes
k_wetdry <- biol$kgw2d / 1000000000
# WARNING only works for CCA because YOY.txt rows 2:end are already in numbers
# and we are merging out the incorrect and irrelevant YOY row 1 (Time=0)
# also hardcoded for sardine example
recnums <- YOY_ss %>%
mutate(yr = as.integer(round(YOY_ss$Time)/365)) %>%
mutate(recnums = SAR.0/k_wetdry/biol$redfieldcn) %>%
filter(yr>0)
file.mort <- file.path(d.name, paste0("output", scenario.name, "Mort.txt"))
mortish <- read.table(file.mort, header = TRUE)
relF_ss <- mortish %>%
select(Time, relF = paste0(code_ss, ".F"))replist$timeseries## Area Yr Era Seas Bio_all Bio_smry SpawnBio Recruit_0 SpawnBio_GP:1
## 948 1 -1 VIRG 1 780801 780801 670770 7515010 670770
## 949 1 0 INIT 1 780801 780801 670770 7515010 670770
## 950 1 1 TIME 1 780801 780801 670770 0 670770
## 951 1 2 TIME 1 780801 780801 670770 7515010 670770
## 952 1 3 TIME 1 780801 780801 670770 7515010 670770
## 953 1 4 TIME 1 780801 780801 670770 7515010 670770
## 954 1 5 TIME 1 780801 780801 670770 7515010 670770
## 955 1 6 TIME 1 780801 780801 670770 7515010 670770
## 956 1 7 TIME 1 780801 780801 670770 7515010 670770
## 957 1 8 TIME 1 780801 780801 670770 7515010 670770
## 958 1 9 TIME 1 780801 780801 670770 7515010 670770
## 959 1 10 TIME 1 780801 780801 670770 7515010 670770
## 960 1 11 TIME 1 780801 780801 670770 7515010 670770
## 961 1 12 TIME 1 780801 780801 670770 7515010 670770
## 962 1 13 TIME 1 780801 780801 670770 7515010 670770
## 963 1 14 TIME 1 780801 780801 670770 7515010 670770
## 964 1 15 TIME 1 780801 780801 670770 7515010 670770
## 965 1 16 TIME 1 780801 780801 670770 7515010 670770
## 966 1 17 TIME 1 780801 780801 670770 7515010 670770
## 967 1 18 TIME 1 780801 780801 670770 7515010 670770
## 968 1 19 TIME 1 780801 780801 670770 7515010 670770
## 969 1 20 TIME 1 780801 780801 670770 7515010 670770
## 970 1 21 TIME 1 780801 780801 670770 7515010 670770
## 971 1 22 TIME 1 780801 780801 670770 7515010 670770
## 972 1 23 TIME 1 780801 780801 670770 7515010 670770
## 973 1 24 TIME 1 780801 780801 670770 7515010 670770
## 974 1 25 TIME 1 780801 780801 670770 7515010 670770
## 975 1 26 TIME 1 780801 780801 670770 7515010 670770
## 976 1 27 TIME 1 780801 780801 670770 7515010 670770
## 977 1 28 TIME 1 780801 780801 670770 7515010 670770
## 978 1 29 TIME 1 780801 780801 670770 7515010 670770
## 979 1 30 TIME 1 780801 780801 670770 7515010 670770
## 980 1 31 TIME 1 756802 756802 665659 524983 665659
## 981 1 32 TIME 1 690362 690362 623378 2827400 623378
## 982 1 33 TIME 1 604603 604603 539935 3988500 539935
## 983 1 34 TIME 1 519518 519518 460652 1469640 460652
## 984 1 35 TIME 1 445503 445503 396419 2466290 396419
## 985 1 36 TIME 1 381139 381139 334165 3049790 334165
## 986 1 37 TIME 1 325942 325942 287224 395109 287224
## 987 1 38 TIME 1 280660 280660 249930 1729700 249930
## 988 1 39 TIME 1 252762 252762 209170 6266950 209170
## 989 1 40 TIME 1 287801 287801 197295 15581500 197295
## 990 1 41 TIME 1 366777 366777 262996 4411850 262996
## 991 1 42 TIME 1 470148 470148 381126 7458430 381126
## 992 1 43 TIME 1 639361 639361 462595 35057900 462595
## 993 1 44 TIME 1 806949 806949 609345 1013420 609345
## 994 1 45 TIME 1 960723 960723 835664 5072260 835664
## 995 1 46 TIME 1 1040790 1040790 897299 18964400 897299
## 996 1 47 TIME 1 1078800 1078800 901135 10013800 901135
## 997 1 48 TIME 1 1105120 1105120 952573 6411250 952573
## 998 1 49 TIME 1 1092950 1092950 963305 7258590 963305
## 999 1 50 TIME 1 1066070 1066070 917430 14845700 917430
## 1000 1 51 TIME 1 1041010 1041010 884436 7938280 884436
## 1001 1 52 TIME 1 1029420 1029420 894723 6379850 894723
## 1002 1 53 TIME 1 1001240 1001240 866814 11899000 866814
## 1003 1 54 TIME 1 1038180 1038180 839932 26633500 839932
## 1004 1 55 TIME 1 1702200 1702200 949069 180363000 949069
## 1005 1 56 TIME 1 2595370 2595370 1646100 22285800 1646100
## 1006 1 57 TIME 1 3734760 3734760 3034160 51408000 3034160
## 1007 1 58 TIME 1 4350270 4350270 3802420 15069600 3802420
## 1008 1 59 TIME 1 4514660 4514660 4038440 21672700 4038440
## 1009 1 60 TIME 1 4377480 4377480 3875360 40500100 3875360
## 1010 1 61 TIME 1 4062990 4062990 3590330 9674330 3590330
## 1011 1 62 TIME 1 3735730 3735730 3329950 25465600 3329950
## 1012 1 63 TIME 1 3360590 3360590 2983620 16571300 2983620
## 1013 1 64 TIME 1 3045410 3045410 2713740 12931200 2713740
## 1014 1 65 TIME 1 2646410 2646410 2385890 5931850 2385890
## 1015 1 66 TIME 1 2300840 2300840 2070260 14242400 2070260
## 1016 1 67 TIME 1 2248440 2248440 1787800 87193300 1787800
## 1017 1 68 TIME 1 2341650 2341650 1832330 3100600 1832330
## 1018 1 69 TIME 1 2508050 2508050 2217860 1860220 2217860
## 1019 1 70 TIME 1 2399640 2399640 2223730 2166430 2223730
## 1020 1 71 TIME 1 2104730 2104730 1941620 3918870 1941620
## 1021 1 72 TIME 1 1760200 1760200 1591560 11511400 1591560
## 1022 1 73 TIME 1 1590550 1590550 1310150 43677100 1310150
## 1023 1 74 TIME 1 1524330 1524330 1242650 622436 1242650
## 1024 1 75 TIME 1 1514370 1514370 1353990 1162800 1353990
## 1025 1 76 TIME 1 1412320 1412320 1301910 3769370 1301910
## 1026 1 77 TIME 1 1215100 1215100 1097240 6691000 1097240
## 1027 1 78 TIME 1 1019400 1019400 910811 3561130 910811
## 1028 1 79 TIME 1 874551 874551 773544 6746820 773544
## 1029 1 80 TIME 1 1107940 1107940 697751 106816000 697751
## 1030 1 81 FORE 1 1529970 1529970 995807 7518830 995807
## SmryBio_SX:1_GP:1 SmryNum_SX:1_GP:1 sel(B):_1 dead(B):_1 retain(B):_1
## 948 780801 22794900 0.00 0.00 0.00
## 949 780801 22794900 0.00 0.00 0.00
## 950 780801 22794900 0.00 0.00 0.00
## 951 780801 22794900 0.00 0.00 0.00
## 952 780801 22794900 0.00 0.00 0.00
## 953 780801 22794900 0.00 0.00 0.00
## 954 780801 22794900 0.00 0.00 0.00
## 955 780801 22794900 0.00 0.00 0.00
## 956 780801 22794900 0.00 0.00 0.00
## 957 780801 22794900 0.00 0.00 0.00
## 958 780801 22794900 0.00 0.00 0.00
## 959 780801 22794900 0.00 0.00 0.00
## 960 780801 22794900 0.00 0.00 0.00
## 961 780801 22794900 0.00 0.00 0.00
## 962 780801 22794900 0.00 0.00 0.00
## 963 780801 22794900 0.00 0.00 0.00
## 964 780801 22794900 0.00 0.00 0.00
## 965 780801 22794900 0.00 0.00 0.00
## 966 780801 22794900 0.00 0.00 0.00
## 967 780801 22794900 0.00 0.00 0.00
## 968 780801 22794900 0.00 0.00 0.00
## 969 780801 22794900 0.00 0.00 0.00
## 970 780801 22794900 0.00 0.00 0.00
## 971 780801 22794900 0.00 0.00 0.00
## 972 780801 22794900 0.00 0.00 0.00
## 973 780801 22794900 0.00 0.00 0.00
## 974 780801 22794900 0.00 0.00 0.00
## 975 780801 22794900 0.00 0.00 0.00
## 976 780801 22794900 0.00 0.00 0.00
## 977 780801 22794900 0.00 0.00 0.00
## 978 780801 22794900 0.00 0.00 0.00
## 979 780801 22794900 0.00 0.00 0.00
## 980 756802 15804800 18726.50 18726.50 18726.50
## 981 690362 13161800 18808.70 18808.70 18808.70
## 982 604603 12608000 16706.50 16706.50 16706.50
## 983 519518 9745830 15619.40 15619.40 15619.40
## 984 445503 8819750 16786.90 16786.90 16786.90
## 985 381139 8785770 14845.90 14845.90 14845.90
## 986 325942 6119520 12645.70 12645.70 12645.70
## 987 280660 5674330 14587.40 14587.40 14587.40
## 988 252762 9919670 12951.50 12951.50 12951.50
## 989 287801 22084600 7868.48 7868.48 7868.48
## 990 366777 19072000 5729.07 5729.07 5729.07
## 991 470148 20100200 8190.15 8190.15 8190.15
## 992 639361 48383500 3143.99 3143.99 3143.99
## 993 806949 33393200 2645.96 2645.96 2645.96
## 994 960723 27392400 4790.62 4790.62 4790.62
## 995 1040790 37251800 6583.58 6583.58 6583.58
## 996 1078800 34899900 5944.72 5944.72 5944.72
## 997 1105120 29710700 7140.06 7140.06 7140.06
## 998 1092950 27070200 8692.41 8692.41 8692.41
## 999 1066070 32879800 9203.40 9203.40 9203.40
## 1000 1041010 29864000 7616.38 7616.38 7616.38
## 1001 1029420 26285900 8152.40 8152.40 8152.40
## 1002 1001240 29408000 8704.40 8704.40 8704.40
## 1003 1038180 46235800 8039.88 8039.88 8039.88
## 1004 1702200 211241000 44707.60 44707.60 44707.60
## 1005 2595370 163057000 11197.40 11197.40 11197.40
## 1006 3734760 160363000 6753.33 6753.33 6753.33
## 1007 4350270 122436000 5095.49 5095.49 5095.49
## 1008 4514660 103661000 4722.43 4722.43 4722.43
## 1009 4377480 109926000 9981.42 9981.42 9981.42
## 1010 4062990 83245000 11827.10 11827.10 11827.10
## 1011 3735730 81115400 23497.30 23497.30 23497.30
## 1012 3360590 70688500 22660.90 22660.90 22660.90
## 1013 3045410 60051000 20145.30 20145.30 20145.30
## 1014 2646410 45958700 18218.10 18218.10 18218.10
## 1015 2300840 44850000 14376.30 14376.30 14376.30
## 1016 2248440 117115000 13144.20 13144.20 13144.20
## 1017 2341650 81459700 5980.30 5980.30 5980.30
## 1018 2508050 56336800 13878.20 13878.20 13878.20
## 1019 2399640 39727400 24180.50 24180.50 24180.50
## 1020 2104730 30281900 32440.20 32440.20 32440.20
## 1021 1760200 31512400 25004.50 25004.50 25004.50
## 1022 1590550 64585300 18469.30 18469.30 18469.30
## 1023 1524330 43720100 11794.40 11794.40 11794.40
## 1024 1514370 30249400 16533.40 16533.40 16533.40
## 1025 1412320 23827500 19437.50 19437.50 19437.50
## 1026 1215100 22463600 26175.80 26175.80 26175.80
## 1027 1019400 18367200 15482.20 15482.20 15482.20
## 1028 874551 18894100 20013.10 20013.10 20013.10
## 1029 1107940 119274000 0.00 0.00 0.00
## 1030 1529970 87470300 0.00 0.00 0.00
## sel(N):_1 dead(N):_1 retain(N):_1 obs_cat:_1 F:_1 SSB_vir_LH
## 948 0 0 0 0 0.00000000 717012
## 949 0 0 0 0 0.00000000 717012
## 950 0 0 0 0 0.00000000 717012
## 951 0 0 0 0 0.00000000 717012
## 952 0 0 0 0 0.00000000 717012
## 953 0 0 0 0 0.00000000 717012
## 954 0 0 0 0 0.00000000 717012
## 955 0 0 0 0 0.00000000 717012
## 956 0 0 0 0 0.00000000 717012
## 957 0 0 0 0 0.00000000 717012
## 958 0 0 0 0 0.00000000 717012
## 959 0 0 0 0 0.00000000 717012
## 960 0 0 0 0 0.00000000 717012
## 961 0 0 0 0 0.00000000 717012
## 962 0 0 0 0 0.00000000 717012
## 963 0 0 0 0 0.00000000 717012
## 964 0 0 0 0 0.00000000 717012
## 965 0 0 0 0 0.00000000 717012
## 966 0 0 0 0 0.00000000 717012
## 967 0 0 0 0 0.00000000 717012
## 968 0 0 0 0 0.00000000 717012
## 969 0 0 0 0 0.00000000 717012
## 970 0 0 0 0 0.00000000 717012
## 971 0 0 0 0 0.00000000 717012
## 972 0 0 0 0 0.00000000 717012
## 973 0 0 0 0 0.00000000 717012
## 974 0 0 0 0 0.00000000 717012
## 975 0 0 0 0 0.00000000 717012
## 976 0 0 0 0 0.00000000 717012
## 977 0 0 0 0 0.00000000 717012
## 978 0 0 0 0 0.00000000 717012
## 979 0 0 0 0 0.00000000 717012
## 980 319865 319865 319865 319865 0.02578960 714612
## 981 249922 249922 249922 249922 0.02976120 669662
## 982 215679 215679 215679 215679 0.03086530 580126
## 983 220781 220781 220781 220781 0.03297180 495119
## 984 216810 216810 216810 216810 0.04228980 426739
## 985 203003 203003 203003 203003 0.04395030 359820
## 986 193817 193817 193817 193817 0.04203750 309181
## 987 185898 185898 185898 185898 0.05879990 269785
## 988 180176 180176 180176 180176 0.06177890 225881
## 989 176713 176713 176713 176713 0.03370460 212053
## 990 175377 175377 175377 175377 0.01471820 281812
## 991 182145 182145 182145 182145 0.01766680 408593
## 992 64694 64694 64694 64694 0.00592886 494962
## 993 78720 78720 78720 78720 0.00297530 651675
## 994 91349 91349 91349 91349 0.00497869 894014
## 995 103954 103954 103954 103954 0.00693209 960260
## 996 116472 116472 116472 116472 0.00572272 964173
## 997 127937 127937 127937 127937 0.00670227 1019380
## 998 137312 137312 137312 137312 0.00846070 1031160
## 999 140571 140571 140571 140571 0.00952496 982226
## 1000 138208 138208 138208 138208 0.00771030 946620
## 1001 136455 136455 136455 136455 0.00837186 957736
## 1002 135849 135849 135849 135849 0.00947842 928029
## 1003 140378 140378 140378 140378 0.00875973 899133
## 1004 1018980 1018980 1018980 1018980 0.04104210 1021030
## 1005 424892 424892 424892 424892 0.00370006 1760660
## 1006 157952 157952 157952 157952 0.00176277 3244280
## 1007 101970 101970 101970 101970 0.00115664 4065330
## 1008 73612 73612 73612 73612 0.00109124 4317630
## 1009 141543 141543 141543 141543 0.00248166 4144230
## 1010 185664 185664 185664 185664 0.00307760 3839810
## 1011 315029 315029 315029 315029 0.00689851 3563600
## 1012 324810 324810 324810 324810 0.00733011 3193200
## 1013 278613 278613 278613 278613 0.00721515 2904300
## 1014 245273 245273 245273 245273 0.00747680 2553540
## 1015 174959 174959 174959 174959 0.00696702 2215540
## 1016 178687 178687 178687 178687 0.00729344 1913340
## 1017 156851 156851 156851 156851 0.00244962 1959450
## 1018 249209 249209 249209 249209 0.00556691 2372960
## 1019 328456 328456 328456 328456 0.01066620 2381260
## 1020 366142 366142 366142 366142 0.01699410 2081360
## 1021 264745 264745 264745 264745 0.01620270 1705870
## 1022 240040 240040 240040 240040 0.01407260 1403720
## 1023 270350 270350 270350 270350 0.00769022 1330020
## 1024 268981 268981 268981 268981 0.01128080 1450060
## 1025 245319 245319 245319 245319 0.01494860 1395130
## 1026 309845 309845 309845 309845 0.02414220 1177600
## 1027 202432 202432 202432 202432 0.01676890 976322
## 1028 254846 254846 254846 254846 0.02581230 830425
## 1029 0 0 0 0 0.00000000 745853
## 1030 0 0 0 NA 0.00000000 1064460
## ABC_buffer
## 948 NA
## 949 NA
## 950 NA
## 951 NA
## 952 NA
## 953 NA
## 954 NA
## 955 NA
## 956 NA
## 957 NA
## 958 NA
## 959 NA
## 960 NA
## 961 NA
## 962 NA
## 963 NA
## 964 NA
## 965 NA
## 966 NA
## 967 NA
## 968 NA
## 969 NA
## 970 NA
## 971 NA
## 972 NA
## 973 NA
## 974 NA
## 975 NA
## 976 NA
## 977 NA
## 978 NA
## 979 NA
## 980 NA
## 981 NA
## 982 NA
## 983 NA
## 984 NA
## 985 NA
## 986 NA
## 987 NA
## 988 NA
## 989 NA
## 990 NA
## 991 NA
## 992 NA
## 993 NA
## 994 NA
## 995 NA
## 996 NA
## 997 NA
## 998 NA
## 999 NA
## 1000 NA
## 1001 NA
## 1002 NA
## 1003 NA
## 1004 NA
## 1005 NA
## 1006 NA
## 1007 NA
## 1008 NA
## 1009 NA
## 1010 NA
## 1011 NA
## 1012 NA
## 1013 NA
## 1014 NA
## 1015 NA
## 1016 NA
## 1017 NA
## 1018 NA
## 1019 NA
## 1020 NA
## 1021 NA
## 1022 NA
## 1023 NA
## 1024 NA
## 1025 NA
## 1026 NA
## 1027 NA
## 1028 NA
## 1029 NA
## 1030 0names(replist$timeseries)[19]<-"Fmort"How does SS3 estimated biomass compare with true?
plotB <-ggplot() +
geom_line(data=trueB_ss, aes(x=time/stepperyr,y=atoutput, color="True B"),
alpha = 10/10) +
geom_point(data=replist$timeseries, aes(x=Yr,y=Bio_all, color="SS3 Est B"),
alpha = 10/10) +
theme_tufte() +
theme(legend.position = "top") +
labs(colour=scenario.name)
plotB +
facet_wrap(~species, scales="free") 
What about recruitment?
plotR <-ggplot() +
geom_line(data=recnums, aes(x=yr,y=recnums, color="True R"),
alpha = 10/10) +
geom_point(data=replist$timeseries, aes(x=Yr,y=(Recruit_0*1000), color="SS3 Est R"),
alpha = 10/10) +
theme_tufte() +
theme(legend.position = "top") +
labs(colour=scenario.name)
plotR 
And F (understanding that F in the mort.txt file from Atlantis is relative, just looking for patterns).
plotF <-ggplot() +
geom_line(data=relF_ss, aes(x=Time/365, y=relF, color="True F"),
alpha = 10/10) +
geom_point(data=replist$timeseries, aes(x=Yr,y=Fmort, color="SS3 Est F"),
alpha = 10/10) +
theme_tufte() +
theme(legend.position = "top") +
labs(colour=scenario.name)
plotF 