import matplotlib.pyplot as plt
import numpy as np

r = 0.8

X = np.random.uniform(size=100)
Y = 0.5 + 1/(0.3+X) + 0.25*np.random.normal(size=100)

XR = np.argsort(np.argsort(X))
YR = np.argsort(np.argsort(Y))

cr = np.cov(XR,YR)
cr = cr[0,1]/np.sqrt(cr[0,0]*cr[1,1])

plt.clf()
plt.figure(figsize=(4,3.8))
plt.axes([0.17,0.12,0.75,0.8])
plt.hold(True)
plt.plot(X, Y, 'o', color='orange')
plt.grid(True)
plt.xlabel("X", size=16)
plt.ylabel("Y", size=16)
plt.title("Spearman correlation=%.2f" % cr)
plt.savefig("spearman_fig4.pdf")
plt.savefig("spearman_fig4.svg")