<strike id="rrvnz"></strike>

    <noframes id="rrvnz">
    <noframes id="rrvnz"><span id="rrvnz"><pre id="rrvnz"></pre></span><address id="rrvnz"></address>
    <address id="rrvnz"><nobr id="rrvnz"></nobr></address><noframes id="rrvnz">
      <em id="rrvnz"></em>

        首 頁IT知識庫翔宇問吧收藏內容
        當前位置:翔宇亭IT樂園IT知識庫R語言

        R語言中進行多重均值檢驗的方法

        減小字體 增大字體 作者:biye5u.com  來源:本站原創  發布時間:2019-03-28 19:58:01

        本站在《使用R語言進行單因素方差分析》一文中介紹了在R語言中使用aov函數進行方差分析的方法,并給出具體的例子講解了該函數的使用方法及其正態性檢驗和方差齊性檢驗的方法。在實際中,在進行單因素方差分析時,可能檢驗結果拒絕了零假設,即因素A的r個水平效應具有顯著的差異,也就是說r個均值之間有顯著區別。但這并不能說明r個均值之間都會存在這個顯著差異,只能說明至少2個總體的均值有顯著差異。那么,到底是哪些水平之間存在差異呢?我們需要對每兩個水平間的均值做一一比較,即多重比較。

        同時比較任意兩個水平均值間有無明顯差異的問題稱為多重比較。通常采用多重T檢驗方法進行多重比較,這種方法本質上就是針對每組數據進行T檢驗,只不過估計方差時利用的是全體數據,因而自由度變大。具體地說,要比較第i個總體與第j個總體的均值是否相同,即檢驗:(薛毅,陳立萍.R語言在統計中的應用)

        H0: μi = μj,H1:μi ≠ μj,i ≠ j,i,j=1,2,...,r

        在R中,使用pairwise.t.test()函數可以幫助我們完成多重T檢驗,該函數的原型如下:

        pairwise.t.test(x, g, p.adjust.method = p.adjust.methods, pool.sd = !paired,
        paired = FALSE, alternative = c("two.sided", "less", "greater"), ...)

        其中,主要參數的含義如下:

        x 是一個向量,代表響應變量(因變量);

        g 代表響應向量的因子或向量;

        p.adjust.method 調整p值所使用的方法;

        paired 邏輯變量,表示是否做成對的t檢驗,默認為FALSE,此變量和pool.sd不能同時為TRUE;

        alternative 字符串,表示備擇假設。取值為two.sided(雙側檢驗),less(單側小于),greater(單側大于),默認為two.sided.

        下面用一個例子來說明本函數的具體使用方法:

        為了減小失業保險支出,促進就業,政府試圖為失業者提供再就業獎勵:如果失業者可以在限定的時間內重新就業,他將可以獲得一定數額的獎金,這個政策是否有效呢?下圖給出了不同獎勵下再就業的天數。

        R語言多重T檢驗的例子

        (本例子來源于網絡)

        編寫程序如下:

        #定義數據
        x<-c(92,100,85,88,89,90,94,80,78,
        86,108,93,88,89,75,78,72,79,
        96,92,90,77,79,71,82,75,81,
        78,75,76,87,73,83,82,68,72)
        A<-factor(rep(1:4,c(9,9,9,9)))

        #定義數據框
        dat<-data.frame(x,A)

        #方差分析
        aovResult<-aov(x~A,data=dat)

        #提取方差分析表
        summary(aovResult)

        #多重比較
        pairwise.t.test(x,A)

        方差分析結果如下圖所示:

         多重比表T檢驗

        由上圖方差分析的結果來看,p值=0.0433<0.05,則可以認為有充分理由認為水平間的效用是不同的。

        那么,不同產生在哪些水平間呢,下面是多重比較的結果:

        多重比較的結果

         從上圖可以看出1和4之間的p值為0.04,則可以認為有充足的理由認為水平1和水平4之間差異顯著,即無獎金組和高獎金組的再就業時間有顯著性差異。而其它組之間沒有差異。

        上圖還給出了p值調整所用的方法:holm修正的方法。

        因為在多重T檢驗時,當因素的水平較多,而檢驗又是同時進行的,則多次重復使用T檢驗會增大犯第一類錯誤的概率,得到的“有顯著差異”的結論不可靠。為了克服這一缺陷,統計學家提出了多種方法來調整P值,這里使用的是holm調整方法。

        -------------------------------------------

        本文為本站原創,如需轉載請注明出處:翔宇亭IT樂園(www.sunyaoshi.com)-R語言中進行多重均值檢驗的方法(http://www.sunyaoshi.com/article/R/2019/6402.html)

        微信搜索“優雅的代碼”關注本站的公眾號,或直接使用微信掃描下面二維碼關注本站公眾號,以獲取最新內容。

        個人成長離不開各位的關注,你的關注就是我繼續前行的動力。

        知識評論評論內容只代表網友觀點,與本站立場無關!

           評論摘要(共 0 條,得分 0 分,平均 0 分) 查看完整評論
        愿您的精彩評論引起共鳴,帶來思考和價值。
        用戶名: 查看更多評論
        分 值:100分 90分 80分 70分 60分 40分 20分
        內 容:
        驗證碼:
        關于本站 | 網站幫助 | 廣告合作 | 網站聲明 | 友情連接 | 網站地圖
        本站部分內容來自互聯網,如有侵權,請來信告之,謝謝!
        Copyright © 2007-2021 biye5u.com. All Rights Reserved.
        特黄特色三级在线观看 - 视频 - 在线观看 - 影视资讯 -心晴网