from scipy.stats import pearsonr
import matplotlib.pyplot as plt

#Create two random variable
X = [1,2,3,4,5,6,7,8,9,10,11,12,13]
Y = [7,8,8,11,10,11,12,15,20,14,16,15,19]

#plot the variables to show linearity
plt.scatter(X,Y)
plt.show()