import numpy as np
from sklearn.metrics import r2_score
from scipy import stats
import statsmodels.api as sm
import math

df1 = np.loadtxt('SAOBS_1981_2007_DD.csv', delimiter=',')
df2 = np.loadtxt('APHRO_1981_2007_DD.csv', delimiter=',')

correlation_r2 = r2_score(df1, df2)