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

熱線電話：13121318867

登錄
首頁精彩閱讀SAS中的協(xié)方差分析
SAS中的協(xié)方差分析
2017-05-29
收藏

SAS中的協(xié)方差分析

所謂的協(xié)方差分析，就是在方差分析的基礎(chǔ)上加上協(xié)變量這一額外因素，而方差分析則只考慮組變量這一因素。協(xié)變量可以有一個，也可以有多個。

在這篇文章中，我只講述單變量的協(xié)方差分析。在醫(yī)學(xué)上通常用來判斷治療前后帶來的差異性結(jié)果是否與治療前的結(jié)果是否存在線性關(guān)系，如果存在線性關(guān)系，則通過線性模型去掉這一因素所帶來的影響。

以人體增重為例，假如在服用某種藥物之前體重為X，服用藥物之后體重為X1，則減肥效果通常是用Y = X - X1來表示。那么X則稱之為協(xié)變量。因?yàn)樵诒容^不同組間的減肥結(jié)果的時候我們需要去除服藥前體重的不同而帶來的誤差，協(xié)方差模型就是通過對Y和X以及組變量CLASS建立一個線性回歸模型，通過模型來求得X的回歸系數(shù)β，然后通過Y-β(X - X平均值)得到調(diào)整之后的Y，通過這一調(diào)整，使得由于不同的療前體重所帶來的誤差被剔除，相當(dāng)于使得大家在療前都處于同一水平上，進(jìn)而可以以調(diào)整后的Y對組間的減肥效果進(jìn)行方差分析。

在SAS里，可以通過各種過程步來進(jìn)行求解，例如reg過程、glm過程。

我就選glm過程，以下述數(shù)據(jù)作為例子簡單寫一下如何通過SAS來進(jìn)行協(xié)方差分析。


例：

比較三種豬飼料A1，A2，A3對豬增重的影響，測得每頭豬的增重（Y）和出生重（X），數(shù)據(jù)列在表4－3中。問三種飼料對豬增重是否有顯著不同的效果？


表4－3    不同飼料對豬增重的影響



首先確定X為出生體重，組變量為class，反應(yīng)變量為Y，建立回歸模型，如下所示：

data tmp;

input x y @@;

class = scan("A1,A2,A3",ceil(_n_/8));

if class ='A1' then do;

k1 =0; k2 =1;

end;

else if class ='A2' then do;

k1 =1; k2 =0;

end;

else do;

k1 =0; k2 =0;

end;

cards;

16 8513 83 11 6512 76 12 80 16 91 14 84 17 90

17 9716 9018 10018 9521 10322 10619 9918 94

22 8924 9120 8323 9525 10027 10230 10532 110

;

run;

ods output ParameterEstimates = stat;

proc glm data = tmp;

model y =  x k1 k2;

run;

ods output close;

得到結(jié)果如下所示：


模型的x的回歸系數(shù)在0.05的水平上是明顯不為0的，因此可以認(rèn)為x與y存在線性關(guān)系，那么接下來就需要去掉這一因素不同水平差異而帶來的變異了。
proc sql noprint;

create table tmp1 as

        select a.*,mean(a.x) as mean_X,b.Estimate as beta, y -  beta * (a.x - calculated                     mean_X) as y1 label = "調(diào)整后的Y"

          from tmp a,stat b

          where b.Parameter = 'x';
quit;

可以看到上述sql過程得到調(diào)整之后的y，剔除了不同水平差異的X之后，那么接下里就可以進(jìn)行方差分析了，這里就不再贅述了。


數(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(), // 加隨機(jī)數(shù)防止緩存 type: "get", dataType: "json", success: function (data) { $('#text').hide(); $('#wait').show(); // 調(diào)用 initGeetest 進(jìn)行初始化 // 參數(shù)1：配置參數(shù) // 參數(shù)2：回調(diào)，回調(diào)的第一個參數(shù)驗(yàn)證碼對象，之后可以使用它調(diào)用相應(yīng)的接口 initGeetest({ // 以下 4 個配置參數(shù)為必須，不能缺少 gt: data.gt, challenge: data.challenge, offline: !data.success, // 表示用戶后臺檢測極驗(yàn)服務(wù)器是否宕機(jī) new_captcha: data.new_captcha, // 用于宕機(jī)時表示是新驗(yàn)證碼的宕機(jī) 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); }