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

熱線電話：13121318867

登錄
首頁精彩閱讀單因素方差分析（aov）-R版本
單因素方差分析（aov）-R版本
2017-11-03
收藏

單因素方差分析（aov）-R版本

R版本的方差分析

#做方差分析有三個假設(shè)，需要提前進行檢驗。1.每個處理效應(yīng)和隨機誤差是可加的。2.正態(tài)獨立性，檢驗誤差應(yīng)該是正態(tài)分布的。3.方差齊次性。水平間的方差應(yīng)該相等。
數(shù)據(jù)導(dǎo)入
    wheat<- data.frame(
    + x=c(5.4,5.3,6.1,5.1,5.6,5.1,4.5,6,6.6,6.5,5.8,6.2,5.6), a=factor(c(rep(1,3),rep(2,4),rep(3,3),rep(4,3)))) #rep：重復(fù)產(chǎn)生數(shù)據(jù)

正態(tài)分布檢驗

    shapiro.test(wheatx[wheata==1])

Shapiro-Wilk normality test
data: wheatx[wheata == 1]
W = 0.84211, p-value = 0.2196

    shapiro.test(wheatx[wheata==2])

Shapiro-Wilk normality test
data: wheatx[wheata == 2]
W = 0.94078, p-value = 0.6591

    shapiro.test(wheatx[wheata==3])

Shapiro-Wilk normality test
data: wheatx[wheata == 3]
W = 0.87097, p-value = 0.2983

    shapiro.test(wheatx[wheata==4])

Shapiro-Wilk normality test

data: wheatx[wheata == 4]
W = 0.96429, p-value = 0.6369

    正態(tài)檢驗，4組數(shù)據(jù)的p值都大于0.05 都不能拒絕原假設(shè)，每組數(shù)據(jù)是正態(tài)的，所以滿足方差分析的一個條件，接下來做方差齊性檢驗

原假設(shè)這個原假設(shè)H0: 各個水平是等方差的

    bartlett.test(wheatx,wheata)#方差齊性檢驗，

Bartlett test of homogeneity of variances

    1
    2

data: wheatxandwheata
Bartlett’s K-squared = 0.44019, df = 3, p-value = 0.9318

    p>0.05，不能拒絕原假設(shè)，具有方差齊性，可以進行方差分析。

方差分析

    wheat.aov <- aov(x~a, data=wheat)
    summary(wheat.aov)
    Df Sum Sq Mean Sq F value Pr(>F)
    a 3 3.002 1.0007 6.523 0.0123 *
    Residuals 9 1.381 0.1534

Signif. codes: 0 ‘’ 0.001 ‘’ 0.01 ‘’ 0.05 ‘.’ 0.1 ‘ ’ 1

    p<0.05 ，拒絕原假設(shè)，認為濃度對wheat有影響。

R代碼

wheat<- data.frame(  
    x=c(5.4,5.3,6.1,5.1,5.6,5.1,4.5,6,6.6,6.5,5.8,6.2,5.6), a=factor(c(rep(1,3),rep(2,4),rep(3,3),rep(4,3))))   


#正態(tài)分布檢驗
shapiro.test(wheat$x[wheat$a==1])
shapiro.test(wheat$x[wheat$a==2])
shapiro.test(wheat$x[wheat$a==3])
shapiro.test(wheat$x[wheat$a==4])
#正態(tài)檢驗，4組數(shù)據(jù)的p值都大于0.05 都不能拒絕原假設(shè)，每組數(shù)據(jù)是正態(tài)的，所以滿足方差分析的一個條件，接下來做方差齊性檢驗     

#原假設(shè)這個原假設(shè)H0: 各個水平是等方差的
bartlett.test(wheat$x, wheat$a)#方差齊性檢驗，
#p>0.05，不能拒絕原假設(shè)，具有方差齊性，可以進行方差分析。

wheat.aov <- aov(x~a, data=wheat)#方差分析
summary(wheat.aov)
#p<0.05 拒絕原假設(shè)，認為濃度對wheat有影響。

數(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)用相應(yīng)的接口 initGeetest({ // 以下 4 個配置參數(shù)為必須，不能缺少 gt: data.gt, challenge: data.challenge, offline: !data.success, // 表示用戶后臺檢測極驗服務(wù)器是否宕機 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); }