1、statasave 命令FileSave As例 1. 表 1.为某一降压药临床试验数据,试从键盘输入 Stata,并保存为 Stata 格式文件。STATA 数据库的维护排序SORT 变量名 1 变量名 2 变量更名rename 原变量名 新变量名STATA 数据库的维护删除变量或记录drop x1 x2 /* 删除变量 x1 和 x2drop x1-x5 /* 删除数据库中介于 x1 和 x5 间的所有变量(包括 x1 和 x5)drop if x0 /* 产生新变量 y,其值为所有 x0 的对数值 log(x),当x0 /* 保留 x0 的所有记录,其余记录删除STATA 数据库的维护替
2、换已存在的变量值replace 变量表达式replace bolck=6 if block=0 /* 将 block0 的数全部替换为 6。replace z=. if z 14.02t = -2.7698 t = -2.7698 t = -2.7698P |t| = 0.0218 P t = 0.9891STATA 命令ttesti #obs #mean #sd #valttesti 10 12.59 1.632619 14.02配对设计 t 检验ttest 变量 1=变量 2STATA 命令ttest x1=x2STATA 结果ttest x1=x2Paired t test-Variab
3、le | Obs Mean Std. Err. Std. Dev. 95% Conf. Interval-+-x1 | 10 12.59 .5162794 1.632619 11.42209 13.75791x2 | 10 13.27 .3415813 1.080175 12.49729 14.04271-+- diff | 10 -.6799999.5204272 1.645735 -1.857288 .4972881-Ho: mean(x1 - x2) = mean(diff) = 0Ha: mean(diff) 0t = -1.3066 t = -1.3066 t = -1.3066P
4、|t| = 0.2237 P t = 0.8881成组设计 t 检验ttest 变量 1=变量 2, unpaired unequalttest 变量,by(分组变量)unequalttesti #obs1 #mean1 #sd1 #obs2 #mean2 #sd2 ,unequal unpaired 表示非配对的,如不选就作配对检验unequal 表示假设两组方差不齐,如不选表示假设两组方差达到齐性例(成组设计)分别测得 14 例老年人煤饼病人及 11 例正常人的尿中 17 酮类固醇排出量(mg/dl)如下,试比较两组的均数有无差别STATA 命令ttest x1=x2, unpairedt
5、test x, by(g)STATA 结果ttest x1=x2,unpTwo-sample t test with equal variances-Variable | Obs Mean Std. Err. Std. Dev. 95% Conf. Interval-+-x1 | 14 4.377857 .3875 1.449892 3.540714 5.215x2 | 11 5.528182 .5232431 1.735401 4.362324 6.69404-+-combined | 25 4.884 .3306453 1.653227 4.201582 5.566418-+-diff |
6、 -1.150325 .636752 -2.467547 .1668972-Degrees of freedom: 23Ho: mean(x1) - mean(x2) = diff = 0Ha: diff 0t = -1.8066 t = -1.8066 t = -1.8066P |t| = 0.0839 P t = 0.9580两组资料间的方差齐性检验sdtest 变量名 1 = 变量名 2sdtest 变量,by(分组变量)sdtesti #obs1 #mean1 #sd1 #obs2 #mean2 #sd2单因素方差分析及方差齐性检验oneway 因变量 分组变量,选择项noanova
7、/* 不打印方差分析表missing /* 将缺省值作为单独的一组tabulate /* 打印各组的基本统计量表 简写:tscheffe /* Scheffe 法 简写:schBonferroni /* Bonferroni 法 简写:bonsidak /* Sidak 法 简写:si各组均数两两比较oneway x group,noanova schComparison of var3 by group(Scheffe)Row Mean-|Col Mean | 1 2-+-2 | -.425| 0.426|3 | -.91 -.485| 0.024 0.330STATA 软件及其应用-III
8、秩和检验和相关与回归分析秩变换配对资料的秩和检验;两组资料的秩和检验;多组资料的秩和检验;直线相关分析;等级相关分析;直线回归分析;秩变换genrank 新变量= 原变量egen 新变量=rank(原变量)配对资料的秩和检验signrank 变量 1 = 变量 2 if 变量 1!=变量 2两组资料的秩和检验ranksum 观察值变量,by(分组变量)两组资料的秩和检验两组资料的秩和检验例 3 用复方猪胆胶囊治疗老年性慢性支气管炎患者 403 例,疗效见第(1) (3)栏。问该药对此两型支气管炎疗效是否相同?两组资料的秩和检验expand franksum x,by(g)多组资料的秩和检验kwallis 观察值变量, by(分组变量)多组资料的秩和检验直线相关和回归分析correlate 变量名pwcorr 变量名 , sig star(#) regress 因变量 自变量predict 新变量 , stdp stdfstdp 计算估计值的标准误 stdf 估计预测值 y 的标准差sig 打印相关系数假设检验之 P 值star(#) 如果相关系数的假设检验之 P 值小于#,则在相关系数旁打印星号