For a statistical experiment I am calculating lots of Gerrity Skill Scores for contingency tables derived from Temperature forecasts/observations.
I tried the package for R, and when I use some relative probabilities as baseline that I have derived from climatological data, I sometimes get Scores greater than 1. I am pretty sure this should not happen.
Here is a script that can be run in batch mode R which should reproduce the error:
I know that the contingency table does not represent the baseline probabilities, but I am still puzzled that this happens.
I also coded the procedure in matlab myself, with the information from the original paper and it reproduces the error.
Is there a condition about the baseline probabilities that I am not aware of?
Any Help is welcome.