Stata软件基本操作和数据分析入门Word文档格式.docx

1、10 12 11 14.7 12 11.4 13 13.8 12.3 12 操作如下:gen d=x1-x2 产生配对差值的变量d sktest d 正态性检验 正态性检验结果如下:. sktest d Skewness/Kurtosis tests for Normality - joint - Variable | Pr（Skewness） Pr（Kurtosis） adj chi2（2） Probchi2 -+- d | 0.279 0.774 1.43 0.4885 正态性检验的无效假设为:资料正态分布 相应的备选假设为:资料非正态分布 ,=0.05，由于正态性检验的P值=0.4018

2、9,，故可以认为资料近似服从正态分布。ttest d=0 配对t检验: H:=0 vs H:0，0d1d=0.05 结果如下:One-sample t test - Variable | Obs Mean Std. Err. Std. Dev. 95% Conf. Interval -+- d | 10 -.6799999 .5204272 1.645735 -1.857288 .4972881 - Degrees of freedom: 9 Ho: mean（d） = 0 Ha: mean 0 t = -1.3066 t = -1.3066 t = -1.3066 P |t| = 0.223

3、7 P t = 0.8881 P值=0.2237,，故认为治疗前后的血红蛋白的平均数差异没有统计学意义。即:没有足够的证据可以认为治疗前后的血红蛋白的总体平均数不同。如果已知差值的样本量，样本均数和样本标准差，可以用立即命令如下（如，已知样本量为10，差值的样本均数为-0.66，差值的标准差为1.65，则输入命令如下:ttesti 样本量 样本均数 样本标准差 0 本例为: ttesti 10 -0.66 1.65 0 得到下列结果如下:. ttesti 10 .66 1.65 0 - | Obs Mean Std. Err. Std. Dev. 95% Conf. Interval -+-

4、x | 10 .66 .5217758 1.65 -.5203389 1.840339 - Degrees of freedom: mean（x） = 0 t = 1.2649 t = 1.2649 t = 1.2649 t = 0.8812 P |t| = 0.2377 P t = 0.1188 结果解释与结论同上述相同。如果对于小样本的情况下，差值不满足正态分布，则用Match-Sign-rank test，操作如下:signrank 差值变量名=0 假如本例不满足正态分布（为了借用上例资料，而假定的，实际上本例满足正态分布）则 H:差值的中位数,0 0（其意义是治疗前的血红蛋白配大于治疗

5、后的血红蛋白的概率,治疗前的血红蛋白小于治疗后的血红蛋白的概率） 差值的中位数0 1本例为 signrank d=0 Wilcoxon signed-rank test sign | obs sum ranks expected -+- positive | 4 18 27 negative | 5 36 27 zero | 1 1 1 -+- all | 10 55 55 unadjusted variance 96.25 adjustment for ties 0.00 adjustment for zeros -0.25 - adjusted variance 96.00 d = 0 z

6、 = -0.919 Prob |z| = 0.3583 P值=0.3583，故没有足够的证据说明两个总体不同。二、平行对照设计的两组资料平均水平统计检验 如果两组资料的方差齐性和相互独立的，并且每组资料服从正态分布（大样本资料可以忽略正态性问题），则用成组t检验，否则可以用成组Wilcoxon秩和检验。例2 为研究噪声对纺织女工子代智能是否有影响，一研究人员在某纺织厂随机抽取接触噪声95dB（A）、接触工龄5年以上的纺织女工及同一单位、条件与接触组相近但不接触噪声的女职工，其子女（学前幼儿）作为研究对象，按韦氏学前儿童智力量表（中国修订版）测定两组幼儿智商，结果如下。问噪声对纺织女工子代智能有

7、无影响,（接触组group=0，不接触组group=1） 资料及其结果如下:group x 0 79 0 93 0 91 0 92 0 94 0 77 0 93 0 74 0 91 0 101 0 83 0 73 0 88 0 102 0 90 0 100 0 81 0 91 0 83 0 106 0 84 0 78 0 87 0 95 0 101 1 101 1 100 1 114 1 86 1 106 1 107 1 107 1 94 1 89 1 104 1 98 1 110 1 89 1 103 1 121 1 94 1 95 1 92 1 109 1 120 1 104 方差齐性检验

8、 H, :vs H:012112=0.1 两组方差齐性的检验命令（仅适合两组方差齐性检验） sdtest x,by（group） Variance ratio test Group | Obs Mean Std. Err. Std. Dev. 95% Conf. Interval -+- 0 | 25 89.08 1.822928 9.11464 85.31766 92.84234 1 | 25 101.52 1.900982 9.504911 97.59657 105.4434 -+- combined | 50 95.3 1.577456 11.1543 92.12998 98.47002

9、 - sd（0） = sd（1） F（24,24） observed = F_obs = 0.920 F（24,24） lower tail = F_L = F_obs = 0.920 F（24,24） upper tail = F_U = 1/F_obs = 1.087 sd（0） sd（1） F_obs = 0.4195 P F_U = 0.8389 P F_obs = 0.5805 P值=0.8389,，因此可以认为两组方差齐性的。正态性检验:资料服从正态分布 vs H:资料偏态分布 01每一组资料正态性检验 sktest x if group=0 x | 0.927 0.326 1.05 0.5926 . sktest x if group=1 Variable | P