99999久久久久久亚洲,欧美人与禽猛交狂配,高清日韩av在线影院,一个人在线高清免费观看,啦啦啦在线视频免费观看www

熱線電話:13121318867

登錄
首頁精彩閱讀R語言實現(xiàn)常用多重比較方法
R語言實現(xiàn)常用多重比較方法
2018-03-03
收藏

R語言實現(xiàn)常用多重比較方法

在單因素方差分析ANOVA中,如果該因素影響比較顯著,那么需要進一步利用多重比較方法比較該因素不同水平的影響,確定不同水平下該因素的影響是否顯著。常見的多重比較方法主要有兩種,LSD法和Tukey HSD法。下面對R語言中,這兩種多重比較方法的實現(xiàn)進行舉例。
前期數(shù)據(jù)如下,影響因素為group,指標為value:
> head(tarD)
                   value  group      sample   time
A0522W11NC1 0.0002053745 normal A0522W11NC1 11week
A0522W11NC2 0.0031773712 normal A0522W11NC2 11week
A0522W11NC3 0.0060378288 normal A0522W11NC3 11week
A0522W11NC4 0.0017626931 normal A0522W11NC4 11week
A0522W11NC5 0.0018035261 normal A0522W11NC5 11week
A0522W11NC6 0.0036690067 normal A0522W11NC6 11week

> tmp <- aov(value ~ group, tarD)
最小顯著差數(shù)檢驗法(LSD法)

> res <- LSD.test(tmp, 'group', p.adj = 'bonferroni')
> print(res$groups)
       trt        means M
1   normal 2.576910e-03 a
2    drug3 7.552555e-04 b
3    drug2 7.269247e-05 b
4 high_fat 6.220610e-05 b
5    drug1 2.954733e-05 b
Tukey氏固定差距檢驗法(Tukey HSD)
> TukeyHSD(tmp)
  Tukey multiple comparisons of means
    95% family-wise confidence level

Fit: aov(formula = value ~ group, data = tarD)

$group
                         diff           lwr          upr     p adj
drug2-drug1      4.314514e-05 -0.0015468705 0.0016331608 0.9999916
drug3-drug1      7.257082e-04 -0.0008643074 0.0023157239 0.6929965
high_fat-drug1   3.265877e-05 -0.0015149488 0.0015802664 0.9999969
normal-drug1     2.547362e-03  0.0009997549 0.0040949700 0.0002613
drug3-drug2      6.825631e-04 -0.0009487586 0.0023138847 0.7563196
high_fat-drug2  -1.048637e-05 -0.0016005020 0.0015795293 1.0000000
normal-drug2     2.504217e-03  0.0009142017 0.0040942330 0.0004945
high_fat-drug3  -6.930494e-04 -0.0022830651 0.0008969662 0.7277757
normal-drug3     1.821654e-03  0.0002316386 0.0034116699 0.0175538
normal-high_fat  2.514704e-03  0.0009670961 0.0040623113 0.0003161
> TukeyHSD(tmp)$group
                         diff           lwr          upr        p adj
drug2-drug1      4.314514e-05 -0.0015468705 0.0016331608 0.9999915820
drug3-drug1      7.257082e-04 -0.0008643074 0.0023157239 0.6929965170
high_fat-drug1   3.265877e-05 -0.0015149488 0.0015802664 0.9999969171
normal-drug1     2.547362e-03  0.0009997549 0.0040949700 0.0002612744
drug3-drug2      6.825631e-04 -0.0009487586 0.0023138847 0.7563195891
high_fat-drug2  -1.048637e-05 -0.0016005020 0.0015795293 0.9999999705
normal-drug2     2.504217e-03  0.0009142017 0.0040942330 0.0004944674
high_fat-drug3  -6.930494e-04 -0.0022830651 0.0008969662 0.7277757202
normal-drug3     1.821654e-03  0.0002316386 0.0034116699 0.0175537862
normal-high_fat  2.514704e-03  0.0009670961 0.0040623113 0.0003161003

數(shù)據(jù)分析咨詢請掃描二維碼

若不方便掃碼,搜微信號:CDAshujufenxi

數(shù)據(jù)分析師資訊
更多

OK
客服在線
立即咨詢
客服在線
立即咨詢
') } function initGt() { var handler = function (captchaObj) { captchaObj.appendTo('#captcha'); captchaObj.onReady(function () { $("#wait").hide(); }).onSuccess(function(){ $('.getcheckcode').removeClass('dis'); $('.getcheckcode').trigger('click'); }); window.captchaObj = captchaObj; }; $('#captcha').show(); $.ajax({ url: "/login/gtstart?t=" + (new Date()).getTime(), // 加隨機數(shù)防止緩存 type: "get", dataType: "json", success: function (data) { $('#text').hide(); $('#wait').show(); // 調(diào)用 initGeetest 進行初始化 // 參數(shù)1:配置參數(shù) // 參數(shù)2:回調(diào),回調(diào)的第一個參數(shù)驗證碼對象,之后可以使用它調(diào)用相應的接口 initGeetest({ // 以下 4 個配置參數(shù)為必須,不能缺少 gt: data.gt, challenge: data.challenge, offline: !data.success, // 表示用戶后臺檢測極驗服務器是否宕機 new_captcha: data.new_captcha, // 用于宕機時表示是新驗證碼的宕機 product: "float", // 產(chǎn)品形式,包括:float,popup width: "280px", https: true // 更多配置參數(shù)說明請參見:http://docs.geetest.com/install/client/web-front/ }, handler); } }); } function codeCutdown() { if(_wait == 0){ //倒計時完成 $(".getcheckcode").removeClass('dis').html("重新獲取"); }else{ $(".getcheckcode").addClass('dis').html("重新獲取("+_wait+"s)"); _wait--; setTimeout(function () { codeCutdown(); },1000); } } function inputValidate(ele,telInput) { var oInput = ele; var inputVal = oInput.val(); var oType = ele.attr('data-type'); var oEtag = $('#etag').val(); var oErr = oInput.closest('.form_box').next('.err_txt'); var empTxt = '請輸入'+oInput.attr('placeholder')+'!'; var errTxt = '請輸入正確的'+oInput.attr('placeholder')+'!'; var pattern; if(inputVal==""){ if(!telInput){ errFun(oErr,empTxt); } return false; }else { switch (oType){ case 'login_mobile': pattern = /^1[3456789]\d{9}$/; if(inputVal.length==11) { $.ajax({ url: '/login/checkmobile', type: "post", dataType: "json", data: { mobile: inputVal, etag: oEtag, page_ur: window.location.href, page_referer: document.referrer }, success: function (data) { } }); } break; case 'login_yzm': pattern = /^\d{6}$/; break; } if(oType=='login_mobile'){ } if(!!validateFun(pattern,inputVal)){ errFun(oErr,'') if(telInput){ $('.getcheckcode').removeClass('dis'); } }else { if(!telInput) { errFun(oErr, errTxt); }else { $('.getcheckcode').addClass('dis'); } return false; } } return true; } function errFun(obj,msg) { obj.html(msg); if(msg==''){ $('.login_submit').removeClass('dis'); }else { $('.login_submit').addClass('dis'); } } function validateFun(pat,val) { return pat.test(val); }