matplotlib

Question

का उपयोग कर एक sccatterplot में एक विश्वास एलिप्स को बनाना matplotlib का उपयोग कर एक sccatterplot में एक विश्वास एलिप्स बनाने के लिए कैसे?

निम्नलिखित कोड स्कैटर प्लॉट बनाने तक काम करता है। फिर, क्या कोई भी डालने से परिचित है स्कैटर प्लॉट पर विश्वास इलिप्स?

import numpy as np 
import matplotlib.pyplot as plt 
x = [5,7,11,15,16,17,18] 
y = [8, 5, 8, 9, 17, 18, 25] 

plt.scatter(x,y) 
plt.show() 

एसएएस से विश्वास एलिप्स के संदर्भ में निम्नलिखित संदर्भ है।

http://support.sas.com/documentation/cdl/en/grstatproc/62603/HTML/default/viewer.htm#a003160800.htm

sas में कोड इस तरह है:

proc sgscatter data=sashelp.iris(where=(species="Versicolor")); 
    title "Versicolor Length and Width"; 
    compare y=(sepalwidth petalwidth) 
      x=(sepallength petallength) 
     /reg ellipse=(type=mean) spacing=4; 
run; 

Answer 1

निम्नलिखित कोड एक एक, दो, तीन और मानक विचलन आकार दीर्घवृत्त ड्रॉ:

x = [5,7,11,15,16,17,18] 
y = [8, 5, 8, 9, 17, 18, 25] 
cov = np.cov(x, y) 
lambda_, v = np.linalg.eig(cov) 
lambda_ = np.sqrt(lambda_) 
from matplotlib.patches import Ellipse 
import matplotlib.pyplot as plt 
ax = plt.subplot(111, aspect='equal') 
for j in xrange(1, 4): 
    ell = Ellipse(xy=(np.mean(x), np.mean(y)), 
        width=lambda_[0]*j*2, height=lambda_[1]*j*2, 
        angle=np.rad2deg(np.arccos(v[0, 0]))) 
    ell.set_facecolor('none') 
    ax.add_artist(ell) 
plt.scatter(x, y) 
plt.show() 

Answer 2

स्वीकृत उत्तर देने के बाद, मैंने पाया कि मैं टी सही ढंग से वृत्त का चतुर्थ भाग का चयन नहीं करता है जब थीटा की गणना, के रूप में यह np.arccos पर निर्भर करता है:

'possible duplicate' और Joe Kington's solution on github पर एक नज़र ले रहा है, मैं इस के लिए नीचे अपने कोड पानी पिलाया:

import numpy as np 
import matplotlib.pyplot as plt 
from matplotlib.patches import Ellipse 

def eigsorted(cov): 
    vals, vecs = np.linalg.eigh(cov) 
    order = vals.argsort()[::-1] 
    return vals[order], vecs[:,order] 

x = [5,7,11,15,16,17,18] 
y = [25, 18, 17, 9, 8, 5, 8] 

nstd = 2 
ax = plt.subplot(111) 

cov = np.cov(x, y) 
vals, vecs = eigsorted(cov) 
theta = np.degrees(np.arctan2(*vecs[:,0][::-1])) 
w, h = 2 * nstd * np.sqrt(vals) 
ell = Ellipse(xy=(np.mean(x), np.mean(y)), 
       width=w, height=h, 
       angle=theta, color='black') 
ell.set_facecolor('none') 
ax.add_artist(ell) 
plt.scatter(x, y) 
plt.show()