《Stata拟合优度检验命令小结》由会员分享，可在线阅读，更多相关《Stata拟合优度检验命令小结（9页珍藏版）》请在金锄头文库上搜索。
1、C:UsersuserDesktopgoodness of fit2goodness of fit2.do 2016年年1月月29日日 0:07 /*医用统计方法 第三版 金丕焕 第十七章第二节 p238 /*例17-4 （样本分布与总体分布比较）研究四环素样品出厂时干燥失重情况是否服从于正态分布 /*mgof 命令 数据格式： 第1列 c：干燥失重分组 第2列 exppro：正态分布的理论百分数，不是累计百分数 第3列 bsnum：实际频数 */ clear set obs 9 input c exppro obsnum 1 1.62 2 2 6.02 10 3 16.25 20 4 26.
2、11 48 5 26.11 49 6 16.25 19 7 6.02 8 8 1.41 3 9 0.19 1 capture end mgof c=exppro fweight=obsnum,mc ksmirnov reps(100) / 结果：Ksmirnov-Smirnov D=0.0449977 Exact P-value=0.49 / 与书中计算结果D=0.0452略有差异 /*ksmirnov 命令 数据格式： 第1列 c：干燥失重分组 第2列 expcpro：正态分布的累计理论百分数 第3列 bsnum：实际频数 */ clear set obs 9 input c expcpro
3、 obsnum 1 0.0162 2 2 0.0764 10 3 0.2389 20 4 0.5 48 5 0.7611 49 6 0.9236 19 7 0.9838 8 8 0.9979 3 9 0.9998 1 capture end expand obsnum -1- C:UsersuserDesktopgoodness of fit2goodness of fit2.do 2016年年1月月29日日 0:07 ksmirnov c=expcpro / 结果：D=0.0452 p-value=0.521 / 与书中计算结果D=0.0452一致 /*例17-5 （样本分布与总体分布比较）
4、研究5年中各月份的脊髓灰质炎发病有无季节差别 /*mgof 命令 数据格式： 第1列 c：月份分组 第2列 exppro：均匀分布的理论百分数，不是累计百分数 第3列 obsnum：实际频数 */ clear set obs 12 input c exppro obsnum 1 8.333 63 2 8.333 78 3 8.333 140 4 8.333 117 5 8.333 105 6 8.333 101 7 8.333 144 8 8.333 127 9 8.333 79 108.333 87 118.333 58 128.333 48 capture end mgof c=exppr
5、o fweight=obsnum,mc ksmirnov reps(100) / 结果：Ksmirnov-Smirnov D=0.096193 Exact P-value=0.0000 / 与书中计算结果D=0.0962一致 /*ksmirnov 命令 数据格式： 第1列 c：月份分组 第2列 expcpro：均匀分布的累计理论百分数 第3列 obsnum：实际频数 */ clear set obs 12 input c expcpro obsnum 1 0.08333 63 2 0.16666 78 3 0.24999 140 4 0.33332 117 5 0.41665 105 6 0.
6、49998 101 7 0.58331 144 -2- C:UsersuserDesktopgoodness of fit2goodness of fit2.do 2016年年1月月29日日 0:07 8 0.66664 127 9 0.74997 79 100.8333 87 110.91663 58 120.99996 48 capture end expand obsnum ksmirnov c=expcpro / 结果：D=0.0962 Exact P-value=0.0000 / 与书中计算结果D=0.0962一致 /*例17-6 （两样本分布比较）研究两种方法抽样调查人口构成是否有
7、统计学意义 /*mgof 命令 数据格式： 第1列 c：年龄组 第2列 pro_bsample：乙样本人口构成，不是累计百分数 第3列 obsnum_asample：甲实际频数（55236） */ clear set obs 17 input c pro_bsample obsnum_asample 1 0.105 5855 2 0.137 7567 3 0.105 4750 4 0.113 5634 5 0.106 4971 6 0.094 4198 7 0.071 3480 8 0.058 3148 9 0.052 2983 100.043 2707 110.041 2541 120.03
8、4 1988 130.021 1657 140.009 1657 150.006 1049 160.003 718 170.002 331 capture end mgof c=pro_bsample fweight=obsnum_asample,mc ksmirnov reps(100) / 结果：Ksmirnov-Smirnov D=0.071996 Exact P-value=0.0000 / 与书中计算结果D=0.072一致 /*mgof 命令 数据格式： 第1列 c：年龄组 第2列 pro_asample：甲样本人口构成，不是累计百分数 第3列 obsnum_bsample：乙样本实
9、际频数（49827） -3- C:UsersuserDesktopgoodness of fit2goodness of fit2.do 2016年年1月月29日日 0:07 */ clear set obs 17 input c pro_asample obsnum_bsample 1 0.106 5232 2 0.137 6826 3 0.086 5232 4 0.102 5630 5 0.09 5282 6 0.076 4684 7 0.063 3538 8 0.057 2890 9 0.054 2591 100.049 2143 110.046 2043 120.036 1694 13
10、0.03 1046 140.03 448 150.019 299 160.013 149 170.006 100 capture end mgof c=pro_asample fweight=obsnum_bsample,mc ksmirnov reps(100) / 结果：Ksmirnov-Smirnov D=0.07201 Exact P-value=0.0000 / 与书中计算结果D=0.072一致 /*ksmirnov 命令 数据格式： 第1列 c：年龄组 第2列 num_sample：甲样本人口实际频数，乙样本人口实际频数 第3列 group：甲、乙样本 */ clear set o
11、bs 34 input c num_sample group 1 58551 2 75671 3 47501 4 56341 5 49711 6 41981 7 34801 8 31481 9 29831 1027071 1125411 1219881 1316571 1416571 -4- C:UsersuserDesktopgoodness of fit2goodness of fit2.do 2016年年1月月29日日 0:07 1510491 16718 1 17331 1 1 52322 2 68262 3 52322 4 56302 5 52822 6 46842 7 35382
12、8 28902 9 25912 1021432 1120432 1216942 1310462 14448 2 15299 2 16149 2 17100 2 capture end expand num_sample ksmirnov c,by(group) / 结果：D=0.0720 Exact P-value=0.0000 / 与书中结算结果D=0.0720一致 /*ksmirnov 命令 手动计算 数据格式： 第1列 c：年龄组 第2列 num_sample：甲样本人口实际频数，乙样本人口实际频数 第3列 group：甲、乙样本 */ clear set obs 17 input c
13、cpro_asample obsnum_bsample 1 58551 2 75671 3 47501 4 56341 5 49711 6 41981 7 34801 8 31481 9 29831 1027071 1125411 1219881 -5- C:UsersuserDesktopgoodness of fit2goodness of fit2.do 2016年年1月月29日日 0:07 1316571 1416571 1510491 16718 1 17331 1 1 52322 2 68262 3 52322 4 56302 5 52822 6 46842 7 35382 8 2
14、8902 9 25912 1021432 1120432 1216942 1310462 14448 2 15299 2 16149 2 17100 2 capture end expand num_sample bysort group:cumul c,generate(proc) equal collapse (first) proc,by(c group) sort group c reshape wide proc, i(c) j(group) gen proclag=proc2-proc1 sum proclag,d di r(max) line proc2 proc1 c,sort
15、 / 结果：D=0.07200599 / 与书中结算结果D=0.0720一致,P=0.05的界值请参考医用统计方法本例公式计算 /*gof chitest chitesti命令 数据格式： 第1列 c：年龄组 第2列 obsnum_asample：乙样本实际频数 第3列 obsnum_bsample：乙样本按照加样本人口构成计算的频数 */ clear set obs 17 input c obsnum_bsample expnum_bsample 1 52325281.662 2 68266826.299 3 52324285.122 4 56305082.354 5 52824484.43
16、 -6- C:UsersuserDesktopgoodness of fit2goodness of fit2.do 2016年年1月月29日日 0:07 6 46843786.852 7 35383139.101 8 28902840.139 9 25912690.658 1021432441.523 1120432292.042 1216941793.772 1310461494.81 14448 1494.81 15299 946.713 16149 647.751 17100 298.962 capture end gof obsnum_bsample expnum_bsample / 结果：pearsons X2=2574.882 P Chi2 = 0.0000 / 与上述其他命令运行结果不一致，因为分组不一样，具体见显示的分组 /*tabulate 命令 数据格式： 第1列 c：年龄组 第2列 n
