【stata系列】——組間係數差異檢驗

2021-02-21 弛悟

        常見的三種組間係數差異檢驗的方法是:引入交叉項(chow檢驗)、基於似無相關模型的檢驗方法(suest)和費舍爾組合檢驗(permutation test)。這三種方法不同於ttest單變量分析;前者是組間係數檢驗,後者是單變量(均值/中位數)差異檢驗。

 

        (1)引入交叉項(chow檢驗):直接在線性回歸中加入交乘項然後檢驗交乘項的係數顯著性即可,這是用得最多最廣泛的一種方法;也叫做大名鼎鼎的chow檢驗,可以通過chowtest 命令快速完成。使用的前提條件是:假設幹擾項同方差、獨立同分布;而且假設控制變量係數在兩組之間無明顯差異。若不滿足前提條件,如變量係數在兩組之間存在差異,或存在異方差的情形,用chow檢驗就會存在問題。變量係數在兩組之間存在差異的解決方法是可以加入更多的交乘項,存在異方差的情形的解決方法是用robust或cluster聚類穩健標準誤來解決。

 

        (2)基於似無相關模型的檢驗方法(suest):在假設兩組樣本的幹擾項相關的前提下進行似無相關模型檢驗。前提條件比chow檢驗更為寬鬆,可以允許兩個幹擾項有不同的分布。在stata中根據suest命令操作即可。需要注意的是,suest適用於截面數據,不適用於面板數據;不能用xtreg個體固定效應模型進行檢驗;除此之外,suest中不能用robust進行異方差修正(這一點無關緊要,因為robust不會影響回歸係數,不加也不影響test結果)。

 

reg y x1 x2 if group==0

est store a

reg y x1 x2 if group==0

est store b

suest a b

// 對兩組樣本進行suest似無相關估計;

test [a_mean]y=[b_mean]y

test [a_mean]x1=[b_mean]x1

test [a_mean]x2=[b_mean]x2

//  檢驗組間係數差異;

 

        (3)費舍爾組合檢驗(permutation test):原理是bootstrap抽樣,檢驗兩組係數差異是否異於零,所以每次結果可能不一樣。具體做法是先算出兩組係數差異,然後通過估計該統計量所處分布來推出經驗p值,判斷組間差異的顯著性(在stata中可通過bdiff命令來實現)。費舍爾檢驗要求兩組的虛擬變量個數一樣(i.ind後各個ind就成了虛擬變量,保證兩組各個ind內個數一樣,要求太苛刻);

 

bdiff, group(xx) model(y x1 x2) reps(1000) detail

// bdiff檢驗;

// group括號內必須是二元虛擬變量,若有多組,需刪除超過0-1之外的其他組;

// model括號內是回歸模型;

// reps括號內是bootstraps抽樣次數,可選擇1000次-5000次;

// detail選項可進一步呈現兩組的實際估計係數。

