统计描述的Stata实现

1、-作者xxxx-日期xxxx统计描述的Stata实现【精品文档】第二章 统计描述的Stata实现本章主要用到的Stata主要命令*描述性统计： summary 变量名 weight,detail频数表: tabulate 变量名 weight计算均数： means 变量名 weight频数分布图： histogram 变量名,bin() start()*Stata的命令可以用前几个字母代替，如summary可以由su代替，tabulate可以由tab代替。下面就直接用缩写，不再赘述。例2-1 某市1982年调查120名20岁男子的身高（cm）资料如下，试编制频数表和频数图。数据格式如下：x12

2、3456789101112131415161718192021222324252627281802930171313233343536371753839404117242434445464748491655051525354555657585960616263646517266676869707172737475767778798081828384858687174888990911739216593949596979899169100101102168103104105170106177107108109110172111112113114115116117118119120 为制作频数表，

3、键入Stata命令：.gen f=int(x-160)/2)*2+160产生用以作频数表的新变量“f”.tab f对变量“f”作频数表“gen”命令产生新变量“f”，将各观察值转换成相应该组的下限值。 int为取整函数，结果为括号内函数值的整数部分，如int(3.24)=3。“160”为第一组的下限，“2”为组距。以第一例观察值0)/2)*2+160=160，则它应归入“160”组。结果如下： f | Freq. Percent Cum.-+-+-hist f,start(160) width(2)作频数图结果如下：例2-2某医生测定230名正常成年男子的空腹血清胰岛素样生长因子-1（IGF-

4、1F）水平，整理后编制为频数分布表（表2-2），请根据该频数分布表作频数图。表2-2 230名正常人空腹血清胰岛素样生长因子-1水平中位数的计算 IGF-1F （1） 人数 （2）频率(%)(3)累计频数（4） 累计频率(%) （5）=（4）/ 50 30 13.04 30 150 71 30.87 101 250 49 21.30 150 350 28 12.17 178 450 14 6.09 192 550 12 5.22 204 650 10 4.35 214 750 8 3.48 222 850 5 2.17 227 9501050 3 1.30 230 合计230Stata数据格式

5、如下： xf150302150713250494350285450146550127650108750898505109503tab x weight=f制作汇总表的频数表在本例中，x的取值为各组的下限，而f是各组的频数，所以需要用Stata的频数选项weight=。结果： x | Freq. Percent Cum.-+- -+-hist x weight=f,start(50) width(100)制作频数图例2-3 随机测量某地10名2030岁健康男性居民血清铁含量（），测量值分别为6.58，7.42，15.32，15.78，17.60，17.98，15.21，17.53，20.11，2

6、2.64，试求其平均血清铁含量。Stata数据为x12345678910su x计算均数标准差结果： Variable | Obs Mean Std. Dev. Min Max-+- x | 10 所以平均平均血清铁含量为15.617。例2-4 计算例2-1的频数表（表2-3）中120名男子的平均身高显然，如果用Stata，对于例2-1的资料不需要事先转化为频数表，然后再计算均数。不妨用例2-2作为例子展示Stata计算频数表资料的均数。su x weight=f计算频数表资料的均数结果： Variable | Obs Weight Mean Std. Dev. Min Max-+- x |

7、10 230 290 223.4123 50 950例2-5 8名麻疹易感儿接种麻疹疫苗3周后，其血凝抑制抗体滴度分别为1:4，1:8，1:16，1:32，1:64，1:128，1:256，1:512。试求其平均抗体滴度。数据格式：x1428316432564612872568512gen logx=log(x)计算x的自然对数su logx计算自然对数的均数，得到disp exp(3.81231)计算均数的反对数，即指数得到几何均数为7另外，Stata还有一个直接计算几何均数的命令：means x结果为： Variable | Type Obs Mean 95% Conf. Interval

8、-+- x | Arithmetic 8 127.5 -20.53203 275.532 | Geometric 8 45.25483 10.94481 187.1206 | Harmonic 8 16.06275 . . -Missing values in confidence interval(s) for harmonic mean indicate that confidence interval is undefined for corresponding variable(s).Consult Reference Manual for details.Means命令计算算术均数、

9、几何均数以及调和均数，其中几何均数为45.25483，和前面的结果相同。例2-6 某地区50名麻疹易感儿童接种麻疹疫苗3周后，测其血凝抑制抗体滴度，如表2-4中第（1）栏和第（2）栏，求平均抗体滴度。表2-4 50名麻疹易感儿童平均抗体滴度计算表 抗体滴度 （1） 人数 （2） 滴度倒数 （3）（4）（5）=（2）×（4）1:4 1 41:8 2 81:16 6 161:3210 321:6416 641:128 81281:256 52561:512 2512 合 计50数据如下：xf14128231664321056416612887256585122gen logx=log(x

10、)计算x的自然对数su logx weight=f计算自然对数的均数，得到disp exp()计算均数的反对数，即指数得到几何均数为也可以用means x weight=f结果：(analytic weights assumed) Variable | Type Obs Mean 95% Conf. Interval-+- x | Arithmetic 8 95.76 -2.30755 193.8275 | Geometric 8 22.59723 145.2029 | Harmonic 8 32.82051 . . -Missing values in confidence interval

11、(s) for harmonic mean indicate that confidence interval is undefined for corresponding variable(s).Consult Reference Manual for details.例2-7 5名成年男子的体重（kg）为60，70，75，80， 90，求中位数。例2-8 某医生测定了6名正常成年男子的空腹血清胰岛素样生长因子-1（IGF-1F）水平为150，170，185，245，265，280，求中位数。注：这两个例子用软件算实在是浪费，可以用su x,detail来处理。对于summary，加上det

12、ail后可以输出百分位数，其中的P50就是中位数。 x- Percentiles Smallest 1% 60 60 5% 60 7010% 60 75 Obs 525% 70 80 Sum of Wgt. 550% 75 Mean 75 Largest Std. Dev. 75% 80 7090% 90 75 Variance 12595% 90 80 Skewness 0例2-9 8名7岁男孩的身高（cm）分别为116，118，119，120，121，123，125，126。请计算这8名7岁男孩身高的第25百分位数和第90百分位数。Stata数据为：x116118119120121123125126Stata命令为：su x,d结果为： x- Percentiles Smallest 1% 116 116 5% 116 11810% 116 119 Obs 8 120 Sum of Wgt. 850% 120.5 Mean 12175% 124 12190% 126 123