from scipy import stats
x = [2019,2020,2021,2022]
y = [161000,182000,257000,282000]
slope,intercept,r,p,std_err = stats.linregress(x,y)
def myfunc(x):
return slope * x + intercept
print(r)
ZnJvbSBzY2lweSBpbXBvcnQgc3RhdHMKCnggPSBbMjAxOSwyMDIwLDIwMjEsMjAyMl0KeSA9IFsxNjEwMDAsMTgyMDAwLDI1NzAwMCwyODIwMDBdCgpzbG9wZSxpbnRlcmNlcHQscixwLHN0ZF9lcnIgPSBzdGF0cy5saW5yZWdyZXNzKHgseSkKCmRlZiBteWZ1bmMoeCk6CiAgIHJldHVybiBzbG9wZSAqIHggKyBpbnRlcmNlcHQKCnByaW50KHIp