STATA最常用命令大全

statasave命令FileSaveAs例1.

表1.為某一降壓藥臨床試驗(yàn)數(shù)據(jù)，試從鍵盤輸入Stata，并保存為Stata格式文件。STATA數(shù)據(jù)庫的維護(hù)排序SORT變量名1變量名2……變量更名rename原變量名新變量名STATA數(shù)據(jù)庫的維護(hù)刪除變量或記錄dropx1x2

/*刪除變量x1和x2dropx1-x5

/*刪除數(shù)據(jù)庫中介于x1和x5間的所有變量(包括x1和x5)dropifx<0

/*刪去x1<0的所有記錄dropin10/12

/*刪去第10~12個(gè)記錄dropifx==.

/*刪去x為缺失值的所有記錄dropifx==.|y==.

/*刪去x或y之一為缺失值的所有記錄dropifx==.&y==.

/*刪去x和y同時(shí)為缺失值的所有記錄drop_all

/*刪掉數(shù)據(jù)庫中所有變量和數(shù)據(jù)STATA的變量賦值用generate產(chǎn)生新變量generate新變量＝表達(dá)式generatebh=_n

/*

將數(shù)據(jù)庫的內(nèi)部編號(hào)賦給變量bh。generategroup=int((_n-1)/5)+1

/*

按當(dāng)前數(shù)據(jù)庫的順序，依次產(chǎn)生5個(gè)1，5個(gè)2，5個(gè)3……。直到數(shù)據(jù)庫結(jié)束。generateblock=mod(_n,6)

/*

按當(dāng)前數(shù)據(jù)庫的順序，依次產(chǎn)生1,2,3,4,5,0。generatey=log(x)ifx>0

/*

產(chǎn)生新變量y，其值為所有x>0的對(duì)數(shù)值log(x)，當(dāng)x<=0時(shí)，用缺失值代替。egen產(chǎn)生新變量setobs12egena=seq()

/*產(chǎn)生1到N的自然數(shù)egenb=seq(),b(3)

/*產(chǎn)生一個(gè)序列，每個(gè)元素重復(fù)＃次egenc=seq(),to(4)

/*產(chǎn)生多個(gè)序列，每個(gè)序列從1到＃egend=seq(),f(4)t(6)/*產(chǎn)生多個(gè)序列，每個(gè)序列從#1到＃2encode

字符變量名，gen(新數(shù)值變量名)作用：將字符型變量轉(zhuǎn)化為數(shù)值變量。STATA數(shù)據(jù)庫的維護(hù)保留變量或記錄keepin10/20

/*保留第10～20個(gè)記錄，其余記錄刪除keepx1-x5

/*保留數(shù)據(jù)庫中介于x1和x5間的所有變量(包括x1和x5)，其余變量刪除keepifx>0

/*保留x>0的所有記錄，其余記錄刪除STATA數(shù)據(jù)庫的維護(hù)替換已存在的變量值replace

變量＝表達(dá)式replace

bolck=6ifblock==0

/*

將block＝0的數(shù)全部替換為6。replace

z=.ifz<0

/*

將所有小于0的z值用缺失值代替。replaceage=25in17/*

將第17條記錄中的變量age替換為25。forvarx*:replaceX=0ifX==.

/*

將所有第一個(gè)字母為x的變量替換為0，如果該變量的值為缺失值縱向連接數(shù)據(jù)庫Ex3-3.dta:

x0

x11.

3550

24502.

2000

24003.

3000

18004.

3950

32005.

3800

3250use"E:\教學(xué)\上機(jī)\ex3-2.dta",clearl

x0

x1

g1.

2450

1450

22.

2100

2400

23.

2300

3800

24.

1590

4200

2appendusingE:\教學(xué)\上機(jī)\ex3-3.dtal

x0

x1

g1.

2450

1450

22.

2100

2400

23.

2300

3800

24.

1590

4200

25.

3550

2450

.6.

2000

2400

.7.

3000

1800

.8.

3950

3200

.9.

3800

3250

.橫向聯(lián)接數(shù)據(jù)庫Ex3-5.dta:

bh

y0

y1

x01.

1

35

79.2

22.

3

45

47.4

83.

4

52

34.6

64.

6

66

28.0

9命令.drop_all.useE:\教學(xué)\MPH上機(jī)\ex3-5.dta.sortbh.save"E:\教學(xué)\MPH上機(jī)\ex3-5.dta",replacefileE:\教學(xué)\MPH上機(jī)\ex3-5.dtasaved.useE:\教學(xué)\MPH上機(jī)\ex3-4.dta.sortbh.mergebhusingE:\教學(xué)\MPH上機(jī)\ex3-5.dta結(jié)果

bh

x0

x1

y0

y1

_merge1.

1

12

24

35

79.2

32.

2

15

26

.

.

13.

3

16

49

45

47.4

34.

4

18

57

52

34.6

35.

5

20

68

.

.

16.

6

9

.

66

28

2列數(shù)據(jù)接龍Stack變量名，into（新變量名）|group(#)[clearwide]示例統(tǒng)計(jì)描述及區(qū)間估計(jì)定量資料的一般描述均數(shù)、標(biāo)準(zhǔn)差、百分位數(shù)、中位數(shù)summarize

[變量名]

[,detail]統(tǒng)計(jì)描述及區(qū)間估計(jì)百分位數(shù)centile

[變量名]

[,centile(#[#...])ccinormalmeansdlevel(#)]統(tǒng)計(jì)描述及區(qū)間估計(jì)定性資料的一般描述按照分類變量給出頻數(shù)和構(gòu)成比tabulate變量名例2.

有三組(group)患者，男女(sex)若干人，sex=1表示男性，sex=0表示女性。測(cè)得其血紅蛋白濃度(x1,%)和紅細(xì)胞計(jì)數(shù)(x2,萬/mm3)，資料存入c:\mydata\ex2.dta。試對(duì)其進(jìn)行描述。見ex5-2.usec:\mydata\ex2.tabgroup.tabsex.tabgroup,sum(x1).tabgroup,sum(x2)統(tǒng)計(jì)描述及區(qū)間估計(jì)可信區(qū)間的估計(jì)ci

變量

[,level(#)binomialpoissonexposure(觀察數(shù)變量)by(分組變量)total]cii

觀察數(shù)均數(shù)標(biāo)準(zhǔn)差[,level(#)]

level(#)

/*指定可信度，缺失時(shí)為95(%)by(分組變量)

/*指定按分組變量分別估計(jì)均數(shù)的可信區(qū)間total

/*指定除按分組變量估計(jì)可信區(qū)間外，還對(duì)整個(gè)數(shù)據(jù)估計(jì)，僅用于指定了by(分組變量)時(shí).usec:\mydata\ex2.sortgroup/*在用by(分組變量)前，必須對(duì)分組變量排序.cix1x2,by(group)STATA的作圖作圖命令GRAPHgraph

[變量名][,圖形類型通用選擇項(xiàng)特殊選擇項(xiàng)]圖形類型histogram/*直方圖，為缺省值。oneway/*一維散點(diǎn)圖twoway/*二維散點(diǎn)圖、線圖matrix/*二維散點(diǎn)圖陣bar/*條圖、百分條圖pie/*圓（餅）圖box/*箱式圖star/*星形圖STATA的作圖作圖命令GRAPH常用選項(xiàng)bin(#)

/*將數(shù)據(jù)分幾組,缺省為5。freq

/*指定縱軸用頻數(shù)表示，否則為頻率。normal/*給直方圖加上相應(yīng)正態(tài)曲線。xlab/ylab/[(#,……,#)]/*指定坐標(biāo)軸的界點(diǎn)。b2/l2[（“字符串”）]/*指定坐標(biāo)軸的副標(biāo)題。STATA的作圖如何利用STATA繪制頻數(shù)分布圖？

例

130名14歲女孩身高資料。grax,bin(10)freqnormalxlab(124,128,132,136,140,144,148,152,156,160,164)ylab(5,10,15,20,25,30,35,40)數(shù)值變量資料的描述均數(shù)、幾何均數(shù)、中位數(shù)、百分位數(shù)極差、四分位數(shù)間距、方差、標(biāo)準(zhǔn)差

變異系數(shù)對(duì)稱分布

均數(shù)±標(biāo)準(zhǔn)差偏態(tài)分布

中位數(shù)±四分位數(shù)間距數(shù)值變量資料的描述means

[變量名]summarize

[變量名]

[,detail]centile[變量名]

[,centile(#)]其他選項(xiàng)detail

/*

詳細(xì)描述，缺失時(shí)為簡(jiǎn)單描述centile(#)

/*

指定需要計(jì)算的百分位數(shù)某市1997年12歲男童120人的身高(cm)資料如下sum

xsumx,dsumxifx<140sumxifx<140,dcentilexcentilex,centile(25,50,75)例有五份血清的抗體效價(jià)為1:10,1:20,1:40,1:80,1:160,描述其抗體滴度的平均水平。means

xSTATA的作圖作圖命令graph

簡(jiǎn)寫gragra

[變量名][,圖形類型通用選擇項(xiàng)特殊選擇項(xiàng)]圖形類型histogram/*直方圖oneway/*一維散點(diǎn)圖twoway/*二維散點(diǎn)圖、線圖matrix/*二維散點(diǎn)圖陣bar

/*條圖、百分條圖pie

/*圓（餅）圖box

/*箱式圖star

/*星形圖直方圖數(shù)值變量資料的統(tǒng)計(jì)分析樣本均數(shù)與總體均數(shù)比較的t檢驗(yàn)配對(duì)設(shè)計(jì)t檢驗(yàn)成組設(shè)計(jì)t檢驗(yàn)方差齊性檢驗(yàn)樣本均數(shù)與總體均數(shù)比較的t檢驗(yàn)ttest

變量名=#valttesti

#obs

#mean

#sd

#val例問題：統(tǒng)計(jì)量與參數(shù)不同的兩種可能其一：抽樣誤差

(偶然的、隨機(jī)的、較小的)其二：本質(zhì)上的差別

(必然的、大于隨機(jī)誤差)例樣本：某醫(yī)生隨機(jī)抽查10名某病患者的血紅蛋白，求得其均數(shù)為12.59(g/dl),標(biāo)準(zhǔn)差為1.632619(g/dl)。問題：該病患者的平均Hb含量是否與正常人的平均Hb含量相同(正常人的平均Hb含量為14.02(g/dl)。STATA命令ttest

變量名=#valttest

x

=14.02STATA結(jié)果ttestx=14.02One-samplettest----------------------------------------------------------------------------Variable|

Obs

Mean

Std.Err.

Std.Dev.

[95%Conf.Interval]---------+------------------------------------------------------------------

x|

10

12.59

.5162794

1.632619

11.42209

13.75791----------------------------------------------------------------------------Degreesoffreedom:9

Ho:mean(x)=14.02Ha:mean<14.02

Ha:mean~=14.02

Ha:mean>14.02

t=

-2.7698

t=

-2.7698

t=

-2.7698P<t=

0.0109

P>|t|=

0.0218

P>t=

0.9891STATA命令ttesti

#obs

#mean

#sd

#val

ttesti

10

12.59

1.632619

14.02配對(duì)設(shè)計(jì)t檢驗(yàn)ttest

變量1=變量2STATA命令ttest

x1=x2STATA結(jié)果ttestx1=x2Pairedttest-------------------------------------------------------------------Variable|

Obs

Mean

Std.Err.

Std.Dev.

[95%Conf.Interval]---------+---------------------------------------------------------

x1|

10

12.59

.5162794

1.632619

11.42209

13.75791

x2|

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)<0

Ha:mean(diff)~=0

Ha:mean(diff)>0t=

-1.3066

t=

-1.3066

t=

-1.3066P<t=0.1119

P>|t|=0.2237

P>t=0.8881成組設(shè)計(jì)t檢驗(yàn)ttest

變量1=變量2,unpaired

[unequal]ttest變量,by（分組變量）[unequal]ttesti

#obs1#mean1#sd1

#obs2#mean2#sd2

[,unequal]

unpaired

表示非配對(duì)的，如不選就作配對(duì)ｔ檢驗(yàn)unequal

表示假設(shè)兩組方差不齊，如不選表示假設(shè)兩組方差達(dá)到齊性例（成組設(shè)計(jì)）分別測(cè)得14例老年人煤餅病人及11例正常人的尿中17酮類固醇排出量(mg/dl)如下，試比較兩組的均數(shù)有無差別STATA命令ttestx1=x2,unpairedttestx,by(g)STATA結(jié)果ttestx1=x2,unpTwo-samplettestwithequalvariances----------------------------------------------------------------------------Variable|

Obs

Mean

Std.Err.

Std.Dev.

[95%Conf.Interval]---------+------------------------------------------------------------------

x1|

14

4.377857

.3875

1.449892

3.540714

5.215

x2|

11

5.528182

.5232431

1.735401

4.362324

6.69404---------+------------------------------------------------------------------combined|

25

4.884

.3306453

1.653227

4.201582

5.566418---------+------------------------------------------------------------------diff|

-1.150325

.636752

-2.467547

.1668972----------------------------------------------------------------------------Degreesoffreedom:23

Ho:mean(x1)-mean(x2)=diff=0

Ha:diff<0

Ha:diff~=0

Ha:diff>0

t=

-1.8066

t=

-1.8066

t=

-1.8066P<t=

0.0420

P>|t|=

0.0839

P>t=

0.9580兩組資料間的方差齊性檢驗(yàn)sdtest變量名1=變量名2sdtest變量,by（分組變量）sdtesti

#obs1#mean1#sd1#obs2

#mean2#sd2單因素方差分析及方差齊性檢驗(yàn)oneway因變量

分組變量，[選擇項(xiàng)]noanova

/*不打印方差分析表missing

/*將缺省值作為單獨(dú)的一組tabulate

/*打印各組的基本統(tǒng)計(jì)量表

簡(jiǎn)寫：tscheffe

/*Scheffe法

簡(jiǎn)寫：schBonferroni

/*Bonferroni法

簡(jiǎn)寫：bonsidak

/*Sidak法

簡(jiǎn)寫：si各組均數(shù)兩兩比較oneway