R验证包

Question

我解决了读取数据的问题，并从 R 的验证包中得到了一些图。现在我想知道如何在 R 中排列数据，即我排列了像 (120, 396) 这样的观察数据，其中 120 是前导时间和 396 是总初始时间。预测数据有一个额外的集合成员维度，如 (120, 396,10)。在将这些数据集放入验证包的验证功能之前，我喜欢这样

obs=round(runif(100))
forecast=runif(100)

我不知道为什么这个过程是必要的以及它做了什么，但在投入验证功能之前是必要的。然后使用我使用的验证功能，

verify(obs,forecast,obs.type="binary",frcst.type="prob")

然后一切正常，我可以将验证功能的结果用于其他事情。但我不知道验证包是如何工作的，我的数据排列是否正确？上述数据集120、396中，哪个维度作为样本维度？

Answer 1

您必须告诉verify 您希望它查看哪些数据以及它是什么类型的数据。因此，您显示的示例行有两个向量，其中 obs 是 0 或 1 并且是二进制变量，而 forecast 在 0 和 1 之间是连续的并且是概率。

您的数据看起来好像在观察和预测中都是连续的，您可能希望将verify 应用于每对或三元组的第二个元素。如果是这样，查看此示例可能会帮助您了解 verify 在这种情况下如何工作。

x <- 120:130
ob <- data.frame(cbind(x, y=3 * x + 40 * rnorm(11) ) )
ob
fc <- data.frame(cbind(x, y=3 * x + 10 * rnorm(11), z = rep(10,11) ) )
fc
V  <- verify(ob$y, fc$y, obs.type = "cont", frcst.type = "cont" )
summary(V)

显然，当您验证自己的观察和预测时，您需要使用自己的变量名称。