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