MODE CSI calculations result in spurious results

[bikegeek]

Describe the Problem

There is an issue with the way the MODE CSI statistics are being calculated (inconsistent with what is documented, and results are not consistent). This was discovered and reported in the METplus Github Discussion:

• https://github.com/dtcenter/METplus/discussions/2496

This affects the R script calculations in METviewer

• https://github.com/dtcenter/METviewer/issues/516

Expected Behavior

Environment

Describe your runtime environment: 1. Machine: (e.g. HPC name, Linux Workstation, Mac Laptop) 2. OS: (e.g. RedHat Linux, MacOS) 3. Software version number(s)

To Reproduce

Describe the steps to reproduce the behavior: 1. Go to '...' 2. Click on '....' 3. Scroll down to '....' 4. See error Post relevant sample data following these instructions: https://dtcenter.org/community-code/model-evaluation-tools-met/met-help-desk#ftp

Relevant Deadlines

List relevant project deadlines here or state NONE.

Funding Source

Define the source of funding and account keys here or state NONE.

Define the Metadata

Assignee

• [ ] Select engineer(s) or no engineer required
• [ ] Select scientist(s) or no scientist required

Labels

• [x] Select component(s)
• [ ] Select priority
• [x] Select requestor(s)

Projects and Milestone

• [ ] Select Organization level Project for support of the current coordinated release
• [ ] Select Repository level Project for development toward the next official release or add alert: NEED CYCLE ASSIGNMENT label
• [ ] Select Milestone as the next bugfix version

Define Related Issue(s)

Consider the impact to the other METplus components.

• [ ] METplus, MET, METdataio, METviewer, METexpress, METcalcpy, METplotpy

Bugfix Checklist

See the METplus Workflow for details.

• [ ] Complete the issue definition above, including the Time Estimate and Funding Source.
• [ ] Fork this repository or create a branch of main_\. Branch name: bugfix_<Issue Number>_main_<Version>_<Description>
• [ ] Fix the bug and test your changes.
• [ ] Add/update log messages for easier debugging.
• [ ] Add/update unit tests.
• [ ] Add/update documentation.
• [ ] Add any new Python packages to the METplus Components Python Requirements table.
• [ ] Push local changes to GitHub.
• [ ] Submit a pull request to merge into main_\. Pull request: bugfix <Issue Number> main_<Version> <Description>
• [ ] Define the pull request metadata, as permissions allow. Select: Reviewer(s) and Development issues Select: Organization level software support Project for the current coordinated release Select: Milestone as the next bugfix version
• [ ] Iterate until the reviewer(s) accept and merge your changes.
• [ ] Delete your fork or branch.
• [ ] Complete the steps above to fix the bug on the develop branch. Branch name: bugfix_<Issue Number>_develop_<Description> Pull request: bugfix <Issue Number> develop <Description> Select: Reviewer(s) and Development issues Select: Repository level development cycle Project for the next official release Select: Milestone as the next official version
• [ ] Close this issue.

[bikegeek]

• the corresponding calcOBJACSI code (in METviewer/R_work/include/agg_statbootstrap.R) is found in the METcalcpy/util/mode_2d_arearat_statistics.py module_
• the corresponding calcOBJCSI code (n METviewer/R_work/include/agg_stat_bootstrap.R) is found in the _METcalcpy/util/mode_2d_ratiostatistics.py module

[bikegeek]

I removed the Y1 and Y2 Bounds values in METviewer

OBJCSI plot with updated equation:

OBJACSI plot with updated equation (no change from plot using the original equation):