Tengo un conjunto de datos 3D ajustado usandoscipy.linalg.lstsq función.
Estaba usando:

# best-fit quadratic curve    
A = np.c_[np.ones(data.shape[0]), data[:,:2], np.prod(data[:,:2], axis=1), data[:,:2]**2]    
C,_,_,_ = scipy.linalg.lstsq(A, data[:,2])    
#evaluating on grid      
Z = np.dot(np.c_[np.ones(XX.shape), XX, YY, XX*YY, XX**2, YY**2], C).reshape(X.shape)

Pero, ¿cómo puedo obtener el valor R ^ 2 de esto para la superficie ajustada?
¿Hay alguna forma de verificar la importancia del resultado del ajuste?

Cualquier idea relacionada con eso sería muy apreciada.
gracias.