|
22 | 22 | # CURRENT SURVEY YEAR BEING TESTED: 2019
|
23 | 23 | ####################################################################################################
|
24 | 24 | # Define current survey year
|
25 |
| -SURVEY_YEAR = 2017 |
| 25 | +SURVEY_YEAR = 2013 |
26 | 26 |
|
27 | 27 | # Initialization configuration
|
28 | 28 | init_config_path = f"C:/Users/Brandyn/Documents/GitHub/echopop/config_files/initialization_\
|
|
65 | 65 | for excl in parameters["exclude_age1"]:
|
66 | 66 | # ---- Transect analysis
|
67 | 67 | survey.transect_analysis(exclude_age1=excl, stratum=stratum, verbose=False)
|
| 68 | + print( |
| 69 | + f"{SURVEY_YEAR} transect analysis [{counter}/8] success!" |
| 70 | + ) |
68 | 71 | # ---- Stratified analysis (transect analysis)
|
69 | 72 | survey.stratified_analysis(bootstrap_ci_method=parameters["bootstrap_ci_method"],
|
70 | 73 | transect_replicates=parameters["transect_replicates"],
|
71 | 74 | verbose=False)
|
| 75 | + print( |
| 76 | + f"{SURVEY_YEAR} stratified analysis (transect) [{counter}/8] success!" |
| 77 | + ) |
72 | 78 | # ---- Fit variogram
|
73 | 79 | survey.fit_variogram(verbose=False)
|
| 80 | + print( |
| 81 | + f"{SURVEY_YEAR} variogram fitting [{counter}/8] success!" |
| 82 | + ) |
74 | 83 | # ---- Iterate across different extrapolation schema
|
75 | 84 | for extrap in parameters["extrapolate"]:
|
76 | 85 | # ---- Kriging analysis (no variogram fitting)
|
77 | 86 | survey.kriging_analysis(extrapolate=extrap, verbose=False)
|
| 87 | + print( |
| 88 | + f"{SURVEY_YEAR} kriging analysis (default variogram) [{counter}/8] success!" |
| 89 | + ) |
78 | 90 | # ---- Apply best-fit variogram
|
79 | 91 | survey.kriging_analysis(best_fit_variogram=True, extrapolate=extrap, verbose=False)
|
| 92 | + print( |
| 93 | + f"{SURVEY_YEAR} kriging analysis (best-fit variogram) [{counter}/8] success!" |
| 94 | + ) |
80 | 95 | # ---- Stratified analysis (kriging analysis)
|
81 | 96 | survey.stratified_analysis(dataset="kriging",
|
82 | 97 | bootstrap_ci_method=parameters["bootstrap_ci_method"],
|
83 | 98 | transect_replicates=parameters["transect_replicates"],
|
84 | 99 | verbose=False)
|
| 100 | + print( |
| 101 | + f"{SURVEY_YEAR} stratified analysis (kriging) [{counter}/8] success!" |
| 102 | + ) |
85 | 103 | # ---- Test reports
|
86 | 104 | survey.generate_reports(reports=["aged_length_haul_counts",
|
87 | 105 | "kriging_input",
|
|
0 commit comments