资源描述:
《《非参数统计》与MATLAB编程第二章描述性统计.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、第二章描述性统计2.1表格法和图形法表2.1灯丝寿命数据10773689776799459985773815465718084886261799863656662798668746182659863716211665886479787977868976748573806878897258698272927888771038863688881647375906289717174707470856165817562947185848363926881627983936165629265646683707081777284675958738
2、3786666947763667568767376907871101784359676171779196756476727774658286797466869689817185995992946268727760878475775145631028567878084936976897559778368726792898296a=Columns1through171077368977679945998577381546571808479986365666279866874618265986371626479787977868976748
3、573806878897258927888771038863688881647375906289717470856165817562947185848363926881936165629265646683707081777284675978666694776366756876737690787110178617177919675647672777465828679746681718599599294626872776087847577518567878084936976897559778368726792Columns18throug
4、h2088626111665886982727174706279835873834359678696894563102898296注:a不能复制到MATLAB中。b=reshape(a,200,1);[min(b)max(b)]ans=43116n=histc(b,[40,50,60,70,80,90,100,110,120])n=21052644522410n/200ans=0.01000.05000.26000.32000.22500.11000.02000.00500bar([40,50,60,70,80,90,100,110,
5、120],n,'histc')表2.2灯丝寿命数据频率分布灯丝寿命(小时)个数频率(%)40—502150—6010560—70522670—80643280—904522.590—1002211100—11042110-12010.5合计200100c=[23224191063211]bar(c)表2.7赔款样本数据直方图x=[12222335789];y=[11234121111];plot(x,y,'*','markersize',30,'linewidth',4)axis([0,10,0,5])图2.8[mean([1,2,2
6、,2,2,3,3,5,7,8,9])median([1,2,2,2,2,3,3,5,7,8,9])]ans=43切尾均值:trimmeantrimmean(x,percent),其中percent为0到100。当x为向量时,切尾均值为:先去掉x的percent/2的最大值和最小值,再求算术平均值。注意:若length(x)×percent/2恰好为整数时,则就x按大小排序后,两端去掉length(x)×percent/2个数,如果不是整数,则去掉length(x)×percent/2四舍五入数(若小数点后正好等于5,则取小的整数,如1
7、.5还是取1)。如:a=[4912.42333.54556677699]回车trimmean(a,50)回车ans=39.4833(25/2)%×length(x)=2.5,两端切掉2个数。即等于:mean([12.42333.5455667])回车ans=39.4833然后当:trimmean(a,51)回车ans=39.3750(51/2)%×length(x)=2.55,即两端切掉3个数。mean([2333.54556])回车ans=39.3750Trimmean也可处理矩阵,不过只是计算矩阵每列的切尾均值。继续上面,如:b=
8、a';c=[1:10]';d=[b,c];trimmean(d,51)回车ans=39.37505.50002.2.2表示离散程度的数值[var(b)std(b)range(b)]ans=145.454812.06057