降水預報R包的驗證

Question

我觀察並模擬了降水。 我想使用verification包來驗證模型的技能。

為此，我必須將pred值轉換為概率預測，並將obs轉換為二進制。 例如， precipitation amount>2mm 。 那就是使用以下腳本

A<- verify(obs, pred, frcst.type = "prob", obs.type = "binary")

summary(A)

我的問題是如何使用上述腳本將pred轉換為概率，將obs轉換為二進制

我的數據長度是122。這是數據。

> Obs

 [1]   1.9  3.7  3.4  1.1  1.0  1.0  0.6  2.8  1.4  1.9  1.6  1.3  3.6  3.6  3.9

 [16]   3.4  2.7  5.9  3.8  1.8  2.2  2.5  2.6  2.4  2.6  1.5  1.0 11.2  2.1  0.6

 [31]  1.8  3.8  4.0  1.3  2.4  2.9  3.9  2.7  1.1  2.6  0.6  2.5  3.1  4.9  3.2

 [46]  2.1  3.6  0.1  2.4  2.3  3.7  0.6  4.0 12.5  4.3  2.3  3.0  3.4  2.3  2.1

 [61]  3.9  2.6  1.9  2.0  3.7  2.2  0.7  0.2  2.5  2.8  4.1  4.2  3.1  5.9  2.8

 [76]  2.1  3.7  1.2  2.3  2.6  6.3  1.9  2.9  6.3  7.1  4.0  5.3  1.1  2.5  3.6

 [91]  5.7  2.9  1.9  0.9  3.0  2.0  1.6  3.3  4.5  4.3  3.3  2.0  0.6  0.9  0.6

[106]  1.1  1.2  2.7  0.5  3.1  1.2  0.1  0.5  1.3  0.7  3.1  2.4  2.7  1.1  1.7

[121]  3.5  2.2

> Pred

  [1]  0.3  0.7 1.6 4.8 2.2 2.4 1.8 3.1 3.4 2.2 3.0 1.0 0.5 1.2 1.4 1.4 0.7 0.9

 [19] 1.6 1.0 0.8 0.6 1.2 1.7 1.9 2.6 2.7 2.7 5.8 9.4 5.3 3.4 4.1 3.9 3.0 2.8

 [37] 3.7 4.6 6.3 3.3 3.0 4.6 3.8 4.1 3.2 3.4 5.3 7.5 5.9 3.4 2.1 2.6 2.7 3.4

 [55] 5.4 9.7 5.9 3.6 2.9 2.6 2.0 2.9 4.2 5.1 4.2 3.0 3.5 5.2 3.7 3.2 3.2 3.7

 [73] 3.1 4.5 4.8 3.8 6.1 5.3 4.8 8.1 5.9 4.8 6.3 6.6 8.4 6.6 4.8 4.3 1.3 1.3

 [91] 3.1 4.5 2.8 3.6 3.6 4.2 3.6 4.1 3.2 3.4 3.8 3.1 2.4 1.6 0.6 1.3 2.1 1.9

[109] 2.0 1.9 1.5 1.1 0.8 0.9 1.7 1.0 0.7 1.1 1.1 1.1 1.6 2.1

Answer 1

您不能僅將選項設置為連續數據嗎？ A <-verify（obs，pred，frcst.type =“ cont”，obs.type =“ cont”）

如果要將obs轉換為二進制，那很容易。 只需將所有大於2mm的值替換為1，並將所有小於2mm的值替換為0。obs [which（obs => 2）] <-1 obs [which（obs <2）] <-0您可以對預測，但是您的預測類型也是二進制。 否則，您將不得不以某種方式調整您的預測模型，即預測降雨的概率大於2mm，而不是預測實際降雨。 您將不得不從確定性模型轉換為概率模型。