我观察并模拟了降水。我想使用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发布于 2014-06-02 19:34:41
您不能只将选项设置为连续数据吗?A<- verify(obs,pred,frcst.type = "cont",obs.type = "cont")
如果您想要将ob转换为二进制,这非常简单。只需将所有大于2 2mm的值替换为1,将所有小于2 2mm的值替换为0。obswhich(obs=>2)<-1 obswhich(obs<2)<-0您可以对预测执行相同的操作,但您的预测类型也是二进制的。否则,您必须调整您的预测模型,使其预测降雨量大于2毫米的概率,而不是预测实际降雨量。您将不得不从确定性模型切换到概率模型。
https://stackoverflow.com/questions/22166096
复制相似问题