权数与SUDAAN简介.ppt

權數與SUDAAN簡介 張新儀 4/14/2001 醫療保健政策研究組 國家衛生研究院 權數 分析調查資料所用的權數是該抽樣單位被抽取之機率的倒數，也是該單位在群體所代表的數目 如果抽樣方法是以不等機率抽樣，在估計母群體的一些統計值，必須加權才能得到不偏估計量 SAS加權方法有二 1.直接乘上權數：每個觀測值乘以Wi，得 到的結果代表整個人群， ?Wi＝總人口數 2.乘上權數再調到觀察人數：每個觀測值乘以Wi /?Wi ?N ，?Wi＝總人口，N＝總樣本數 甚麼是SUDAAN？ SUDAAN的全名是SURVEY DATA ANALYSIS，顧名思義，這是專為分析抽樣調查資料的軟體。 是Research Triangle Institute所發展出來的軟體，原來是SAS下的一模組(Module)，現可獨立執行 網址 /units/shsp/sud1.cfm 為甚麼要用SUDAAN？ 一般統計分析都假設是簡單逢機抽樣（SRS），而抽樣調查常常不是SRS，也就是說每個人（家戶）被抽到的機率不是相等的，他們的代表性也不是相等的，所以要以特別方式處理。 抽樣調查資料的特性 a). 每個人（家戶）被抽到的機率不相等 b). 分層：同一層內同質性高 c). 聚集現象 d). 拒訪及其他 SUDAAN主要處理的有 A).權數：用來反應每個抽樣單位被抽到的機率， 也可 以 說是反應該單位在母群體中的代表性； B).抽樣設計：影響變異數（標準差）的估算， 進而影響統計檢定。 抽樣調查資料用一般假設SRS的軟體分析的問題 a). 點估計可用一般軟體加權後得到； b). 該估計值的標準誤常被低估； c). 信賴區間太窄； d). 統計檢定總是顯著。 加權與不加權的統計檢定 ＊所有百分比都經過加權 1. 加權到總人口數後之?2 ＝79441.8, p＝0.001 2. 加權到總人口數後再調整回樣本數， ?2 ＝17.1, p=0.146 3. 加權且考慮抽樣方法(design effect)， ?2 ＝108.0, p=0.0001 BMI (Mean? s.e. ) 營養調查的SUDAAN程式（MUST FOLLOW!） proc sort data=nahsit; by stra;run; /* 按地區層排序 */; data city; /*輸入每一層的鄉鎮市區數 */; input stra $ citycnt; cards; 1 18 2 30 3 19 4 6 5 23 6 64 7 196 ； 營養調查的SUDAAN程式（MUST FOLLOW!） data t; merge nahsit city; by stra; /* 合併地區層和該層之鄉鎮市區數*/; strata=stra+0; /* 將文字型變項改為數值型 */; cityn=city+0; sexn=sex+0; run; proc sort; by strata citycnt; run; /* 排序(SUDAAN要求一定要按層排序)*/; PROC CORSSTAB data=t filetype=sas DESIGN=WOR; /*開始SUDAAN /; NEST STRATA CITYCNT; TOTCNT CITYCNT _MINUS1_; WEIGHT ewt; SUDAAN語法可分為四大部分 1.程序Procedures: 如PROC descriptive, PROC crosstab etc… 2.宣告抽樣設計等相關事項： 如DESIGN=, WEIGHT, NEST, TOTCNT 3.和計算有關的部分： 如SUBGROUP, LEVELS, RECODE, SUBPOPN 4.輸出格式(output)： 如TITLE, FOOTNOTE, SETENV, PRINT, OUTPUT, FORMAT等。 SUDAAN要求的資料格式 SUDAAN 可讀ASCII, SAS, SPSS 等資料檔 所有的變數必須是數值型 資料必