相關焦點

  • 擴增子統計繪圖6韋恩圖:比較組間共有和特有OTU或分類單元
    擴增子系列課程也分為三部曲:第一部《擴增子圖表解讀》:加速大家對同行文章的解讀能力。第二部《擴增子分析解讀》:學習數據分析的基本思路和流程。第三部《擴增子統計繪圖》:即是對結果進行可視和統計檢驗,達到出版級的圖表結果。
  • Stata:合成控制法 (Synthetic Control Method)操作及應用
    為此,Abadie et al. (2010)進行了一系列的安慰劑檢驗,依次將donor pool中的每個州作為假想的處理地區(假設也在1988年通過控煙法),而將加州作為控制地區對待,然後使用合成控制法估計其「控煙效應」,也稱為「安慰劑效應」。通過這一系列的安慰劑檢驗,即可得到安慰劑效應的分布,並將加州的處理效應與之對比。
  • 基尼係數
    Income & Wealth Gini,貧富差異界線
  • 委託檢驗、型式檢驗、監督檢驗的區別
    監督檢驗:是政府行為,一般由各級質監系統對市場上銷售的產品進行抽樣並送檢,委託人是抽樣的質監局,檢測結果也是僅對樣品負責,檢測費用由政府專項財政經費承擔。 委託檢驗:檢驗結果僅對送樣的樣品負責,屬於商業民事行為,而且是按照委託人的要求進行檢驗,不具有權威性;檢驗費用是由委託人承擔。
  • Excel函數應用篇:函數計算年金終值係數-附年金終值係數表
    一、年金終值係數指固定的間隔時間相等的期間(如以年為單位)分期支付(存入)1元金額,經過若干年後按複利計算的累計本利之和。
  • 微言041:基尼係數
    習近平總書記在二十國集團領導人第十一次峰會的開幕致辭中表示,現在世界的基尼係數已經達到0.7左右,超過了公認的0.6危險線
  • 互助問答第200期:中介效應的stata操作問題
    stata做中介效應的相關命令較多,諸如medeff,paramed,medsem,rmpw和vaf (分解)、ml_mediation (跨層的中介效應)、sgmediation
  • 常用建築材料的太陽輻射反射係數
    實際在綠色建築設計和評價中,往往存在疑問的是:什麼樣的面層或鋪裝材料,才能滿足太陽輻射反射係數不小於0.4的要求?遺憾的是綠建評價標準、評價細則、評價技術圖集均未能給出常用材料的太陽輻射反射係數清單。我們嘗試在一些官方資料中尋找依據,有如下線索可供綠色建築設計和評價參考:1.
  • 小經驗:用Excel表格繪製檢驗科室內質控圖(附視頻)
    用Excel表格繪製檢驗科室內質控圖,是在LIS系統不支持質控數據管理,或者基層醫院缺乏資金不能購買專業質控軟體的前提下,做出的折中方案。此處以Microsoft Excel軟體講述如何用Excel表格處理質控中的平均數、標準差SD、變異係數CV值及繪製質控圖。適合零基礎,簡單易學可行,如功能過於複雜反而不易操作上手。
  • 【秘訣】品管如何做好來料檢驗(IQC)工作?
    抽樣嚴謹抽樣檢驗是一個統計學概率存在風險,這種風險可能把不良品漏檢而放到生產線上去,也可能把合格品判退,從而間接造成本公司的損失,所以在抽樣時,一定要儘可能讓樣本反映母體的品質狀況。在抽樣計劃中,根據現行的國際或國家標準,可以控制在5%以下的風險率,但在具體抽取樣本的過程中,可能因個人偏好而造成母體的隨機性不夠,所以在抽樣方法上要特別注意嚴謹性。
  • 【技巧】如何做好IQC來料檢驗工作?一步步交給你
    1、抽樣嚴謹抽樣檢驗是一個統計學概率存在風險,這種風險可能把不良品漏檢而放到生產線上去,也可能把合格品判退,從而間接造成本公司的損失,所以在抽樣時,一定要儘可能讓樣本反映母體的品質狀況。在抽樣計劃中,根據現行的國際或國家標準,可以控制在5%以下的風險率,但在具體抽取樣本的過程中,可能因個人偏好而造成母體的隨機性不夠,所以在抽樣方法上要特別注意嚴謹性。
  • 豬肉檢驗章也分顏色 紅章藍章哪個好?
    市民範女士致電本報熱線,向記者訴說著她的疑惑,她在(浙江嘉興)秀洲區某農貿市場買豬肉時發現,有的豬肉蓋著藍色的檢驗合格章,有的豬肉蓋著紅色檢驗合格章
  • 【可靠度】新版可靠性規範對荷載分項係數調整的影響
    《建築結構可靠性設計統一標準》GB50068-2018將於2019年4月1日起實施,最大的變化是取消了舊版規範中恆載控制、活載控制工況的分項係數區分,統一為恆荷載分項係數1.3、活荷載分項係數1.5。
  • 華語【話題】在奧瑞岡辯論賽裡面,證據檢驗是如何進行的,有什麼原則?
    相當於證據檢驗=提問證人=質詢對方辯手。      (只說證據檢驗部分)而在這個過程中,要注意的是,提問和回答同法庭辯論中一樣並不是臨場的,律師和證人事先都會做極其充分的準備(圍繞證詞、證據。雖然證人同樣不能撒謊,但有無準備區別甚大,同律師),在這個過程中,質詢者的真正目標往往不是很多人所興奮的「直接拆掉某份數據」,而是在評委(法官)面前揭示出這份證據證明論點的真正效力,以及這個論點與對方本身立場的差異(距)。       一場辯論賽是由很多論點組成的,證明這些論點的叫論據,而這個證明的行為叫論證。
  • 檢驗批容量計算
    1、查專業規範確定計算單位:一般按立方、平方、構件數、批次、房間數2、查專業規範確定有、無限值:有限值:有具體量化限制的(如:多少為一個檢驗批
  • 案例分享丨口罩產品的出廠檢驗和型式檢驗分別是哪些?
    請聯繫CIO小助手(微信號:13332878028)本欄目將定期發布口罩產品的出廠檢驗和型式檢驗分別是哪些
  • 中國性別薪酬差異報告:一線城市男女差異最小「同工不同酬」問題依然嚴重
    BOSS直聘發布的《2016年中國性別薪酬差異報告》顯示,中國女性平均月薪僅為男性的77%。
  • 汽油顏色有差異是有貓膩?
    相同標號的油品為什麼顏色有明顯差異?這讓人不禁心生疑竇,「兩桶油」是不是有誰在玩貓膩?針對這一疑問,中國石油大學化工學院王延臻教授在之前接受媒體採訪時給出了以下解釋:出現顏色不同有多方面的原因,一個比方說它存儲系統受到汙染的話,顏色會變深,就是它的油罐或者管道,當深色物質汙染的時候會出現這種情況。第二個就是本身 汽油被氧化後,顏色也是變深。
  • 【檢驗動態】中檢集團檢驗公司與坦尚尼亞標準局PVoC合作項目正式啟動
    2016年12月1日,經坦尚尼亞標準局同意,中檢集團檢驗公司在中國大陸、香港、澳門和臺灣四大區域正式啟動PVoCO業務。今年9月,中檢集團領導、檢驗公司總經理孔祥月和坦尚尼亞標準局(簡稱TBS)局長在京籤署了「CCIC與TBS就關於PVoC項目合作的協議書」,TBS於9月30日在其官網就雙方的合作發布了正式公告。