例2.5软件实现（spss/sas/stata/r）

例2-5 69例类风湿关节炎（RA）患者血清EBV-VCA-lgG抗体滴度的分布如下，求其平均抗体滴度。
抗体滴度 1:10 1:20 1:40 1:80 1:160 1:320 1:640 1:1280
人数 4 3 10 10 11 15 14 2

SPSS代码及结果

先点击转换–计算变量 再点击分析—描述统计—描述

COMPUTE lgx = LG10(x).
EXECUTE.
DESCRIPTIVES
  VARIABLES=lgx /SAVE
  /STATISTICS=MEAN STDDEV MIN MAX KURTOSIS SKEWNESS.

SAS代码及结果

data ex2_5;
  input x f @@;
  y=log10(x);
cards;
10      4
20      3
40     10
80     10
160    11
320    15
640    14
1280    2
;
proc means noprint;
  var y;
  freq f;
  output out=b
  mean=logmean;
run;
data c;
  set b;
  g=10**logmean;
proc print data=c;
  var g;
run;

**程序说明：**第一步创建数据集，有三个变量x,f和y，x是抗体滴度的倒数，f是频数，y是x的对数。第二步用means过程计算y的描述性统计量，将计算得到的均数输出到数据集b中，用变量名logmean表示，noprint表示不在output窗口显示means过程的结果。第三步新建数据集c，新产生g变量，为变量logmean的反对数，第四步将数据集c的结果显示在output窗口中。

运行结果：

STATA代码及结果

ameans x  /*计算变量x的算术均数，几何均数和调和均数*/

结果：

R代码及结果

R可实现的包较多，此处不做过多描述。

library(readxl)
stat.desc(x)