两个群体的比较-Stata教学课件

两个群体的比较-Stata教学因為中文相容性問題有一些亂碼，辨識不易以j2、j3為例j2問受訪者拾.2.通常您平均每週大約花多少時間做家務工作？_ 小時j3問受訪者拾.3.通常您的配偶平均每週大約花多少時間做家務工作？_小時Ming-chi Chen2社會統計我們的資料裡有變數標籤，但是因為相容性的關係會有亂碼查看是否有亂碼？Data-data editor在j2這個變數名稱上click一下，下面一整欄的數值都反白了滑鼠右鍵-variable-properties-label出現的中文是通常您平均牢週大約花多少時間做家務工作把亂碼改好也將j3變數標籤的亂碼改好Ming-chi Chen3社會統計查看變數有無異常值關掉Data editor視窗用box plot來看有無極端值Graphics-easy graphs-box plot-main-在variable的空格裡鍵入j2Ming-chi Chen4社會統計用box plot來看有無極端值Ming-chi Chen5社會統計同樣方法也可以查看j3的極端值也可以直接在指令欄Ming-chi Chen6社會統計這就是指令欄這就是指令欄Ming-chi Chen7社會統計在指令欄裡直接鍵入Graph box j2然後按enterMing-chi Chen8社會統計Summarize varname,detail指令欄鍵入summarize j2,detail或statistics-summaries,tables,&tests-summary statistics-summary statisticsMing-chi Chen9社會統計.通常您平均每週大約花多少時間做家務工作？-Percentiles Smallest 1%0 0 5%0 010%0 0 Obs 192425%2 0 Sum of Wgt.192475%20 998太愛做家事了吧！高得不合理Ming-chi Chen10社會統計Recode極端值我們到85q1_format.txt去看，發現J2 J3 996不知道 998不適用 999拒答所以要把995以上定義為system missingRecode j2 995/max=.這裡的句點.就是Stata系統定義的缺失值。Ming-chi Chen11社會統計.summarize j2,detail 通常您平均每週大約花多少時間做家務工作？-Percentiles Smallest 1%0 0 5%0 010%0 0 Obs 184925%2 0 Sum of Wgt.184975%15 105一週只有168小時，所以應該合理換算，以一天16小時算，一週112小時12.inspect j2j2:通常您平均每週大約花多少時間做家務工作 Number of Observations-Total Integers Nonintegers|#Negative -|#Zero 305 305 -|#Positive 1544 1544 -|#-|#Total 1849 1849 -|#.Missing 75+-0 168 1924 (47 unique values)用inspect來看大致分佈以及缺失個案數Data-describe data-inspect variables13Recode j2 168=112Ming-chi Chen14社會統計.inspect j2j2:通常您平均每週大約花多少時間做家務工作 Number of Observations-Total Integers Nonintegers|#Negative -|#Zero 305 305 -|#Positive 1544 1544 -|#-|#Total 1849 1849 -|#.Missing 75+-0 112 1924 (46 unique values)15.sum j2,detail 通常您平均每週大約花多少時間做家務工作？-Percentiles Smallest 1%0 0 5%0 010%0 0 Obs 184925%2 0 Sum of Wgt.184975%15 10516.inspect j3j3:通常您的配偶平均每週大約花多少時間做家 Number of Observations-Total Integers Nonintegers|#Negative -|#Zero 263 263 -|#Positive 1661 1661 -|#-|#Total 1924 1924 -|#.#Missing -+-0 999 1924 (54 unique values)17.summarize j3,detail通常您的配偶平均每週大約花多少時間做家務工作？Percentiles Smallest1%0 05%0 010%0 0 Obs 192425%4 0 Sum of Wgt.192475%996 99818Missing value&recodeRecode j3 990/max=.Recode j3 168=112Ming-chi Chen19社會統計.recode j3 168=112(j3:4 changes made).inspect j3j3:通常您的配偶平均每週大約花多少時間做家 Number of Observations-Total Integers Nonintegers|#Negative -|#Zero 263 263 -|#Positive 1144 1144 -|#-|#Total 1407 1407 -|#.Missing 517+-0 150 1924 (50 unique values)20.summarize j3,detail通常您的配偶平均每週大約花多少時間做家務工作？Percentiles Smallest1%0 05%0 010%0 0 Obs 140725%2 0 Sum of Wgt.140775%21 11221Recode j3 112/max=112Tabulate j3Ming-chi Chen22社會統計-+-Ming-chi Chen23社會統計來看看男女的差別A1.這題是性別，男是1，女是2。Data-data editor-找的A1這個變數-滑鼠右鍵Variable-properties-label改成性別Value label-define/modify-define-label name輸入gender-OK-value鍵入1-text鍵入男-OKvalue鍵入1-text鍵入男-OK-cancel-close-value label選擇gender-OK關掉Data editor視窗Ming-chi Chen24社會統計男女的家務分擔是否有不同？Statistics-Summaries,tables,&tests-tables-One/Two-way table of summary statistics自變數依變數Ming-chi Chen25社會統計差別很大嗎？|Summary of|通常您平均每週大約花多少時間做家務工作|性別|Mean Std.Dev.Freq.-+-男|6.0485537 10.23684 968 女|18.330306 16.287017 881-+-Total|11.900487 14.791877 1849Ming-chi Chen26社會統計母體變異數未知但已知相等Statistics-Summaries,tables,&tests-Classical tests of hypotheses-Group mean comparison tests依變數自變數信賴水準Ming-chi Chen27社會統計.ttest j2,by(a1)level(99)Two-sample t test with equal variances-Group|Obs Mean Std.Err.Std.Dev.99%Conf.Interval-+-+-+-Ho:diff=0 degrees of freedom=1847 Ha:diff 028母體變異數未知但已知不相等以上的方法是假設母體變異數未知但已知相等。不管樣本大小，統計軟體一般用t檢定那如果母體變異數未知但已知不相等，怎麼辦？Ming-chi Chen29社會統計母體變異數未知但已知不相等Statistics-Summaries,tables,&tests-Classical tests of hypotheses-Group mean comparison tests變異數不相等自由度需要比較複雜，由Welch提出的運算方式Ming-chi Chen30社會統計男女性負擔家務工作時數的差異，在母體變異數未知但已知不等的情況下.ttest j2,by(a1)unequal welch level(99)Two-sample t test with unequal variances-Group|Obs Mean Std.Err.Std.Dev.99%Conf.Interval-+-+-+-Ha:diff 0Ming-chi Chen31社會統計變異數相等與否的Levene檢定Statistics-Summaries,tables,&tests-Classical tests of hypotheses-Group variance comparison tests依變數自變數Ming-chi Chen32社會統計變異數相等與否的Levene檢定.sdtest j2,by(a1)level(99)Variance ratio test-Group|Obs Mean Std.Err.Std.Dev.99%Conf.Interval-+-+-Ho:ratio=1 degrees of freedom=967,880 Ha:ratio 1sd(男)/sd(女)不等於一，p值顯示可以拒斥變異數相等的虛無假設Ming-chi Chen33社會統計根據Levene檢定的結果，選擇變異數不相等的假設比較正確。也就是男性分擔家務的時數顯著地少於女性。Ming-chi Chen34社會統計已婚未婚者的家務工作負擔的比較A5為受訪者的婚姻狀況1為未婚，2為已婚，3為其他已婚者家務負擔比較大嗎？Ming-chi Chen35社會統計已婚未婚者的家務工作負擔的比較仿照男女的比較得到如下的錯誤回報.ttest j2,by(a5)level(99)more than 2 groups found,only 2 allowedr(420);這是因為a5這個變數有三個變數值：未婚、已婚和其他要用條件是來限制，僅比較未婚者和已婚者Ming-chi Chen36社會統計Statistics-Summaries,tables,&tests-Classical tests of hypotheses-Group mean comparison testsMing-chi Chen37社會統計變異數相等.ttest j2 if a5!=3,by(a5)level(99)Two-sample t test with equal variances-Group|Obs Mean Std.Err.Std.Dev.99%Conf.Interval-+-+-+-Ho:diff=0 degrees of freedom=1835 Ha:diff 0Ming-chi Chen38社會統計變異數不相等.ttest j2 if a5!=3,by(a5)unequal welch level(99)Two-sample t test with unequal variances-Group|Obs Mean Std.Err.Std.Dev.99%Conf.Interval-+-+-+-Ha:diff 0Ming-chi Chen39社會統計Levene檢定.sdtest j2 if a5!=3,by(a5)level(99)Variance ratio test-Group|Obs Mean Std.Err.Std.Dev.99%Conf.Interval-+-+-Ho:ratio=1 degrees of freedom=305,1530 Ha:ratio 1無法拒斥變異數相等的虛無假設Ming-chi Chen40社會統計兩層群體的比較已婚男女間，未婚男女間是否有差異？婚姻是否不利於女性（至少就花在家務勞動上的時間而言）？Ming-chi Chen41社會統計變異數相等Statistics-Summaries,tables,&tests-Classical tests of hypotheses-Group mean comparison testsMing-chi Chen42社會統計.by a5,sort:ttest j2 if a5!=3,by(a1)level(99)-a5=未婚Two-sample t test with equal variances-Group|Obs Mean Std.Err.Std.Dev.99%Conf.Interval-+-+-+-Ho:diff=0 degrees of freedom=304 Ha:diff 0多重比較變異數相等43多重比較變異數相等-a5=已婚Two-sample t test with equal variances-Group|Obs Mean Std.Err.Std.Dev.99%Conf.Interval-+-+-+-Ho:diff=0 degrees of freedom=1529 Ha:diff 044多重比較變異數不相等.by a5,sort:ttest j2 if a5!=3,by(a1)unequal welch level(99)-a5=未婚Two-sample t test with unequal variances-Group|Obs Mean Std.Err.Std.Dev.99%Conf.Interval-+-+-+-Ha:diff 045多重比較變異數不相等-a5=已婚Two-sample t test with unequal variances-Group|Obs Mean Std.Err.Std.Dev.99%Conf.Interval-+-+-+-Ha:diff 046多層次比較變異數相等檢定.by a5,sort:sdtest j2 if a5!=3,by(a1)level(99)-a5=未婚Variance ratio test-Group|Obs Mean Std.Err.Std.Dev.99%Conf.Interval-+-+-Ho:ratio=1 degrees of freedom=176,128 Ha:ratio 147多層次比較變異數相等檢定-a5=已婚Variance ratio test-Group|Obs Mean Std.Err.Std.Dev.99%Conf.Interval-+-+-Ho:ratio=1 degrees of freedom=783,746 Ha:ratio 148Box Plot箱型圖的比較Ming-chi Chen49社會統計單身男性和已婚男性是否有差別？單身女性和已婚女性是否有差別？Ming-chi Chen50社會統計配對樣本結婚對女性不利？前例的分析中，我們比較已婚者與未婚者從事家務時間的差異，由此來推論婚前婚後可能產生的變化。但婚前組與婚後組是由不同受訪者所構成的獨立樣本。如果是否結婚與某些人格特質有關，則我們不知道是因為婚姻本身造成行為上的改變，還是具有某種行為傾向的人比較容易選擇婚姻。即我們的分析可能隱藏自我選擇 self-selection的問題Ming-chi Chen51社會統計配對樣本為了證明婚姻對從事家務時間的影響不是來自於自我選擇，更好的分析樣本為長期追蹤資料（longitudinal data），即能追蹤同一個受訪者，在婚前及婚後所產生行為上的變化。但這種樣本資料的蒐集十分費時費力。Ming-chi Chen52社會統計配對樣本夫妻之間從事家務的時間是否有顯著的差異？我們可以用兩種方式來分析：將已婚男性與已婚女性當作兩獨立樣本，比較所有先生的平均值與太太的平均值是否有差異？Ming-chi Chen53社會統計配對樣本但夫妻從事家務的時間不是獨立事件，先生多分擔，太太自然可以少做。因此應該比較同一家庭中，夫與妻從事家務的時間是否有差異，而不是比較所有的夫的平均值與所有妻的平均值。Ming-chi Chen54社會統計Statistics-Summaries,tables,&tests-Classical tests of hypotheses-Mean comparison tests,paired data1st-2ndMing-chi Chen55社會統計夫妻之間的家務分工Paired t test-Variable|Obs Mean Std.Err.Std.Dev.99%Conf.Interval-+-+-diff|1380 -1.517391 .6578304 24.43732 -3.214199 .1794161-Ho:mean(diff)=0 degrees of freedom=1379 Ha:mean(diff)0配偶間相減，但是是妻減夫還是夫減妻？僅知夫妻間有差異比配偶少，且達顯著水準Ming-chi Chen56社會統計配對樣本如果要比較先生與太太從事家務時間的多寡，則應該如何分析？男女分開分析Ming-chi Chen57社會統計產生新的變數並定義其計算式Generate h_work=(j3 j2)Replace h_work=(j2 j3)if a1=2Ming-chi Chen58社會統計One sample mean comparison testStatistics-Summaries,tables,&tests-Classical tests of hypotheses-one sample mean comparison testMing-chi Chen59社會統計.ttest h_work=0,level(99)One-sample t test-Variable|Obs Mean Std.Err.Std.Dev.99%Conf.Interval-+-Ho:mean=0 degrees of freedom=1379 Ha:mean 0已婚女性的負擔Ming-chi Chen60社會統計質化變數（比例）的比較K1問如果母親外出工作，對還沒上學的小孩比較不好。1非常贊成，2贊成，3不贊成，4非常不贊成，5無意見，6不知道，7不瞭解題意，9拒答，0未答Recode k1(1 2=1)(3 4=0)(else=.)把這個依變數變成1和0兩個數值而已。Ming-chi Chen61社會統計Statistics-Summaries,tables,&tests-Classical tests of hypotheses-Group proportion testMing-chi Chen62社會統計.prtest k1,by(a1)level(99)Two-sample test of proportion 男:Number of obs=935 女:Number of obs=861-Variable|Mean Std.Err.z P|z|99%Conf.Interval-+-男|.7754011 .0136478 .7402468 .8105554 女|.7584204 .0145876 .7208453 .7959956-+-diff|.0169806 .0199765 -.0344753 .0684366-Ho:diff=0 Ha:diff 0Ming-chi Chen63社會統計