numpy.linalg.lstsqからp値を取得するには？

Question

私はnumpy.linalg.lstsqから得たフィットのp値を計算したいと思います。ここでおもちゃの例：

import numpy as np 
x = np.array([[ 58295.62187335],[ 45420.95483714],[ 3398.64920064],[ 977.22166306],[ 5515.32801851],[ 14184.57621022],[ 16027.2803392 ],[ 15313.01865824],[ 6443.2448182 ]]) 
y = np.array([ 143547.79123381, 22996.69597427, 2591.56411049, 661.93115277, 8826.96549102, 17735.13549851, 11629.13003263, 14438.33177173, 6997.89334741]) 

a, res, rank, s = np.linalg.lstsq(x, y) 

前の質問（get the R^2 value from scipy.linalg.lstsq）から私は、しかし、私はまた、p値を計算したいと思い、R²を得るようになった知っています。

ありがとうございます。

Answer 1

あなたはマイナー補正（だと思う） `R、P = pearsonr（X、* yを）` scipy.stats

import numpy as np 
from scipy.stats import pearsonr 
x = np.array([ 58295.62187335, 45420.95483714, 3398.64920064, 977.22166306, 5515.32801851, 14184.57621022, 16027.2803392 , 15313.01865824, 6443.2448182 ]) 
y = np.array([ 143547.79123381, 22996.69597427, 2591.56411049, 661.93115277, 8826.96549102, 17735.13549851, 11629.13003263, 14438.33177173, 6997.89334741]) 

r, p = pearsonr(x,y)