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＆lt; - 验证（obs，pred，frcst.type =＆＃34; cont＆＃34;，obs.type =＆＃34; cont＆＃34;）

如果你想将obs转换为二进制，那很容易。只需将所有大于2毫米的值替换为1，将所有小于2毫米的值替换为0。 OBS [其中（OBS =→2）]＆LT; -1 OBS [其中（OBS＆2）]＆LT; -0 您可以对预测执行相同的操作，但是您的预测类型也是二进制的。否则，您必须以预测降雨概率大于2毫米的方式调整预测模型，而不是预测实际降雨量。您必须从确定性模型切换到概率模型。

降水预报R包的验证

1 个答案: