醫(yī)學(xué)統(tǒng)計(jì)學(xué)實(shí)習(xí)課：STATA軟件的應(yīng)用(1)

1、 醫(yī)學(xué)統(tǒng)計(jì)學(xué)實(shí)習(xí)課STATA軟件的應(yīng)用(1)Stata 簡(jiǎn)介 由美國(guó)計(jì)算機(jī)資源中心（Computer Resource Center）研制；現(xiàn)在為Stata公司的產(chǎn)品；/目前最高版本13.0；操作靈活、簡(jiǎn)單、易學(xué)易用，同時(shí)具有數(shù)據(jù)管理軟件、統(tǒng)計(jì)分析軟件、繪圖軟件、矩陣計(jì)算軟件和程序語言的特點(diǎn) 。Stata 用戶界面結(jié)果窗口變量名窗口命令回顧窗口命令窗口菜單及工具欄Stata 語法格式命令是由關(guān)鍵詞 參數(shù) 選項(xiàng)等構(gòu)成的字符串 *所有命令、函數(shù)、變量名等都區(qū)分大小寫 特殊選項(xiàng) 關(guān)鍵詞 命令參數(shù) ,命令選項(xiàng) 1.關(guān)鍵詞：指明執(zhí)行什么命令，一條命令中必須出現(xiàn)。如排序，sort 2.命令參數(shù): 指明對(duì)哪

2、個(gè)變量進(jìn)行操作。Stata 變量STATA的變量名：由英文字母，數(shù)字，下劃線組成大小寫有區(qū)別 Aa關(guān)鍵字或系統(tǒng)變量不得用作變量名。系統(tǒng)變量 _N（樣本量） _n（內(nèi)部編號(hào)） _pi（圓周率）Stata 數(shù)據(jù)輸入與儲(chǔ)存數(shù)據(jù)輸入直接從鍵盤輸入；使用數(shù)據(jù)編輯窗口輸入；打開已有數(shù)據(jù)文件和拷貝；粘貼方式交互數(shù)據(jù) 。1. 直接從鍵盤輸入命令：input 變量名列表數(shù)據(jù)輸入以“end”結(jié)束不可重復(fù)輸入變量.input x x 1. 1 2. 2 3. 3 4. 4 5. 5 6. end . input x y x y 1. 1 2 2. 3 4 3. 5 6 4. 4 . 5. . 8 6. 8 7 7.

3、 end. input x y x y 1. 1 2 . 6. 8 7 7. end. input x zx already definedr(110);2. 使用數(shù)據(jù)編輯窗口輸入進(jìn)入數(shù)據(jù)編輯窗口直接輸入“edit”在工具條上點(diǎn)擊雙擊每一列的標(biāo)題可以更改變量名每一行代表一條記錄每一列代表一個(gè)變量3. 打開已有數(shù)據(jù)文件4. 與表格式文件交互粘貼與Excel，Word的表格等直接通過復(fù)制粘貼的方式交互數(shù)據(jù)。 在數(shù)據(jù)量不大時(shí)，這種方式操作極為方便。Stata 數(shù)據(jù)輸入與儲(chǔ)存數(shù)據(jù)存儲(chǔ)save命令 save “文件路徑” , replace 文件后綴是“.dta”save c:mydatadata1 ,

4、replaceFile Save As菜單方式命令方式Stata 數(shù)據(jù)打開數(shù)據(jù)打開use命令 use “文件路徑”文件后綴是“.dta”use c:data1File openStata 最基本操作list (縮寫：l)變量列表，列出內(nèi)存中數(shù)據(jù)的內(nèi)容display 表達(dá)式 (縮寫：di) 直接輸出表達(dá)式的計(jì)算結(jié)果計(jì)算di 2+3di sqrt(25) Stata 運(yùn)算Stata的加、減、乘、除及乘方運(yùn)算符依次是+、*、/ 和 。Stata的關(guān)系運(yùn)算符有：=(等于)、(大于)、=(大于等于)、(小于)、0 【產(chǎn)生新變量y，其值為所有x0的對(duì)數(shù)值log(x)，當(dāng)x=0時(shí)，用缺失值代替】Attent

5、ion！Stata 變量賦值egen 新變量=fcn（參數(shù)）,選擇項(xiàng) fcn是功能項(xiàng)： seq(), fill(), rsum(), rmean(), rmin(), rmax(), rfirst(), rlast()例：set obs 12 /*產(chǎn)生一個(gè)空數(shù)據(jù)庫，記錄數(shù)為12 egen a=seq() /*產(chǎn)生1_n的自然數(shù) egen b=seq(),from(1) to(3) block(3) /*重復(fù)3次 egen c=fill(2 5) /*產(chǎn)生等差數(shù)列 egen d=rsum(a b c)/*多個(gè)變量對(duì)應(yīng)觀察值的計(jì)算數(shù)據(jù)庫的維護(hù) 數(shù)據(jù)庫的排序 刪除變量或記錄 保留變量或記錄 替換變量

6、值 變量更名 展開和壓縮數(shù)據(jù) 數(shù)據(jù)結(jié)構(gòu)轉(zhuǎn)換數(shù)據(jù)庫的排序例：ex1.dtasort 變量清單gsort +/-變量清單 sort x 對(duì)變量x中數(shù)值從小到大進(jìn)行排列 gsort x 對(duì)變量x中數(shù)值從小到大進(jìn)行排列 gsort -x 對(duì)變量x中數(shù)值從大到小進(jìn)行排列 Stata 數(shù)據(jù)庫的刪除操作刪除變量或記錄drop x1 x2 /* 刪除變量x1和x2drop x1-x5 /* 刪除介于x1和x5的所有變量(包括x1和x5)drop if x0 /* 刪去x0 /* 保留x0的所有記錄，其余記錄刪除Stata 數(shù)據(jù)庫的替換替換已存在的變量值replace 變量表達(dá)式replace bolck=6

7、if block=0 /* 將block0的數(shù)全部替換為6。replace z=. if z0 /* 將所有小于0的z值用缺失值代替。replace age = 25 in 17/* 將第17條記錄中的變量age替換為25。Stata 數(shù)據(jù)庫的更名變量更名 rename 原變量名 新變量名 rename var1 x 【把var1 更名為 x 】展開和壓縮數(shù)據(jù)expand f /*展開數(shù)據(jù)*/contract a b /*壓縮數(shù)據(jù)*/例： ex2.dtaStata 數(shù)據(jù)庫結(jié)構(gòu)維護(hù)行列互換 xpose, clearvarname format(%fmt) 選項(xiàng)： varname /*將原數(shù)據(jù)庫中

8、變量名放在一個(gè)名為 _varname的變量中 format(%fmt) /*轉(zhuǎn)換同時(shí)，將所有變量轉(zhuǎn)換為規(guī)定格式%fmt 例：ex7.dta xpose ,clear format(%9.1f) xpose ,clear varnameStata 數(shù)據(jù)庫結(jié)構(gòu)維護(hù)列數(shù)據(jù)接龍 stack 變量名，into（新變量名）|group（#）clear wide 例：ex6.dta stack a b c ,into(x) stack a b c d e f, group(3)group(#)表示將原數(shù)據(jù)庫中的變量轉(zhuǎn)為幾個(gè)組練習(xí)輸入數(shù)據(jù)包含變量(exercise.dta)將變量更名為x,y將數(shù)據(jù)按y由大到小

9、排序產(chǎn)生新變量id，使之等于每個(gè)觀察對(duì)象的當(dāng)前順序產(chǎn)生新變量z，當(dāng)y取值小于5時(shí)使z=0，當(dāng)y大于等于5時(shí)，使z=1；刪除y取值缺失的記錄。答案rename var1 x rename var2 ygsort -ygen id=_ngen z=1 replace z=0 if y5drop if y=.原始數(shù)據(jù)某地100名成年男子紅細(xì)胞計(jì)數(shù)（1012個(gè)/升 ）資料如下:4.095.335.624.635.184.275.073.603.315.324.884.314.125.334.404.793.925.464.815.094.204.133.944.415.264.665.295.235.

10、583.534.544.684.484.404.764.814.574.973.945.484.275.105.785.123.604.014.755.806.015.505.364.184.334.844.744.604.764.584.344.724.813.844.174.853.294.914.454.434.994.494.355.265.045.384.935.414.523.864.994.244.504.924.135.055.145.055.174.555.425.704.676.184.375.404.154.084.714.124.794.89頻數(shù)圖的制作例 ex3.dt

11、a use d:ex3.dta“ histogram x ,bin(10) freq xlab(3.2(0.3)6.2) ylab(0(5)20) norm l1(人數(shù)) b1(紅細(xì)胞)bin(#) 指定直方圖中的分組數(shù)，缺省值為5freq 指定以頻數(shù)為縱軸畫圖，缺省時(shí)為以頻率為縱軸xlab ylab 指定坐標(biāo)軸刻度b1()l1() 在圖形下方添加標(biāo)題在圖形左側(cè)添加標(biāo)題normal 在直方圖上加上正態(tài)分布曲線graph7 x ,bin(10) freq xlab(3.2(0.3)6.2) ylab(0(5)20) norm l1(人數(shù)) b1(紅細(xì)胞) 練習(xí)（ex5.dta） 某市1982年1

12、10名7歲男童的身高（cm）試畫出頻數(shù)分布圖 數(shù)據(jù)集（ex5.dta）112.4119.3120.1116.4117.0110.8125.1115.4122.5121.5118.3116.3126.2124.4112.7114.8123.2114.8126.1113.2119.8123.5110.3122.8124.4122.0113.0116.3120.0117.7121.5116.8132.5117.4114.9121.8119.7125.6118.4123.8130.5128.0118.1122.4118.4117.2114.3123.0121.0119.5129.1120.0119.7120.6120.7112.3129.3114.7120.8119.6125.8114.9122.9118.9123.1120.1124.1120.7120.7120.2113.0116.3126.6120.0127.1122.8120.4124.1116.8112.2118.0118.5117.1119.8108.2114.3117.9114.3122.7122.1111.5125.2125.2120.8113.2120.5119.1116.0116.21