आर^2

Question

पर Numpy Lstsq अवशिष्ट मूल्य कनवर्ट करना मैं नीचे के रूप में कम से कम वर्ग प्रतिगमन कर रहा हूं (univariate)। मैं आर^2 के संदर्भ में परिणाम का महत्व व्यक्त करना चाहता हूं। Numpy unscaled अवशिष्ट का एक मूल्य देता है, यह सामान्य करने का एक समझदार तरीका क्या होगा।

field_clean,back_clean = rid_zeros(backscatter,field_data) 
num_vals = len(field_clean) 
x = field_clean[:,row:row+1] 
y = 10*log10(back_clean) 

A = hstack([x, ones((num_vals,1))]) 
soln = lstsq(A, y) 
m, c = soln [0] 
residues = soln [1] 

print residues 

Answer 1

देखें http://en.wikipedia.org/wiki/Coefficient_of_determination

आपका आर 2 मूल्य =

1 - residual/sum((y - y.mean())**2) 

जो एक उदाहरण के रूप

1 - residual/(n * y.var()) 

के बराबर है:

import numpy as np 

# Make some data... 
n = 10 
x = np.arange(n) 
y = 3 * x + 5 + np.random.random(n) 

# Note that polyfit is an easier way to do this... 
# It would just be "model, resid = np.polyfit(x,y,1,full=True)[:2]" 
A = np.vstack((x, np.ones(n))).T 
model, resid = np.linalg.lstsq(A, y)[:2] 

r2 = 1 - resid/(y.size * y.var()) 
print r2