Tom, you are wildly off topic as usual.
If you want to show the world that you have reproduced the negative cosine prediction of the singlet state, but you leave off the vertical axis of the plot in order to hide the fact that you did not get perfect anti-correlation at equal settings but only correlation -0.98, you are a cheat.
Take a look at
http://rpubs.com/gill1109/JCS2
All the formulas in that code are taken from Joy's own description of how to simulate his model.
We are not talking about Christian's original S^3 model (whatever that is) but a derived version of it (Joy's derivation, not mine!) which can be expressed in S^2 terms. Implemented in Java code written by Chantal Roth, authorized by Joy; and another Javascript version written by Daniel Sabsay, authorized by Joy, proudly exhibited on his web page. My version is a port of Chantal's to R, and it has her blessing, at least.
The formulas in all these simulation experiments are taken from Joy's paper "Whither..." http://arxiv.org/abs/1301.1653 (appendix on simulation models). Cab is defined in formula (A27), Na and Nb are given by the expressions just after (A28) and (A29). The procedure generating an outcome +/-1 for the product of Alice and Bob's measurement, or 0 for "no state", with the help of an auxiliary randomization, is taken from Roth and Sabsay's code, but appears to me also to be a true reflection of Joy's own description. The numerical values of the "phase shifts" are taken from Joy's paper. Joy agrees that these are not quite spot on, yet.
Please note: I am not saying that the ommission of the vertical axis was done deliberately and with the purpose to deceive. I am not accusing any particular person of cheating.