如何獲得Stata中兩個比例的差異的置信區間

例如

webuse highschool 
tab race sex, col chi2 


    1=white, | 
    2=black, | 1=male, 2=female 
    3=other |  male  female |  Total 
-----------+----------------------+---------- 
    White |  1,702  1,850 |  3,552 
      |  87.82  86.73 |  87.25 
-----------+----------------------+---------- 
    Black |  201  249 |  450 
      |  10.37  11.67 |  11.05 
-----------+----------------------+---------- 
    Other |  35   34 |  69 
      |  1.81  1.59 |  1.69 
-----------+----------------------+---------- 
    Total |  1,938  2,133 |  4,071 
      | 100.00  100.00 | 100.00 

     Pearson chi2(2) = 1.9652 Pr = 0.374

中的白色人種誰是男性和女性的比例不同的是87.82 - 86.73 = 1.09，我想這種差異的置信區間。

來源

2013-05-16 P Sellaz

看看我的回答，請訪問http一個相關的問題：//stats.stackexchange .com/a/49877/7071 –

您可以使用即時形式的-prtest-代替，那就是-prtesti-。

的缺點是，你必須輸入計數和手動比例：

你的榜樣：

prtesti 1702 0.8782 1850 0.8673 

Two-sample test of proportions      x: Number of obs =  1702 
                y: Number of obs =  1850 
------------------------------------------------------------------------------ 
    Variable |  Mean Std. Err.  z P>|z|  [95% Conf. Interval] 
-------------+---------------------------------------------------------------- 
      x |  .8782 .0079276      .8626622 .8937378 
      y |  .8673 .0078874      .851841  .882759 
-------------+---------------------------------------------------------------- 
     diff |  .0109 .0111829      -.0110181 .0328181 
      | under Ho: .0112015  0.97 0.331 
------------------------------------------------------------------------------ 
     diff = prop(x) - prop(y)         z = 0.9731 
    Ho: diff = 0 

    Ha: diff < 0     Ha: diff != 0     Ha: diff > 0 
Pr(Z < z) = 0.8347   Pr(|Z| < |z|) = 0.3305   Pr(Z > z) = 0.1653

來源

2013-05-20 07:41:33

prtest命令是你所需要的：prtest sex, by(race)。你的變量不應該包含兩個以上的組。

webuse highschool 
tab race sex, col chi2 

// dummies 
gen is_black = (race == 2) if race < 3 
gen is_female = (sex == 2) if !mi(sex) 

// proportions test 
prtest is_female, by(is_black)

來源

2013-05-17 17:26:10

謝謝，我給出的例子不是我的實際應用程序 - 我的數據差異遠大於％。但''by'選項似乎不允許'prop' –

錯誤是我的 - 我指的是'prtest'命令，其操作方式與'ttest'命令相似，但比例相同。我完全重寫了我的答案。 –

再次感謝，但我所有的變量/做/採取超過2個層次，所以我得到'XXX帶3個值，而不是2' –

如何獲得Stata中兩個比例的差異的置信區間

回答

相關問題