Wie bekomme ich p-Wert von numpy.linalg.lstsq?

Question

Ich möchte den p-Wert der Anpassung berechnen, die ich von numpy.linalg.lstsq erhalten habe. Hier ist ein Spielzeug Beispiel:

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) 

aus vorheriger Frage (get the R^2 value from scipy.linalg.lstsq) Ich weiß, bekommen R² zu bekommen, aber ich würde auch den p-Wert berechnen möge.

Vielen Dank im Voraus.

Answer 1

Sie 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)