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