Jump to content
The Weather Forums

Recommended Posts

Hello

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:

 

library('verification')
A=matrix(c(2394,481,0,0,0,0,4718,1399,0,0,0,0,1740,1755,0,0,0,0,3202,594,0,0,0,0,1717),nrow=5)
pb=c(0.0632,0.2692,0.1461,0.2338,0.2877)
multi.cont(A,pb)

 

 

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.

  • Like 1
Link to post
Share on other sites
  • 4 weeks later...

Hello

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:

 

library('verification')

A=matrix(c(2394,481,0,0,0,0,4718,1399,0,0,0,0,1740,1755,0,0,0,0,3202,594,0,0,0,0,1717),nrow=5)

pb=c(0.0632,0.2692,0.1461,0.2338,0.2877)

multi.cont(A,pb)

 

 

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.

I work in matlab, send me that code if you want and I can look. 

  • Like 1
Link to post
Share on other sites

Join the conversation

You can post now and register later. If you have an account, sign in now to post with your account.

Guest
Reply to this topic...

×   Pasted as rich text.   Paste as plain text instead

  Only 75 emoji are allowed.

×   Your link has been automatically embedded.   Display as a link instead

×   Your previous content has been restored.   Clear editor

×   You cannot paste images directly. Upload or insert images from URL.

Loading...
×
×
  • Create New...