jdiixjxjs

(1) # coding: utf-8

(2) import scipy.stats

(3) Q1=scipy.stats.norm(0,1).ppf(0.25)

(4) Q3=scipy.stats.norm(0,1).ppf(0.75)

(5) Upperfence=scipy.stats.norm(0,1).cdf(Q3+1.5*(Q3-Q1))

(6) Lowerfence=scipy.stats.norm(0,1).cdf(Q1-1.5*(Q3-Q1))

(7) probUL=round(Upperfence-Lowerfence,4)

(8) probOutliers=1-probUL

(9) print (u'Q1-μ =%.4fσ ,Q3-μ=%.4fσ '%(Q1,Q3))

(10) print (u'IQR=Q3-Q1=%.4f σ '%(Q3-Q1))

(11) print (u'Q3+1.5×IQR-μ =%.4fσ '%(Q3+1.5*(Q3-Q1)))

(12) print (u'Q1-1.5×IQR-μ =%.4fσ '%(Q1-1.5*(Q3-Q1)))

(13) print (u'P(Q1-1.5×IQR<x<Q3+1.5×IQR)=%.4f'%(probUL))

(14) print (u'在上下限之外的概率=%.4f%%'%(100*probOutliers))