我想知道是否有人可以帮助我使用 R 验证包进行一些验证分析。我有两组降水数据，一组是在 20x20x100（纬度 x 经度 x 天）矩阵上的观察，另一组是 20x20x100x5 上的模型结果，其中最后一个维度是集合成员。换句话说，该模型在相同的 100 天内运行了 5 次。
我的目标是根据这些数据绘制可靠性图。数据以英寸/天为单位。
我的主要问题是：1）如何根据这个数据集得到预测概率？换句话说，如何从英寸/天到二进制或分类预测？
2）我想首先我需要选择一些参考值（三分位数？），我将使用它们来获得我的二进制/分类结果。但我不知道如何做到这一点。
任何帮助将不胜感激。谢谢你。亚历克斯·K。
您可以通过首先定义一个阈值来更改为二进制 [0,1]，例如平均降水量 (x1)，然后确定预测值是超过 [1] 还是不超过 [0] 阈值。有了预测类别 (0,1) 后，您现在可以根据观察到的值对其进行测试，以查看哪些同意或不同意...例如，您可以去
预测 <- ifelse(observed$value==x1,0, 1 & ifelse(observed$value>x1,1,0)) 观察到 <- ifelse(observed$value<=x1,0,1)
接下来，您可以使用验证包 Ver<-verify(forecast,observed, frcst.type = "binary", obs.type = "binary") summary(Ver) 中的 verify 函数评估预测和观察之间的一致性
这并不详尽，我希望它有所帮助。理查德：悉尼大学