add missing dimension for evalp output shape

1548ec44 · Thibault Hallouin · e0722aa7 · 1548ec44
Commit 1548ec44 authored 2 years ago by Thibault Hallouin
Hide whitespace changes
Inline Side-by-side

Showing

with 1 addition and 1 deletion
+1 -1
--- a/evalhyd-python.ipynb
+++ b/evalhyd-python.ipynb
@@ -94,7 +94,7 @@
    "import evalhyd\n",
    "\n",
    "# compute Brier score\n",
-    "#   > shape: {sites: 1, lead times: 1, subsets: 1, samples: 1}\n",
+    "#   > shape: {sites: 1, lead times: 1, subsets: 1, samples: 1, thresholds: 2}\n",
    "evalhyd.evalp(obs, prd, [\"BS\"], thr, events=\"high\")"
   ],
   "metadata": {

 %% Cell type:markdown id: tags:

 ### Deterministic evaluation

 %% Cell type:code id: tags:

 ``` python
 import numpy

 # define streamflow observations
 #   > shape: {1, time: 4}
 obs = numpy.array(
    [[4.7, 4.3, 5.5, 2.7]]
 )

 # define streamflow predictions
 #   > shape: {series: 1, time:4}
 prd = numpy.array(
    [[5.3, 4.2, 5.7, 2.3]]
 )
 ```

 %% Cell type:code id: tags:

 ``` python
 import evalhyd

 # compute Nash-Sutcliffe efficiency
 #   > shape: {series: 1, subsets: 1, samples: 1}
 evalhyd.evald(obs, prd, ["NSE"])
 ```

 %% Cell type:markdown id: tags:

 ### Probabilistic evaluation

 %% Cell type:code id: tags:

 ``` python
 import numpy

 # define streamflow observations
 #   > shape: {sites: 1, time: 5}
 obs = numpy.array(
    [[4.7, 4.3, 5.5, 2.7, 4.1]]
 )

 # define streamflow predictions
 #   > shape: {sites: 1, lead times: 1, members: 3, time: 5}
 prd = numpy.array(
    [[[[5.3, 4.2, 5.7, 2.3, 3.1],
       [4.3, 4.2, 4.7, 4.3, 3.3],
       [5.3, 5.2, 5.7, 2.3, 3.9]]]]
 )

 # define streamflow thresholds
 #   > shape: {sites: 1, thresholds: 2}
 thr = numpy.array([[4., 5.]])
 ```

 %% Cell type:code id: tags:

 ``` python
 import evalhyd

 # compute Brier score
-#   > shape: {sites: 1, lead times: 1, subsets: 1, samples: 1}
+#   > shape: {sites: 1, lead times: 1, subsets: 1, samples: 1, thresholds: 2}
 evalhyd.evalp(obs, prd, ["BS"], thr, events="high")
 ```