《非参数统计》与MATLAB编程 第二章 描述性统计

第二章描述性统计

2．1 表格法和图形法

表2.1 灯丝寿命数据

107 73 68 97 76 79 94 59 98 57 73 81 54 65 71 80 84 88 62 61 79 98 63 65 66 62 79 86 68 74 61 82 65 98 63 71 62 116 65 88 64 79 78 79 77 86 89 76 74 85 73 80 68 78 89 72 58 69 82 72 92 78 88 77 103 88 63 68 88 81 64 73 75 90 62 89 71 71 74 70 74 70 85 61 65 81 75 62 94 71 85 84 83 63 92 68 81 62 79 83 93 61 65 62 92 65 64 66 83 70 70 81 77 72 84 67 59 58 73 83 78 66 66 94 77 63 66 75 68 76 73 76 90 78 71 101 78 43 59 67 61 71 77 91 96 75 64 76 72 77 74 65 82 86 79 74 66 86 96 89 81 71 85 99 59 92 94 62 68 72 77 60 87 84 75 77 51 45 63 102 85 67 87 80 84 93 69 76 89 75 59 77 83 68 72 67 92 89 82 96

a =

Columns 1 through 17

107 73 68 97 76 79 94 59 98 57 73

81 54 65 71 80 84

79 98 63 65 66 62 79 86 68 74 61

82 65 98 63 71 62

64 79 78 79 77 86 89 76 74 85 73

80 68 78 89 72 58

92 78 88 77 103 88 63 68 88 81 64

73 75 90 62 89 71

74 70 85 61 65 81 75 62 94 71 85

84 83 63 92 68 81

93 61 65 62 92 65 64 66 83 70 70

81 77 72 84 67 59

78 66 66 94 77 63 66 75 68 76 73

76 90 78 71 101 78

61 71 77 91 96 75 64 76 72 77 74 65 82 86 79 74 66

81 71 85 99 59 92 94 62 68 72 77 60 87 84 75 77 51

85 67 87 80 84 93 69 76 89 75 59 77 83 68 72 67 92

Columns 18 through 20

88 62 61

116 65 88

69 82 72

71 74 70

62 79 83

58 73 83

43 59 67

86 96 89

45 63 102

89 82 96

注：a不能复制到MATLAB中。

b=reshape(a,200,1);

[min(b) max(b)]

ans =

43 116

n=histc(b,[40,50,60,70,80,90,100,110,120])

n =

2

10

52

64

45

22

4

1

n/200

ans =

0.0100

0.0500

0.2600

0.3200

0.2250

0.1100

0.0200

0.0050

bar([40,50,60,70,80,90,100,110,120],n,'histc')

表2.2 灯丝寿命数据频率分布

灯丝寿命（小时）个数频率（%）

40—50 2 1

50—60 10 5

60—70 52 26

70—80 64 32

80—90 45 22.5

90—100 22 11

100—110 4 2

110-120 1 0.5

合计200 100

c=[2 32 24 19 10 6 3 2 1 1]

bar(c)

表2.7 赔款样本数据直方图x=[1 2 2 2 2 3 3 5 7 8 9];y=[1 1 2 3 4 1 2 1 1 1 1];

plot(x,y,'*','markersize',30,'linewidth',4)

axis([0,10,0,5])

图2.8

[mean([1,2,2,2,2,3,3,5,7,8,9]) median([1,2,2,2,2,3,3,5,7,8,9])] ans =

4 3

切尾均值：trimmean

trimmean(x,percent) ，其中percent为0到100。当x为向量时，切尾均值为：先去掉x的percent/2的最大值和最小值，再求算术平均值。注意：若length(x)×percent/2恰好为整数时，则就x按大小排序后，两端去掉length(x)×percent/2个数，如果不是整数，则去掉length(x)×percent/2四舍五入数（若小数点后正好等于5，则取小的整数，如1.5还是取1）。

如：a=[4 9 12.4 23 33.5 45 56 67 76 99] 回车

trimmean(a,50) 回车

ans =

39.4833

（25/2）%×length(x)=2.5,两端切掉2个数。即等于：

mean([12.4 23 33.5 45 56 67]) 回车

ans =

39.4833

然后当：

trimmean(a,51) 回车

ans =

39.3750

（51/2）%×length(x)=2.55,即两端切掉3个数。

mean([23 33.5 45 56]) 回车

ans =

39.3750

Trimmean 也可处理矩阵，不过只是计算矩阵每列的切尾均值。继续上面，如：b=a';c=[1:10]';d=[b,c];trimmean(d,51) 回车